裂隙岩体渗透性影响因素的研究

采用UDEC离散单元法中关于裂隙岩体开挖模拟及水力全耦合分析模型,分析裂隙岩体洞室开挖后,因围岩应力与水力耦合作用导致裂隙隙宽变化及渗流变化的过程。为了更直观地了解耦合作用对裂隙岩体渗透特性的影响,以隧洞开挖为例,用开挖后隧洞内总涌水量来表征岩体的渗透特性。利用数值试验的方法,研究了块体边界大小、初始应力比、裂隙隙宽和裂隙夹角对开挖后隧洞内涌水量变化的影响,进而可以看出它们对裂隙岩体渗透性的影响。并得出如下结论：随着块体尺寸和初始应力比的增大,隧洞内总涌水量减少;随初始隙宽的增大涌水量增加并当达到某一固定值时保持不变;隧洞涌水量在θ2/θ1=3.5,其中θ1=30°,即两组节理的夹角为75°处达到最大。     

基于SEM的可变形块体离散元法研究

针对边坡工程中岩土体连续-非连续渐进破坏的特点,提出一种新的变形体离散元方法（DEM）。与传统有限单元法（FEM）不同,弹簧元法（SEM）通过构建一组广义弹簧系统描述单元的力学行为。弹簧元法中的一个广义弹簧可以具有多个方向的刚度系数,确定广义弹簧系统的构造形式及其各刚度系数表达式是弹簧元法的核心。以三角形单元为例,介绍平面弹簧元的基本理论。对任何二维正交广义弹簧系统,通过定义广义弹簧变形与单元应变之间的关系,直接对比单元的应变能与弹簧系统的弹性势能即可得到广义弹簧刚度系数的表达形式。定义泊松刚度系数和纯剪刚度系数两个系统参数,描述正交广义弹簧之间的联系。对任意泊松比的材料,该方法都可准确地描述泊松效应的影响,计算结果与传统有限元法一致。该方法不需要求得有限元单元刚度矩阵的具体形式,具有直接方便、物理意义明确的优点,应用该方法给出任意4节点单元弹簧系统的构造形式及其各刚度系数的表达式。基于SEM的可变形块体离散元法,用弹簧元中的广义弹簧求解块体变形,用离散元中的接触弹簧计算块体间作用力,在单元节点的控制方程中实现弹簧元-离散元耦合计算,通过接触弹簧的状态实现材料由连续到非连续的破坏过程。在基于连续介质离散元法（CDEM）程序的基础上实现弹簧元-离散元耦合程序,应用耦合程序计算均质土坡在重力作用下的弹塑性变形和基覆边坡在重力作用下的破坏,初步证明该方法用于边坡变形渐进破坏分析的可行性。     

基于混合离散的砌石挡土墙边坡极限承载力下限分析

将极限分析下限法理论、混合数值离散思想和线性规划结合起来研究砌石挡土墙边坡的极限承载力。采用三角形有限单元离散土体来模拟土体的连续介质力学特性,构建土体静力许可应力场的约束条件,采用块体单元离散砌石体来模拟砌石体的非连续介质力学特性,构建砌石挡土墙的静力许可应力场的约束条件;同时建立有限元单元和块体单元交界面的约束条件;然后以超载系数为目标函数建立求解砌石挡土墙极限承载力的下限法线性规划模型,并使用内点算法进行最优化求解,获得边坡的极限荷载（或安全系数）和对应的应力场。通过3个算例的分析验证了所提方法的正确性。所提方法是将混合数值离散思想引入极限分析领域的一次尝试。     

基于SPH-DEM流-固耦合算法的滑坡涌浪模拟

采用光滑粒子流体动力学(SPH)方法和离散单元法(DEM),并基于达西渗透试验原理,提出一种SPH-DEM耦合算法用于处理宏观尺度下离散体与流体之间的相互作用问题,并基于该流-固耦合方法采用FORTRAN语言建立滑坡涌浪数学模型。模拟了块体滑坡问题,并与试验数据进行对比以验证该耦合方法的有效性。模拟分析了离散体滑坡产生涌浪以及涌浪传播的过程,分析比较了数值模拟得到的最大涌浪高度与经验公式结果间的关系。研究表明:采用该耦合方法模拟水下块体滑坡问题得到的结果与试验数据吻合良好;模拟离散体滑坡涌浪问题时,较清晰地反映了滑坡体对水的排挤到涌浪的产生过程、水在滑坡体中的渗透过程和滑坡体与水耦合作用发生变形的过程;不同公式分析方法得到的最大涌浪高度之间均有一定差别,算例中滑坡体受重力作用沿倾斜坡面滑动入水过程更接近于潘家铮法的垂直运动模式,因此,模拟结果与之最为接近。     

基于刚体平动运动单元的上限有限元研究

刚性块体极限分析上限法常应用于岩土工程稳定性研究,然而应用时需假定刚性块体破坏模式并递推繁琐的几何关系。为此,提出一种适应性更广的基于非线性规划模型的刚体平动运动单元上限有限元法,并解决了其优化模型初始值的确定问题。通过引入有限单元思想,将计算区域离散成刚体单元,同时以单元速度和节点坐标作为决策变量,由上限定理建立非线性规划模型获得上限解。利用编制的上限有限元程序进行边坡和浅埋隧道稳定性算例验证,表明运动单元上限有限元法能调整速度间断线至较优方位,所得破坏模式特征鲜明,上限解精度高,可广泛应用于边坡、隧道等稳定性分析研究。     

基于有限元法对深基坑围护桩水平位移的分析

利用有限单元法,通过有限元程序来进行模拟计算。将所研究问题假定为平面应变问题;岩土体为理想弹塑性模型,采用平面四节点等参单元,遵循适于岩土的Drucker—Prager屈服准则;桩与岩土体之间通过设置接触单元来模拟基坑围护桩桩与岩土体的相互作用,重点分析建筑物与基坑的不同距离、基坑围护桩的刚度变化和围护桩底端土体物理力学参数变化时基坑开挖对围护桩变形的影响。     

滑坡运动过程仿真分析

滑坡是一个动态过程，滑坡体的运动是一个集滑动，转动，拉张等运动方式的复杂运动过程，传统的极限平衡和计算和有限元分析均无法描述滑坡的运动学特点和运动过程。非连续变形分析（DDA）是最近发展起来的一种新的离散数值分析方法。该方法基于块体的运动学理论及数值分析，可以开展块体的静力和动力学计算。应用非连续变形分析方法对长江三峡区新滩滑坡的运动全过程进行了数值模拟研究，模拟方案充分依据该滑坡的地质，地形特征，按不同岩土体和地质结构面类型进行块体单元的划分，共划分成504个块体单元。模拟结果表明，新滩滑坡是以斜坡中部姜家坡一带的局部破坏为其运动的开始阶段，并进一步牵引上部滑体和推动下部滑体。代表性块体单元的位移变化曲线和滑动速度变化曲线反映了滑动过程中滑坡体块体系统的变形是非连续的，各处块体的动态形态各异，从而很好地再现了新滩滑坡的整个动态过程，揭示了滑坡的运动机制。     

基于离散单元法的滑坡堆积及其涌浪计算

采用离散单元法对水库库岸边坡的动态破坏过程进行了研究,给出了计算流程并开发了相应的软件。介绍了离散单元法的基本原理,分析了常用的计算滑坡速度方法的不足及用离散单元法的优势,基于离散单元法提供的块体运动信息并结合波动方程模拟了滑坡体滑入水库时所激起的涌浪。此外,采用离散单元法研究了滑坡的堆积形状。通过算例验证了方法的可行性与程序的正确性。算例结果表明：块体在坡面上运动时,用离散单元法计算得到的块体的速度和能量法与潘家铮法一致,且更方便、更高效;由于涌浪的叠加,多个块体所激起的涌浪高度比单个块体的涌浪高一些,这表明在进行涌浪分析时,考虑多块体的涌浪叠加效应是至关重要的,滑坡的堆积形状为滑坡的灾害评估提供了依据。通过自行开发的离散单元法软件模拟了块体的滑动和涌浪的产生及其传播过程,其模拟结果为研究边坡的滑动破坏的启动、破坏过程及滑坡的风险评价提供了新的思路。     

基于PFC~(2D)的贵阳西二环滑坡破坏过程模拟

正近年来,随着计算机技术的发展,数值模拟方法在滑坡领域得到广泛应用。目前的数值模拟方法中,主要包括确定性分析方法和非确定性分析方法两类,而确定性分析方法又可分为连续介质分析与非连续介质分析方法。其中,连续介质数值分析方法有有限单元法、边界元法、有限差分法等,非连续介质分析方法有块体离散元法、颗粒离散元法、关键块体理论、不连续变形分析(DDA法)等[1]-[2]。颗粒流法(PFC)作为离散单元法的一种,是基于离散元的方法来模拟圆盘或球颗粒介质的运动及其相互作用,可以真实地模拟滑坡的失稳破坏过程。目前已广泛应用于滑坡破坏运动分析之中[3]-[5]。     

模拟岩石破裂过程的块体单元离散弹簧模型

在变形体离散元的基础上建立块体单元离散弹簧模型,并应用于岩石破裂过程的数值模拟研究。该模型以连续介质力学理论为基础,将块体单元离散为具有明确物理意义的弹簧系统,通过对弹簧系统的能量泛函求变分获得各弹簧的刚度系数,进而可以直接利用弹簧刚度求解单元的变形和应力,提高计算效率。以重力作用下的岩质边坡计算为例,通过与传统的有限元进行对比,验证该模型弹性计算结果的正确性。在该基础上,引入Mohr-Coulomb与最大拉应力的复合破坏准则,判断单元的破坏状态及破裂方向。当单元的内部破坏面确定后,则通过块体切割的方式实现单元破坏,并建立单元边界和单元内部的双重破裂机制,实现块体由连续到非连续的破裂过程,进而显示的模拟裂纹的形成和扩展。最后,以巴西圆盘劈裂、单轴压缩破裂以及三点弯曲梁等典型算例验证该方法,结果表明该方法可以较好地模拟拉伸、压剪等应力状态下裂纹的形成和扩展,从而可模拟岩石介质由连续到非连续的破裂过程。     

偏心荷载作用下软弱下卧层承载力验算方法研究

以均布荷载作用下软弱下卧层竖向附加应力的简化计算方法为基础,导出了偏心荷载作用下软弱下卧层竖向附加应力的简化计算公式,并提出偏心荷载作用下软弱下卧层承载力验算的基本原则。     

高州大井岩组变质重结晶作用和变形作用的关系

高州大井岩组泥质岩中变斑晶种类繁多,变质重结晶作用分成四个阶段：Ｍ１发生于变形前埋藏过程中,Ｍ２、Ｍ３属进变质阶段,Ｍ４为退变阶段。变斑晶的生长主要集中在Ｄ１幕后Ｄ２幕前的一个相对“静态”的时期内,变形与变质重结晶作用存在一定的时间间隙。     

承受负摩擦力的桩基沉降计算的迭代法

结合负摩擦力产生的机理,认为软土地区的桩基的负摩擦力问题主要是沉降问题,而不是强度问题。利用半理论半经验的桩基沉降计算方法［１］ ,对承受负摩擦力的桩基的中性点及桩基的沉降、桩的分担比进行简单的迭代计算,由此可确定承受负摩擦力的桩基的沉降,计算结果与实测结果的对比表明,该法简单可行。     

祁连山地区1310年以来湿润指数及其年际变幅的变化与突变分析

本文通过用树木年轮资料重建的祁连山地区5～7月份1310年以来的湿润指数序列，建立了一个反映该地区湿润指数年际变幅序列，对这两个序列分别进行了等级分类和干湿、强弱的时段分析，并用最大熵谱分析法对这两个序列的不同时段分别进行了周期分析。采用HK突变检验方法对这两个序列分别进行了突变分析，发现祁连山地区的湿润指数及其年际变幅存在明显的突变年份。     

The Vendian succession of northeastern Spitsbergen: Petrogenesis of a dolomite-tillite association

The sedimentary record of late Precambrian time is magnificently displayed in the highland snowfields of northeastern Spitsbergen (Svalbard). Vendian strata are represented essentially by the Polarisbreen Group which consists mostly of dolostone and includes two dolomitic glacial units. The oldest sediments in the Polarisbreen Group compose the Elbobreen Formation (c. 400 m), which is divided into four laterally-persistent members. The Lower Carbonate Member (E1, 125 m) contains a distinctive basal dark-grey limestone (with microspar-filled synaeresis cracks) suggested to be of lagoonal origin and associated with minor dolostone, shale and chert. Higher parts of the member are dominantly dolostone, partly stromatolitic, with some shale and sandstone; shallow subtidal to intertidal deposition is indicated by the dominance of intraclastic lithologies and relics of anhydrite. Penecontemporaneous dolomite is partially overprinted by microsparry dolomite, thought to be of groundwater origin.The redefined Petrovbreen Member (E2) consists of diamictite and other detrital dolostone. Pronounced thickness variations (2–40 m) are thought to be original depositional features. The member represents the deposits of a short glacial period in which the following depositional processes are inferred: lodgement (massive diamictite), subaqueous meltout (massive and bedded diamictite), ice-rafting (lithologies bearing dropstones, and possibly also diamictite), redeposition by sediment gravity flows (some diamictite and conglomerate; rhythmite and shale), current winnowing (thin tabular conglomerate), subaerial or subaqueous meltwater action (channelled conglomerate and sandstone), periglacial shrinkage (diamictite wedge-fillings).The MacDonaldryggen Member (E3, 230 m) is a monotonous succession of shaly dolostone of lagoonal origin. It grades up into the Slangen Member (E4, 25 m) which consists of subtidal to intertidal dolarenite with anhydrite relics succeeded by fenestral dolostone that was fractured and cemented by saline groundwaters in an emergent environment.The Wilsonbreen Formation (160 m) represents a return to glacial deposition, but this time longer-lasting and with substantial extra-basinal material represented. The Gropbreen Member (W1, 28–107 m) and the Ormen Member (W3, 44–139 m) consist dominantly of dolomitic diamictite with subordinate conglomerate and sandstone and are separated by a Middle Carbonate Member (W2, 3–30 m) which contains distinctive rhythmitic and stromatolitic limestone as well as sandstone. The same depositional processes can be recognised as in the Petrovbreen Member, but the Wilsonbreen Formation is overall of somewhat more continental aspect (lower proportion of rhythmite and dropstone structures). In addition there are: basal breccia and crack-fillings formed by frost-shattering of the underlying cemented dolostone, tabular sandstone thought to be formed by wave reworking of outwash, a striated (terrestrial) cobble pavement, supraglacially-derived breccia horizons, periglacial wedges filled by sand and the W2 assemblage of possible lacustrine origin.The Dracoisen Formation (525 m) represents an abrupt return to non-glacial conditions. An upward-deepening wave-dominated succession of pure dolostone (D1, 20 m) and impure dolostone (D2, 105 m) is succeeded by offshore black shale (D3, 150 m) and then by a very-shallow water succession of evaporite lacustrine aspect with a dolostone containing evaporite relics (D5, 10 m) separating dolomitic sandstone and shale (D4, 80 m and D6, 150 m). The contact with the transgressive Cambrian sandstones above is poorly exposed.Deposition of the succession dominantly under marine conditions is inferred, but it is difficult to rule out a lacustrine environment at any particular horizon. This dolomite—tillite association can be explained by penecontemporaneous (and minor secondary) dolomite formation in marginal environments (with warm climatic indicators at some levels) being sharply interrupted, because of rapid climatic changes, by glacial sediments containing abundant detrital dolomite. Since the latter sediments make up only 17% of the 1080m-thick succession, glacial conditions only occupied a small proportion of Vendian time.     

Garnet-forming reactions and recrystallization in high-grade mylonite zones, MacRobertson Land, east Antarctica

Proterozoic granulite facies gneisses in MacRobertson Land, east Antarctica, are cut by numerous D5 mylonite-ultramylonite zones of probable Cambrian age. In garnet-absent mafic two-pyroxene gneisses and garnet-bearing charnockitic orthogneisses, the mylonite-ultramylonite zones are characterized by the growth of garnet at the expense of ilmenite, pyroxene and plagioclase. Textures within each mylonite zone can vary from protomylonitic to ultramylonitic. A range of mineral textures involving M5 garnet is developed corresponding to variations in deformation intensity. In protomylonites, garnet occurs as coronas on orthopyroxene-plagioclase and ilmenite-plagioclase boundaries, and as overgrowths on earlier garnet. In ultramylonites, fine-grained orthopyroxene-plagioclase-garnet ± quartz ± clinopyroxene intergrowths and poikilitic garnet are common. Garnet growth in all shear zones is accompanied by shifts in the compositions of neoblastic minerals occurring with garnet, consistent with local chemical equilibrium having been attained during recrystallization. Mylonitization is inferred to have occurred at P ∼ 6.5 kbar. Temperature estimates for M5 vary between 550 and 797 C, which may reflect variations and uncertainties associated with the calibrations used and/or partial re-equilibration during cooling. The presence of post-tectonic, coronate garnet in some mylonite zones indicates that garnet continued to form exclusively in the mylonite zones after movement had ceased and is interpreted to reflect the effects of localized strain heating.     

Ridgelet变换在地震数据处理中的应用

最近，Stanford大学E J Candes和D L Donoho教授提出了一种新的多尺度变换-Ridgelet变换，它特别适合于具有直线或超平面奇性的高维信号的描述，而且具有较高的逼近精度。将它应用于地震处理中，即可获得较以前其它方法无法达到的精度和效果。     

Book reviews

Book review in this article
Cross-bedding, Bedforms and Paleocurrents, by D. M. Rubin
Beach and Nearshore Sediments and Processes, ed. By R. A. Davis Jr.,
Sedimentary Processes on the Amazon Continental Shelf, ed. by C. A. Nittrouer and D. J. DeMaster
Fjords: Processes and Products, by J. P. M. Syvitski, D. C. Burrell, and J. M. Skei
Alluvial Soils, ed. by J. Gerrard,
Electron Micrographs (TEM, SEM) of Clays and Clay Minerals, by K.-H. Henning and M. Störr
The Origins of Angiosperms and their Biological consequences, ed. by E. M. Friis, W. G. Chaloner & P. R. Crane
The Motion of Allochthonous Terranes Across the North Pacific Basin, by M. G. Debiche, A. Cox and D. Engebretson
Geomorphic Systems of North America, ed. by W. L. Graf
Approaches to Interpretation of Sedimentary Environments, ed. by D. J. Cant and F. J. Hein
Carbonate Depositional Environments: Modern and Ancient. Part I. Reefs.
European Dinantian Environments, ed. by J. Miller, A.E. Adams, and V.P. Wright
Mesozoic and Cenozoic Oceans, ed. by K. J. Hsü
Marine Minerals: Advances in Research and Resource Assessment, ed. by P.G. Teleki, M.R. Dobson, J.R. Moore and U. von Stackelberg.
Sedimentation and Mineral Deposits in the Southwestern Pacific Ocean, ed. by D.S. Cronan
A Practical Approach to Sedimentology, by Ray Lindholm
Clastic Particles, Scanning Electron Microscopy and Shape Analysis of Sedimentary and Volcanic Claste, ed. by John R. Marshall     

Evolution of fluid pathways of Charroux-Civray tonalite (part II): Numerical study of microcrack networks

Transport properties in two tonalite rock samples are studied using image analysis and S.E.M. data. Microcracks are considered to represent the most important small scale fluid vector of the rock. These microcracks can be caused by thermal or mechanical stresses history, or decompression. The proposed model of microporosity is geometrically constrained by the spatial distribution of mineral grains. It was obtained at the grain scale. Image analysis of tonalite half core provides a 2D representation of primary mineral species distribution, which defines the grain-boundary network. S.E.M. stereological counting provides the main features of intergranular microcracks at the scale of the pores: probabilities of opening and mean apertures. These data allow one to reconstruct, via Monte-Carlo simulations, a model of two open microvoids type in 2D space, according to the probability of opening of each microcrack types. The proposed model is constructed with two crack family, with aperture lower than 0.1 μm and with higher aperture. Geometry of the percolation backbone is further determined from these simulations: tortuosity of the percolation backbone is equal to (1.48)² and (1.6)² respectively for unaltered and fairly altered tonalite facies. In terms of path types, composition of the percolation backbone shows also strong differences according to the facies.     

Book reviews

Book review in this article
Iron-Formation: Facts and Problems, by A. F. Trendall and R. C. Morris.
Palaeomagnetism—Principles and Applications in Geology, Geophysics and Archaeology, by D. H. Tarling.
Dolomites of the Monterey Formation and other Organic-rich Units, ed. by R. E. Garrison, M. Kastner and D. H. Zenger.
Lake Sedimentology, by L. Håkanson and M. Jansson.