CN109653740A - 一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 - Google Patents
一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 Download PDFInfo
- Publication number
- CN109653740A CN109653740A CN201811601635.1A CN201811601635A CN109653740A CN 109653740 A CN109653740 A CN 109653740A CN 201811601635 A CN201811601635 A CN 201811601635A CN 109653740 A CN109653740 A CN 109653740A
- Authority
- CN
- China
- Prior art keywords
- acid fracturing
- supporting
- self
- carbonate
- rock
- 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
- 239000002253 acid Substances 0.000 title claims abstract description 195
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000011435 rock Substances 0.000 claims abstract description 103
- 230000000694 effects Effects 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 27
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 22
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 230000007797 corrosion Effects 0.000 claims abstract description 9
- 238000005260 corrosion Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000035699 permeability Effects 0.000 claims description 19
- 239000007788 liquid Substances 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 claims description 2
- 208000010392 Bone Fractures Diseases 0.000 description 65
- 206010017076 Fracture Diseases 0.000 description 65
- 239000000243 solution Substances 0.000 description 12
- 239000011148 porous material Substances 0.000 description 6
- 125000005587 carbonate group Chemical group 0.000 description 5
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005530 etching Methods 0.000 description 3
- 238000004215 lattice model Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 235000002639 sodium chloride Nutrition 0.000 description 2
- 239000011780 sodium chloride Substances 0.000 description 2
- 239000000758 substrate Substances 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
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 235000013267 Pinus ponderosa Nutrition 0.000 description 1
- 241000555277 Pinus ponderosa Species 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 229920001169 thermoplastic Polymers 0.000 description 1
- 239000004416 thermosoftening plastic Substances 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
- 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
-
- 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
- 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
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
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)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，包括以下步骤：收集碳酸盐岩酸压自支撑裂缝基础参数；采取岩样模拟酸岩反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比；建立酸压自支撑裂缝应力敏感下的缝宽变化方程；建立考虑碳酸盐岩酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度预测模型；建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果。该预测模型更加接近实际地层的情况，进而有助于酸压自支撑裂缝渗透率预测的准确性。
Description
技术领域
本发明涉及石油与天然气勘探与开发领域，具体涉及到一种碳酸盐岩酸压自支撑裂缝渗透率预测方法。
背景技术
我国碳酸盐岩油气资源量丰富，其中海相碳酸盐岩油气资源量约385×108t油当量、湖相碳酸盐岩油气资源量约45×108t油当量，平均探明率仅为11％左右；碳酸盐岩资源目前已经成为我国勘探开发的重要领域。国内外勘探开发经验表明，酸压是非均质碳酸盐岩油气藏有效动用的关键技术。酸压是通过酸液对裂缝表面进行非均匀刻蚀，形成具有一定导流能力的酸蚀裂缝。酸蚀裂缝后形成的渗透率好坏是评价酸压成功与否的重要指标，有效预测酸压裂缝渗透率对于改善酸压效果具有重要意义。
目前针对酸压裂缝渗透率的预测多采用立方定律进行计算，根据立方定律，酸压裂缝渗透率只与初始裂缝宽度有关。然而，碳酸盐岩酸压自支撑裂缝的导流能力具有以下特点：酸压过程中酸液对岩石的溶蚀将显著降低岩石强度，岩石强度对酸压裂缝渗透率的影响不能忽略；当酸压裂缝较小时，酸压裂缝中气体滑脱效应不容忽视，且滑脱效应随着裂缝变小而增强；气体在微纳米尺度酸压微裂缝中流动时粘度将显著影响裂缝渗透率；酸压井在生产过程中，孔隙压力下降导致储层有效应力增加，使得裂缝渗透率显著降低。目前的酸压裂缝渗透率预测方法无法计量上述因素的影响。
发明内容
本发明的目的是提供一种碳酸盐岩酸压自支撑裂缝渗透率预测方法。
为达上述目的，本发明的一个实施例中提供了一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，其能够有助于提高酸压自支撑裂缝渗透率预测的准确性。
本发明实施例的具体技术方案是：
一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，其包括以下步骤：
收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数；
采取井下储层岩样，选用酸压用酸液类型，模拟地层温度环境下酸压过程中酸与岩石的反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比；
根据获取的杨氏模量和泊松比，基于酸压自支撑裂缝压缩性、碳酸盐岩基质压缩性对酸压自支撑裂缝宽度的影响，建立酸压自支撑裂缝应力敏感下的缝宽变化方程；
基于酸压自支撑裂缝应力敏感下的缝宽变化方程，建立考虑碳酸盐岩酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度预测模型；
建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果。
在一个优选的实施方式中，在所述步骤收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数中，收集参数具体为：碳酸盐岩初始地层压力、当前地层压力、地层温度、气体摩尔质量、气体分子密度、气体分子碰撞直径、理想气体粘度、酸压自支撑裂缝孔隙度、酸压自支撑裂缝迂曲度基本参数。
在一个优选的实施方式中，在所述采取井下储层岩样，选用酸压用酸液类型，模拟地层温度环境下酸压过程中酸与岩石的反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比步骤中，包括步骤：
将标准岩样放入岩心夹持器中，施加围压，在剪切速率一定下以恒定排量注入酸液；
将经过酸液作用后的标准岩样开展三轴岩石力学测试实验，通过应力-应变曲线获取酸-岩作用前后的岩石弹性模量和泊松比。
在一个优选的实施方式中，在所述步骤根据获取的杨氏模量和泊松比，基于酸压自支撑裂缝压缩性、碳酸盐岩基质压缩性对酸压自支撑裂缝宽度的影响，建立酸压自支撑裂缝应力敏感下的缝宽变化方程中，包括：
基于酸压自支撑裂缝压缩性的影响，建立裂缝宽度变化方程；
基于碳酸盐岩基质压缩性的影响，建立裂缝宽度变化方程；
叠加所述酸压自支撑裂缝压缩性影响、所述碳酸盐岩基质压缩性影响、所述应力敏感影响，得到应力敏感下酸压自支撑裂缝宽度总变化方程。
在一个优选的实施方式中，在步骤建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果中，
碳酸盐岩酸压自支撑裂缝渗透率预测表达式为：
Kt＝Kvs+Kk
其中，Kt为考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率，m2；Kvs为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝滑脱流动渗透率，m2；Kk为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝克努森扩散渗透率，m2。
在一个优选的实施方式中，所述考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝滑脱流动渗透率计算表达式为：
在一个优选的实施方式中，所述考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝克努森扩散渗透率计算表达式为：
综上所述，本发明具有以下优点：改进的技术方案是充分考虑了酸压自支撑裂缝岩石强度、酸压自支撑裂缝中气体流动同时存在连续流以及滑脱流效应、酸压自支撑裂缝微纳米尺度受限空间内气体流动粘度变化，以及酸压井生产过程中有效应力增加导致微裂缝闭合对酸压裂缝渗透率的综合影响，建立了酸压自支撑裂缝渗透率预测模型，该预测模型更加接近实际地层的情况，进而有助于酸压自支撑裂缝渗透率预测的准确性。
附图说明
图1为本发明酸压自支撑裂缝立方网格模型示意图；
图2为本发明理想气体与真实气体黏度在不同微裂缝宽度下随压力的变化情况；
图3为本发明碳酸盐岩在酸压前后岩石的抗压强度变化；
图4为本发明酸压自支撑裂缝在不同杨氏模量下渗透率随压力的变化情况；
图5为本发明酸压自支撑裂在不同泊松比下渗透率随压力的变化情况；
图6为本发明酸压自支撑裂缝在理想气体粘度和真实气体粘度下的渗透率变化情况；
图7为本发明酸压自支撑裂缝在不同压缩系数下的渗透率随压力变化情况。
具体实施方式
结合附图和本发明具体实施方式的描述，能够更加清楚地了解本发明的细节。但是，在此描述的本发明的具体实施方式，仅用于解释本发明的目的，而不能以任何方式理解成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变形，这些都应被视为属于本发明的范围。
为了提高酸压自支撑裂缝渗透率预测的准确性，在本申请中申请人提出了一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，该方法包括以下步骤：
步骤1、收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数。
步骤1中收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数中，收集参数具体为：碳酸盐岩初始地层压力、当前地层压力、地层温度、气体摩尔质量、气体分子密度、气体分子碰撞直径、理想气体粘度、酸压自支撑裂缝孔隙度、酸压自支撑裂缝迂曲度基本参数。
步骤2、采取井下储层岩样，选用酸压用酸液类型，模拟地层温度环境下酸压过程中酸与岩石的反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比。
步骤2中包括以下内容：
步骤2.1、模拟储层环境的酸-岩反应过程
将井下储层岩样加工成直径为2.54cm，长度为5.0cm的标准岩样放入岩芯夹持器中，施加围压6.9MPa，在剪切速率170s-1下以恒定排量3ml/min注入酸液，模拟酸液在裂缝壁面动态流动、反应、滤失过程；其中实验温度设置为地层温度90℃，注入酸液时间为30min。
步骤2.2获取酸-岩反应后的岩石力学参数
采用高温高压三轴岩石力学测试系统测试岩石的三轴力学参数，首先将经过酸液作用后的岩样用热塑性塑料套住，然后将装有试件的保护筒放入压力室中，按0.05MPa/s的加载速率用油向岩样加围压及孔压，然后按0.5MPa/s加轴向应变使岩样破裂，通过应力-应变曲线计算酸-岩作用前后的岩石弹性模量和泊松比。
步骤3、根据获取的杨氏模量和泊松比，基于酸压自支撑裂缝压缩性、碳酸盐岩基质压缩性对酸压自支撑裂缝宽度的影响，建立酸压自支撑裂缝应力敏感下的缝宽变化方程。
步骤3中包括以下内容：
步骤3.1、酸压自支撑裂缝压缩性的影响：
为了研究地层压力变化引起的应力敏感对碳酸盐岩自支撑裂缝宽度和碳酸盐岩基质压缩等方面的影响，将碳酸盐岩储层进行网格化处理，得到若干个基质裂缝立方网格(见图1(i))；在网格中，基质单元是长度为a的正方体，基质之间是宽度为b的自支撑裂缝，最终得到理想化的碳酸盐岩储层立方网格模型。假设碳酸盐岩基质中孔隙极小，储层渗流空间全由自支撑裂缝提供，这里重点针对基质之间的自支撑裂缝进行讨论。该模型中，碳酸盐岩自支撑裂缝宽度b、自支撑裂缝高度a、自支撑裂缝孔隙度φ满足下式关系：
式中，φ为碳酸盐岩自支撑裂缝孔隙度，无量纲；a为碳酸盐岩自支撑裂缝高度，m；b为碳酸盐岩自支撑裂缝宽度，m；
根据定义，地层净压力由下式表示：
σ＝pob-pp (2)
式中，σ为净压力，MPa；pob为碳酸盐岩上覆岩石压力，MPa；pp为碳酸盐岩地层孔隙压力，MPa；
当孔隙压力和上覆岩石压力变化时，
Δσ＝(pob-pob0)-(pp-pp0) (3)
式中，Δσ为净压力改变量，MPa；pob0为碳酸盐岩初始上覆岩石压力，MPa；pp0为碳酸盐岩地层初始孔隙压力，MPa；
酸压自支撑裂缝压缩系数可表示为：
式中，cf为酸压自支撑裂缝压缩系数，MPa-1；Δφ为酸压自支撑裂缝孔隙度变化量，无量纲；φ0为酸压自支撑裂缝初始孔隙度，无量纲；
在地层上覆岩石压力不变的情况下，酸压自支撑裂缝压缩系数也可表达为：
式中，c0为酸压自支撑裂缝初始压缩系数，Pa-1；δ为酸压自支撑裂缝压缩系数变化速率，Pa-1；
结合式(1)、式(3)、式(4)，酸压自支撑裂缝的压缩性引起的裂缝宽度变化为：
Δbf＝-b0cf[(pob-pob0)-(pp-pp0)] (6)
式中，Δbf为酸压自支撑裂缝的压缩性引起的裂缝宽度变化，m；b0为酸压自支撑裂缝初始宽度，m；
步骤3.2、碳酸盐岩基质压缩性的影响：
碳酸盐岩体积模量由下式表示：
式中，K为碳酸盐岩体积模量，MPa；Δphydrostatic为地层静水压力，MPa；ΔV为碳酸盐岩体积改变量，m3；V0为碳酸盐岩初始体积，m3；
由于地层孔隙压力变化对碳酸盐岩基质压缩的影响远大于上覆岩石压力对基质压缩的影响；因此，式(7)可表达为：
式中，V为碳酸盐岩体积，m3；
碳酸盐岩体积机械应变可表达为：
式中，-εvm为碳酸盐岩体积机械应变，无量纲；
碳酸盐岩体积模量由岩石泊松比和杨氏模量表示：
式中，E为碳酸盐岩杨氏模量，MPa；υ为碳酸盐岩泊松比，无量纲；
由式(9)、式(10)得：
碳酸盐岩线性应变为体积机械应变的三分之一，故其线性应变表示为
式中，εlm为碳酸盐岩线性应变，无量纲；
碳酸盐岩线性应变定义为：
式中，Δam为由碳酸盐岩基质压缩性引起的基质长度变化量，m；a0为基质初始长度，m；
碳酸盐岩基质长度变化量与酸压自支撑裂缝宽度变化量相等，但变化趋势相反：
Δam＝-Δbm (14)
式中，Δbm为由碳酸盐岩基质压缩性引起酸压自支撑裂缝宽度变化量，m；
由式(12)、式(13)、式(14)得出在地层压力变化下，由碳酸盐岩基质压缩性引起的酸压自支撑裂缝宽度变化Δbm表示为：
步骤3.3、地层应力敏感下酸压自支撑裂缝宽度总变化方程：
综合式(6)、(15)得由压力变化导致的酸压自支撑裂缝宽度总变化量Δbt：
式中，Δbt为由压力变化导致的酸压自支撑裂缝宽度总变化量，m。
步骤4、基于酸压自支撑裂缝应力敏感下的缝宽变化方程，建立考虑碳酸盐岩酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度预测模型。
对微纳米受限空间内的气体，由于稀薄性的增加使得其表现出来的粘度和宏观气体的粘度有很大不同，其具体表现为有效粘度。
气体平均分子自由程表达为：
式中，λ为气体平均分子自由程，m；kB为玻尔兹曼常数，J/K；T为地层温度，K；Γ为气体分子碰撞直径，m；p为地层压力，Pa；
酸压自支撑裂缝中的气体克努森数为气体平均分子自由程与自支撑裂缝宽度的比值：
式中，Kn为气体克努森数，无因次；
考虑酸压自支撑裂缝宽度变化的气体克努森数表达为：
式中，Knb为考虑酸压自支撑裂缝宽度变化的气体克努森数，无因次；
计算气体有效粘度的表达式为：
式中，μeff0为气体有效粘度，Pa·s；μ为理想气体粘度，Pa·s；
考虑酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度表达式为：
式中，μeff为考虑酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度，Pa·s。
步骤5、建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果。
在立方网格模型基础上(见图1(i))，通过进一步简化，得到缝长、缝高均为a，缝宽为b的长方体酸压自支撑裂缝简化模型(见图1(ii))；考虑酸压自支撑裂缝有不同形状，定义长方体自支撑裂缝简化模型的截面几何参数为：
式中，ζ为酸压自支撑裂缝纵横比，无因次；
考虑缝宽变化下酸压自支撑裂缝纵横比为：
式中，ζb为考虑缝宽变化的酸压自支撑裂缝纵横比，无因次；当气体克努森数为10-3<Knb<10-1，气体分子间碰撞和气体与酸压自支撑裂缝壁面的碰撞均不可忽略，气体分子在壁面发生滑脱，通过修正滑脱边界条件，则气体滑脱流动可表达为：
式中，Jvs为考虑气体粘度变化和应力敏感效应的酸压自支撑裂缝气体滑脱流动流量，kg·s-1；A(ζb)为滑脱流动的裂缝截面形状因子，无因次；M为气体摩尔质量，kg·mol-1；τ为自支撑裂缝迂曲度，无因次；R为普适气体常数，J·mol-1K-1；α为稀有效应系数，无因次；m为气体滑脱常数，无因次；
式中，tanh(x)为双曲正切函数，i为计数符号；
当克努森数Knb≥10，气体分子与壁面碰撞占主导地位，为克努森扩散，则气体扩散量用克努森方程表达为：
式中，Jk为考虑气体粘度变化和应力敏感效应的酸压自支撑裂缝气体克努森流动流量，kg·s-1。
式中，B(ζb)为克努森流动的裂缝形状因子，无因次；
酸压自支撑裂缝中气体分子之间碰撞频率为：
式中，fm为气体分子之间碰撞频率，s-1；v为气体平均热运动速度，m·s-1；n为单位体积气体分子数，m-3；
酸压自支撑裂缝中气体分子与壁面碰撞频率为：
式中，fw为气体分子与自支撑裂缝壁面碰撞频率，s-1；
滑脱流动权重系数为气体分子之间碰撞频率占总碰撞频率的比值：
式中，wvs为气体分子之间碰撞频率占总碰撞频率的比值，无因次；
克努森流动权重系数为气体分子与自支撑壁面碰撞频率占总碰撞频率的比值：
式中，wk为气体分子与壁面碰撞频率占总碰撞频率的比值，无因次；
各传输机理对碳酸盐岩酸压自支撑裂缝渗透率的贡献可分别表达为：
式中，Kvs为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝滑脱流动渗透率，m2；Kk为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝克努森扩散渗透率，m2；Vstd为标准状态下气体摩尔体积，m3·mo1-1；
考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率可表达为：
Kt＝Kvs+Kk (34)
式中，Kt为考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率，m2。
计算实例
采用的模拟基础参数如表1所示。
表1模拟基础数据表
(2)计算结果
现有技术在计算自支撑裂缝渗透率时，未考虑岩石力学参数变化、真实气体粘度变化及应力敏感效应对渗透率的影响。而本发明则能综合考虑以上因素，得出更为符合实际情况的结果。
由图2可知，随着地层压力降低，气体粘度变化逐渐增大，在低地层压下(小于5MPa)，气体粘度迅速减小，酸压自支撑裂缝初始宽度越小，粘度随压力变化越明显。
由图3可以看出，酸蚀前后岩芯的应力应变特征发生一定变化，由于酸液的溶蚀使岩芯孔隙空间增大，曲线初期的压实致密过程持续时间更长。同时可以看出酸蚀后岩芯表现出更强的塑性，发生形变的趋势增加。根据应力应变曲线，酸处理前后的杨氏模量分别为3.0×1010Pa和2.0×1010Pa，泊松比分别为0.25和0.35。酸处理前后岩石抗压强度从169MPa变为155MPa。
酸压过程中酸液对岩石的溶蚀将显著降低岩石杨氏模量，同时增大岩石泊松比，进而对酸压自支撑裂缝渗透率造成不同的影响。如果不考虑岩石力学参数的变化，将使得渗透率预测结果产生较大误差。由图4和图5可知，碳酸盐岩力学参数(泊松比和杨氏模量)对酸压自支撑裂缝渗透率有较大的影响，在中等地层压力(5MPa)附近，岩石力学参数对自支撑裂缝渗透率影响最为显著。
由图6可知，地层压力较高时，由于气体粘度变化较小，因而两种情况下(考虑粘度变化和不考虑粘度变化)酸压自支撑裂缝渗透率无明显差别。在低地层压力下(小于5MPa)，由于气体粘度急剧减小，所以两种情况下的自支撑裂缝渗透率呈现明显差别，这种趋势在自支撑裂缝初始宽度较小的情况下更加明显。
由图7可知，自支撑裂缝压缩性对裂缝渗透率影响很大，在低地层压力下(小于5MPa)，不同裂缝压缩性对裂缝渗透率影响逐渐减小。其它条件相同时，自支撑裂缝压缩性越大(c0越大、δ越小)，裂缝越容易被压缩，裂缝渗透率的值越小。总体而言，酸压自支撑裂缝渗透率与其裂缝压缩性呈负相关。
本发明具有以下优点：改进的技术方案是充分考虑了酸压自支撑裂缝岩石强度、酸压自支撑裂缝中气体流动同时存在连续流以及滑脱流效应、酸压自支撑裂缝微纳米尺度受限空间内气体流动粘度变化，以及酸压井生产过程中有效应力增加导致微裂缝闭合对酸压裂缝渗透率的综合影响，建立了酸压自支撑裂缝渗透率预测模型，该预测模型更加接近实际地层的情况，进而有助于酸压自支撑裂缝渗透率预测的准确性。
虽然结合附图对本发明的具体实施方式进行了详细地描述，但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内，本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。
Claims (7)
1.一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，其包括以下步骤：
收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数；
采取井下储层岩样，选用酸压用酸液类型，模拟地层温度环境下酸压过程中酸与岩石的反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比；
根据获取的杨氏模量和泊松比，基于酸压自支撑裂缝压缩性、碳酸盐岩基质压缩性对酸压自支撑裂缝宽度的影响，建立酸压自支撑裂缝应力敏感下的缝宽变化方程；
基于酸压自支撑裂缝应力敏感下的缝宽变化方程，建立考虑碳酸盐岩酸压自支撑裂缝应力敏感效应和受限空间效应的气体有效粘度预测模型；
建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果。
2.如权利要求1所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，在所述步骤收集碳酸盐岩酸压自支撑裂缝初始尺寸、地层压力参数、气体性质参数中，收集参数具体为：
碳酸盐岩初始地层压力、当前地层压力、地层温度、气体摩尔质量、气体分子密度、气体分子碰撞直径、理想气体粘度、酸压自支撑裂缝孔隙度、酸压自支撑裂缝迂曲度。
3.如权利要求1所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，在所述采取井下储层岩样，选用酸压用酸液类型，模拟地层温度环境下酸压过程中酸与岩石的反应过程，获取酸溶蚀碳酸盐岩后的杨氏模量和泊松比步骤中，包括步骤：
将标准岩样放入岩心夹持器中，施加围压，在剪切速率一定下以恒定排量注入酸液；
将经过酸液作用后的标准岩样开展三轴岩石力学测试实验，通过应力-应变曲线获取酸-岩作用前后的岩石弹性模量和泊松比。
4.如权利要求1所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，在所述步骤根据获取的杨氏模量和泊松比，基于酸压自支撑裂缝压缩性、碳酸盐岩基质压缩性对酸压自支撑裂缝宽度的影响，建立酸压自支撑裂缝应力敏感下的缝宽变化方程中，包括：
基于酸压自支撑裂缝压缩性的影响，建立裂缝宽度变化方程；
基于碳酸盐岩基质压缩性的影响，建立裂缝宽度变化方程；
叠加所述酸压自支撑裂缝压缩性影响、所述碳酸盐岩基质压缩性影响、所述应力敏感影响，得到应力敏感下酸压自支撑裂缝宽度总变化方程。
5.如权利要求1所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，在步骤建立考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率动态预测模型，获得碳酸盐岩酸压自支撑裂缝渗透率预测结果中，
碳酸盐岩酸压自支撑裂缝渗透率预测表达式为：
Kt＝Kvs+Kk
其中，Kt为考虑应力敏感效应以及气体在自支撑裂缝微纳米尺度受限空间内粘度变化影响的碳酸盐岩酸压自支撑裂缝渗透率，m2；Kvs为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝滑脱流动渗透率，m2；Kk为考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝克努森扩散渗透率，m2。
6.如权利要求5所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，所述考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝滑脱流动渗透率计算表达式为：
7.如权利要求5所述的一种碳酸盐岩酸压自支撑裂缝渗透率预测方法，所述考虑应力敏感效应以及气体粘度变化的碳酸盐岩酸压自支撑裂缝克努森扩散渗透率计算表达式为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811601635.1A CN109653740B (zh) | 2018-12-26 | 2018-12-26 | 一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811601635.1A CN109653740B (zh) | 2018-12-26 | 2018-12-26 | 一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109653740A true CN109653740A (zh) | 2019-04-19 |
| CN109653740B CN109653740B (zh) | 2020-11-10 |
Family
ID=66116663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811601635.1A Active CN109653740B (zh) | 2018-12-26 | 2018-12-26 | 一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109653740B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
| CN111255445A (zh) * | 2020-01-17 | 2020-06-09 | 西南石油大学 | 获取评价模型的方法与系统及评价方法与评价系统 |
| CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
| CN112347411A (zh) * | 2020-10-30 | 2021-02-09 | 西南石油大学 | 一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法 |
| CN114117791A (zh) * | 2021-11-26 | 2022-03-01 | 西安石油大学 | 一种碳酸盐岩酸化压裂数值模拟方法 |
| CN114166726A (zh) * | 2022-02-14 | 2022-03-11 | 西南石油大学 | 一种岩心渗透率张量敏感性损伤测试设备与评价方法 |
| CN114441423A (zh) * | 2020-10-20 | 2022-05-06 | 中国石油化工股份有限公司 | 一种预测碳酸盐岩自支撑酸蚀裂缝导流能力的方法 |
| CN115032368A (zh) * | 2022-06-07 | 2022-09-09 | 西南石油大学 | 一种压裂裂缝自支撑导流能力全过程评价方法 |
| CN117077575A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102748017A (zh) * | 2012-07-11 | 2012-10-24 | 西南石油大学 | 用于油气田开发的自支撑裂缝测试分析装置及方法 |
| CN105041288A (zh) * | 2015-07-13 | 2015-11-11 | 中国石油大学(北京) | 一种碳酸盐岩油气藏裂缝转向酸压方法 |
| CN106290784A (zh) * | 2016-08-02 | 2017-01-04 | 中国石油大学(华东) | 一种酸岩反应单岩板裂缝模拟装置及其工作方法 |
| CN107066769A (zh) * | 2017-06-09 | 2017-08-18 | 西南石油大学 | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 |
| CN107191173A (zh) * | 2016-03-14 | 2017-09-22 | 中国石油化工股份有限公司 | 一种碳酸盐岩的酸压效果的评价方法 |
| CN108204231A (zh) * | 2016-12-20 | 2018-06-26 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层体积酸压导流能力的评价方法 |
| CN108343433A (zh) * | 2018-02-28 | 2018-07-31 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
| CN108661619A (zh) * | 2018-05-10 | 2018-10-16 | 西南石油大学 | 一种确定碳酸盐岩油气藏酸蚀填砂裂缝铺砂浓度的方法 |
| CN108979612A (zh) * | 2018-08-09 | 2018-12-11 | 西南石油大学 | 一种致密油气藏压裂酸化复杂裂缝流动能力优化方法 |
| CN109033674A (zh) * | 2018-08-09 | 2018-12-18 | 西南石油大学 | 一种靶向压裂酸化井裂缝参数优化方法 |
-
2018
- 2018-12-26 CN CN201811601635.1A patent/CN109653740B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102748017A (zh) * | 2012-07-11 | 2012-10-24 | 西南石油大学 | 用于油气田开发的自支撑裂缝测试分析装置及方法 |
| CN105041288A (zh) * | 2015-07-13 | 2015-11-11 | 中国石油大学(北京) | 一种碳酸盐岩油气藏裂缝转向酸压方法 |
| CN107191173A (zh) * | 2016-03-14 | 2017-09-22 | 中国石油化工股份有限公司 | 一种碳酸盐岩的酸压效果的评价方法 |
| CN106290784A (zh) * | 2016-08-02 | 2017-01-04 | 中国石油大学(华东) | 一种酸岩反应单岩板裂缝模拟装置及其工作方法 |
| CN108204231A (zh) * | 2016-12-20 | 2018-06-26 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层体积酸压导流能力的评价方法 |
| CN107066769A (zh) * | 2017-06-09 | 2017-08-18 | 西南石油大学 | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 |
| CN108343433A (zh) * | 2018-02-28 | 2018-07-31 | 西南石油大学 | 用于计算页岩微裂缝缝宽变化下气体质量传输参数的方法 |
| CN108661619A (zh) * | 2018-05-10 | 2018-10-16 | 西南石油大学 | 一种确定碳酸盐岩油气藏酸蚀填砂裂缝铺砂浓度的方法 |
| CN108979612A (zh) * | 2018-08-09 | 2018-12-11 | 西南石油大学 | 一种致密油气藏压裂酸化复杂裂缝流动能力优化方法 |
| CN109033674A (zh) * | 2018-08-09 | 2018-12-18 | 西南石油大学 | 一种靶向压裂酸化井裂缝参数优化方法 |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110609974A (zh) * | 2019-09-09 | 2019-12-24 | 西南石油大学 | 一种考虑蚓孔扩展的酸压裂缝动态滤失计算方法 |
| CN111255445A (zh) * | 2020-01-17 | 2020-06-09 | 西南石油大学 | 获取评价模型的方法与系统及评价方法与评价系统 |
| CN111255445B (zh) * | 2020-01-17 | 2022-06-17 | 西南石油大学 | 获取评价模型的方法与系统及评价方法与评价系统 |
| CN111963158B (zh) * | 2020-08-12 | 2022-04-01 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
| CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
| CN114441423A (zh) * | 2020-10-20 | 2022-05-06 | 中国石油化工股份有限公司 | 一种预测碳酸盐岩自支撑酸蚀裂缝导流能力的方法 |
| CN112347411A (zh) * | 2020-10-30 | 2021-02-09 | 西南石油大学 | 一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法 |
| CN114117791A (zh) * | 2021-11-26 | 2022-03-01 | 西安石油大学 | 一种碳酸盐岩酸化压裂数值模拟方法 |
| CN114166726A (zh) * | 2022-02-14 | 2022-03-11 | 西南石油大学 | 一种岩心渗透率张量敏感性损伤测试设备与评价方法 |
| US11604132B1 (en) | 2022-02-14 | 2023-03-14 | Southwest Petroleum University | Testing device and evaluation method for sensitivity damage of core permeability tensor |
| CN115032368A (zh) * | 2022-06-07 | 2022-09-09 | 西南石油大学 | 一种压裂裂缝自支撑导流能力全过程评价方法 |
| CN115032368B (zh) * | 2022-06-07 | 2023-08-04 | 西南石油大学 | 一种压裂裂缝自支撑导流能力全过程评价方法 |
| CN117077575A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置 |
| CN117077575B (zh) * | 2023-10-16 | 2024-02-23 | 西安石油大学 | 一种碳酸盐岩油藏酸压酸液滞留时间计算方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109653740B (zh) | 2020-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109653740A (zh) | 一种碳酸盐岩酸压自支撑裂缝渗透率预测方法 | |
| Tan et al. | Hydraulic fracture vertical propagation behavior in transversely isotropic layered shale formation with transition zone using XFEM-based CZM method | |
| Macini et al. | Laboratory measurements of non-Darcy flow coefficients in natural and artificial unconsolidated porous media | |
| Lacy et al. | Fracture width and embedment testing in soft reservoir sandstone | |
| Senetakis et al. | The small-strain shear modulus and damping ratio of quartz and volcanic sands | |
| CN106501144B (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| Hu et al. | Mechanisms of fracturing fluid spontaneous imbibition behavior in shale reservoir: A review | |
| Nabawy et al. | Pore-throat characterization in highly porous and permeable sandstones | |
| CN108590601B (zh) | 一种注水扩容膨胀施工参数优选的实验方法 | |
| CN107806339B (zh) | 一种压裂裂缝导流能力实验方法 | |
| CN108386176B (zh) | 一种天然裂缝与人工裂缝延伸规律物模试验方法 | |
| Uchida et al. | Geomechanical study of the Mallik gas hydrate production field trials | |
| CN108399270A (zh) | 一种确定页岩地层中各向异性泥质比例的方法 | |
| CN106290714A (zh) | 一种裂缝性储层物理模拟方法 | |
| Zhang et al. | A study on effect of seepage direction on permeability stress test | |
| Xu et al. | Reducing residual oil saturation: underlying mechanism of imbibition in oil recovery enhancement of tight reservoir | |
| CN109459795A (zh) | 一种孔隙性砂岩断裂带封闭性评价方法 | |
| CN110159251A (zh) | 一种试验预测页岩水力压裂效果的方法 | |
| CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
| CN110295894A (zh) | 一种建立水平井产能预测模型的方法 | |
| CN106545335B (zh) | 致密油藏衰竭式开发合理生产压差的确定方法 | |
| Gao et al. | Effect of pressure pulse stimulation on imbibition displacement within a tight sandstone reservoir with local variations in porosity | |
| Shen et al. | An experimental investigation of the nonlinear gas flow and stress‐dependent permeability of shale fractures | |
| Cruichshank et al. | Anchoring methodologies for pore-scale network models: application to relative permeability and capillary pressure prediction | |
| CN106442269A (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 |