CN106295235B - 一种页岩地层天然裂缝剪切滑移量的计算方法 - Google Patents
一种页岩地层天然裂缝剪切滑移量的计算方法 Download PDFInfo
- Publication number
- CN106295235B CN106295235B CN201610806787.XA CN201610806787A CN106295235B CN 106295235 B CN106295235 B CN 106295235B CN 201610806787 A CN201610806787 A CN 201610806787A CN 106295235 B CN106295235 B CN 106295235B
- Authority
- CN
- China
- Prior art keywords
- shear
- fracture
- stress
- mpa
- natural
- 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.)
- Active
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种页岩地层天然裂缝剪切滑移量的计算方法，依次包括以下步骤：(A)计算天然裂缝所受的正应力σn和剪应力τ；(B)基于裂缝中流体的压力以及步骤(A)的计算结果，判断天然裂缝的开启方式；(C)将缝长为a的天然裂缝划分为长度相等的N个单元体，每个单元体长度为(D)基于步骤(B)的判断结果和步骤(C)的单元体划分方式，采用位移不连续方法计算天然裂缝的剪切滑移量。本发明根据天然裂缝的基础参数以及其受力状态判断天然裂缝的开启方式，在给定的天然裂缝长度、角度、摩擦系数以及页岩地层岩石力学参数的前提下，采用位移不连续方法计算天然裂缝的剪切滑移量，原理可靠，操作简单，为页岩地层水力压裂施工参数优化提供了重要的指导依据。
Description
技术领域
本发明涉及石油工程领域水力压裂过程中一种页岩地层天然裂缝剪切滑移量的计算方法。
背景技术
水力压裂技术是油气藏增产改造的重要措施。水力压裂是利用地面高压泵组，以超过地层吸收能力的排量将压裂液泵入地层来产生裂缝，然后继续注入带有支撑剂(砂粒)的压裂液，使裂缝继续延伸并在其中充填支撑剂，当压裂液返排后，在地层压力作用下，支撑剂在裂缝中起到支撑裂缝的作用，阻止裂缝闭合，从而在地层中形成具有一定长度、允许流体流动的填砂裂缝。
清水压裂是水力压裂的一种形式，被广泛应用于页岩油气藏的增产改造中。它的特点是不加入支撑剂，仅通过泵入低粘度压裂液，在页岩地层中产生人工裂缝，并沟通天然裂缝，以形成复杂的裂缝网络体系，提高流体的流动能力。通常情况下，页岩地层中原有的天然裂缝处于闭合状态，无法为流体提供流动通道。在清水压裂过程中，当人工裂缝与天然裂缝相交时，会促使天然裂缝产生剪切滑移，同时由于裂缝表面凹凸不平，即使天然裂缝中未充填支撑剂，裂缝表面凸起也可以相互支撑，使天然裂缝在闭合压力的作用下保持一定的开启程度，以此为流体提供流动通道，达到改善油气流动条件和油气井增产的目的。
裂缝剪切滑移量是指裂缝两个表面相对位移的大小，为了使裂缝表面的凸起相互支撑，天然裂缝必须具备一定的剪切滑移量，否则，裂缝的两个粗糙表面在闭合压力的作用下将会完全啮合而无法为流体提供流动通道。现有的研究结果表明，天然裂缝的剪切滑移量影响了水力压裂的增产效果，因此，准确地计算剪切滑移量对压裂施工参数优化具有重要的指导意义。
天然裂缝的剪切滑移量与裂缝的长度、角度以及裂缝面的摩擦系数有关。裂缝角度是指沿顺时针方向，裂缝与水平最大主应力的夹角。摩擦系数指裂缝滑动时，摩擦力与裂缝面所受的正压力的比值。剪切滑移量计算的目的即是在给定的裂缝长度、角度以及裂缝面摩擦系数的基础上，计算天然裂缝不同位置处的剪切滑移量。
位移不连续方法是岩体工程研究中常用的数值计算方法，由于其直接以裂缝面上的相对位移为未知量，因此在求解裂缝面受力大小和位移量时更加简单、方便。
目前国内外已有的研究成果大多集中于天然裂缝剪切滑移机理的研究，而对剪切滑移量的研究仍停留在定性的认识上，缺乏定量化的计算方法，难以满足生产的需求。
发明内容
本发明的目的在于提供一种页岩地层天然裂缝剪切滑移量的计算方法，该方法原理可靠，操作简单，能够在给定天然裂缝和岩石基础参数的前提下，计算裂缝不同位置处的剪切滑移量。
为达到以上技术目的，本发明提供以下技术方案。
本发明根据天然裂缝的基础参数以及其受力状态判断天然裂缝的开启方式，以此为基础，在给定的天然裂缝长度、角度、摩擦系数以及页岩地层岩石力学参数的前提下，采用位移不连续方法计算天然裂缝的剪切滑移量。
一种页岩地层天然裂缝剪切滑移量的计算方法，依次包括以下步骤：
(A)计算天然裂缝所受的正应力σn和剪应力τ：
式中：σn—裂缝壁面所受的正应力，MPa；
τ—裂缝所受的剪应力，MPa；
σH—水平最大主应力，MPa；
σh—水平最小主应力，MPa；
θ—裂缝的角度，rad。
(B)基于裂缝中流体的压力以及步骤(A)的计算结果，判断天然裂缝的开启方式：
(1)天然裂缝法向张开，判断依据为：此时，天然裂缝所受的剪应力小于裂缝的剪切强度，但缝内流体压力大于裂缝所受的正应力，裂缝发生张性破坏。
(2)天然裂缝剪切破坏，判断依据为：此时，天然裂缝所受的剪应力大于裂缝的剪切强度，裂缝发生剪切破坏。
其中：σn—裂缝壁面所受的正应力，MPa；
P—裂缝内流体压力，MPa；
τ0—岩石内聚力，MPa；
Kf—裂缝面的摩擦系数，无因次；
τ—裂缝所受的剪应力，MPa。
(C)将缝长为a的天然裂缝划分为长度相等的N个单元体，每个单元体长度为
(D)基于步骤(B)的判断结果和步骤(C)的单元体划分方式，采用位移不连续方法计算天然裂缝的剪切滑移量：
(1)当天然裂缝法向张开时，任意一个单元体i所受应力与其剪切滑移量满足以下方程：
(2)当天然裂缝剪切破坏时，任意一个单元体i所受应力与其剪切滑移量满足以下方程：
其中，为应力边界影响系数，其表达式为：
式中：σn—裂缝壁面所受的正应力，MPa；
P—裂缝内流体压力，MPa；
τ—裂缝所受的剪应力，MPa；
—单元体i所受的正应力，MPa；
—单元体i所受的剪应力，MPa；
—第j个单元体的法向位移量，m；
—第j个单元体的剪切滑移量，m；
Di n—第i个单元体的法向位移量，m；
N—单元体的个数；
G—岩石的剪切应变模量，MPa；
E—岩石的杨氏模量，MPa；
μ—泊松比，无因次；
a—天然裂缝长度，m；
i—第i个单元体；
j—第j个单元体。
联立求解N个单元体的应力—滑移量方程组，即可得到每个单元体的剪切滑移量，最终得出天然裂缝的剪切滑移量。
所述步骤(A)中，采用了二维线弹性理论计算裂缝面所受的正应力σn和剪应力τ(参考文献：周健，陈勉，金衍，等.压裂中天然裂缝剪切破坏机制研究[J].岩石力学与工程学报，2008，27：2637～2641)。
对步骤(B)中天然裂缝开启方式的判断依据做如下说明：
当天然裂缝所受的剪应力超过剪切强度时，裂缝发生剪切破坏而产生剪切裂缝，此时，裂缝内流体压力满足：
τ＞τ0+Kf(σn-P)
可得：
当天然裂缝发生张性破坏时，裂缝内流体压力需大于裂缝所受的正压力，以支撑裂缝壁面，防止裂缝闭合。同时，裂缝所受的剪应力应小于裂缝的剪切强度，因此，裂缝内流体压力满足：
式中：σn—裂缝壁面所受的正应力，MPa；
P—裂缝内流体压力，MPa；
τ0—岩石内聚力，MPa；
Kf—裂缝面的摩擦系数，无因次；
τ—裂缝所受的剪应力，MPa。
所述步骤(D)中，天然裂缝剪切滑移量的推导过程如下：
(1)由Green函数可得天然裂缝边界上任意一点的位移为：
(2)在步骤(1)的基础上，根据弹性力学理论，任一单元体i的剪应力和法应力可以由单元体j的位移不连续量计算得到：
(3)根据天然裂缝的破坏形式，确定裂缝中每一个单元体的边界条件。
对于法向张开的天然裂缝：
对于剪切破坏的天然裂缝：
(4)将边界条件带入步骤(2)中，即可得到各单元体的剪切滑移量。
式中：—天然裂缝边界上任意一点的位移，m；
A—与该点位置有关的系数；
σi s—第i个单元体所受的剪应力，MPa；
σi n—第i个单元体所受的正应力，MPa；
N—单元体的个数；
—第j个单元体的法向位移量，m；
—第j个单元体的剪切滑移量，m；
Di n—第i个单元体的法向位移量，m；
σn—裂缝壁面所受的正应力，MPa；
P—裂缝内流体压力，MPa；
τ—裂缝所受的剪应力，MPa；
为应力边界影响系数，当天然裂缝未发生弯曲时，均等于零。
与现有技术相比，本发明的有益效果是：采用该方法可以较为精确地计算页岩地层天然裂缝的剪切滑移量，该方法计算方式简单，为页岩地层水力压裂施工参数优化提供了重要的指导依据。
附图说明
图1是天然裂缝剪切滑移量的计算结果。
具体实施方式
下面以某页岩地层实际参数为例，结合附图，对本发明的步骤进行详细说明。
页岩地层各项参数表
一种页岩地层天然裂缝剪切滑移量的计算方法，依次包括以下步骤：
(A)计算天然裂缝所受的正应力σn和剪应力τ：
式中：σn—裂缝壁面所受的正应力，MPa；
τ—裂缝所受的剪应力，MPa；
σH—水平最大主应力，MPa；
σh—水平最小主应力，MPa；
θ—裂缝的角度，rad。
各项参数的计算结果
(B)基于裂缝中流体的压力以及步骤(A)的计算结果，可判断该天然裂缝满足：因此，该裂缝为剪切破坏。
(C)将该天然裂缝划分为长度相等的150个单元体，每个单元体长度为0.1m。
(D)基于步骤(B)的判断结果和步骤(C)的单元体划分方式，采用位移不连续方法计算天然裂缝的剪切滑移量：
其中，为应力边界影响系数，其表达式为：
式中：σn—裂缝壁面所受的正应力，MPa；
P—裂缝内流体压力，MPa；
τ—裂缝所受的剪应力，MPa；
—单元体i所受的正应力，MPa；
—单元体i所受的剪应力，MPa；
—第j个单元体的法向位移量，m；
—第j个单元体的剪切滑移量，m；
Di n—第i个单元体的法向位移量，m；
N—单元体的个数；
G—岩石的剪切应变模量，MPa；
E—岩石的杨氏模量，MPa；
μ—泊松比，无因次；
a—天然裂缝长度，m；
i—第i个单元体；
j—第j个单元体。
联立求解150个单元体的应力—滑移量方程，即可得到每个单元体的剪切滑移量，最终得出天然裂缝的剪切滑移量(图1)。
Claims (2)
1.一种页岩地层天然裂缝剪切滑移量的计算方法，依次包括以下步骤：
(A)计算天然裂缝壁面所受的正应力σn和剪应力τ；
(B)基于裂缝内流体压力以及步骤(A)的计算结果，判断天然裂缝的开启方式：
时，天然裂缝法向张开；
时，天然裂缝剪切破坏；
其中：σn—天然裂缝壁面所受的正应力，MPa，
P—裂缝内流体压力，MPa，
τ0—岩石内聚力，MPa，
Kf—裂缝面的摩擦系数，无因次，
τ—裂缝所受的剪应力，MPa；
(C)将天然裂缝长度为a的天然裂缝划分为长度相等的N个单元体，每个单元体长度为
(D)基于步骤(B)的判断结果和步骤(C)的单元体划分方式，采用位移不连续方法计算天然裂缝的剪切滑移量：
当天然裂缝法向张开时，任意一个单元体i所受应力与其剪切滑移量满足以下方程：
其中i＝1，2，3…N；j＝1，2，3…N
当天然裂缝剪切破坏时，任意一个单元体i所受应力与其剪切滑移量满足以下方程：
其中i＝1，2，3…N；j＝1，2，3…N
其中，为应力边界影响系数，其表达式为：
式中：σn—天然裂缝壁面所受的正应力，MPa，
P—裂缝内流体压力，MPa，
τ—裂缝所受的剪应力，MPa，
—单元体i所受的正应力，MPa，
—单元体i所受的剪应力，MPa，
—第j个单元体的法向位移量，m，
—第j个单元体的剪切滑移量，m，
Di n—第i个单元体的法向位移量，m，
N—单元体的个数，
G—岩石的剪切应变模量，MPa，
E—岩石的杨氏模量，MPa，
μ—泊松比，无因次，
a—天然裂缝长度，m，
i—第i个单元体，
j—第j个单元体；
联立求解N个单元体的应力—滑移量方程组，即可得到每个单元体的剪切滑移量，最终得出天然裂缝的剪切滑移量。
2.如权利要求1所述的一种页岩地层天然裂缝剪切滑移量的计算方法，其特征在于，所述步骤(A)计算天然裂缝壁面所受的正应力σn和剪应力τ，过程如下：
式中：σH—水平最大主应力，MPa，
σh—水平最小主应力，MPa，
θ—裂缝的角度，rad。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610806787.XA CN106295235B (zh) | 2016-09-07 | 2016-09-07 | 一种页岩地层天然裂缝剪切滑移量的计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610806787.XA CN106295235B (zh) | 2016-09-07 | 2016-09-07 | 一种页岩地层天然裂缝剪切滑移量的计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106295235A CN106295235A (zh) | 2017-01-04 |
| CN106295235B true CN106295235B (zh) | 2018-10-30 |
Family
ID=57709811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610806787.XA Active CN106295235B (zh) | 2016-09-07 | 2016-09-07 | 一种页岩地层天然裂缝剪切滑移量的计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106295235B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110872943B (zh) * | 2018-08-30 | 2021-07-13 | 中国石油化工股份有限公司 | 一种确定地层射孔方向的方法 |
| CN110017135B (zh) * | 2019-02-15 | 2022-05-20 | 西南石油大学 | 一种裂缝性地层井壁裂缝扩展压力预测方法 |
| CN110593842B (zh) * | 2019-10-22 | 2020-05-05 | 西南石油大学 | 实验确定页岩储层水力压裂自支撑裂缝导流能力的方法 |
| CN112033812B (zh) * | 2020-08-12 | 2022-11-08 | 成都北方石油勘探开发技术有限公司 | 一种水力剪切压裂剪胀导流能力的测试方法及系统 |
| CN114662316B (zh) * | 2022-03-25 | 2026-01-09 | 中海石油（中国）有限公司天津分公司 | 一种计算空间三维天然裂缝壁面内摩擦系数的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1916359A (zh) * | 2005-11-28 | 2007-02-21 | 长庆石油勘探局 | 一种实现重复压裂造新缝的方法 |
| CN105114065A (zh) * | 2015-08-06 | 2015-12-02 | 中国石油天然气股份有限公司 | 水力压裂裂缝延伸过程的模拟方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10267131B2 (en) * | 2012-08-13 | 2019-04-23 | Schlumberger Technology Corporation | Competition between transverse and axial hydraulic fractures in horizontal well |
-
2016
- 2016-09-07 CN CN201610806787.XA patent/CN106295235B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1916359A (zh) * | 2005-11-28 | 2007-02-21 | 长庆石油勘探局 | 一种实现重复压裂造新缝的方法 |
| CN105114065A (zh) * | 2015-08-06 | 2015-12-02 | 中国石油天然气股份有限公司 | 水力压裂裂缝延伸过程的模拟方法 |
Non-Patent Citations (1)
| Title |
|---|
| 压裂中天然裂缝剪切破坏机制研究;周健等;《岩石力学与工程学报》;20080630;第27卷;第2637-2641页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106295235A (zh) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106295235B (zh) | 一种页岩地层天然裂缝剪切滑移量的计算方法 | |
| Frash et al. | Laboratory simulation of binary and triple well EGS in large granite blocks using AE events for drilling guidance | |
| Lian et al. | A study on casing deformation failure during multi-stage hydraulic fracturing for the stimulated reservoir volume of horizontal shale wells | |
| CN110017135B (zh) | 一种裂缝性地层井壁裂缝扩展压力预测方法 | |
| CN107609258B (zh) | 一种页岩重复压裂转向裂缝起裂压力的计算方法 | |
| CN109025942A (zh) | 一种致密气藏斜井压裂不规则多裂缝的产量计算方法 | |
| CN109374867A (zh) | 一种基于离散元的砂砾岩水力压裂模拟方法 | |
| CN103590824A (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
| CN107437127B (zh) | 一种油井停喷地层压力预测方法 | |
| CN106649963B (zh) | 体积压裂复杂缝网平均裂缝长度和等效裂缝条数确定方法 | |
| CN103161435A (zh) | 一种稠油热采直井试井解释方法 | |
| CN108829945A (zh) | 一种基于压裂施工压力实时评估增产改造体积方法 | |
| CN107480383A (zh) | 一种通过测压数据监测注水动态裂缝的方法 | |
| CN105257278B (zh) | 一种支撑剂嵌入深度的获取方法 | |
| CN110348031A (zh) | 水平井压裂近井筒裂缝扭曲形态数值模拟方法 | |
| CN111222281A (zh) | 一种气藏型储气库注采管柱冲蚀失效风险测定方法 | |
| CN104963675B (zh) | 一种预测水平井暂堵酸化酸液推进距离的方法 | |
| CN106295234A (zh) | 一种页岩地层水力压裂裂缝剪切滑移量的计算方法 | |
| CN116988770A (zh) | 一种油田分层注水方法、装置、电子设备及存储介质 | |
| Frash et al. | Laboratory-scale-model testing of well stimulation by use of mechanical-impulse hydraulic fracturing | |
| CN105317430B (zh) | 测定各向异性地层水平井井壁围岩周向应力的方法 | |
| CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
| CN111914494B (zh) | 一种致密储层水平井压后压裂液返排率预测方法及系统 | |
| CN106442173A (zh) | 一种闭合应力作用下剪切裂缝开度的预测方法 | |
| CN113901678B (zh) | 一种工程注浆用水泥浆凝结质量评价装置及评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |