Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method          下载免费PDF全文

A. Verde  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(6):865-886

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Displacement velocity effects on rock fracture shear strengths     

《Journal of Structural Geology》2016

Triaxial shear tests are performed to assess the effects of displacement velocity and confining pressure on shear strengths and dilations of tension-induced fractures and smooth saw-cut surfaces prepared in granite, sandstone and marl specimens. A polyaxial load frame is used to apply confining pressures between 1 and 18 MPa with displacement velocities ranging from 1.15 × 10⁻⁵ to 1.15 × 10⁻² mm/s. The results indicate that the shearing resistances of smooth saw-cut surfaces tend to be independent of the displacement velocity and confining pressure. Under each confinement the peak and residual shear strengths and dilation rates of rough fractures increase with displacement velocities. The sheared-off areas increase when the confining pressure increases, and the displacement rate decreases. The velocity-dependent shear strengths tend to act more under high confining pressures for the rough fractures in strong rock (granite) than for the smoother fractures in weaker rocks (sandstone and marl). An empirical criterion that explicitly incorporates the effects of shear velocity is proposed to describe the peak and residual shear strengths. The criterion fits well to the test results for the three tested rocks.  相似文献   

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation     

《Computers and Geotechnics》2016

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.  相似文献   

渝东南下古生界页岩构造裂缝形成及分布控制因素          下载免费PDF全文

岳锋  程礼军  焦伟伟  王飞 《地质科学》2016,(4):1090-1100

天然裂缝对页岩气储集和渗流有重要影响,但目前对页岩构造裂缝研究不够深入。本文通过露头、岩心裂缝观察和分析,探讨了渝东南下古生界页岩构造裂缝类型、形成机理及分布控制因素。基于裂缝特征及力学成因,将剪切裂缝分为高角度剪切裂缝、倾斜滑脱裂缝和水平滑脱裂缝。倾斜滑脱裂缝是在上覆岩层重力和水平构造应力共同作用下沿应力集中的软弱面发生剪切滑动形成,水平滑脱裂缝是在构造挤压应力作用条件下主要沿页理面方向的剪切或层间滑动形成;页理发育程度及岩层曲率是控制水平滑脱裂缝形成的关键因素,岩石矿物组成、构造作用及岩层厚度是控制其它类型构造裂缝形成和分布的主要因素,页岩岩石力学层的划分需要根据裂缝特征、岩性变化特征、岩石力学参数及沉积界面等综合确定,岩层厚度与层间构造裂缝密度呈负相关关系。  相似文献   

北大巴山闹阳坪锌萤石矿区断裂特征与找矿预测     

张文高  杨兴科  韩鹏飞  樊培贺  何虎军  Moussounda Kounga Claude 《大地构造与成矿学》2016,(2):323-334

闹阳坪锌萤石矿床位于北大巴山逆冲推覆褶皱带平利隆起东侧,为该区首次发现的受构造裂隙控制的气水–热液型锌萤石矿床。本文从矿区断裂入手,通过对该区矿床地质、断裂特征、矿体特征的研究,总结了该区断裂的演化序列及其对成矿的控制机制与规律。认为NW-SE向断裂组是矿区的主要控矿构造,矿区断裂变形发育演化序列为晚印支–早燕山期形成近EW向断裂F1,早燕山期形成NW-SE向断裂F7和NE-SW向断裂F4,随后的右行走滑作用叠加有张扭性应力,形成有利于成矿热液运移充填成矿的张扭性右行走滑断裂。并明确了成矿期应变椭球体,恢复了成矿期主压应力方向为NNW向(340°~350°)。在此基础上,预测平面上闹阳坪矿区F7与F8断裂之间为下一步找矿重点地段,剖面上K3萤石矿体下伏北东方向可能存在隐伏矿体。  相似文献   

江西千子岭铅锌矿床地质特征简介   总被引：1，自引：0，他引：1  

何发林  周先军  李宗朋  罗振辉  李淑琴 《矿产与地质》2016,(5):768-773

江西千子岭铅锌矿床是钦杭成矿带江西段核部新发现的一处隐伏矿床,成矿定位于信江中新生代断陷盆地北缘大断裂;铅锌矿体赋存在EW向硅化破碎带中,矿石为块状热液硅质岩-闪锌矿-方铅矿(黄铜矿)-萤石组合,容矿热液硅质岩从地表到深部表现为无矿化→铅锌矿化→铅锌铜矿化的垂向矿化分带特征。EW向容矿破碎带卷入上白垩系紫红色砂岩,推断铅锌(铜)矿化属喜山运动的产物。  相似文献   

甘肃党河南山地区石块地金矿地质地球化学特征     

贾宝林  董泰山  金辉 《甘肃地质》2016,25(3):58-63

石块地金矿位于党河南山金成矿带上,是由断裂构造控制的岩浆热液型金矿。通过1∶5万水系沉积物测量,在本区圈定了As-3综合异常,并进一步利用1∶1万土壤测量进行解剖,发现了石块地金异常区,经过后期工程揭露,发现了有价值的金矿体7条。综合研究表明,勘查区金矿体赋存岩性均为花岗斑岩,并严格受NW—SE向断裂构造的控制。主矿体的延伸方向与土壤测量圈定的Au异常长轴展布方向一致,充分证明了在该区运用化探手段找矿的重要指导作用。  相似文献   

三维高密度电阻率法在堤坝工程中的模拟试验     

张汝凯  张晓培  吕守航  陈兵兵  王杰 《世界地质》2016,35(4):1133-1137

为检测堤坝中裂隙、洞穴及渗漏层等隐患的存在,利用三维高密度电阻率法对数字模型和物理模型建立的裂隙隐患进行探测,结果表明室内建立的物理模型在立体图和分层切片图中裂隙的形态有些变化,但其分布范围、位置与相对各自模型中实际的比例尺寸相差不大,说明数值模拟得出的结论可以为实际工程进行理论指导,具有一定的参考价值。  相似文献   

构造物理模拟和PIV有限应变分析对构造裂缝预测的启示   总被引：1，自引：0，他引：1       下载免费PDF全文

沈礼  贾东  尹宏伟  魏东涛  陈竹新  孙闯  崔键 《高校地质学报》2016,22(1):171

油气储层中构造裂缝发育与有限应变状态关系密切,为了探索有限应变分析与构造裂缝预测的新技术方法,此次研究设计完成了一组单侧挤压收敛模型的物理模拟实验,并引入粒子图像测速(PIV,Particle Image Velocimetry)技术对实验过程进行了定量化分析。实验模型在垂向上为含粘性层的多层结构,实验结果形成了一个肉眼可见的箱状褶皱。通过PIV技术可以获取实验模型变形演化过程中各阶段的位移场数据,计算出各阶段的增量应变,实现从初始状态到褶皱形成之后整个变形过程的有限应变分析,探讨构造裂缝成因机制和分布规律,进行定量化裂缝预测。挤压变形过程初期,应变分布范围很广,有限应变较弱(约4%~8%),在挤压方向上的线应变表现为弱压应变,在垂向上的线应变表现为弱张应变,这种现象是褶皱和断层产生前平行层缩短和层增厚的纯剪变形结果,也是区域型张裂缝和剪裂缝形成的主要机制。褶皱和断层即将发育之时至发育之后,应变局限在断层发育的剪切带及附近区域,有限应变表现为较强(达20%)的剪切应变和剪切张应变,是断层面附近简单剪切变形作用的结果,也是局部型剪裂缝和张剪裂缝形成的主要机制。  相似文献   

Stability and capsizing analysis of nonlinear ship rolling in wind and stochastic beam seas     

《Applied Ocean Research》2016

Considering the actual seaway condition, stability and capsizing of nonlinear ship rolling system in stochastic beam seas is of significant importance for voyage safety. Safe zone are defined in the phase space plan of the unperturbed Hamilton system to qualitatively distinguish ship motions as capsize and noncapsize. Capsize events are defined by solutions passing out of the safe zone. The probability of such an occurrence is studied by virtue of the random Melnikov function and the concept of phase space flux. In this paper, besides conventional wave excitation, the effect of wind load is also taken into account. The introduction of wind load will lead to asymmetry, in other words, it transforms the symmetric heteroclinic orbits into asymmetric homoclinic orbits. For asymmetric dynamical system, the orbital analytic solutions and its power spectrum are not readily available, and the technique of discrete time Fourier transformation (DTFT) is used. In the end, as verification of theoretical critical significant wave height, capsizing probability contour diagram is generated by means of numerical simulation. The contour diagram shows that these analytical methods provide reliable and predictive results about the likelihood of a vessel capsizing in a given seaway condition.  相似文献