层理对岩石抗压强度影响研究
- 格式:pdf
- 大小:1.94 MB
- 文档页数:2
Chenmical Intermediate当代化工研究2016·0539技术应用与研究层理对岩石抗压强度影响研究OO穆景福(陕西延长石油(集团)有限责任公司研究院OO陕西OO710075)摘要:页岩作为页岩气的主要聚集岩层,由于其具有薄页状的层理,使得页岩与其它常规储层岩石力学性质显著不同,具有明显的各向异性特征。
依然假定页岩为各向同性岩体,必然会产生不可避免的错误,在现场工作中造成不必要的损失。
本文对不同层理倾角页岩进行单轴和三轴压缩实验获得力学参数,研究页岩的各向异性强度特征以及变形各向异性特征。
由于离散元理论适用于不连续岩石的特点,利用离散元软件PFC2D对含有层理页岩进行数值模拟研究,实现对页岩力学参数的确定,从而研究不同角度下含有单一层理和含有一组层理页岩的抗压强度。
关键词:页岩;各向异性;颗粒流软件中图分类号:TU452 文献标识码:AExploration of the anisotropic mechanical of shaleMu Jingfu(Research institute, Shaanxi Yanchang Oil(Group) Co., Ltd, Shanxi , 710075)Abstract :Shale which is the set of rock of shale gas is obviously anisotropic different from conventional reservoir rock, since it has thin pageof bedding. It must cause unnecessary loss in the actually work if still assumes that shale is isotropic. In this paper the anisotropy of shale strength and deformation is studied base on the mechanical of rock of different angle of bedding which are obtained through the uniaxial and triaxial compression test. The anisotropy of shale of single bedding and a group of bedding is studied. The mechanical of shale is simulated by discrete element software PFC2D.Keyword :Shale ;Anisotropy ;PFC引言自然界岩体中存在着大量的节理、层理等不连续面,具有明显的各向异性。
层理弱面的存在对岩体的强度以及变形规律产生明显的影响,依然假设岩体为连续介质,将会产生不可避免的错误。
当前,许多学者开展了大量关于节理、层理弱面对岩体力学性质影响的试验研究。
但是由于岩石的不均匀性、复杂性以及技术原因,室内试验获得力学参数费时、费力而且花费高。
如若将室内试验和数值模拟试验结合,将会更好的研究和了解层理面的存在对岩体力学性质的影响。
对于含有层理的不连续体,不宜采用连续介质模型或 “等效模型”进行分析,采用离散元软件PFC2D进行仿真模拟,能更准确的分析含有层理岩石试样的力学参数。
一、室内试验1.试样制备及实验过程岩石单轴压缩实验是指,是岩石试样在无侧压的情况下,在轴向压力作用下发生变形直至破坏的试验。
页岩层理发育明显,所以本次试验所用岩石试样取自志留系龙马溪组中厚层状黑色页岩,碳质含量高,试验仪器为全伺服控制岩石力学实验系统。
由于取自露头的页岩形状很不规则,不能直接用于实验。
实验前必须对页岩进行加工,由于取芯设备、时间和技术原因,只能分别与层理方向为0度,45度,90度角度取芯并编号,岩石试样为圆柱型,直径为25mm,长度为50mm。
轴向和径向应变仪夹持并固定在试件上,试件放入三轴应力试验仪器,无围压,轴向加载速度为0.35mm/min,通过数据采集系统获得试验过程中的轴向应力,从而得到该试件的抗压强度。
2.试验结果岩石是矿物颗粒集合体,具有明显的非均质性,同时由于岩石在长期复杂的地质作用产生裂隙、空洞等缺陷,所以试样的力学参数具有明显的差异性,离散度大。
本次试验三个层理倾角的取样不少于5块,以减少试样及试验产生的误差。
层理倾角为0度岩心抗压强度分布为101.1MPa~170.4Mpa,平均值137.0Mpa;层理倾角为45度岩心抗压强度分布为71.3MPa~106.9MPa,平均值90.4Mpa;层理倾角为90度岩心抗压强度分布为94.9Mpa~158.6MPa,平均值132.4MPa。
随着倾角的增大,试样的抗压强度先减小后增大层理倾角不同,岩石的抗压强度有明显的差异。
二、层理试件数值试验模拟PFC2D程序作为一种特殊的离散单元法,克服了传统的连续介质力学模型的宏观连续性假设,无需事先指定岩石或节理的本构关系, 可以解决基于连续介质力学建立本构模型难以处理的复杂岩石断裂力学问题。
PFC 模型是直接破坏型模型, 其强度并不按指定的应力—应变关系发生变化, 而是随着内部微观结构发生改变。
1.试验模型建立PFC2D模型建立的依据是所取岩心试件室内试验结果获得的。
根据室内试验获得的0度、45度、90度岩石试件抗压强度、弹性模量和泊松比的平均值,通过不断调试程2016·0540技术应用与研究Chenmical Intermediate当代化工研究序,获得与室内试验相符合的数值模型的颗粒微观力学参数。
数值试验选择点接触本构模型,模型为长方形,长度为100mm,宽度为50mm,最小颗粒半径0.3mm,颗粒半径比1.66,颗粒密度2600kg/m3,力学参数与实际基本一致。
建立层理面倾角为0度、15度、30度、45度、60度、75度、90度等七种不同倾角的模型,微观力学参数分别为法向粘结力为70Pa,剪切粘结力为420Pa,摩擦系数为0.1。
层理面位于模型的中心位置,含有一组层理面模型层理面间距为12mm。
2.数值试验过程含有层理模型生成以后,进行双轴压缩模拟环境的建立。
模型周围的四个墙单元来对模型施加力或者位移从而完成实验模拟。
对于双轴实验,上端墙单元与下端墙单元作为加载板,通过伺服控制机制,侧面两个墙单元施加速度,从而保持恒定的围压。
关闭侧面墙单元的伺服控制,删除侧面的两个墙单元,就可以对模型进行单轴压缩实验。
上方墙单元与下方墙单元分别以速度Vp 彼此移动靠近从而产生压力进行双轴压缩实验。
如果速度设定为一个单一的速度,如此大的加速度产生的施加给试件的惯性力将会对试件产生破坏。
为了避免惯性力的产生,设定一个小的加速步S p 在Np步循环以后达到最终速度Vp。
对双轴压缩环境设定结束后,即可进行单轴和双轴压缩试验数值模拟。
3.单轴压缩数值试验结果分析如上所述建立模型,进行单轴压缩实验,获取倾角45度试件数值模拟应力应变曲线,并与实验应力应变曲线对比。
红色虚线为数值试验模拟结果,蓝色实线为室内试验获得的数据,基于PFC2D的数值试验能很好的模拟含有层理岩石的室内压缩试验,所获得的数值模拟力学参数与真实试验获得的力学参数基本一致。
通过PFC2D程序内置FISH语言,编制程序,获得含有单一层理和一组(三条)层理在7种层理面倾角情况下的抗压强度。
对比室内试验数据与数值模拟数据,数据差异控制在6%以内,差异较小。
对于含有单一层理模型抗压强度曲线呈现“单肩型”,符合Jeager提出的单一弱面理论,在0度到30度时,抗压强度基本相等变化不明显;30度到60度时,抗压强度明显减小,并在60时取得最小值;60度到90度时,抗压强度逐渐增大,90度与0度抗压强度基本相等。
对于含有一组层理模型抗压强度曲线不再为“单肩 型”,在0度到30度时,抗压强度略微减小;30度到60度时,抗压强度明显减小;在60度到75度之间取得最小值;60度到90度时,抗压强度逐渐增大;90度与0度抗压强度基本相等。
含有一组层理模型抗压强度小于同一倾角的含有单一层理模型的抗压强度,层理面的数量影响模型的抗压强度,层理面增多,抗压强度减小。
4.双轴压缩数值试验结果分析修改目标围压值为10M P a,20M P a,30M P a, 40MPa,50MPa,运行程序,模拟在多种围压情况下,数值试件的抗压强度的变化。
随着围压增大,试件的抗压强度明显增大,最低抗压强度值依然在60度倾角试件获得,曲线逐渐趋于平滑,呈现出“双肩型”,各向异性性质逐渐不明显,当围压达到单轴抗压强度一半时,试件中存在的弱面即不在发生作用。
对含有一组三条平行层理试件分别施加10M P a、 20MPa、30MPa、40MPa和50MPa围压进行双轴试验数值模拟。
在相同围压下,抗压强度曲线随着试件倾角的增大表现为先减小后增大,在60度时取得最小值。
在围压下,由于含有多条层理弱面,弱面效应对于低角度和高角度倾角也变得明显,含有一组层理试件曲线不再与含有单一层理的 “双肩型”曲线相同。
三、结论通过对含有层理页岩进行室内力学参数试验以及含有单一层理和一组层理试件模型的数值试验,得出以下主要结论:1.0度、45度、90度三种倾角试件室内单轴试验表明,随着角度增大,单轴抗压强度先减小后增大,层理弱面效应明显。
2.颗粒流程序能很好的模拟含有层理岩石的数值模拟试验。
单一层理与一组(三条)层理模型抗压强度倾角效应曲线不同;含有单一层理模型抗压强度倾角效应曲线为“单肩型”,而含有一组(三条)层理模型抗压强度倾角效应曲线,表现为先缓慢减小到急剧减小然后急剧增大的趋势。
3.在施加围压情况下,含有单一层理模型抗压强度倾角效应曲线为“双肩型”,一组层理模型表现为先减小后增大的趋势。
参考文献[1]张鸿翔,页岩气:全球油气资源开发的新亮点——我国页岩气开发的现状与关键问题.中国科学院院刊,2010(04). [2]周健,池毓蔚,池永等.砂土双轴试验的颗粒流模拟[J].岩土工程学报,2011,23(6):701—704.[3]Cundall P A.Practice Flow Code in 2 Dimensions [M].Minnesota: Itasca Consulting Group,Inc,1996.[4]Cundall P A. Practice Flow Code in 2 Dimensions Fish in pfc[M].Minnesota: Itasca Consulting Group, Inc,1996.[5]Cundall P A. PFC2D User’s Manual [M].Minnesota: Itasca Consulting Group,Inc,1996.[6]Jeager J C.Shear failure of anisotropic rocks[J] .Geol.Mang.1960,97,65-72.[7]Hoek E. Brittle failure of rock. In Rock Mechanics Engineering Practics[J],G Stagg and O.C.Zienkiewicz.eds.,Wiley,London.1968,99.•【作者简介】穆景福(1985~),男,陕西延长石油(集团)有限责任公司研究院,研究方向:储层改造。