CN110397430B - 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 - Google Patents
一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 Download PDFInfo
- Publication number
- CN110397430B CN110397430B CN201910611601.9A CN201910611601A CN110397430B CN 110397430 B CN110397430 B CN 110397430B CN 201910611601 A CN201910611601 A CN 201910611601A CN 110397430 B CN110397430 B CN 110397430B
- Authority
- CN
- China
- Prior art keywords
- acid
- fracture
- crack
- seam
- width
- 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
- 239000002253 acid Substances 0.000 title claims abstract description 219
- 239000007788 liquid Substances 0.000 title claims abstract description 58
- 239000011435 rock Substances 0.000 title claims abstract description 53
- 238000009826 distribution Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 42
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 20
- 238000002347 injection Methods 0.000 claims abstract description 40
- 239000007924 injection Substances 0.000 claims abstract description 40
- 230000008569 process Effects 0.000 claims abstract description 25
- 239000000243 solution Substances 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000035699 permeability Effects 0.000 claims description 17
- 238000010276 construction Methods 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 13
- 238000005530 etching Methods 0.000 claims description 12
- 239000011148 porous material Substances 0.000 claims description 12
- 241000361919 Metaphire sieboldi Species 0.000 claims description 9
- 241001233061 earthworms Species 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 230000036632 reaction speed Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 10
- 230000000694 effects Effects 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 6
- 208000010392 Bone Fractures Diseases 0.000 description 118
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 239000007789 gas Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000007797 corrosion Effects 0.000 description 3
- 238000005260 corrosion Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000020477 pH reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 208000003044 Closed Fractures Diseases 0.000 description 1
- 102000020897 Formins Human genes 0.000 description 1
- 108091022623 Formins Proteins 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 229910000019 calcium carbonate Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
Images
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
- 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
-
- 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
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
本发明涉及一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，包括：将前置液酸压过程分为压裂造缝时段和注酸时段，采用二维PKN裂缝延伸扩展模型模拟压裂造缝时段裂缝延伸扩展过程，得到压裂造缝时段t1结束时的裂缝长度
裂缝宽度
和缝内压力p0；在注酸时段t2内，将t2等分为n个时间步，分别得到第一个时间步结束时缝长方向的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力p1；同样得到注酸时段t2内任意一个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力pj(j＝1,2,3,···,n)；计算前置液酸压裂缝闭合后的裂缝导流能力分布wkf。本发明能够准确预测酸蚀裂缝沿缝长方向上的导流能力分布，优化碳酸盐岩前置液酸压设计方案，提升碳酸盐岩开发效果。
Description
技术领域
本发明涉及油气田开发酸化压裂储层改造领域，具体涉及一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法。
技术背景
碳酸盐岩储层是一种重要的油气资源储集体，具有巨大的开发潜力。碳酸盐岩储层主要以石灰岩、白云岩为主，埋深一般超过5000m(齐宝权,赵佐安,贺洪举,等.应用测井储层品质综合评价指数来判别储层的有效性——以四川盆地西部地区二叠系栖霞组为例[J].天然气工业,2018,38(02):25-32)。这类储层具有较低的孔隙度、渗透率以及强非均质性，在自然条件下很难进行油气资源开采。因此，碳酸盐岩储层开发主要是让岩石成分中的碳酸钙与酸发生反应来达到改善储层的目的，可以通过基质酸化或者酸压裂的方式来实现这一目的，从而改善储层流通性能，前者主要形成具有高渗透性的酸蚀蚓孔通道，后者主要形成具有高导流能力的人工酸蚀裂缝。
酸压是指在高于地层破裂压力条件下，将酸液注入地层，在地层中形成人工裂缝，同时酸液与裂缝壁面岩石发生反应，非均匀刻蚀裂缝壁面，形成沟槽状、凹凸不平的刻蚀或者形成酸蚀蚓孔，施工结束后裂缝不完全闭合，最终形成具有一定宽度和较高导流能力的人工裂缝，从而实现油气井的增产。
为了得到高导流能力酸蚀裂缝，国内外逐渐发展形成了前置液酸压工艺。前置液酸压工艺首先需要向储层注入不与地层发生化学反应的压裂液(前置液)来压开地层形成人工裂缝，然后再向地层注入酸液，对人工裂缝壁面产生非均匀溶蚀，使其无法完全闭合，从而获得高导流能力酸蚀裂缝。前置液酸压后的裂缝导流能力的大小直接可以反应储层增产改造的效果好坏，因此预测酸压裂缝导流能力是一项非常重要的工作(M.J Economides,A.D Hill,D.Zhu.Petroleum production systems[J].Journal of Petroleum Scienceand Engineering,1996,15:2-4)。
前置液酸压裂缝的导流能力大小与地层闭合应力、地层非均质性、岩石酸溶蚀体积以及酸蚀的不规则程度等均有关系。目前，常用室内岩心实验法(苟波,郭建春,陈迟,等.实验确定碳酸盐岩油气藏酸压裂缝导流能力分布的方法,CN 106522935 A.2017)和基于各地碳酸盐岩储层改造资料和相关实验的导流能力经验公式计算方法来获取裂缝导流能力。室内岩心实验法获取裂缝导流能力最为直观、简便，但是缺点是难以实现地层真实环境下的裂缝导流能力预测，并且受限于实验岩心样品的尺寸大小，从而无法得到沿裂缝长度方向上的真实导流能力分布。而根据经验公式计算得到的裂缝导流能力在应用上更为广泛，也更加容易实现，但是往往存在一定的误差。为了得到更加准确的导流能力计算结果，需要综合考虑多种因素对裂缝导流能力的影响。
发明内容
本发明目的在于提供一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，该方法考虑了地层渗透率非均质性、酸岩反应和酸蚀蚓孔滤失对裂缝导流能力的影响，能够准确预测前置液酸压形成的酸蚀裂缝沿缝长方向上的导流能力分布，优化碳酸盐岩前置液酸压设计方案，提升碳酸盐岩开发效果。
为达到以上技术目的，本发明采用以下技术方案。
一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，依次包括以下步骤：
(1)、将前置液酸压过程分为压裂造缝时段和注酸时段，施工时间分别为t1和t2，采用二维PKN裂缝延伸扩展模型模拟压裂造缝时段裂缝延伸扩展过程，通过求解PKN模型得到压裂造缝时段t1结束时的裂缝长度
裂缝宽度
和缝内压力p0；
(2)、在注酸时段t2内，将t2等分为n个时间步，等时间步长为Δt2，在注酸开始时，以步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件，利用裂缝内酸岩反应平衡方程，计算得到第一个时间步结束时缝长方向的缝内酸液浓度
(3)、将步骤(2)中得到的缝长方向的缝内酸液浓度
代入考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程，得到第一个时间步结束时缝长方向的酸蚀裂缝宽度
(4)、将步骤(3)中得到的酸蚀裂缝宽度
代入缝内酸液流动平衡方程，得到第一个时间步结束时的缝内酸液流动压力p1；
(5)、将步骤(3)、(4)中分别得到的酸蚀裂缝宽度
和缝内酸液流动压力p1作为第二个时间步开始时的初始条件，重复(2)、(3)、(4)，得到第二个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力p2。同样地，可以得到注酸时段t2内任意一个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力pj，j为第j个时间步(j＝1,2,3,···,n)，注酸时段t2结束时的缝内酸液浓度、酸蚀裂缝宽度和缝内酸液流动压力分别为
和pn；
(6)、将步骤(5)得到的酸蚀裂缝宽度
代入裂缝导流能力计算公式，得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
利用注酸时段结束时的裂缝导流能力分布
计算得到前置液酸压裂缝闭合后的裂缝导流能力分布wkf。
在本发明中，所述步骤(1)中将前置液酸压过程分为压裂造缝时段和注酸时段，施工时间分别为t1和t2，采用二维PKN裂缝延伸扩展模型模拟前置液酸压裂缝延伸扩展过程，通过求解PKN模型得到压裂造缝时段t1结束时的裂缝长度
裂缝宽度
和缝内压力p0，过程如下：
1)运用二维PKN裂缝延伸扩展模型，通过下式计算压裂造缝时段t1结束时的裂缝长度
(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京：石油工业出版社,1998)：
式中：
——压裂造缝时段结束时的裂缝长度，m；
Q1——前置液的施工排量，m3/min；
Hf——裂缝高度，m；
C1——前置液的滤失系数，m/min0.5。
2)通过下式计算该时段结束时的缝口处的裂缝宽度(李颖川.采油工程[M].石油工业出版社,2009)：
式中：
——缝口处的裂缝宽度，m；
x——以缝口为原点的缝内任意横截面处的位置坐标，m；
μ1——前置液的粘度，mPa·s；
E——杨氏模量，MPa；
ν——泊松比，无量纲。
通过下式计算压裂造缝时段t1结束时缝长方向上的裂缝宽度
(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998)：
式中：
——造缝时段结束时的裂缝宽度，m。
3)通过下式计算压裂造缝时段t1结束时的缝内压力p0(L.Dang.C.Zhou.Simulation of effective fracture length of prepad acidfracturing considering multiple leak-off effect[J].Natural Gas Industry B,2019,6:64-70)：
式中：p0——造缝时段结束时的缝内压力，MPa；
q——缝内任意横截面处的流量，m3/min。
缝内任意横截面处的流量计算如下(王鸿勋，张士诚.水力压裂设计数值计算方法[M].北京：石油工业出版社，1998)：
在本发明中，所述步骤(2)过程如下：在注酸时段t2内，将t2等分为n个时间步，等时间步长为Δt2，在注酸开始时，以步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件，利用裂缝内酸岩反应平衡方程，计算得到第一个时间步结束时缝长方向的缝内酸液浓度
Dang.C.Zhou.Simulation of effective fracture length of prepad acidfracturing considering multiple leak-off effect[J].Natural Gas Industry B,2019,6:64-70)：
式中：
——第一个时间步结束时的缝内酸液浓度，kmol/m3；
CW——裂缝壁面处的酸液浓度分布，kmol/m3，一般情况下视为0，即完全反应；
kg——酸岩反应速度，m/min；
t2——酸液注入施工时间，min；
μ2——酸液粘度，mPa·s；
k——储层基质渗透率，10-3um2；
pe——储层压力，MPa；
d——酸蚀蚓孔长度，m。
在本发明中，所述步骤(3)过程如下：将步骤(2)中得到的缝长方向的缝内酸液浓度
代入考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程，得到第一个时间步结束时缝长方向的酸蚀裂缝宽度
(赵立强,缪尉杰,罗志锋.闭合酸蚀裂缝导流能力模拟研究[J].油气藏评价与开发,2019,9(02):25-32)：
式中：
——第一个时间步结束时的酸蚀裂缝宽度，m；
β——酸溶解岩石的能力，kg/kmol；
ρr——储层岩石密度，kg/m3；
η——滤失酸液中发生酸岩反应的酸液百分比，一般情况下视为0，即忽略滤失酸液的酸岩反应。
在本发明中，所述步骤(4)过程如下：将步骤(3)中得到的酸蚀裂缝宽度
代入缝内酸液流动平衡方程，得到第一个时间步结束时的缝内酸液流动压力p1(胡晋阳.裂缝性碳酸盐岩储层水平井酸化温度场模型研究[D].西南石油大学,2016)：
式中：p1——缝内酸液流动压力，MPa。
在本发明中，所述步骤(5)过程如下：将步骤(3)、(4)中分别得到的酸蚀裂缝宽度
和缝内酸液流动压力p1作为第二个时间步开始时的初始条件，重复(2)、(3)、(4)，得到第二个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力p2。同样地，可以得到注酸时段t2内任意一个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力pj，j为第j个时间步(j＝1,2,3,···,n)，注酸时段t2结束时的缝内酸液浓度、酸蚀裂缝宽度和缝内酸液流动压力分别为
和pn(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998)。
在本发明中，所述步骤(6)中将步骤(5)得到的酸蚀裂缝宽度
代入裂缝导流能力计算公式，得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
利用注酸时段结束时的裂缝导流能力分布
可以计算得到整个前置液酸压裂缝闭合后的裂缝导流能力分布wkf，过程如下：
1)在前置液酸压施工过程中，产生的人工裂缝尚未开始闭合，可以认为整个前置液酸压过程中闭合应力为零，求取岩石闭合应力为零时的酸蚀裂缝导流能力(J.Deng,J.Mou.A New Acid-Fracture Conductivity Model Based on the SpatialDistributions of Formation Properties[C].SPE International Symposium andExhibition on Formation Damage Control,2010)，即得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
式中：
——闭合应力为零时的导流能力，10-3um2·cm；
C2——酸蚀蚓孔滤失系数，m/min0.5；
λD——以缝口为原点的缝内任意横截面处的无因次位置坐标，无量纲；
σD——地层渗透率非均匀程度，无量纲；
σ(lnk)——地层渗透率的标准差；
a1,a2,a3——误差系数，a1＝1.82，a2＝3.25，a3＝0.12；
erf(y)——误差函数，根据y值大小在误差补偿函数表中得到相应值。
2)在前置液酸压施工结束后，裂缝在地层闭合应力作用下开始逐渐闭合，将注酸时段结束时的导流能力
代入下式，计算得到整个前置液酸压裂缝闭合后的裂缝导流能力分布wkf(J.Deng,J.Mou.A New Correlation of Acid Fracture ConductivitySubject to Closure Stress[C].SPE Hydraulic Fracturing Technology Conferenceand Exhibition,2011)：
式中：wkf——前置液酸压裂缝闭合后的裂缝导流能力，10-3um2·cm；
σC——地层闭合应力，MPa；
E——杨氏模量，MPa。
与现有技术相比，本发明的有益效果在于：
(1)缝长方向上的酸蚀缝宽求解基于经典二维PKN裂缝延伸扩展过程、酸岩反应过程和酸蚀蚓孔滤失过程，能够准确获取注酸过程任意时间的裂缝宽度分布；
(2)考虑了缝内酸液流动过程压力的重新分布，能够实时更新缝内酸液流速，实现任意时间缝内酸液流动反应过程的酸液浓度的准确计算；
(3)考虑了碳酸盐岩储层渗透率的非均质性对酸蚀裂缝导流能力的影响，实现了沿缝长方向上裂缝导流能力分布的准确预测。
本发明克服了现有方法未能综合考虑酸液流动、酸岩反应、酸蚀蚓孔滤失、储层渗透率非均质性等多种因素对前置液酸压裂缝导流能力的影响，对于优化前置液酸压施工方案设计，提升前置液酸压改造效果具有重要的指导意义。
附图说明
图1是本发明压裂造缝时段结束时的缝长-缝宽分布图。
图2是本发明压裂造缝时段结束时的缝内压降图。
图3是本发明注酸时段结束时的缝内酸液分布图。
图4是本发明注酸时段结束时的裂缝导流能力分布图。
具体实施方式
下面结合附图和具体实例对本发明作进一步说明。
各种参数如表1所示：
表1实施例中计算导流能力的各类参数
一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，依次包括以下步骤：
(1)将前置液酸压过程按照施工时间分为压裂造缝时段和注酸时段，对应的施工时间分别为t1和t2，采用二维PKN裂缝延伸扩展模型模拟前置液酸压裂缝延伸扩展过程，通过求解PKN模型得到压裂造缝时段t1结束时的裂缝长度
裂缝宽度
和缝内压力p0(计算结果分别如图1和图2所示)：
1)运用二维PKN裂缝延伸扩展模型，通过下式计算压裂造缝时段t1结束时的裂缝长度
(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京：石油工业出版社,1998)：
式中：
——压裂造缝时段结束时的裂缝长度，m；
Q1——前置液的施工排量，m3/min；
Hf——裂缝高度，m；
C1——前置液的滤失系数，m/min0.5。
2)然后通过计算该时段结束时的缝口处的裂缝宽度(李颖川.采油工程[M].石油工业出版社,2009)：
式中：
——缝口处的裂缝宽度，m；
x——以缝口为原点的缝内任意横截面处的位置坐标，m；
μ1——前置液的粘度，mPa·s；
E——杨氏模量，MPa；
ν——泊松比，无量纲。
最后通过下式计算压裂造缝时段t1结束时缝长方向上的裂缝宽度
(王鸿勋,张士诚.水力压裂设计数值计算方法[M].北京:石油工业出版社,1998)：
式中：
——造缝时段结束时的裂缝宽度，m。
3)通过下式计算压裂造缝时段t1结束时的缝内压力p0(L.Dang.C.Zhou.Simulation of effective fracture length of prepad acidfracturing considering multiple leak-off effect[J].Natural Gas Industry B,2019,6:64-70)：
式中：p0——造缝时段结束时的缝内压力，MPa；
q——缝内任意横截面处的流量，m3/min。
缝内任意横截面处的流量计算如下(王鸿勋，张士诚.水力压裂设计数值计算方法[M].北京：石油工业出版社，1998)：
(2)在注酸时段t2内，将t2等分为n个时间步，等时间步长为Δt2，在注酸开始时，以步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件，利用裂缝内酸岩反应平衡方程，计算得到第一个时间步结束时的缝长方向的缝内酸液浓度
1)求解缝长方向的缝内酸液浓度
需要将裂缝内酸岩反应平衡方程进行离散化，酸岩反应平衡方程为：
式中：Cf——缝内酸液浓度，kmol/m3；
wf——裂缝宽度，m；
CW——裂缝壁面处的酸液浓度分布，kmol/m3，一般视为0，即完全反应；
vx——沿缝长方向酸液流速，m3/min；
vl——裂缝壁面酸蚀蚓孔滤失速度，m3/min；
kg——酸岩反应速度，m/min。
其中，流体在裂缝中流动的运动方程和滤失方程分别为：
式中：p——缝内压力，MPa；
μ2——酸液粘度，mPa·s；
k——储层基质渗透率，10-3um2；
pe——储层压力，MPa；
d——酸蚀蚓孔长度，m。
离散后的裂缝内酸岩反应平衡方程为：
其中：
式中：
——在第n个时间步结束时刻，第i个、第i+1个、第i-1个缝长单元所对应的酸液浓度，kmol/m3；
Δx——缝长单元长度，m；
Δt2——时间步长，min；
ki——第i个缝长单元所对应的储层基质渗透率，10-3um2。
计算过程所需要的酸液浓度边界条件为：
式中：C0——缝口处酸液浓度，kmol/m3，保持恒定不变。
2)将步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件和上述给定的酸液浓度边界条件一起代入离散后的裂缝内酸岩反应平衡方程，使用MATLAB编程计算，可得到第一个时间步结束时缝长方向的缝内酸液浓度
(3)将步骤(2)中得到的缝长方向的缝内酸液浓度
代入考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程，得到第一个时间步结束时的缝长方向的酸蚀裂缝宽度
1)求解缝长方向的酸蚀裂缝宽度
需要将考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程进行离散化，酸岩反应平衡方程为：
式中：β——酸溶解岩石的能力，kg/kmol；
ρr——储层岩石密度，kg/m3；
η——滤失酸液中发生酸岩反应的酸液百分比，一般情况下视为0，即忽略滤失酸液的酸岩反应。
离散后的裂缝宽度变化方程为：
2)将步骤(2)中得到的缝长方向的缝内酸液浓度
代入离散后的裂缝宽度变化方程，使用MATLAB编程计算，可得到第一个时间步结束时的缝长方向的酸蚀裂缝宽度
(4)将步骤(3)中得到的酸蚀裂缝宽度
代入缝内酸液流动平衡方程，得到第一个时间步结束时的缝内酸液流动压力p1：
1)求解缝长方向的酸蚀裂缝宽度
需要将缝内酸液流动平衡方程进行离散化，缝内酸液流动平衡方程为：
其中，酸液在裂缝中流动的运动方程和滤失方程分别为：
离散后的缝内酸液流动平衡方程为：
其中：
计算过程所需要的缝内压力边界条件为：
2)将步骤(3)中得到的酸蚀裂缝宽度
和上述给定的缝内压力边界条件一起代入离散后的缝内酸液流动平衡方程，使用MATLAB编程计算，可得到第一个时间步结束时的缝内酸液流动压力p1。
(5)将步骤(3)、(4)中分别得到的第一个时间步结束时的酸蚀裂缝宽度
和缝内酸液流动压力p1作为第二个时间步开始时的初始条件，重复(2)、(3)、(4)，得到第二个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力p2。同样地，可以得到注酸时段t2内任意一个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力pj，j为第j个时间步(j＝1,2,3,···,n)，注酸时段t2结束时的缝内酸液浓度、酸蚀裂缝宽度和缝内酸液流动压力分别为
和pn。
(6)将步骤(5)得到的酸蚀裂缝宽度
代入裂缝导流能力计算公式，得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
利用注酸时段结束时的裂缝导流能力分布
计算得到整个前置液酸压裂缝闭合后的裂缝导流能力分布wkf，过程如下：
1)在前置液酸压施工过程中，产生的人工裂缝尚未开始闭合，可以认为整个前置液酸压过程中闭合应力为零，求取岩石闭合应力为零时的酸蚀裂缝导流能力(J.Deng,J.Mou.A New Acid-Fracture Conductivity Model Based on the SpatialDistributions of Formation Properties[C].SPE International Symposium andExhibition on Formation Damage Control,2010)：
式中：
——闭合应力为零时的导流能力，10-3um2·cm；
C2——酸蚀蚓孔滤失系数，m/min0.5；
λD——以缝口为原点的缝内任意横截面处的无因次位置坐标，无量纲；
σD——地层渗透率非均匀程度，无量纲；
σ(lnk)——地层渗透率的标准差；
a1,a2,a3——误差系数，a1＝1.82，a2＝3.25，a3＝0.12；
erf(y)——误差函数，根据y值大小在误差补偿函数表中得到相应值。
2)在前置液酸压施工结束后，裂缝在地层闭合应力作用下开始逐渐闭合，将注酸时段结束时(此时闭合应力为零)的导流能力
代入下式，最终计算得到整个前置液酸压裂缝闭合后的裂缝导流能力分布(J.Deng,J.Mou.A New Correlation of AcidFracture Conductivity Subject to Closure Stress[C].SPE Hydraulic FracturingTechnology Conference and Exhibition,2011)：
式中：wkf——前置液酸压裂缝闭合后的裂缝导流能力，10-3um2·cm；
σC——地层闭合应力，MPa；
E——杨氏模量，MPa。
Claims (6)
1.一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，依次包括以下步骤：
(1)、将前置液酸压过程分为压裂造缝时段和注酸时段，施工时间分别为t1和t2，采用二维PKN裂缝延伸扩展模型模拟压裂造缝时段裂缝延伸扩展过程，通过求解PKN模型得到压裂造缝时段t1结束时的裂缝长度
裂缝宽度
和缝内压力p0；
(2)、在注酸时段t2内，将t2等分为n个时间步，等时间步长为Δt2，在注酸开始时，以步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件，利用裂缝内酸岩反应平衡方程，计算得到第一个时间步结束时缝长方向的缝内酸液浓度
(3)、将步骤(2)中得到的缝长方向的缝内酸液浓度
代入考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程，得到第一个时间步结束时缝长方向的酸蚀裂缝宽度
(4)、将步骤(3)中得到的酸蚀裂缝宽度
代入缝内酸液流动平衡方程，得到第一个时间步结束时的缝内酸液流动压力p1；
(5)、将步骤(3)、(4)中分别得到的酸蚀裂缝宽度
和缝内酸液流动压力p1作为第二个时间步开始时的初始条件，重复(2)、(3)、(4)，得到第二个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力p2，同样得到注酸时段t2内任意一个时间步结束时的缝内酸液浓度
酸蚀裂缝宽度
和缝内酸液流动压力pj，j为第j个时间步(j＝1，2，3，…，n)；
(6)、将步骤(5)得到的酸蚀裂缝宽度
代入裂缝导流能力计算公式，得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
利用注酸时段结束时的裂缝导流能力分布
计算得到前置液酸压裂缝闭合后的裂缝导流能力分布wkf。
2.如权利要求1所述的一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，其特征在于，所述步骤(1)过程如下：
通过下式计算压裂造缝时段t1结束时的裂缝长度
式中：Q1——前置液的施工排量，m3/min；
Hf——裂缝高度，m；
C1——前置液的滤失系数，m/min0.5；
通过下式计算压裂造缝时段t1结束时缝长方向上的裂缝宽度
式中：
——缝口处的裂缝宽度，m；
x——以缝口为原点的缝内任意横截面处的位置坐标，m；
μ1——前置液的粘度，mPa·s；
E——杨氏模量，MPa；
ν——泊松比，无量纲；
通过下式计算压裂造缝时段t1结束时的缝内压力p0：
式中：q——缝内任意横截面处的流量，m3/min。
3.如权利要求1所述的一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，其特征在于，所述步骤(2)过程如下：在注酸时段t2内，将t2等分为n个时间步，等时间步长为Δt2，在注酸开始时，以步骤(1)中得到的裂缝宽度
和缝内压力p0作为初始条件，利用裂缝内酸岩反应平衡方程，计算得到第一个时间步结束时缝长方向的缝内酸液浓度
式中：CW——裂缝壁面处的酸液浓度分布，kmol/m3；
kg——酸岩反应速度，m/min；
μ2——酸液粘度，mPa·s；
k——储层基质渗透率，10-3um2；
pe——储层压力，MPa；
d——酸蚀蚓孔长度，m。
4.如权利要求1所述的一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，其特征在于，所述步骤(3)过程如下：将步骤(2)中得到的缝长方向的缝内酸液浓度
代入考虑缝内酸蚀蚓孔滤失和酸岩反应的裂缝宽度变化方程，得到第一个时间步结束时缝长方向的酸蚀裂缝宽度
式中：β——酸溶解岩石的能力，kg/kmol；
ρr——储层岩石密度，kg/m3；
η——滤失酸液中发生酸岩反应的酸液百分比。
5.如权利要求1所述的一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，其特征在于，所述步骤(4)过程如下：将步骤(3)中得到的酸蚀裂缝宽度
代入缝内酸液流动平衡方程，得到第一个时间步结束时的缝内酸液流动压力p1：
式中：
——第一个时间步结束时的缝长方向酸液流速，m3/min；
μ2——酸液粘度，mPa·s。
6.如权利要求1所述的一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法，其特征在于，所述步骤(6)过程如下：
1)求取岩石闭合应力为零时的酸蚀裂缝导流能力，即得到注酸时段t2内任意一个时间步结束时的缝长方向的裂缝导流能力分布
式中：
——裂缝平均宽度，m；
C2——酸蚀蚓孔滤失系数，m/min0.5；
λD——以缝口为原点的缝内任意横截面处的无因次位置坐标，无量纲；
σD——地层渗透率非均匀程度，无量纲；
σ(lnk)——地层渗透率的标准差；
a1,a2,a3——误差系数，a1＝1.82，a2＝3.25，a3＝0.12；
erf(y)——误差函数，根据y值大小在误差补偿函数表中得到相应值；
2)将注酸时段结束时的导流能力分布
代入下式，计算得到前置液酸压裂缝闭合后的裂缝导流能力分布wkf：
式中：σC——地层闭合应力，MPa；
E——杨氏模量，MPa。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910611601.9A CN110397430B (zh) | 2019-07-08 | 2019-07-08 | 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910611601.9A CN110397430B (zh) | 2019-07-08 | 2019-07-08 | 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110397430A CN110397430A (zh) | 2019-11-01 |
| CN110397430B true CN110397430B (zh) | 2020-03-27 |
Family
ID=68324013
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910611601.9A Active CN110397430B (zh) | 2019-07-08 | 2019-07-08 | 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110397430B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111197476B (zh) * | 2020-01-08 | 2020-10-09 | 西南石油大学 | 一种考虑复杂滤失介质的酸压裂缝体刻蚀形态计算方法 |
| CN113671148B (zh) * | 2020-05-14 | 2024-05-28 | 中国石油天然气股份有限公司 | 一种酸蚀裂缝有效长度的实验测试方法 |
| CN112345436B (zh) * | 2020-09-22 | 2023-07-25 | 中国石油天然气股份有限公司 | 酸蚀裂缝导流能力分布测试方法 |
| CN112287533B (zh) * | 2020-10-21 | 2022-02-01 | 西南石油大学 | 一种裂缝性碳酸盐岩酸蚀裂缝导流能力计算方法 |
| CN112347411B (zh) * | 2020-10-30 | 2022-02-01 | 西南石油大学 | 一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法 |
| CN112360448B (zh) * | 2020-11-23 | 2021-06-18 | 西南石油大学 | 一种利用水力裂缝蠕变扩展确定压后焖井时间的方法 |
| CN112949050B (zh) * | 2021-02-09 | 2022-11-04 | 中国石油天然气股份有限公司 | 酸蚀参数的确定方法、装置、计算机设备及存储介质 |
| CN116838309B (zh) * | 2023-06-15 | 2024-03-19 | 长江大学 | 一种碳酸盐岩储层酸压裂缝有效长度测定方法 |
| CN117307152B (zh) * | 2023-11-28 | 2024-02-09 | 四川省能源地质调查研究所 | 页岩气井压裂裂缝导流能力动态评估方法及其装置 |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150107835A1 (en) * | 2013-10-17 | 2015-04-23 | Schlumberger Technology Corporation | Well treatment with shapeshifting particles |
| WO2016072877A1 (en) * | 2014-11-06 | 2016-05-12 | Schlumberger Canada Limited | Fractures treatment |
| US9938818B2 (en) * | 2015-01-14 | 2018-04-10 | Chevron U.S.A. Inc. | Hydraulic fracture permeability characterization from electromagnetic measurements |
| CN105134159B (zh) * | 2015-07-02 | 2016-03-23 | 西南石油大学 | 一种适用于缝洞型碳酸盐岩储层的高效酸压方法 |
| US20170174980A1 (en) * | 2015-12-17 | 2017-06-22 | Schlumberger Technology Corporation | Bio-fiber treatment fluid |
| CN106649963B (zh) * | 2016-10-14 | 2019-10-01 | 东北石油大学 | 体积压裂复杂缝网平均裂缝长度和等效裂缝条数确定方法 |
| CN106522935B (zh) * | 2016-12-26 | 2019-04-26 | 西南石油大学 | 实验确定碳酸盐岩油气藏酸压裂缝导流能力分布的方法 |
| CN106593429B (zh) * | 2017-01-09 | 2018-05-01 | 西南石油大学 | 一种致密储层体积裂缝导流能力测试方法 |
| CN107191169B (zh) * | 2017-06-26 | 2019-03-15 | 中国石油天然气股份有限公司 | 用于低产低效老油井的重复压裂方法 |
| CN107545113B (zh) * | 2017-09-08 | 2020-02-07 | 西南石油大学 | 非常规油气藏水力压裂复杂缝网形成过程模拟方法 |
| CN108442911A (zh) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | 一种页岩气水平井重复压裂水力裂缝参数优化设计方法 |
| CN108691526B (zh) * | 2018-05-10 | 2019-03-15 | 西南石油大学 | 实验确定碳酸盐岩油气藏粗糙裂缝酸刻蚀导流能力的方法 |
| CN108979612B (zh) * | 2018-08-09 | 2020-07-07 | 西南石油大学 | 一种致密油气藏压裂酸化复杂裂缝流动能力优化方法 |
| CN109681179B (zh) * | 2018-12-25 | 2021-01-22 | 成都理工大学 | 一种形成差异化刻蚀的酸压工艺方法 |
-
2019
- 2019-07-08 CN CN201910611601.9A patent/CN110397430B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110397430A (zh) | 2019-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110397430B (zh) | 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 | |
| Williams et al. | Design of acid fracturing treatments | |
| CN109751029B (zh) | 一种深层页岩气压裂的方法 | |
| CN108071378B (zh) | 一种适用于致密油气藏的co2泡沫压裂方法 | |
| Dang et al. | Simulation of effective fracture length of prepad acid fracturing considering multiple leak-off effect | |
| EP2820238B1 (en) | Methods, devices, and kits relating to the implementation of gelled fluids having a transient gel microstructure | |
| CN109710965B (zh) | 一种水平井人工压裂裂缝参数有效性的评价方法 | |
| Ketova et al. | Testing of preformed particles polymer gel technology on core filtration models to limit water inflows | |
| Du et al. | Experimental study on the autogenic acid fluid system of a high-temperature carbonate reservoir by acid fracturing | |
| CN109796951B (zh) | 近井地带解堵剂及其制备方法和应用方法 | |
| CN113987965B (zh) | 一种暂堵转向裂缝的预测方法及装置 | |
| CN110714755B (zh) | 水驱油藏剩余油二次富集速度快速预测方法 | |
| CN107760293B (zh) | 一种压裂预前置液处理方法 | |
| RU2597305C1 (ru) | Способ разработки нефтяной залежи в карбонатных коллекторах | |
| CN104712299B (zh) | 适合气井控水增气压裂的设计方法 | |
| Xiang-E | Development and application of the water control and profile modification technology in China oil fields | |
| CN108590571A (zh) | 一种低渗裂缝型见水油藏对于注水井调剖剂用量确定方法 | |
| CN113530534B (zh) | 储层改造方法及装置 | |
| CN112211609B (zh) | 一种深层碳酸盐岩油气藏重复酸压方法 | |
| CN113468476B (zh) | 一种聚合物驱增油量预测方法及装置 | |
| CN111852433B (zh) | 一种提高白云岩储层多尺度裂缝导流能力的方法 | |
| CN111075414B (zh) | 一种浅层中高渗超稠油油藏有效动用工艺 | |
| Zhang et al. | Process, mechanism and impacts of scale formation in alkaline flooding by a variable porosity and permeability model | |
| Jiang et al. | Steam and Gas Push (SAGP)-4; Recent Theoretical Developments and Laboratory Results Using Layered Models | |
| Liu et al. | Development and field application of strongly resilient temporary plugging diversion agent for fracturing |
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 |