深部隧道围岩分区破裂颗粒流模拟研究

深部岩体处于"三高"环境下,表现出不同于浅埋岩体的性质,其变形破裂规律更为复杂(分区破裂、片帮、塑性流动、岩爆等),为了准确描绘深部隧道围岩变形破裂规律,采用PFC从微观角度研究深部岩体的宏观响应。研究发现:随着隧道埋深增加,压力增大,由浅部围岩表面塑性破坏变为深部围岩破裂扩展,破裂区域呈交替分区破裂向深部发展,破裂区的间距与岩性和埋深有关;从横断面看,拱腰先出现破裂,然后拱脚出现破裂,最后贯通形成破裂区;若围岩表面施加外力,破裂区域减小,分区向深部移动,因此预应力锚杆有效地改善了围岩承载特性;研究结果与模型试验吻合,符合深部岩体卸荷作用下的变形破坏规律。结论可为深部岩体工程设计施工提供参考。     

Dam-reservoir-massed foundation system and travelling wave along reservoir bottom

In the present paper, effects of non-uniform excitation due to spatially variation of seismic waves under the reservoir bottom on linear and nonlinear responses of arch dams are studied. Foundation is assumed to be massed and infinite elements are utilized to model semi-infinite medium via the far-end boundary of the foundation FE model. A continuum crack propagation model based on the smeared crack approach in 3D space is introduced. Reservoir's water is assumed compressible and the coupled system is solved using the staggered method. As a case study, Amir-Kabir double curvature arch dam in Iran is selected to investigate seismic behavior of the system. Two cases are analyzed in which wave travelling velocities are taken as 650 m/s and infinity. According to the results, non-uniform excitation leads to reduction in seismic response of the system and, in addition, frequency content of the response is different completely in comparison with the system under the uniform excitation.     

Investigating landslide‐related cracks along the edge of two loess platforms in northwest China

Cracks are widely developed along the edge of loess platforms in northwest China. Field surveys reveal that these cracks can be grouped into shallow and deeply penetrating ones. The former occur at a small distance from the platform edge, normally penetrate into the top unsaturated loess with the penetration depth being controlled by the joints in loess. The latter penetrate deeper into the saturated loess farther away from the platform edge. These cracks control the inflow and drainage of irrigation water. The shallow penetrating crack can fail as a slide or fall with a volume of up to hundreds of cubic meters. The deeply penetrating crack can fail as a flow‐like landslide with a volume of thousands of cubic meters or more. A full‐scale field test simulating irrigation on the platform surface was conducted. The two types of crack can be interconnected so that the water applied in the test finally flowed into the deep crack and was discharged from the platform. Analysis of soil stress states and the results of field test show that the deeply‐penetrating cracks could have both positive as well as negative effects on slope stability. On the one hand, water can flow more freely in the cracks, and the loess could be saturated and trigger a landslide. On the other hand, the water can drain more easily along the crack and slope stability could be enhanced as the groundwater level is suppressed. Copyright © 2012 John Wiley & Sons, Ltd.     

Nonlinear seismic response analysis of arch dam-foundation systems- part I dam-foundation rock interaction

The arch dam–foundation rock dynamic interaction and the nonlinear opening and closing effects of contact joints on arch dam are important to the seismic response analysis of arch dams. Up to date, there is not yet a reasonable and rigorous procedure including the two factors in seismic response analysis. The methods for the analysis of arch dam–foundation rock dynamic interaction in frequency domain are not suitable to the problem with nonlinear behaviors, in this paper, so an analysis method in time domain is proposed by combining the explicit finite element method and the transmitting boundary, and the dynamic relaxation technique is adopted to obtain the initial static response for dynamic analysis. Moreover, the influence of arch dam–foundation dynamic interaction with energy dispersion on seismic response of designed Xiaowan arch dam in China is studied by comparing the results of the proposed method and the conventional method with the massless foundation, and the local material nonlinear and nonhomogeneous behaviors of foundation rock are also considered. The reservoir water effect is assumed as Westergaard added mass model in calculation. The influence of the closing–opening effects of contact joints of arch dam on the seismic response will be studied in another paper.     

Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.     

Seismic response analysis of arch dam-water-rock foundation systems

The effect of water compressibility on the seismic responses of arch dams is not well understood. In this paper, a numerical model is developed with rigorous representation of the dynamic interaction between arch dam-water-rock foundation. The model is applied to the seismic response a nalysis of an arch dam with a height of 292m designed to aseismic intensity of IX. It is shown that consideration of the water compressibility clearly decreases the stress responses at key positions of the dam, while the added mass model gives a conservative estimate.     

Earthquake response of concrete arch dams: a plastic–damage approach

There are several alternatives to evaluate seismic damage‐cracking behavior of concrete arch dams, among which damage theory is the most popular. A more recent option introduced for this purpose is plastic–damage (PD) approach. In this study, a special finite element program coded in 3‐D space is developed on the basis of a well‐established PD model successfully applied to gravity dams in 2‐D plane stress state. The model originally proposed by Lee and Fenves in 1998 relies on isotropic damaged elasticity in combination with isotropic tensile and compressive plasticity to capture inelastic behaviors of concrete in cyclic or dynamic loadings. The present implementation is based on the rate‐dependent version of the model, including large crack opening/closing possibilities. Moreover, with utilizing the Hilber–Hughes–Taylor time integration scheme, an incremental–iterative solution strategy is detailed for the coupled dam–reservoir equations while the damage–dependent damping stress is included. The program is initially validated, and then, it is employed for the main analyses of the Koyna gravity dam in a 3‐D modeling as well as a typical concrete arch dam. The former is a major verification for the further examination on the arch dam. The application of the PD model to an arch dam is more challenging because the governing stress condition is multiaxial, causing shear damage to become more important than uniaxial states dominated in gravity dams. In fact, the softening and strength loss in compression for the damaged regions under multiaxial cyclic loadings affect its seismic safety. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of Injected Water on Hydraulic Fracturing Deduced from Acoustic Emission Monitoring

—In order to investigate the effects of injected water in hydraulic fracturing, experiments were conducted on cubic granite specimens, comparing fracturings induced by conventional water injection with those induced by pressurization of a urethane sleeve, thereby realizing "hydraulic fracturing" without the use of fracturing fluid. In both experiments, a shear type mechanism was found to be dominant in fault plane solutions of AE events. However, in the case of water injection, cracks extended rapidly with large drops in hole water pressure and bursts of AE, whereas in pressurization by the urethane sleeve, cracks extended stepwise with no such large drops in hole pressure and no bursts of AE. The difference in crack extension in the two experiments can be analyzed by comparing relations between crack length and stress intensity factor of mode I at a crack tip. The observation and analysis indicate that existence of fracturing fluid like water helps initiated cracks to extend rapidly and widely in hydraulic fracturing in actual HDR fields. Received September 12, 1996, accepted January 24, 1997     

库水模型及库水深度对拱坝动力响应的影响

基于接触非线性有限元模型,以锦屏一级拱坝为例,库水分别采用附加质量模型、可压缩流体有限元模型、不可压缩流体有限元模型计算了正常蓄水位及运行低水位时坝体的动力响应,结果表明:库水模型对拱坝动力响应有较大影响,随库水深度的增大,各模型计算结果差异增大;相比于流体可压缩模型,采用不可压缩流体模型所得动力响应普遍偏大;运行低水位工况,由于静水压力减小导致拱效应减弱,从而降低了拱坝的整体性,因此运行低水位工况各缝开度普遍高于正常蓄水位工况,且其拉应力范围较大,因此,运行低水位工况将对抗震设计起控制作用。     

Experimental study of dynamic damage of an arch dam

Seismic responses of a 292‐m high arch dam were studied by experiment on a shaking table. The model system included the arch dam with contraction joints, a part of a reservoir, and a partial foundation with a topographic feature near the dam. Potential rock wedges on the abutments and the mechanical properties including uplift on the kinematic planes were carefully simulated. A damping boundary consisting of a viscous liquid was introduced to simulate the effect of dynamic energy emission to far field, which made the dynamic interaction between the dam and the foundation be adequately represented in the model test of an arch dam system. Dynamic responses of the arch dam system under a sequence of seismic loadings in increasing strength were examined. Eleven cracks or overstresses on the model dam due to the earthquake excitations were observed, and consequently, its natural frequency dropped by about 14%, but the model dam was stable under the hydrostatic load of the impounded water after the test. Copyright © 2006 John Wiley & Sons, Ltd.     

Distinguishing oil and water layers by interpreting acoustic logging data with changing well diameters

During surveys, water layers may interfere with the detection of oil layers. In order to distinguish between oil and water layers, research on the properties of well diameters and oil and water layers and their relation to acoustic logging rules is essential. Using Hudson's crack theory, we simulated oil and water layers with different well diameters or crack parameters (angle and number density). We found that when the well radius increases from 0.03 m to 0.05 m, the variation ratio of compressional wave amplitude for the oil layer is less than that for the water layer. The difference of Stoneley wave amplitude between the crack parameters (angle and number density) is greater in the case of the water layer than in the case of the oil layer. The response sensitivity of wave energy is greater for the water layer than that for the oil layer. When the well radius increases from 0.05 m to 0.14 m, the maximum excitation intensity for oil layer is greater than that for the water layer. We conclude that the propagation of an elastic wave is affected by medium composition and well diameter, and the influence has certain regularity. These results can guide further reservoir logging field exploration work.     

一种缝洞型碳酸盐岩储层胶结指数m计算新方法

碳酸盐岩储集空间类型多样、孔隙结构复杂,具有很强的非均质性,使得传统的阿尔奇公式应用效果不佳,主要原因之一在于胶结指数m的无法准确确定.基于多孔介质理论,在深入分析不同孔隙空间对储层导电贡献外,进一步考虑了孔隙形态对导电性能的影响,尤其是裂缝倾角,进而给出了全新的计算胶结指数m的方法,并深入剖析组成碳酸盐岩复杂孔隙空间的各部分对胶结指数m的影响.认为对于物性较差的缝洞型储层,裂缝倾角对胶结指数m的影响很大,当裂缝角度较小时,裂缝会使胶结指数m变小,而当裂缝角度较大时,裂缝反而会使胶结指数m变大;当储层物性很好时,裂缝倾角的影响则可以忽略.实际应用表明,利用基于裂缝倾角的多孔介质模型得到的胶结指数m所计算的含水饱和度精度相比传统方法有较大提高,且该方法适用于任何类型的储层,具有很好的应用性.     

有缝拱坝-地基系统非线性地震波动反应分析方法

拱坝坝体中缝界面在地震作用下的开、合模拟对拱坝地震反应有重要影响。目前在拱坝地震反应分析中应用的一些模型,如Ｆｅｎｖｅｓ模型、Ｄｏｗｌｉｎｇ的三参数模型,存在对缝界面的接触状态模拟精度不高的问题,这对于准确评价工程上关心的缝界面最大张开率是不足的。此外,坝基中存在的断层、节理、裂隙待软弱夹层的非线性力学特性和无限地基能量辐射作用对拱坝地震反应也有重要影响。在拱坝地震反应分析中同时考虑这三种影响因素     

A comparative study of the different procedures for seismic cracking analysis of concrete dams

Different procedures are compared for the three-dimensional seismic cracking analysis of gravity and arch dams during strong earthquakes. The fracture procedures include the extended finite element method with cohesive constitutive relations, crack band finite element method with plastic-damage relations, and the finite element Drucker−Prager elasto-plastic model. These procedures are used to analyze the nonlinear dynamic response of Koyna dam to the 1967 Koyna earthquake and the seismic cracking of the Dagangshan arch dam subjected to design earthquake. The cracking process and profiles of the two dams using the three different procedures are compared. The applicability and the suitability of the three procedures for seismic cracking analysis of gravity and arch dams are discussed.     

Seismic stability assessment of an arch dam-foundation system

A seismic stability assessment of arch dam-foundation systems is presented using a comprehensive approach,in which the main factors that significantly influence the seismic response of an arch dam-foundation system are considered.A large scale finite element model with over 1 million degrees of freedom is constructed for the Baihetan arch dam(289 m high),which is under construction in the Southwest of China.In particular,the complicated geological conditions with faults intersecting interlayer shear weakness zones at the dam base and the dam abutment resisting force body is modeled in the analysis.Three performance indices are adopted to assess the seismic stability of the arch dam.The results demonstrate that the opening of the joints of the Baihetan arch dam is small and the water stop installed between the joints would not be torn during a design earthquake.The yielding formed in the interface between the dam and foundation does not reach the grouting curtain that would remain in an elastic state after an earthquake.The yielding zones occurring on the upper portion of the dam faces extend 1/8 thickness of block section into the dam body and thus cantilever blocks need not be concerned with sliding stability.The faults and interlayer shear weakness zones in the near field foundation exhibit severe yielding,and a potential sliding surface is penetrated.Although the factor of safety against sliding of the surface fluctuates with a decreased trend during an earthquake,the minimum instantaneous value reaches 1.02 and is still larger than 1.0.Therefore,a conclusion is drawn that the Baihetan arch dam-foundation system will remain stable under the design earthquake.     

大跨度拱桥地震动输入模式研究

地震动输入是大跨度桥梁地震反应分析的重要一环。从明确大跨度拱桥的临界跨度入手,探讨了大跨度拱桥地震动输入模式中的行波效应、三向输入及其地震动选取问题。基于某大跨度钢管混凝土拱桥,建立了有限元分析模型,考查了行波效应及三向地震动输入对拱关键截面内力、减隔震支座位移及粘滞阻尼器冲程的影响。结果表明:给出的大跨度拱桥临界跨度确定方法合理,行波效应对拱顶轴力有重要影响,不考虑三向地震动同时作用会明显低估大跨度拱桥的地震反应。确定大跨度拱桥地震动输入模式时需考虑行波效应与三向地震动同时输入。     

Earthquake response of arch dams to spatially varying ground motion

The response of two arch dams to spatially varying ground motions recorded during earthquakes is computed by a recently developed linear analysis procedure, which includes dam–water–foundation rock interaction effects and recognizes the semi‐unbounded extent of the rock and impounded water domains. By comparing the computed and recorded responses, several issues that arise in analysis of arch dams are investigated. It is also demonstrated that spatial variations in ground motion, typically ignored in engineering practice, can have profound influence on the earthquake‐induced stresses in the dam. This influence obviously depends on the degree to which ground motion varies spatially along the dam–rock interface. Thus, for the same dam, this influence could differ from one earthquake to the next, depending on the epicenter location and the focal depth of the earthquake relative to the dam site. Copyright © 2009 John Wiley & Sons, Ltd.     

Crack healing in rocksalt via diffusion in adsorbed aqueous films: Microphysical modelling versus experiments

Microcracks within the excavation damaged or disturbed zone (EDZ) in a salt-based radioactive waste repository (or an energy storage facility) can heal/seal by mechanical closure driven by compaction creep, by surface-energy-driven processes like diffusive mass transfer, and by recrystallization. It follows that permeability evolution in the excavation damaged zone around a backfilled or plugged cavity will in the short term be dominated by mechanical closure of the cracks, while in the longer term diffusive mass transfer effects are expected to become more important. This paper describes a contribution to assessing the integrity of radioactive waste repositories sited in rocksalt formations by developing a microphysical model for single crack healing in rocksalt. More specifically, single crack healing models for cracks containing a thin adsorbed water film are developed. These microphysical models are compared with single crack healing experiments, which conclusively demonstrate diffusion controlled healing. Calibration of unknown model parameters, related to crack surface diffusivity, against the experimental data enable crack healing rates under repository conditions to be estimated. The results show that after the stress re-equilibration that follows repository sealing, crack disconnection can be expected on a timescale of a few years at laboratory humidity levels. However, much longer times are needed under very dry conditions where adsorbed aqueous films are very thin.     

近断层地震动对上承式拱桥动力响应的影响

为了揭示近断层地震作用下上承式钢筋混凝土拱桥的动力响应特点,以西南山区某上承式拱桥为背景,用OpenSEES平台建立了全桥非线性动力分析模型,探讨了近断层地震动的输入方式、脉冲效应和竖向地震动等三个关键因素对桥梁动力响应的影响规律。研究结果表明:地震输入方式对拱圈地震响应的影响较小,但对拱上立柱地震响应的影响很大,尤其是拱顶附近的短立柱,在抗震分析中,建议偏安全地采用三向地震输入方式;脉冲效应对拱桥地震响应的影响非常大,会导致拱圈、拱上立柱和桥面板地震响应大幅增加,桥面板残余平面转角甚至增大6倍以上;竖向地震动对拱圈轴力和面内弯矩、拱上立柱纵向弯矩和剪力的影响很大,拱顶处的面内弯矩放大倍数最大可达2.95,总体来说,采用规范所建议的方法考虑竖向地震是偏保守的。     

Quantifying Damage, Saturation and Anisotropy in Cracked Rocks by Inverting Elastic Wave Velocities

Crack damage results in a decrease of elastic wave velocities and in the development of anisotropy. Using non-interactive crack effective medium theory as a fundamental tool, we calculate dry and wet elastic properties of cracked rocks in terms of a crack density tensor, average crack aspect ratio and mean crack fabric orientation from the solid grains and fluid elastic properties. Using this same tool, we show that both the anisotropy and shear-wave splitting of elastic waves can be derived. Two simple crack distributions are considered for which the predicted anisotropy depends strongly on the saturation, reaching up to 60% in the dry case. Comparison with experimental data on two granites, a basalt and a marble, shows that the range of validity of the non-interactive effective medium theory model extends to a total crack density of approximately 0.5, considering symmetries up to orthorhombic. In the isotropic case, Kachanov's (1994) non-interactive effective medium model was used in order to invert elastic wave velocities and infer both crack density and aspect ratio evolutions. Inversions are stable and give coherent results in terms of crack density and aperture evolution. Crack density variations can be interpreted in terms of crack growth and/or changes of the crack surface contact areas as cracks are being closed or opened respectively. More importantly, the recovered evolution of aspect ratio shows an exponentially decreasing aspect ratio (and therefore aperture) with pressure, which has broader geophysical implications, in particular on fluid flow. The recovered evolution of aspect ratio is also consistent with current mechanical theories of crack closure. In the anisotropic cases—both transverse isotropic and orthorhombic symmetries were considered—anisotropy and saturation patterns were well reproduced by the modelling, and mean crack fabric orientations we recovered are consistent with in situ geophysical imaging. Our results point out that: (1) It is possible to predict damage, anisotropy and saturation in terms of a crack density tensor and mean crack aspect ratio and orientation; (2) using well constrained wave velocity data, it is possible to extrapolate the contemporaneous evolution of crack density, anisotropy and saturation using wave velocity inversion as a tool; 3) using such an inversion tool opens the door in linking elastic properties, variations to permeability.