CN112177598B - 一种考虑压裂液压缩性的地层起裂压力预测方法 - Google Patents
一种考虑压裂液压缩性的地层起裂压力预测方法 Download PDFInfo
- Publication number
- CN112177598B CN112177598B CN201910584484.1A CN201910584484A CN112177598B CN 112177598 B CN112177598 B CN 112177598B CN 201910584484 A CN201910584484 A CN 201910584484A CN 112177598 B CN112177598 B CN 112177598B
- Authority
- CN
- China
- Prior art keywords
- pressure
- shaft
- time
- fracturing fluid
- wellhead
- 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
- 238000000034 method Methods 0.000 title claims abstract description 87
- 239000012530 fluid Substances 0.000 title claims abstract description 66
- 238000005336 cracking Methods 0.000 title claims description 57
- 239000004568 cement Substances 0.000 claims abstract description 38
- 230000008859 change Effects 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 32
- 239000011435 rock Substances 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 22
- 238000010168 coupling process Methods 0.000 claims abstract description 22
- 238000005859 coupling reaction Methods 0.000 claims abstract description 22
- 239000007787 solid Substances 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000010276 construction Methods 0.000 claims description 32
- 239000011148 porous material Substances 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 32
- 230000000977 initiatory effect Effects 0.000 abstract description 20
- 238000013461 design Methods 0.000 abstract description 9
- 238000011161 development Methods 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000006073 displacement reaction Methods 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
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)
- Geophysics (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明提供了一种考虑压裂液压缩性的地层起裂压力预测方法，属于页岩气开发技术领域。该方法包括：S1：建立井筒压裂液流动的离散化井筒模型，获得井底压力随时间变化的规律和井口压力随时间变化的规律；S2：建立射孔条件下的井筒‑水泥环‑近井壁地层三维流固耦合力学模型，获得井筒憋压过程中近井壁应力场的分布规律和渗流场的分布规律；S3：反演地层起裂时间及井口起裂压力。本发明建立了憋压过程中井筒压力场及速度场的计算模型，在此基础上，结合岩石力学等相关知识，完善了页岩气井水力压裂起裂压力及起裂时间的计算方法，在一定程度上为深入认识页岩气井水力压裂的起裂机理及开展压裂设计提供了进一步的理论基础及设计依据。
Description
技术领域
本发明属于页岩气开发技术领域，具体涉及一种考虑压裂液压缩性的地层起裂压力预测方法。
背景技术
我国页岩气资源丰富。经初步评价，我国陆域页岩气地质资源潜力为134.42万亿立方米，可采资源潜力为25.08万亿立方米(不含青藏区)。但是页岩基质渗透率低，低孔低渗，需通过水力压裂作业才可以获得工业性气流。目前，我国常规压裂方式为水平分段多簇射孔压裂，该方法由于操作简单，工艺相对比较成熟，目前已广泛应用于我国页岩气资源开发中。
但是，在深层、超深层页岩气储层水力压裂作业中发现，部分地区由于埋藏较深，岩石致密，闭合压力高，导致破裂压力极高，施工极为困难，给页岩气储层改造带来了极大的困难。因此，需要从页岩气储层的物理及力学性质出发，探讨高破裂压力的形成原因及降低手段。页岩气井水力压裂裂缝起裂的预测模型，显得尤为重要。
针对此，不少学者对起裂压力的预测模型展开了研究。研究方法主要基于解析法、有限元法、边界元法等，例如，中国专利公开文献CN107609258B公开了一种页岩重复压裂转向裂缝起裂压力的计算方法，其技术方案与本发明适用范围存在区别；中国专利公开文献CN106555575A公开了一种深层页岩的地层起裂压力的计算方法，其在计算中未考虑流体压裂液、粘度、降阻率等对起裂压力的影响，在预测起裂时间及起裂压力上具有一定的技术局限性；中国专利公开文献CN107587867B公开了一种提升页岩缝网复杂度的重复压裂工艺设计方法，其技术方案与本发明适用范围存在区别。
另外，大部分模型均将孔眼及井筒内压力的加载方式简化为恒定压力，实际上，在压裂初期，井筒内压力并不是立刻升高，而是逐步上升，其增压速率受压裂液排量、粘度及压缩性影响，并直接影响页岩气井的起裂压力及起裂时间。因此，现有的页岩气井起裂压力预测模型难以直接考虑排量、粘度等压裂液参数对起裂压力的影响，故在预测起裂压力及时间上具有一定的技术局限性，有必要建立一个考虑井筒压裂液压缩性的起裂压力计算模型。
发明内容
本发明的目的在于解决上述现有技术中存在的难题，提供一种考虑压裂液压缩性的地层起裂压力预测方法，实现对页岩气井水力压裂的起裂压力及起裂时间的预测。
本发明是通过以下技术方案实现的：
本发明方法包括以下步骤：
S1：建立井筒压裂液流动的离散化井筒模型，获得井底压力随时间变化的规律和井口压力随时间变化的规律；
S2：建立射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型，获得井筒憋压过程中近井壁应力场的分布规律和渗流场的分布规律；
S3：反演地层起裂时间及井口起裂压力。
所述步骤S1中的建立井筒压裂液流动的离散化井筒模型的操作包括：
沿套管轴向将套管等距离划分为n个井筒微元段，每个井筒微元段长度为Δzi；第一个井筒微元段对应井口，第n个井筒微元段对应井底；
建立井筒压裂液流动的离散化井筒模型如下：
其中，为t时刻第i个井筒微元段的压裂液质量，采用下式计算得到：
P为压力，v为平均流速，表示第i个井筒微元段和第i+1个井筒微元段相交处截面，/>表示第i个井筒微元段和第i-1个井筒微元段相交处截面，Δt为单位时间步长，α为井筒微元段与垂向的夹角，d为套管内壁直径，g为重力加速度，ρ0为标准大气压下的压裂液密度，p0为标准大气压，C为压裂液压缩系数；τw为井筒壁面摩擦阻力。
所述步骤S1中的获得井底压力随时间变化的规律和井口压力随时间变化的规律的操作包括：
(1)在单位时间步长下求解所述井筒压裂液流动的离散化井筒模型，获得该单位时间步长下的井底压力和井口压力；所述井底压力对应求解得到的第n个井筒微元段在该单位时间步长下的压力pn，所述井口压力对应求解得到的第1个井筒微元段在该单位时间步长下的压力p1；
(2)判断pn是否大于等于最大水平主应力的Y倍，如果否，则
t＝t+Δt,然后返回步骤(1)，如果是，则转入步骤(3)；
(3)输出井底压力随时间变化的规律和井口压力随时间变化的规律：所述井底压力随时间变化的规律对应pn在各个单位时间步长下的数值，所述井口压力随时间变化的规律对应p1在各个单位时间步长下的数值。
所述步骤(1)中是采用牛顿迭代法在各个单位时间步长下求解所述井筒压裂液流动的离散化井筒模型；
所述步骤(2)中的Y为1.5。
所述步骤S2的操作包括：
步骤201：获取套管、水泥环和井眼尺寸，结合射孔施工参数，利用有限元法建立射孔条件下的套管-水泥环-近井壁地层三维流固耦合力学模型；
步骤202：对所述射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型中的套管、水泥环赋予材料属性；对所述射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型中的页岩赋予物理力学属性；
步骤203：对所述射孔条件下的套管-水泥环-近井壁地层三维流固耦合力学模型施加远场边界载荷；
步骤204：对所述射孔条件下套管-水泥环-近井壁地层三维流固耦合力学模型的套管内壁及孔眼处施加液压载荷，所述液压载荷的大小为步骤S1中得到的井底压力；
步骤205：计算得到憋压过程中近井壁页岩的应力场的分布规律及渗流场的分布规律。
所述步骤201中采用三维六面体缩减积分孔压单元对所述套管-水泥环-近井壁地层三维力学模型进行网格划分，并在孔眼处对网格进行加密。
所述步骤202中的物理力学属性包括：渗透率、孔隙度、弹性模量、泊松比。
所述步骤S3的操作包括：
步骤301：输入步骤S1获得的井底压力随时间变化的规律；根据步骤S2得到的应力场的分布规律计算得到不同施工时间下的近井壁岩石的最大主应力；
步骤302：根据所述不同施工时间下的近井壁岩石的最大主应力反演地层起裂时间；
步骤303：根据所述地层起裂时间和步骤S1得到的井口压力随时间变化的规律反演井口起裂压力。
所述步骤302的操作包括：
以所述不同施工时间下的近井壁岩石的最大主应力作为插值的基础数据，对该基础数据进行插值，得到井壁页岩的最大主应力达到岩石抗拉强度时的施工时间，该施工时间即为地层起裂时间。
所述步骤303的操作包括：
将步骤S1得到的井口压力随时间变化的规律中的各个井口压力的值作为插值的基础数据，对该基础数据进行插值，然后计算施工时间为地层起裂时间时的井口压力的大小，该井口压力即为井口起裂压力。
与现有技术相比，本发明的有益效果是：本发明基于流体力学相关知识，建立了憋压过程中井筒压力场及速度场的计算模型，在此基础上，结合岩石力学等相关知识，完善了页岩气井水力压裂起裂压力及起裂时间的计算方法，在一定程度上为深入认识页岩气井水力压裂的起裂机理及开展压裂设计提供了进一步的理论基础及设计依据。
附图说明
图1为本发明的基本流程图；
图2为本发明的计算憋压过程中井筒压力场的基本流程图；
图3为本发明的计算近井壁应力场及渗流场的基本流程图；
图4为本发明的计算起裂压力及时间的基本流程图；
图5为实施例中的井口压力计算结果和实测值的对比；
图6为射孔条件下套管-水泥环-地层三维有限元模型；
图7为近井壁页岩最大主应力随憋压时间的变化曲线图。
具体实施方式
下面结合附图对本发明作进一步详细描述：
如图1到图4所示，本发明方法包括：
S1：建立井筒压裂液流动的离散化井筒模型，获得井底压力随时间变化的规律和井口压力随时间变化的规律：获取施工排量、开井压力，结合压裂液降阻率、粘度、压缩性等参数，建立憋压过程中井筒压裂液流动的数学模型，求解井底及井口压力随时间的变化规律；
S2：建立射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型，获得井筒憋压过程中近井壁应力场的分布规律和渗流场的分布规律：基于渗流力学及岩石力学相关知识，结合射孔参数及井筒、水泥环、地层的力学参数，建立井筒-水泥环-近井壁地层三维流固耦合力学模型，计算井筒憋压过程中近井壁应力场及渗流场的分布规律；
S3：反演地层起裂时间及井口起裂压力：结合计算结果，分析页岩最大主应力在井筒憋压过程中的变化规律，结合最大水平主应力准则，反演地层起裂时间及井口起裂压力。
优选地，所述步骤S1的建立憋压过程中井筒压裂液流动的数学模型，具体包括如下内容：
步骤101：计算参数准备：获取压裂液的压缩系数、粘度、密度、降阻率以及施工排量、开井压力、井眼轨迹、射孔位置；
步骤102：进行数值模拟：基于流体力学相关知识，建立压裂憋压过程中井筒流体流动物理模型，采用数值模拟方法进行憋压过程中井筒压力场、速度场的数值模拟，各个时间步内井筒中压裂液的流体压力分布采用牛顿迭代法求解，在此基础上，求解流速场分布，并进入下一时间步计算，直至井底压力达到最大水平主应力的1.5倍以上；具体如下：
采用数学模型计算憋压过程中井筒流体压力场，憋压过程中井筒压裂液流动的数学模型为公式(1)到(4)：
流体压缩性方程：
式中，C为压裂液压缩系数；ρ为压裂液密度；p为开井压力(可通过邻井获得)。
可化为：
式中：ρ0为标准大气压下的压裂液密度，p0为标准大气压；
质量守恒方程：
式中，z为管长方向位置；u为压裂液在管内的流速。
井筒压降方程：
式中，d为套管内壁直径，α为井筒单元与垂向的夹角，g为重力加速度，f为摩擦系数，β为降阻率。
所述的井筒流体流动模型需要预先通过以下方式以建立井筒压裂液流动的离散化井筒模型(该模型是指公式(5)和(6))：
沿套管轴向将套管等距离划分为n个井筒微元段，每个井筒微元段长度为Δzi。
以单个井筒微元段为研究对象，则质量守恒方程可离散为：
式中，为t时刻第i个井筒微元段的压裂液质量；
根据式(2),压裂液密度为关于压力的函数，因此公式(5)中包含了压力。为第i个井筒微元段和第i+1个井筒微元段相交处截面上的平均流速，/>表示第i个井筒微元段和第i-1个井筒微元段相交处截面上的平均流速，Δt为单位时间步长。
同时，井筒压降方程可离散为：
其中，τw表示井筒壁面摩擦阻力。
值得注意的是，由于本模型考虑了流体的压缩性，根据流体压缩性方程，井筒微元段的压裂液密度为：
进一步，需要定义憋压过程中井筒压裂液流动模型的边界条件。首先，定义开始压裂时的井口压力，然后，定义憋压过程中井口的注入排量，该排量可以为定值或关于时间的函数。
进一步，将井筒共划分为n个井筒微元段，每个井筒微元段包括流速及压力两个未知数，即共2n个未知数。每个井筒微元段可以建立连续性方程及井筒压降方程2个方程(即公式(5)和(6))，共可建立2n个方程。因此，方程数个数与未知数个数相同，可以求解。鉴于连续性方程和井筒压降方程均为全隐式的非线性方程，所以采用牛顿迭代法求解。
进一步，计算憋压过程中不同时间下井底压力(对应最后一个井筒微元段的压力值pn)的大小，当井底压力的大小达到最大水平主应力(可通过测井或岩心声发射实验获取)的1.5倍以上，则停止计算。
步骤103：输出计算结果：根据需求，输出计算起裂压力及起裂时间所需要的参数，即：井底压力、井口压力随时间变化的规律，步骤102求解2n个方程得到的是n个井筒微元段在不同时间下的流速及压力，所述井底压力随时间变化的规律对应求解得到的最后一个井筒微元段(即第n个井筒微元段)的压力pn在各个单位时间步长下的数值，所述井口压力随时间变化的规律对应求解得到的第一个井筒微元段的压力p1在各个单位时间步长下的数值。
所述步骤S2的建立套管-水泥环-近井壁地层三维流固耦合力学模型，具体包括如下内容：
步骤201：获取套管、水泥环和井眼尺寸，结合射孔施工参数，基于有限元法建立射孔条件下套管-水泥环-近井壁地层三维流固耦合力学模型；模型划分网格为三维六面体缩减积分孔压单元，为保证计算结果的准确性，在孔眼处可局部加密网格；
步骤202：对模型中的套管及水泥环赋予材料属性：输入套管、水泥环的材料参数；对模型中的页岩赋予物理力学属性：通过测井资料及室内岩心实验获得页岩的物理及力学性质，输入页岩的物理、力学参数；所述页岩的物理及力学参数包括但不限于渗透率、孔隙度、弹性模量、泊松比等参数；页岩的力学本构模型(即弹塑性力学方程)可以根据页岩的力学性质调整，包括但不限于横观各向同性弹性力学模型、基于D-P准则的弹塑性力学模型；
步骤203：对射孔条件下套管-水泥环-近井壁地层三维流固耦合力学模型施加远场边界载荷：在垂直方向上注入垂直主应力，在最大水平主应力方向输入最大水平主应力,在最小水平主应力方向输入最小水平主应力，在岩石外边界施加孔隙压力边界；页岩的边界条件主要包括地层初始应力场及渗流场，在模型外边界施加位移约束条件及孔隙压力边界条件。
步骤204：对射孔条件下套管-水泥环-近井壁地层三维力学模型的套管内壁及孔眼处施加液压载荷，液压载荷的大小为步骤103中的输出结果，即与施工时间相关的井底压力变化值；
步骤205：计算憋压过程中近井壁页岩的应力场及渗流场分布：可以通过离散元法、有限元法或边界元法计算。以有限元法为例，基于岩石力学平衡力学方程、渗流力学方程建立页岩储层流固耦合力学方程，并对力学方程进行离散，得到耦合方程的有限元计算格式，并通过数值方法求解。
优选地，所述步骤S3的计算起裂压力及起裂时间，具体包括如下内容：
步骤301：在步骤103和步骤205的基础上，分析及确定憋压过程中不同施工时间下井底压力(即步骤103得到的井底压力随时间变化的规律)和近井壁岩石的最大主应力。步骤205得到的是近井壁页岩的应力场，通过步骤205计算出的应力场即可采用常规的力学公式计算出最大主应力。
步骤302：基于最大水平主应力准则，根据不同施工时间下近井壁岩石的最大主应力的计算结果反演起裂时间，具体如下：通过步骤S2获取不同施工时间下的页岩最大主应力；以计算结果作为插值的基础数据，对数据进行插值，计算井壁页岩最大主应力达到岩石抗拉强度时的施工时间。根据最大水平主应力准则，当岩石的最大主应力达到岩石的抗拉强度时，岩石开始破坏，因此，该施工时间为起裂时间。
步骤303：根据不同施工时间下的井口压力计算结果(即步骤S1得到的井口压力随时间变化的规律)反演起裂压力，具体如下：通过步骤S1获取不同施工时间下的井口压力；以计算结果作为插值的基础数据，对数据进行插值，计算施工时间为起裂时间时井口压力的大小，即为井口起裂压力。
本发明的实施例为川东南某页岩气井第1段测试压裂。如图1所示，具体包括以下步骤：
S1：建立井筒压裂液流动的离散化井筒模型，获得井底压力随时间变化的规律和井口压力随时间变化的规律：获取施工排量、开井压力，结合压裂液降阻率、粘度、压缩性等参数，建立憋压过程中井筒压裂液流动的数学模型，求解井底及井口压力随时间的变化规律；
S2：建立射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型，获得井筒憋压过程中近井壁应力场的分布规律和渗流场的分布规律：基于渗流力学及岩石力学相关知识，结合射孔参数及井筒、水泥环、地层的力学参数，建立井筒-水泥环-近井壁地层三维流固耦合力学模型，计算井筒憋压过程中近井壁应力场及渗流场的分布规律；
S3：反演地层起裂时间及井口起裂压力：结合计算结果，分析页岩最大主应力在井筒憋压过程中的变化规律，结合最大水平主应力准则，反演地层起裂时间及井口起裂压力。
步骤S1如图2所示，具体包括：
步骤101：根据压裂设计，本段所设计使用的压裂液为滑溜水，其粘度、密度、降阻率通过室内实验可以获取。为方便验证本模型的准确性，计算中施工排量、开井压力等参数来自于现场施工数据；
步骤102：将井筒离散为80个井筒微元段，单个井筒微元段的长度为70m，单位时间步长为1s，即：n＝80；δz＝80；δt＝1；
构造各井筒微元段的连续性方程及井筒压降方程为：
对式(8)、式(9)进行求解的步骤包括：
由式(7)可知，未知数为关于压力的函数，结合式(9)，式(8)可化为关于p的非线性方程，即：
R(p)＝0 (10)
式中，R＝[R1,R2,…,Rn]，p＝[p1,p2,…,pn]
在各时间步长内，以上一步计算的压力分布为初始值，代入式(9)中，通过二分法、弦割法等数值方法从井口微元段开始求解井筒内压裂液的流速分布，代入式(10)中等号左边求取R。同时，采用差分法对非线性方程内的各井筒微元段压力求导，得：
由此，可通过牛顿-迭代法更新井筒内压力分布，即：
p*＝Δp+p0 (13)
式中，p0为试算值，p*为更新结果。
当计算误差小于设定误差时，则停止计算，进入下一步长的计算；若不是，则继续更新。
步骤103：计算结果如图5所示，由图可以看出，在地层破裂前，计算的井口压力与实际施工值较为接近，说明所建立的模型准确性及可靠性较好，可用于起裂压力的计算中。输出憋压过程中不同时间的井底压力及井口压力。
步骤S2如图4所示，包括：
步骤201中，根据井身结构及压裂设计参数，建立射孔条件下套管-水泥环-地层三维有限元模型，如图6所示。
步骤202：基于套管、水泥环的材料参数，对力学模型中的套管及水泥环赋予材料属性；在测井资料及室内岩心实验的基础上，获取页岩的物理及力学性质，对页岩赋予物理力学属性；
步骤203：对射孔条件下套管-水泥环-近井壁地层三维力学模型施加孔隙压力边界条件即位移边界条件，孔隙压力来源于测井资料及钻井数据；施加了初始地应力场，地应力场来源于室内声发射实验及测井相关数据；
步骤204：对射孔条件下套管-水泥环-近井壁地层三维力学模型的套管内壁及孔眼处施加液压载荷，液压载荷的大小为步骤103中的输出结果，即憋压过程中随时间变化的井底压力；
步骤205：计算憋压过程中近井壁页岩的应力场及渗流场分布。
所示步骤S3如图5所示，包括：
步骤301：在步骤103和步骤205的基础上，分析及确定憋压过程中不同施工时间下井底压力和近井壁岩石的最大主应力；该条件下近井壁页岩的最大主应力随时间的变化曲线如图7所示。
步骤302：由计算结果提取用于反演起裂时间的最大水平主应力-憋压时间的数据。根据所取岩心的室内实验结果，页岩抗拉强度为8.23MPa。结合三次样条插值，当最大水平主应力达到8.23MPa，起裂时间为147.0941s。
步骤303：由计算结果提取用于反演起裂压力的憋压时间-井口压力的数据。根据上一步的反演结果，起裂时间为147.0941s。结合三次样条插值，井口起裂压力为81.9403MPa。
根据测试压裂的实际施工结果，起裂时间为149s，起裂压力为81.42MPa。本发明方法计算的结果在起裂压力上与实际计算结果误差为0.64％，在起裂时间上与实际计算结果误差为1.28％。说明本方法的计算方法及所建立的模型较为可靠，可准确地反映页岩气井水力压裂的裂缝起裂机理，并为页岩气井的水力压裂设计提供理论依据。
本发明的方法完善了页岩气井压裂裂缝起裂压力及起裂时间的预测方法，为页岩气井水力压裂设计提供了理论依据。
上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。
Claims (8)
1.一种考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述方法包括以下步骤：
S1：建立井筒压裂液流动的离散化井筒模型，获得井底压力随时间变化的规律和井口压力随时间变化的规律；
S2：建立射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型，获得井筒憋压过程中近井壁应力场的分布规律和渗流场的分布规律；
S3：反演地层起裂时间及井口起裂压力；
所述步骤S1中的建立井筒压裂液流动的离散化井筒模型的操作包括：
沿套管轴向将套管等距离划分为n个井筒微元段，每个井筒微元段长度为Δzi；第一个井筒微元段对应井口，第n个井筒微元段对应井底；
建立井筒压裂液流动的离散化井筒模型如下：
其中，为t时刻第i个井筒微元段的压裂液质量，采用下式计算得到：
P为压力，v为平均流速，表示第i个井筒微元段和第i+1个井筒微元段相交处截面，表示第i个井筒微元段和第i-1个井筒微元段相交处截面，Δt为单位时间步长，α为井筒微元段与垂向的夹角，d为套管内壁直径，g为重力加速度，ρ0为标准大气压下的压裂液密度，p0为标准大气压，C为压裂液压缩系数；τw为井筒壁面摩擦阻力；
所述步骤S2的操作包括：
步骤201：获取套管、水泥环和井眼尺寸，结合射孔施工参数，利用有限元法建立射孔条件下的套管-水泥环-近井壁地层三维流固耦合力学模型；
步骤202：对所述射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型中的套管、水泥环赋予材料属性；对所述射孔条件下的井筒-水泥环-近井壁地层三维流固耦合力学模型中的页岩赋予物理力学属性；
步骤203：对所述射孔条件下的套管-水泥环-近井壁地层三维流固耦合力学模型施加远场边界载荷；
步骤204：对所述射孔条件下套管-水泥环-近井壁地层三维流固耦合力学模型的套管内壁及孔眼处施加液压载荷，所述液压载荷的大小为步骤S1中得到的井底压力；
步骤205：计算得到憋压过程中近井壁页岩的应力场的分布规律及渗流场的分布规律。
2.根据权利要求1所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤S1中的获得井底压力随时间变化的规律和井口压力随时间变化的规律的操作包括：
(1)在单位时间步长下求解所述井筒压裂液流动的离散化井筒模型，获得该单位时间步长下的井底压力和井口压力；所述井底压力对应求解得到的第n个井筒微元段在该单位时间步长下的压力pn，所述井口压力对应求解得到的第1个井筒微元段在该单位时间步长下的压力p1；
(2)判断pn是否大于等于最大水平主应力的Y倍，如果否，则t＝t+Δt,然后返回步骤(1)，如果是，则转入步骤(3)；
(3)输出井底压力随时间变化的规律和井口压力随时间变化的规律：所述井底压力随时间变化的规律对应pn在各个单位时间步长下的数值，所述井口压力随时间变化的规律对应p1在各个单位时间步长下的数值。
3.根据权利要求2所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤(1)中是采用牛顿迭代法在各个单位时间步长下求解所述井筒压裂液流动的离散化井筒模型；
所述步骤(2)中的Y为1.5。
4.根据权利要求1所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤201中采用三维六面体缩减积分孔压单元对所述套管-水泥环-近井壁地层三维力学模型进行网格划分，并在孔眼处对网格进行加密。
5.根据权利要求1所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤202中的物理力学属性包括：渗透率、孔隙度、弹性模量、泊松比。
6.根据权利要求1所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤S3的操作包括：
步骤301：输入步骤S1获得的井底压力随时间变化的规律；根据步骤S2得到的应力场的分布规律计算得到不同施工时间下的近井壁岩石的最大主应力；
步骤302：根据所述不同施工时间下的近井壁岩石的最大主应力反演地层起裂时间；
步骤303：根据所述地层起裂时间和步骤S1得到的井口压力随时间变化的规律反演井口起裂压力。
7.根据权利要求6所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤302的操作包括：
以所述不同施工时间下的近井壁岩石的最大主应力作为插值的基础数据，对该基础数据进行插值，得到井壁页岩的最大主应力达到岩石抗拉强度时的施工时间，该施工时间即为地层起裂时间。
8.根据权利要求6所述的考虑压裂液压缩性的地层起裂压力预测方法，其特征在于：所述步骤303的操作包括：
将步骤S1得到的井口压力随时间变化的规律中的各个井口压力的值作为插值的基础数据，对该基础数据进行插值，然后计算施工时间为地层起裂时间时的井口压力的大小，该井口压力即为井口起裂压力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910584484.1A CN112177598B (zh) | 2019-07-01 | 2019-07-01 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910584484.1A CN112177598B (zh) | 2019-07-01 | 2019-07-01 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112177598A CN112177598A (zh) | 2021-01-05 |
| CN112177598B true CN112177598B (zh) | 2024-03-08 |
Family
ID=73914790
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910584484.1A Active CN112177598B (zh) | 2019-07-01 | 2019-07-01 | 一种考虑压裂液压缩性的地层起裂压力预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112177598B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113109162B (zh) * | 2021-04-13 | 2022-02-01 | 西南石油大学 | 一种基于热流固耦合的岩石起裂压力计算方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7066019B1 (en) * | 1999-06-23 | 2006-06-27 | Schlumberger Technology Corporation | Cavity stability prediction method for wellbores |
| CN103299031A (zh) * | 2010-12-22 | 2013-09-11 | 莫里斯·B·杜西奥尔特 | 用于提高页岩资源产量的多阶段裂缝灌注法 |
| CN105201479A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
| CN106555575A (zh) * | 2015-09-24 | 2017-04-05 | 中国石油化工股份有限公司 | 一种深层页岩的地层起裂压力的计算方法 |
| CN106640061A (zh) * | 2016-12-06 | 2017-05-10 | 西南石油大学 | 一种井筒与地层裂缝耦合流动模拟实验装置及方法 |
| WO2017192635A1 (en) * | 2016-05-03 | 2017-11-09 | Schlumberger Technology Corporation | Methods and systems for analysis of hydraulically-fractured reservoirs |
| CN108798654A (zh) * | 2018-04-27 | 2018-11-13 | 中国石油天然气股份有限公司 | 确定页岩气井的井底压力与时间的对应关系的方法及装置 |
| CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
| CN109209350A (zh) * | 2018-09-12 | 2019-01-15 | 西南石油大学 | 一种煤层斜井套管射孔完井破裂压力的预测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2854371C (en) * | 2011-11-04 | 2019-12-24 | Schlumberger Canada Limited | Modeling of interaction of hydraulic fractures in complex fracture networks |
-
2019
- 2019-07-01 CN CN201910584484.1A patent/CN112177598B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7066019B1 (en) * | 1999-06-23 | 2006-06-27 | Schlumberger Technology Corporation | Cavity stability prediction method for wellbores |
| CN103299031A (zh) * | 2010-12-22 | 2013-09-11 | 莫里斯·B·杜西奥尔特 | 用于提高页岩资源产量的多阶段裂缝灌注法 |
| CN106555575A (zh) * | 2015-09-24 | 2017-04-05 | 中国石油化工股份有限公司 | 一种深层页岩的地层起裂压力的计算方法 |
| CN105201479A (zh) * | 2015-10-09 | 2015-12-30 | 西南石油大学 | 一种页岩储层水平井分段压裂射孔簇参数优化设计方法 |
| WO2017192635A1 (en) * | 2016-05-03 | 2017-11-09 | Schlumberger Technology Corporation | Methods and systems for analysis of hydraulically-fractured reservoirs |
| CN106640061A (zh) * | 2016-12-06 | 2017-05-10 | 西南石油大学 | 一种井筒与地层裂缝耦合流动模拟实验装置及方法 |
| CN108798654A (zh) * | 2018-04-27 | 2018-11-13 | 中国石油天然气股份有限公司 | 确定页岩气井的井底压力与时间的对应关系的方法及装置 |
| CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
| CN109209350A (zh) * | 2018-09-12 | 2019-01-15 | 西南石油大学 | 一种煤层斜井套管射孔完井破裂压力的预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112177598A (zh) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104895550B (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
| CN109374867B (zh) | 一种基于离散元的砂砾岩水力压裂模拟方法 | |
| CN111236908A (zh) | 一种适用在低渗透致密气藏中的多段压裂水平井产能预测模型及产能敏感性分析的方法 | |
| CN103590824B (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
| CN107545113A (zh) | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 | |
| CN102418511B (zh) | 低渗透储层井下关井水井压降试井分析方法 | |
| US20060015310A1 (en) | Method for simulation modeling of well fracturing | |
| CN110984973A (zh) | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 | |
| CN108843303B (zh) | 一种基于泥岩蠕变模型的油水井套管损坏预测方法 | |
| Michael | Hydraulic Fracturing Optimization: Experimental Investigation of Multiple Fracture Growth Homogeneity via Perforation Cluster Distribution | |
| Sobhaniaragh et al. | Numerical simulation of MZF design with non-planar hydraulic fracturing from multi-lateral horizontal wells | |
| Liu et al. | Dynamic behavior of coalbed methane flow along the annulus of single-phase production | |
| CN109577945A (zh) | 一种低渗-超低渗油藏窜流通道判别的实验装置与方法 | |
| CN112177598B (zh) | 一种考虑压裂液压缩性的地层起裂压力预测方法 | |
| Wang et al. | Investigation of the interface cracks on the cement sheath stress in shale gas wells during hydraulic fracturing | |
| Wang et al. | Impact of fluid compressibility for plane strain hydraulic fractures | |
| Yan et al. | Experimental study of hydraulic fracturing for unconsolidated reservoirs | |
| CN115510778A (zh) | 一种陆相页岩储层无限级压裂工艺优化方法及*** | |
| Shan et al. | A new finite element method to predict the fracture initiation pressure | |
| CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
| WO2021072108A1 (en) | Determination of a surface leak rate in an injection well | |
| Roustaei | Yield stress fluid flows in uneven geometries: applications to the oil & gas industry | |
| Xu et al. | A 3D tortuous fracture network construction approach to analyze proppant distribution in post-fractured shale | |
| CN112035993A (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
| CN111950111A (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |