CN102608021A - 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 - Google Patents
一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 Download PDFInfo
- Publication number
- CN102608021A CN102608021A CN201110385585XA CN201110385585A CN102608021A CN 102608021 A CN102608021 A CN 102608021A CN 201110385585X A CN201110385585X A CN 201110385585XA CN 201110385585 A CN201110385585 A CN 201110385585A CN 102608021 A CN102608021 A CN 102608021A
- Authority
- CN
- China
- Prior art keywords
- sand
- fracture
- thin interbed
- shale
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明是一种砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法,包括施工排量和砂比。对砂泥岩薄互层储层的压裂来说,由于岩石塑性强,断裂韧性高,裂缝启裂困难,加之裂缝形态复杂,砂堵事故高发。从压裂设计角度来说,优化排量和砂比是控制砂堵的有效手段,但对砂泥岩薄互层储层尚缺少定量的确定方法。高含泥岩石的应力应变关系是非线性弹塑性,无法用解析方法表达其扩展压力。本发明将高含泥地层视为弹塑性非线性材料,通过力学建模和有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力,由缝内净压力与排量的关系确定排量;将计算砂泥岩地层裂缝宽度剖面看作平面应变问题,得到裂缝扩展形态,用砂泥岩缝宽比来计算临界砂比。
Description
技术领域
本发明是关于一种石油工程领域砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法,主要包括施工排量和砂比,以保证压裂施工过程不发生砂堵事故,同时该发明也适用于其他类型的高含泥塑性储层。
背景技术
砂堵是指加砂压裂过程砂粒聚集在射孔附近或靠近井筒的裂缝中,堵塞压裂液和砂粒的输运路径,造成泵压急剧上升、无法加砂等事故。砂堵对压裂施工造成重大影响,甚至导致压裂失败。影响砂堵的因素较复杂,主要包括地层因素、压裂液因素、工程因素、设计因素等,如何有效避免砂堵的产生一直是压裂设计和施工作业中感到棘手的难题。对超低渗致密储层的压裂来说,储层埋藏深,成岩作用强,且泥质含量高,导致岩石塑性强,断裂韧性高,裂缝启裂困难。若砂泥岩形成薄互层,导致裂缝形态进一步复杂,砂堵事故高发。
砂泥岩薄互层储层压裂涉及高含泥地层中的裂缝扩展。对于高含泥岩石,由于其应力应变关系是非线性的弹塑性,无法用解析方法表达其扩展压力。从压裂设计角度来说,优化排量和砂比是控制砂堵的有效手段,但对砂泥岩薄互层储层尚缺少定量的确定方法。本发明基于高含泥弹塑性地层的力学建模,通过有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力,进而确定施工排量;将计算砂泥岩地层裂缝宽度剖面看作平面应变问题,得到裂缝扩展形态,进而确定临界砂比。
发明内容
本发明的目的是提供砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法。该方法将高含泥地层视为弹塑性非线性材料,通过力学建模和有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力,由缝内净压力与排量的关系确定排量;将计算砂泥岩地层裂缝宽度剖面看作平面应变问题,得到裂缝扩展形态,用砂泥岩缝宽比来计算临界砂比。
本发明目的是这样实现的,将井眼周边的应力状态视为平面应变,根据对称性,建立水力裂缝的半边模型,如图1所示(裂缝端部局部放大如图2所示)。模型尺寸2m*2m,由于对称性,模型左侧固定,模型右侧施加最大有效地应力,上下施加最小有效地应力,井筒及压裂裂缝内为压裂液。为考虑裂缝扩展到不同长度时的缝宽分布情况,模型尺寸采用无量纲化形式,其中初始裂缝长度1000单位。裂缝扩展准则是:对砂岩来说是弹性材料,采用断裂韧性表征,对泥岩来说是塑性材料,采用抗拉强度表征。
分别计算砂岩及泥岩中裂缝扩展所需的缝内净压力。计算砂岩时,材料视为弹性,只需弹性模量及泊松比;计算泥岩时,材料视为非线性弹塑性,除弹性模量、泊松比,还需提供实验测定的应力应变曲线。计算结果显示泥岩扩展压力较砂岩高,说明砂泥岩薄互层储层的扩展压力由泥岩扩展压力决定。由于缝内压裂液的压力与体积、排量成比例,与均匀
考虑水力裂缝的长度较高度、宽度远远大,此时可以将计算裂缝宽度剖面看作平面应变问题,建立如图3所示的力学模型。模型尺寸为5m*5m,沿垂直方向划分3层泥岩夹2层砂岩。施加地应力条件和上覆压力条件,井筒压力应高于孔隙压力与泥岩扩展压力之和,以确保裂缝能够延伸。边界条件为模型的上下端及右端固定,缝内压力作用在左端。
缝内净压力的一般模式为三次多项式
Px(x,y)=A10+A10x+A01y+A20x2+A11xy+A01y2-
A30x3+A21x2y+A12xy2+A33y3 (1)
如果砂、泥层内最小地应力均匀,可以忽略奇函数项,则缝内压力可以进一步简化为
Pn(x,y)=A00+A20x2+A02y2 (2)
附图说明
以下附图仅旨在于对本发明做示意性说明和解释,并不限定本发明的范围。其中,
图1:水力裂缝的半边模型
图2:裂缝端部局部放大
图3:砂泥岩互层中的裂缝形状计算模型
图4:砂泥岩互层裂缝形状
图5:岩石非线性弹塑性应力应变关系
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现举例说明本发明的具体实施方式。
1.由岩石力学实验确定砂岩弹性模量、泊松比、断裂韧性及泥岩的弹性模量、泊松比、抗拉强度、弹塑性应力应变关系。例如某砂岩弹性模量E=23GPa,泊松比v=0.18,抗拉强度St=4.7MPa,断裂韧性KIc=19.5MPa·mm1/2;泥岩弹性模量E=11GPa,泊松比v=0.34,抗拉强度St=10.5MPa,应力应变关系如图5所示。
2.查询或测量目标地层的地应力状态,例如两个水平主应力为46MPa、37MPa,上覆压力60MPa,孔隙压力28MPa。计算出最小有效地应力为37-28=9MPa,最大有效地应力为46-28=18MPa,上覆有效压力60-28=32MPa。
3.建立如图1的水力裂缝的半边有限元数值模型。分别计算砂岩及泥岩中裂缝扩展所需的缝内净压力,此处分别是9MPa、11MPa。
5.确定缝内压力的垂向分布。我们在砂岩中心位置处射孔,压力项A00等于井筒有效压力(井筒压力49-孔隙压力28=21MPa),A20=0,泥岩中心处的净压力为泥岩层的扩展压力(11MPa),由此得到压力分布Pn(y)=21-9.33y2。
6.将计算裂缝宽度剖面看作平面应变问题,建立如图3所示的力学模型,模型尺寸为5m*5m,沿垂直方向划分3层泥岩夹2层砂岩,每层厚度1m,可以在通用有限元软件中完成。计算得到裂缝宽度剖面(图4)。此处泥岩层段平均缝宽1.4mm,相同缝内压力条件下的砂岩平均缝宽3.6mm。
Claims (3)
1.砂泥岩薄互层储层压裂裂缝扩展压力计算方法,其特征在于:砂泥岩薄互层中裂缝扩展压力由泥岩扩展压力确定,由于泥岩是塑性材料,扩展准则采用抗拉强度表征,而不能采用断裂韧性表征;根据对称性建立水力裂缝的半边有限元数值模型,模型左侧固定,模型右侧施加最大地应力,上下施加最小地应力,井筒及压裂裂缝内为压裂液;为考虑裂缝扩展到不同长度时的缝宽分布情况,计算模型的尺寸采用无量纲化形式,其中初始裂缝长度1000单位。
3.砂泥岩薄互层储层压裂临界砂比的确定方法,其特征在于:砂岩/泥岩薄互层中不可避免出现裂缝宽度分布的不均匀,泥岩层段裂缝宽度很窄,导致裂缝形态呈现类似“葫芦串”的形状,加砂过程中砂粒容易在喉结处聚集堵塞造成砂堵,据此提出砂比上限可用砂泥岩缝宽比来计算,即
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110385585.XA CN102608021B (zh) | 2011-11-29 | 2011-11-29 | 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110385585.XA CN102608021B (zh) | 2011-11-29 | 2011-11-29 | 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102608021A true CN102608021A (zh) | 2012-07-25 |
CN102608021B CN102608021B (zh) | 2014-04-02 |
Family
ID=46525579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110385585.XA Expired - Fee Related CN102608021B (zh) | 2011-11-29 | 2011-11-29 | 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102608021B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
CN106761647A (zh) * | 2017-01-13 | 2017-05-31 | 黄仲尧 | 一种估算页岩储层压后平面改造面积的方法 |
CN106951659A (zh) * | 2017-04-05 | 2017-07-14 | 中国石油大学(北京) | 一种高效实现煤层及煤层顶部砂岩合压的射孔方案 |
CN105373648B (zh) * | 2015-10-09 | 2018-10-12 | 西南石油大学 | 一种砂泥岩互层型沉积体储层构型的建模方法 |
CN111368463A (zh) * | 2020-04-08 | 2020-07-03 | 广州海洋地质调查局 | 一种水平井穿层压裂设计方法 |
CN111396014A (zh) * | 2020-03-16 | 2020-07-10 | 中国石油大学(北京) | 薄互层储层改造方法、装置和设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100999990A (zh) * | 2006-11-10 | 2007-07-18 | 大庆油田有限责任公司 | 低渗透油层垂直缝定向射孔压裂方法 |
CN101968348A (zh) * | 2010-09-07 | 2011-02-09 | 中国石油大学(北京) | 对压裂裂缝进行可视化监测的方法 |
CN102220865A (zh) * | 2011-05-13 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种检测灰岩地层孔隙压力的方法 |
-
2011
- 2011-11-29 CN CN201110385585.XA patent/CN102608021B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100999990A (zh) * | 2006-11-10 | 2007-07-18 | 大庆油田有限责任公司 | 低渗透油层垂直缝定向射孔压裂方法 |
CN101968348A (zh) * | 2010-09-07 | 2011-02-09 | 中国石油大学(北京) | 对压裂裂缝进行可视化监测的方法 |
CN102220865A (zh) * | 2011-05-13 | 2011-10-19 | 中国石油天然气股份有限公司 | 一种检测灰岩地层孔隙压力的方法 |
Non-Patent Citations (1)
Title |
---|
董长银等: "砾石充填防砂工艺参数优化设计", 《中国石油大学学报(自然科学版)》, vol. 30, no. 5, 31 October 2006 (2006-10-31), pages 57 - 61 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105373648B (zh) * | 2015-10-09 | 2018-10-12 | 西南石油大学 | 一种砂泥岩互层型沉积体储层构型的建模方法 |
CN105550441A (zh) * | 2015-12-12 | 2016-05-04 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
CN105550441B (zh) * | 2015-12-12 | 2018-05-18 | 山东科技大学 | 一种基于连续介质的工程岩体破裂劣化数值模拟方法 |
CN106761647A (zh) * | 2017-01-13 | 2017-05-31 | 黄仲尧 | 一种估算页岩储层压后平面改造面积的方法 |
CN106761647B (zh) * | 2017-01-13 | 2020-08-14 | 中国石油化工股份有限公司 | 一种估算页岩储层压后平面改造面积的方法 |
CN106951659A (zh) * | 2017-04-05 | 2017-07-14 | 中国石油大学(北京) | 一种高效实现煤层及煤层顶部砂岩合压的射孔方案 |
CN111396014A (zh) * | 2020-03-16 | 2020-07-10 | 中国石油大学(北京) | 薄互层储层改造方法、装置和设备 |
CN111368463A (zh) * | 2020-04-08 | 2020-07-03 | 广州海洋地质调查局 | 一种水平井穿层压裂设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102608021B (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102608021A (zh) | 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 | |
CN105019877B (zh) | 一种页岩水平井压裂破碎区体积的计算方法 | |
CN109374867B (zh) | 一种基于离散元的砂砾岩水力压裂模拟方法 | |
Castro et al. | Consolidation around stone columns. Influence of column deformation | |
Zheng et al. | Numerical investigation on the hydraulic fracture propagation based on combined finite-discrete element method | |
Wang et al. | A fully 3-D, multi-phase, poro-elasto-plastic model for sand production | |
Wang et al. | A 3-D poro-elasto-plastic model for sand production around open-hole and cased & perforated wellbores | |
WO2014176338A3 (en) | Compositions and methods for use of proppant surface chemistry to improve proppant consolidation and flowback control | |
Shahsavari et al. | Mine backfill porewater pressure dissipation: numerical predictions and field measurements | |
Saberhosseini et al. | A new geomechanical approach to investigate the role of in-situ stresses and pore pressure on hydraulic fracture pressure profile in vertical and horizontal oil wells | |
Jeng et al. | Three-dimensional consolidation of a porous unsaturated seabed under rubble mound breakwater | |
Ye et al. | Consolidation and dynamics of 3D unsaturated porous seabed under rigid caisson breakwater loaded by hydrostatic pressure and wave | |
CN116306385A (zh) | 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质 | |
Karev et al. | Modelling of mechanical and filtration processes near the well with regard to anisotropy | |
CN104632173B (zh) | 非天然裂缝致密储层缝网压裂选层方法 | |
Zhou et al. | Rock dilation and its effect on fracture transmissivity | |
Jahani et al. | Prediction of rock strength with matrix acidizing stimulation and induced wormhole by computational simulation methods | |
Guoying | Method to predict proppant flowback in fractured gas well | |
CN104879111A (zh) | 一种实现压裂裂缝参数优化的方法和装置 | |
Nassir et al. | Permeability enhancement of oil sands | |
Kowalska | Influence of the ratio between dilatancy angle and internal friction angle on stress distribution behind a gravity retaining wall | |
Ann et al. | Finite element analysis of a soil nailed slope-some recent experience | |
CN203960807U (zh) | 适于天然主、支流河道交汇处的透水耐冲分流坝 | |
Wan et al. | Prediction of volumetric sand production using a coupled geomechanics-hydrodynamic erosion model | |
Ran et al. | Research on borehole stability of shale based on seepage-stress-damage coupling model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140402 Termination date: 20141129 |
|
EXPY | Termination of patent right or utility model |