Simulation of mixed-mode fracture using SPH particles with an embedded fracture process zone

This paper focuses on the modelling of mixed-mode fracture using the conventional smoothed particle hydrodynamics (SPH) method and a mixed-mode cohesive fracture law embedded in the particles. The combination of conventional SPH and a mixed-mode cohesive model allows capturing fracture and separation under various loading conditions efficiently. The key advantage of this framework is its capability to represent complex fracture geometries by a set of cracked SPH particles, each of which can possess its own mixed-mode cohesive fracture with arbitrary orientations. Therefore, this can naturally capture complex fracture patterns without any predefined fracture topologies. Because a characteristic length scale related to the size of the fracture process zone is incorporated in the constitutive formulation, the proposed approach is independent from the spatial discretisation of the computational domain (or mesh independent). Furthermore, the anisotropic fracture responses of materials can be naturally captured thanks to the orientation of the fracture process zone embedded at the particle level. The performance of the proposed approach demonstrates its potentials in modelling mixed-mode fracture of rocks and similar quasi-brittle materials.     

Thermo-hydro-mechanical modeling of fracturing porous media with two-phase fluid flow using X-FEM technique

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media.     

Fast multipole displacement discontinuity method (FM‐DDM) for geomechanics reservoir simulations

The displacement discontinuity method (DDM) is frequently used in geothermal and petroleum applications for modeling the behavior of fractures in linear‐elastic rocks. The DDM requires O(N²) memory and O(N³) floating point operations (where N is the number of unknowns) to construct the coefficient matrix and solve the linear system of equations by direct methods. Therefore, the conventional implementation of the DDM is not computationally efficient for very large systems of cracks, often limiting its application to small‐scale problems. This work presents an approach for solving large‐scale fracture problems using the fast multipole method (FMM). The approach uses both the DDM and a kernel‐independent version of the FMM along with a preconditioned generalized minimal residual algorithm to accelerate the solution of linear systems of equations using desktop computers. Using the fundamental solutions for constant displacement discontinuity in a two‐dimensional elastic medium, several numerical examples involving fracture networks representing fractured reservoirs are treated. Numerical results show good agreement with analytical solutions and demonstrate the efficiency of the FMM implementation of the DDM for large‐scale simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Full 3-D numerical modeling of borehole electric image logging and the evaluation model of fracture

A full 3-D finite element method numerical modeling program is written based on the principle and technical specification of borehole electric image well logging tool. The response of well logging is computed in the formation media model with a single fracture. The effect of changing fracture aperture and resistivity ratio to the logging response is discussed. The identification ability for two parallel fractures is also present. A quantitative evaluation formula of fracture aperture from borehole electric image logging data is set up. A case study of the model well is done to verify the accuracy of the for-mula. The result indicates that the formula is more accurate than the foreign one.     

南黄海新构造运动

定义了南黄海新构造运动的概念;提出了南黄海新构造运动的研究方法;阐明了南黄海新构造运动的机理,即受滨西太平洋构造域控制,是在老构造基础上的继承和发展;描述了南黄海新构造运动的基本特征;总结了南黄海新构造运动的研究意义。     

冲绳海槽天然气水合物与地质构造的关系

海底天然气水合物大多与通过切穿沉积盖层的断裂的上升烃类流体相关,这些高渗透带包括底辟和泥火山等侵入构造,所以海底断裂、底辟和泥火山等构造周围可能赋存天然气水合物;其次,高沉积速率和巨厚沉积层可使有机质迅速掩埋而保存起来,为天然气水合物的生成提供充足物源,因此,邻近陆坡河谷口的海底沉积扇也是天然气水合物赋存的有利地区;另外,由于陆坡区的水合物沉积层比海盆更容易受外界温压变化的影响发生失稳分解,诱发海底滑坡,所以滑坡与天然气水合物赋存之间的关系也非常密切。冲绳海槽邻近海域具有覆水深、沉积厚度大、沉积速率高和有机质丰富等有利赋存条件,目前的研究已经在该海域发现了天然气水合物赋存的地球物理证据BSR,因此,在现有研究基础上开展断裂、泥火山、海底扇、海底滑坡等与天然气水合物相关的构造研究,可以深入了解天然气水合物在不同地质构造中的分布特征与演化,为更精确地评估其资源潜力提供参考。     

土钉支护稳定分析的模拟退火-随机投点耦合算法

充分地考虑土钉支护的特殊性,建立了一个比较精确的基于一般破裂面的土钉支护稳定性分析方法。对土钉支护内部稳定分析和外部稳定分析分别建立其优化数值分析模型;然后导入模拟退火算法和随机投点法,结合两者的特长建立了一种可有效求解该模型的耦合算法,并编制了相应的计算程序;通过工程实例证明了基于一般破裂面的土钉支护稳定性分析方法的有效性,并得出了若干对土钉支护稳定分析具有参考意义的结论。     

冻土冲击韧度与凿碎比能关系的试验研究

冻土的冲击韧度与试块破坏断口的平整度反映了冻土的冲击特性。冻土凿碎比能是表示凿碎单位体积冻土所消耗的功,它反映了冲击破碎冻土的能力。以冻结粘土为试样,利用摆锤式冲击试验仪进行不同温度、不同冲击能量的带缺口的冲击实验,利用凿测器进行了冻土凿碎比能试验。实验结果表明,冻结粘土的冲击韧度随温度降低而增大,随冲击能量的增加而增大。随温度的降低,试样的破坏断口的凹凸不平度增大,其范围在±3mm之间。随着温度的降低,冲击凿入深度越来越浅,冻土的硬度越来越大,钻进程度越来越难。通过回归分析,凿碎比能随着冲击韧度的增加而增大,且在冲击韧度一定的情况下,随着凿入次数的增加,凿碎比能增加的值越来越小,凿碎比能与冲击韧度的相关性较好,表明用冻土的冲击韧度来反映凿碎比能是可行的。     

裂隙岩体渗流概念模型研究

裂隙岩体中的渗流和传统的多孔介质渗流在机理上存在本质的差别,这种差别主要表现为裂隙岩体在各种尺度上存在的非均质性。模拟裂隙岩体渗流的主要困难在于描述这种非均质性。目前的概念模型,包括等效连续体模型、离散裂隙网络模型和混合模型使用了不同的技术来预测裂隙岩体中的渗流。这些模型基于不同的假设和概念框架,有着各自不同的优缺点。在实际应用时,应当根据研究域的具体特点和所要解决的问题的要求对其选择,此外,还讨论了单裂隙的概念模型。     

有限宽切槽对CCNBD断裂试样应力强度因子的影响

国际岩石力学学会（ISRM）在1995年提出一种新型的岩石断裂韧度试样--人字形切槽巴西圆盘试样（cracked chevron notched Brazilian disc--CCNBD）,对该试样的一个重要力学参数即最小无量纲应力强度因子的标定,以前的分析和计算都没有考虑切槽宽度的影响。然而试样切槽的宽度受切割刀具厚度所限,不能为零。当试样较小时,切槽宽度则相对较大。通过三维边界元计算分析表明,切槽宽度越大,无量纲应力强度因子的标定值就越大;对于ISRM推荐的CCNBD标准试样,得出其最小无量纲因子值为0.954,这比ISRM给出对应值0.84要大13.6 %。同时,小裂纹应力强度因子曲线的变化趋势也发生了质的变化,这可能会导致实验的失败。推荐最小无量纲应力强度因子的标定采用考虑切槽的三维分析。