非饱和土壤水分运动参数的确定——以河南驻马店亚粘土为例

土壤水分运动参数的选取和确定是研究包气带中水分贮存和运动的基础.本文以河南驻马店亚粘土为例,进行了非饱和土壤水分运动参数的测定.通过将实测数据与经验模型拟合,将实验结果、计算结果和拟合结果进行了对比分析,发现所用的经验模型与实际数据吻合较好.研究结果表明选用的经验模型是可靠的,从而确定了土壤水分运动最优参数.     

伊宁吐拉苏火山盆地金矿成矿构造系统与远景评估

吐拉苏火山盆地金矿成矿构造系统由2种类型的成矿构造组成，一种是火山机构中的高角度塌陷式断裂；另一种是爆发．沉积相中的低角度水压式断裂．阿希金矿受控于弧形塌陷式断裂，其北东段坍塌强度大于南东段，最大张裂空间控制的矿体群向NE侧伏．阿希火山机构东侧的阿恰勒河组下可能隐伏着与西侧对应的高角度弧形塌陷式成矿断裂．爆发．沉积相中硅化岩型金矿受控于水压式断裂，不透水凝灰岩层圈闭富水粗火山碎屑岩层，组成压力仓构造并发生水力压裂作用．由水力压裂作用形成的层控水压式成矿断裂出现的几率，远大于后期由非成矿断裂切割抬升它们出现的几率．构造解析表明吐拉苏火山盆地2类成矿构造控制的金矿具有巨大的地质找矿潜力．     

Fully coupled hydraulic fracture simulation using the improved element partition method

The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

柴达木盆地东部石炭系页岩三维裂缝网络模型研究

裂缝既是页岩气的聚集空间,同时也是页岩气的运移通道, 对页岩气成藏起着重要控制作用。裂缝的产状、密度、组合特征等很大程度上决定了页岩气是否具有开发价值。该研究对柴达木盆地东部的13个露头的岩石裂缝进行了野外观测统计,实测裂缝1 000多条,选取了三个具有代表性的页岩露头进行统计分析建模。通过野外观测和室内统计分析发现,工作区页岩裂缝发育有三组(某些地段一组不发育),各组裂缝的方向满足Fisher分布,迹线长度满足对数正态分布。确定了各组裂缝方向、长度的均值和方差,以及各组裂缝的密度,建立了三维裂缝网络模型,绘出了模拟的三维裂缝网络图。裂缝三维模型为圆盘模型,即假设单个裂缝三维空间上具有圆盘形状。建模过程采用了逆建模方法,其基本原则是通过拟合实测的二维露头观测数据,优化裂缝三维大小和密度参数,使模型能够准确再现野外所观测到的实际现象,包括露头面上不同裂缝的数量、密度和长度。该研究对柴达木盆地东部石炭系页岩裂缝描述达到了定量化,为研究页岩中油气运移提供了很好的基础数据。     

Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling

The initiation and propagation of directional hydraulic fracturing (DHF) was investigated based on true tri-axial experiment and finite element modeling. The influences of notch angle, notch length and injection rate on the DHF were investigated. The initiation and propagation of DHF was modeled by a 3D nonlinear finite element method. A comparison between experimental investigation and numerical modeling results indicates that there is a good correlation between unbalanced force (UF) and fracturing. UF can be used to predict the hydraulic fracture initiation and propagation.     

结点虚流量法理论基础阐释及改进算法

含自由面的无压渗流问题本质上是一类非线性自由边值问题,固定网格的结点虚流量法在全域范围内不断扣除虚域流量贡献,从而使该问题得到求解。它具有网格依赖性小、出逸点收敛快等优点,但其内在理论基础尚未被完全揭示。通过引入互补型约束条件建立了结点虚流量法和Signorini型变分不等式提法的等价性桥梁,在此基础上引入过渡区放大系数 对自由面判别准则进行优化,并以砂槽模型试验为例进行验证。对比结果表明,优化后算法数值稳定性更好,计算结果与试验数据吻合度更高。研究成果为超大规模网格的渗控结构优化设计提供了有效分析手段。     

帽峰山地区典型小流域主要森林类型土壤渗透性

采集帽峰山地区典型小流域主要林地类型(青冈林、毛竹林、凤尾竹林、桉树林)原状土样进行土壤剖面(0~20 cm、20~40 cm、40~60 cm)渗透系数及其相关因子的实验研究,结果表明:各种林地类型土壤渗透系数在剖面上均表现出从上到下递减的特点,毛竹、青冈、凤尾竹、桉树林的土壤剖面渗透系数变化范围分别为13.64~4.91、10.49~3.63、8.16~2.96、7.92~3.71 m/d。由于20~40 cm、40~60 cm土层渗透系数差异性相对较小,土壤渗透性剖面可简化为上、下两个亚层,上部亚层渗透系数是下部亚层的1.89~2.42倍。各种林地类型土壤渗透系数与孔隙度、根密度呈正相关,与干容重呈负相关,与根径关系相关性方向具有不确定性。不同林地类型土壤渗透系数与干容重、根密度、根径之间的相关程度存在较大差异,但均表现出与孔隙度之间存在良好相关性。     

Numerical modeling of fluid‐driven fractures in cohesive poroelastoplastic continuum

In this article, we investigate the main parameters that influence the propagation of a fluid‐driven fracture in a poroelastoplastic continuum. These parameters include the cohesive zone, the stress anisotropy, and the pore pressure field. The fracture is driven in a permeable porous domain that corresponds to weak formation by pumping of an incompressible viscous fluid at the fracture inlet under plane strain conditions. Rock deformation is modeled with the Mohr–Coulomb yield criterion with associative flow rule. Fluid flow in the fracture is modeled by the lubrication theory. The movement of the pore fluid in the surrounding medium is assumed to obey the Darcy law and is of the same nature as the fracturing fluid. The cohesive zone approach is used as the fracture propagation criterion. The problem is modeled numerically with the finite element method to obtain the solution for the fracture length, the fracture opening, and the propagation pressure as a function of the time and distance from the pumping inlet. It is demonstrated that the plastic yielding that is associated with the rock dilation in an elastoplastic saturated porous continuum is significantly affected by the cohesive zone characteristics, the stress anisotropy, and the pore pressure field. These influences result in larger fracture profiles and propagation pressures due to the larger plastic zones that are developing during the fracture propagation. Furthermore, it is also found that the diffusion process that is a major mechanism in hydraulic fracture operations influences further the obtained results on the fracture dimensions, plastic yielding, and fluid pressures. These findings may explain partially the discrepancies in net pressures between field measurements and conventional model predictions. Copyright © 2012 John Wiley & Sons, Ltd.     

Structural evolution of a salt-cored,domed, reactivated fault complex,Jebel Madar,Oman

Fault and fracture patterns associated with domal structures are observed around the world both at outcrop and in the subsurface. However, the structural evolution of domes and the influence of previous fault and fracture sets are not yet fully understood. This study provides the first structural evolution analysis of Jebel Madar, a domal outcrop in the Adam Foothills of Northern Oman, and evaluates the role of multiple local tectonic events and the associated salt diapirism on fracture and fault distribution pattern development of a salt-cored domal outcrop.Analyses at Jebel Madar suggest that three local tectonic events with different stress regimes that are tentatively linked to three regional tectonic events resulted in a salt-cored, domed, reactivated fault complex: 1) initial dome-formation and NE–SW oriented mode 1 opening fractures and subsequent grabens; 2) E–W oriented dextral strike-slip faulting; and 3) reactivation and inversion of faults, and final dome formation. Salt emplacement is associated with the first and last tectonic events.This integrated study highlights the complicated evolution of fault and fracture distribution patterns at a salt-cored dome. Caution is therefore required in interpretation of similar structural patterns on top of and around salt-cored domes.