CN102608021B - 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 - Google Patents

一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 Download PDF

Info

Publication number
CN102608021B
CN102608021B CN201110385585.XA CN201110385585A CN102608021B CN 102608021 B CN102608021 B CN 102608021B CN 201110385585 A CN201110385585 A CN 201110385585A CN 102608021 B CN102608021 B CN 102608021B
Authority
CN
China
Prior art keywords
sand
fracture
shale
pressure
mud
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.)
Expired - Fee Related
Application number
CN201110385585.XA
Other languages
English (en)
Other versions
CN102608021A (zh
Inventor
赵海峰
陈勉
金衍
张广清
侯冰
张伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201110385585.XA priority Critical patent/CN102608021B/zh
Publication of CN102608021A publication Critical patent/CN102608021A/zh
Application granted granted Critical
Publication of CN102608021B publication Critical patent/CN102608021B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明是一种砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法，包括施工排量和砂比。对砂泥岩薄互层储层的压裂来说，由于岩石塑性强，断裂韧性高，裂缝启裂困难，加之裂缝形态复杂，砂堵事故高发。从压裂设计角度来说，优化排量和砂比是控制砂堵的有效手段，但对砂泥岩薄互层储层尚缺少定量的确定方法。高含泥岩石的应力应变关系是非线性弹塑性，无法用解析方法表达其扩展压力。本发明将高含泥地层视为弹塑性非线性材料，通过力学建模和有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力，由缝内净压力与排量的关系确定排量；将计算砂泥岩地层裂缝宽度剖面看作平面应变问题，得到裂缝扩展形态，用砂泥岩缝宽比来计算临界砂比。

Description

一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法
技术领域
本发明是关于一种石油工程领域砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法，主要包括施工排量和砂比，以保证压裂施工过程不发生砂堵事故，同时该发明也适用于其他类型的高含泥塑性储层。
背景技术
砂堵是指加砂压裂过程砂粒聚集在射孔附近或靠近井筒的裂缝中，堵塞压裂液和砂粒的输运路径，造成泵压急剧上升、无法加砂等事故。砂堵对压裂施工造成重大影响，甚至导致压裂失败。影响砂堵的因素较复杂，主要包括地层因素、压裂液因素、工程因素、设计因素等，如何有效避免砂堵的产生一直是压裂设计和施工作业中感到棘手的难题。对超低渗致密储层的压裂来说，储层埋藏深，成岩作用强，且泥质含量高，导致岩石塑性强，断裂韧性高，裂缝启裂困难。若砂泥岩形成薄互层，导致裂缝形态进一步复杂，砂堵事故高发。
砂泥岩薄互层储层压裂涉及高含泥地层中的裂缝扩展。对于高含泥岩石，由于其应力应变关系是非线性的弹塑性，无法用解析方法表达其扩展压力。从压裂设计角度来说，优化排量和砂比是控制砂堵的有效手段，但对砂泥岩薄互层储层尚缺少定量的确定方法。本发明基于高含泥弹塑性地层的力学建模，通过有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力，进而确定施工排量；将计算砂泥岩地层裂缝宽度剖面看作平面应变问题，得到裂缝扩展形态，进而确定临界砂比。
发明内容
本发明的目的是提供砂泥岩薄互层储层压裂改造防砂堵工艺参数的确定方法。该方法将高含泥地层视为弹塑性非线性材料，通过力学建模和有限元数值计算确定砂泥岩薄互层储层裂缝扩展压力，由缝内净压力与排量的关系确定排量；将计算砂泥岩地层裂缝宽度剖面看作平面应变问题，得到裂缝扩展形态，用砂泥岩缝宽比来计算临界砂比。
本发明目的是这样实现的，将井眼周边的应力状态视为平面应变，根据对称性，建立水力裂缝的半边模型，如图1所示(裂缝端部局部放大如图2所示)。模型尺寸2m*2m，由于对称性，模型左侧固定，模型右侧施加最大有效地应力，上下施加最小有效地应力，井筒及压裂裂缝内为压裂液。为考虑裂缝扩展到不同长度时的缝宽分布情况，模型尺寸采用无量纲化形式，其中初始裂缝长度1000单位。裂缝扩展准则是：对砂岩来说是弹性材料，采用断裂韧性表征，对泥岩来说是塑性材料，采用抗拉强度表征。
分别计算砂岩及泥岩中裂缝扩展所需的缝内净压力。计算砂岩时，材料视为弹性，只需弹性模量及泊松比；计算泥岩时，材料视为非线性弹塑性，除弹性模量、泊松比，还需提供实验测定的应力应变曲线。计算结果显示泥岩扩展压力较砂岩高，说明砂泥岩薄互层储层的扩展压力由泥岩扩展压力决定。由于缝内压裂液的压力与体积、排量成比例，与均匀
Figure BSA00000623092000021
Figure BSA00000623092000022
考虑水力裂缝的长度较高度、宽度远远大，此时可以将计算裂缝宽度剖面看作平面应变问题，建立如图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)
计算所得的裂缝宽度剖面如图4所示，图中标识出砂岩、泥岩。裂缝形状呈“S”状，砂岩层裂缝往地层内凹，泥岩层裂缝往地层内凸起，导致泥岩层段裂缝宽度很窄，平均缝宽1.4mm。相同缝内压力条件下的砂岩对应平均缝宽3.6mm。砂泥岩互层中泥岩段缝宽较砂岩平均缝宽降低61％。泥岩层的窄缝宽是导致砂堵的机理所在。据此提出砂泥岩薄互层储层压裂防砂堵砂比上限可用砂泥岩缝宽比来计算，即
Figure BSA00000623092000023
附图说明
以下附图仅旨在于对本发明做示意性说明和解释，并不限定本发明的范围。其中，
图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。
4.确定排量：
Figure BSA00000623092000031
此处约为23％。
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。
7.
Figure BSA00000623092000032

Claims (2)

1.砂泥岩薄互层储层压裂施工排量的确定方法，其特征在于：高含泥地层的扩展压力较高，
维持裂缝扩展需提供高的缝内净压力，需要提高排量以维持较高缝内净压力，
Figure FSB0000115200540000011
Figure FSB0000115200540000012
2.砂泥岩薄互层储层压裂临界砂比的确定方法，其特征在于：砂泥岩薄互层中不可避免出现裂缝宽度分布的不均匀，泥岩层段裂缝宽度很窄，导致裂缝形态呈现类似“葫芦串”的形状，加砂过程中砂粒容易在喉结处聚集堵塞造成砂堵，据此提出砂比上限用砂泥岩缝宽比来计算，
Figure FSB0000115200540000013
CN201110385585.XA 2011-11-29 2011-11-29 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 Expired - Fee Related CN102608021B (zh)

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 CN102608021A (zh) 2012-07-25
CN102608021B true 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)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105373648B (zh) * 2015-10-09 2018-10-12 西南石油大学 一种砂泥岩互层型沉积体储层构型的建模方法
CN105550441B (zh) * 2015-12-12 2018-05-18 山东科技大学 一种基于连续介质的工程岩体破裂劣化数值模拟方法
CN106761647B (zh) * 2017-01-13 2020-08-14 中国石油化工股份有限公司 一种估算页岩储层压后平面改造面积的方法
CN106951659A (zh) * 2017-04-05 2017-07-14 中国石油大学(北京) 一种高效实现煤层及煤层顶部砂岩合压的射孔方案
CN111396014B (zh) * 2020-03-16 2021-07-23 中国石油大学(北京) 薄互层储层改造方法、装置和设备
CN111368463A (zh) * 2020-04-08 2020-07-03 广州海洋地质调查局 一种水平井穿层压裂设计方法

Citations (3)

* Cited by examiner, † Cited by third party
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 中国石油天然气股份有限公司 一种检测灰岩地层孔隙压力的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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 (2)

* Cited by examiner, † Cited by third party
Title
砾石充填防砂工艺参数优化设计;董长银等;《中国石油大学学报(自然科学版)》;20061031;第30卷(第5期);57-61 *
董长银等.砾石充填防砂工艺参数优化设计.《中国石油大学学报(自然科学版)》.2006,第30卷(第5期),57-61.

Also Published As

Publication number Publication date
CN102608021A (zh) 2012-07-25

Similar Documents

Publication Publication Date Title
CN102608021B (zh) 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法
Guo et al. Investigation of production-induced stress changes for infill-well stimulation in Eagle Ford Shale
CN108319756B (zh) 一种致密储层体积压裂缝网扩展模拟及表征方法
CN105019877B (zh) 一种页岩水平井压裂破碎区体积的计算方法
Lian et al. A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells
CN109374867B (zh) 一种基于离散元的砂砾岩水力压裂模拟方法
CN109025942A (zh) 一种致密气藏斜井压裂不规则多裂缝的产量计算方法
Zheng et al. Modeling fracture closure with proppant settling and embedment during shut-in and production
CN105277671B (zh) 一种确定页岩地层脆性指数的方法
Wu et al. Influence of rock heterogeneity on hydraulic fracturing: A parametric study using the combined finite-discrete element method
CN101806224A (zh) 井下煤层气抽采产能预测方法
Li et al. Stress redistribution and fracture propagation during restimulation of gas shale reservoirs
Zheng et al. Numerical investigation on the hydraulic fracture propagation based on combined finite-discrete element method
WO2017031857A1 (zh) 一种构建双溶腔盐穴储库地面沉降预测模型的装置
Wang et al. A 3-D poro-elasto-plastic model for sand production around open-hole and cased & perforated wellbores
JP2017025654A (ja) 掘削壁面の安定性評価方法
Jeng et al. Three-dimensional consolidation of a porous unsaturated seabed under rubble mound breakwater
Yu et al. Impact of asymmetric stimulated rock volume on casing deformation in multi-stage fracturing; a case study
CN116717224A (zh) 一种针对低渗致密储层复杂缝网的压裂产能预测方法
CN116562189A (zh) 一种针对动态裂缝漏失的堵漏颗粒材料的优选方法、系统及存储介质
Zhou et al. Rock dilation and its effect on fracture transmissivity
CN112487594B (zh) 油藏水体倍数计算方法及装置
CN108959764A (zh) 低渗油藏水驱扫油面积系数评价方法及系统
Molina et al. Analysis of Well-Completion performance and production optimization of a gas well using computational fluid dynamics
Weidong et al. Advances in research of hydraulic fractures in unconsolidated sands

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