Flow Simulations To Evaluate Upscaling of Permeability

We study upscaling of the permeability for porous media flow on a grid with one million blocks. The purpose is to illustrate how flow simulations can be used to evaluate upscaling methods.     

棒式光纤传感器在地下洞群模型试验中的应用

地质力学模型试验是在复杂条件下进行大型地下洞室围岩稳定性分析的重要研究手段,能够真实地反映地质构造和工程结构的空间关系,准确模拟施工过程对厂房稳定性的影响以及岩体的破坏特征。双江口水电站地下洞群模型试验中采用预留孔径灌胶方式进行高精度的棒式光纤光栅（FBG）位移传感器的埋设,以便对洞室周边的位移进行量测,同时应用FLAC3D进行数值模拟比较。结果表明：棒式光纤光栅传感器量测方法解决了模型试验中内部变形小、量测精度要求高的问题,与数值模拟结果相吻合;试验中传感器埋设的方法是成功的,无论是开挖阶段还是超载阶段传感器均能有效地监测模型内部变形,揭示地下洞室群模型在开挖阶段和超载阶段的变形机制和破坏特征。     

Robust numerical implementation of a 3D rate-dependent model for reservoir geomechanical simulations

A 3D elasto-plastic rate-dependent model for rock mechanics is formulated and implemented into a Finite Element (FE) numerical code. The model is based on the approach proposed by Vermeer and Neher (A soft soil model that accounts for creep. In: Proceedings of the International Symposium “Beyond 2000 in Computational Geotechnics,” pages 249-261, 1999). An original strain-driven algorithm with an Inexact Newton iterative scheme is used to compute the state variables for a given strain increment.The model is validated against laboratory measurements, checked on a simplified test case, and used to simulate land subsidence due to groundwater and hydrocarbon production. The numerical results prove computationally effective and robust, thus allowing for the use of the model on real complex geological settings.     

层状岩体三维可视化构模与数值模拟的集成研究

由于数据结构的差异,用于岩土工程三维可视化的网格与用于数值模拟的计算网格之间存在本质差异,如何通过两种网格的共享实现三维地质模型与数值模拟之间的集成是岩土工程建模的重要内容。在对层状岩体三维可视化网格与数值模拟网格的特点进行剖析的基础上,提出了基于松散模式的三维规则格网与FLAC3D基本元素之间的转化方法,从理论上实现了“可视”与“可算”的结合。将其应用于我国某海底隧道的选线研究中,证明了该方法的合理性和有效性。     

洞群模型试验中微型多点位移计的设计及应用

地质力学模型试验是研究大型地下洞室群围岩稳定性状况的重要方法,但是,许多制约因素使得这类试验难以达到预期效果,其中包括试验过程中物理量的量测问题。以双江口水电站地下洞群稳定性为研究对象,开展了模型试验为准确量测三维应力状态下洞群模型试验中围岩内部的微小位移,设计、制作和安装了高精度微型多点位移计,获得了围岩在开挖全过程中的位移变化规律曲线,同时,还开展了与模型试验相对应的数值模拟。将试验结果与数值模拟结果进行了对比分析。研究结果表明,设计成功的微型多点位移计能够广泛应用在各类岩土工程模型试验之中。     

盐岩油气储库介质地质力学模型相似材料的研制

以江苏金坛盐岩地下油气储库为示范工程,以精铁粉、重晶石粉、石英砂、一级松香和医用酒精为原料,根据相似原理以及金坛盐岩油气储库地层中盐岩、泥岩和盐岩夹层的实测物理力学参数,开展了各种配比条件下材料的单轴压缩、三轴压缩、直剪和巴西劈裂试验,通过原岩和模型材料的重度、泊松比、弹性模量、抗压强度、抗拉强度、黏聚力、内摩擦角等基本物理力学参数以及相应变形特征和强度破坏特性的对比分析,研制出基本满足相似条件要求的盐岩油气储库介质（盐岩、泥岩和盐岩夹层）的地质力学模型相似材料,为开展盐岩地下油气储库运营的变形与失效破坏地质力学模型试验提供了可靠的材料保证。     

CO2地质封存力学问题的数值模拟方法综述

随着社会的发展,人类对化石能源的依赖导致大量的CO2排入大气层,从而引起全球变暖、海平面上升等一系列全球性气候问题。CO2地质封存是当前CO2减排最有前景的技术,但大量的CO2注入地层易诱发相应的力学问题：地表隆起变形、盖层完整性受损、断层活化等。因此,为了降低CO2地质封存带来的安全风险,理论分析、数值模拟和响应面方法等手段被应用于此类力学问题的分析中。由于数值模拟方法能够解决大尺度范围内复杂几何模型的多场耦合问题,数值模拟成为当前在CO2地质封存力学领域中应用最广泛的方法。因此,对CO2地质封存若干力学问题的数值模拟方法进行了全面的综述。首先,简要介绍了多孔介质的温度-渗流-力学-化学（THMC）多场耦合原理,并对数值模拟解决多场耦合问题的方法进行了归类。然后,详细总结了数值模拟在解决CO2地质封存力学问题方面的国内研究进展。最后,讨论了数值模拟方法在此类力学问题方面的应用缺陷,并提出了若干建议。     

Use of Border Regions for Improved Permeability Upscaling

A procedure for the improved calculation of upscaled grid block permeability tensors on Cartesian grids is described and applied. The method entails the use of a border region of fine-scale cells surrounding the coarse block for which the upscaled permeability is to be computed. The implementation allows for the use of full-tensor permeability fields on the fine and coarse scales. Either periodic or pressure–no flow boundary conditions are imposed over the extended local domain (target block plus border regions) though averaged quantities, used to compute the upscaled permeability tensor, are computed only over the target block region. Flow and transport results using this procedure are compared to those from standard methods for different types of geological and simulation models. Improvement using the new approach is consistently observed for the cases considered, though the degree of improvement varies for different models and flow quantities.     

Assessment of geomechanical characterization of a rock mass using a seismic geophysical technique

Summary In this article the use of the seismic properties of a rock mass to assess the geomechanical behaviour of the rock mass is examined and illustrated. The correlation between the values of seismic velocity and some basic parameters of rock mechanics, such as RQD and fracture frequency, have been established. Moreover, the relationships between the engineering rock-mass classification systems and seismic quality index (QI) have also been derived. The relationships have been derived from an investigation for a proposed dam in Wadi Mujib, Jordan. The ground conditions comprise colluvial and fluvial superficial deposits overlying limestones, mudstones, marls and shales. It is suggested that the relationships obtained in this study could be applied in areas that have similar lithological and structural characteristics.     

Thermodynamic properties of MgSiO3 perovskite derived from large scale molecular dynamics simulations

A molecular dynamics simulation study of MgSiO₃ has been performed using a large sample containing 4096 unit cells. Thermodynamic properties have been extracted using a semiclassical approximation to the correct quantum mechanical treatment, using the calculated density of states and the quantum harmonic formalism for thermodynamic functions. Simulations performed at different temperatures and volumes have given an estimate of the relative contributions due to thermal expansion (quasi-harmonic effects) and direct anharmonic interactions. Comparison of results for mean square atomic displacements with results on smaller samples have shown the limitations of smaller sample sizes.     

Numerical Calculation of Equivalent Cell Permeability Tensors for General Quadrilateral Control Volumes

A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques.     

内蒙古高家梁煤矿巷道加固及围岩支护结构的数值模拟分析

高家梁煤矿位于内蒙古鄂尔多斯市东胜区东南8km处,矿井设计规模年产600×104吨。建井期间在埋深50～150m的运输大巷中部,遇到一种低强度软岩,巷道掘进与施工过程中冒顶、塌方时有发生,常规锚杆支护失效,对施工安全和进度造成很大影响。基于巷道地质情况与软岩特征,根据巷道围岩松动圈理论及围岩与支护作用关系,重点分析了软岩巷道的变形、顶板失稳机理,确定了运输大巷的工程加固方法。经过ADINA数值模拟验算以及后期的运输大巷变形监测曲线都证明加固方案满足运输大巷对支护刚度和强度的要求,同时还将加固方案进行了经济评价,得出了这种加固方案对于高家梁矿的特殊软岩巷道段加固在技术上可行,经济上可以接受的结论。     

裂隙岩体宏观力学参数的二维数值模拟

岩体通常含许多不同尺度的裂隙,因而力学性质非常复杂。由于现场测试及常规实验室试验常受各种条件的限制,因而其不能准确反映岩体的宏观物理力学性质。本文基于有限元分析软件对二维情况下的完整岩体和含不同倾角的软弱夹层的岩体在单轴和双轴压缩下进行了数值模拟,得到岩体变形的应力-应变曲线。该类曲线能反映裂隙对岩体力学参数的影响,对于指导工程实践有一定的意义,并且为获取岩体力学参数提供了新的途径。     

双江口工程地下洞室群围岩损伤及稳定性研究

对四川省双江口水电站工程的地下洞群的节理岩体围岩稳定性进行研究,分别采用了弹塑性模型、节理岩体断裂损伤模型进行数值分析。在弹塑性模型中,采用锚固等效法来模拟锚杆及锚索的支护效果。断裂损伤模型以损伤及断裂力学原理为基础,运用损伤本构模型和损伤演化方程在FLAC3D软件的框架下开发出子程序,考虑了岩体中节理裂隙对洞室围岩稳定性的影响,进而得到用上述的模型和计算程序分析出的围岩塑性区、损伤演化区和位移规律。将两种模型的计算结果与相应工程模型试验的结果进行了比较分析。研究结果表明,损伤模型更接近工程实际情况,与模型试验结果更为接近     

锦屏电站坝厂区卸荷岩体力学参数研究

选择具有代表性的区块,根据结构面的性状和规模拟定模型层次和尺寸,建立岩体数值仿真模型,通过计算机数值仿真,采用卸荷非线性力学理论对锦屏电站坝厂区岩体进行卸荷岩体宏观力学特性及参数的数值分析;同时,用RMR法和BP神经网络法对岩体宏观力学特性及参数进行分析,综合确定岩体力学参数,为设计部门提供设计参考依据.     

Upscaling of Hooke''s Law for Imperfectly Layered Rocks

In this paper, we consider the upscaling of Hooke's law and its parameters on the fine scale, to a similar law with upscaled parameters on a larger scale. It is assumed that the fine scale material properties of the rock are imperfectly layered. In the governing equations, the deviations from perfect layering introduce a small parameter that can be used in perturbation series expansions for the stress, the strain, and the displacement. In the approximation of order zero the upscaled compliance matrix contains the well-known Backus parameters; this approximation holds exactly for a perfect layering. However, many natural rock types are imperfectly layered and in that case the approximation of order zero may not be sufficiently accurate. Therefore, we consider also the first order corrections. The derivation and results are presented both for the most general case and for the much simpler case in which the fine scale Poisson ratio may be assumed constant. From thermodynamic principles, it follows that the compliance tensor is symmetric on the fine scale. However, it is shown that the argument for symmetry cannot be extended to upscaled rigidities. One of the most important conclusions is that upscaled compliance tensors are nonsymmetric when there are trends in the deviations from perfect layering.     

基于有限元法的层状岩体破裂规律探讨

层状岩体的非均质性及各向异性导致其破裂方式及规律与均质岩体有显著不同。对层状岩体分别进行不同方式的单轴、双轴、三轴试验, 分析应力-应变曲线特征; 再利用ANSYS有限元软件进行数值模拟, 观察应力、应变在岩体上的分布, 通过曲线和图件的对比分析, 并结合岩石破裂理论, 总结不同应力状态下层状岩体的破裂方式、顺序及规律; 最后以富台地区为例, 对分析结果进行验证。研究结果表明, 不同受力方式对层状岩体破裂的影响体现在施加的载荷及约束与层面的方位。当应力方向与岩层面平行时, 强度大的石灰岩岩体发生集中应力, 首先破裂; 而应力与岩层面垂直时, 强度小的泥岩岩体首先破裂。岩石试验、数值模拟结果以及实例均成功验证了这个规律。      

岩体卸荷力学特性的数值仿真研究

岩体卸荷力学性质主要表现为非线性和各向异性等特征,与常规岩体试验结果和假定存在着差别,采用数值仿真的方法对卸荷岩体力学特性进行研究,阐述卸荷岩体力学特性研究的重要性;给出了数值仿真分析研究的岩体加、卸载条件;利用大型有限元软件ADINA建立数值仿真模型,对模型进行计算分析,并和已有试验成果进行对比,得出了卸荷岩体力学特性的有关规律。     

An explicitly coupled hydro‐geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks

Modeling hydraulic fracturing in the presence of a natural fracture network is a challenging task, owing to the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and existing natural interfaces. Understanding these complex interactions through numerical modeling is critical to the design of optimum stimulation strategies. In this paper, we present an explicitly integrated, fully coupled discrete‐finite element approach for the simulation of hydraulic fracturing in arbitrary fracture networks. The individual physical processes involved in hydraulic fracturing are identified and addressed as separate modules: a finite element approach for geomechanics in the rock matrix, a finite volume approach for resolving hydrodynamics, a geomechanical joint model for interfacial resolution, and an adaptive remeshing module. The model is verified against the Khristianovich–Geertsma–DeKlerk closed‐form solution for the propagation of a single hydraulic fracture and validated against laboratory testing results on the interaction between a propagating hydraulic fracture and an existing fracture. Preliminary results of simulating hydraulic fracturing in a natural fracture system consisting of multiple fractures are also presented. Copyright © 2012 John Wiley & Sons, Ltd.