CN108590640A - 一种复杂裂缝网络渗透率计算方法 - Google Patents
一种复杂裂缝网络渗透率计算方法 Download PDFInfo
- Publication number
- CN108590640A CN108590640A CN201810312664.XA CN201810312664A CN108590640A CN 108590640 A CN108590640 A CN 108590640A CN 201810312664 A CN201810312664 A CN 201810312664A CN 108590640 A CN108590640 A CN 108590640A
- Authority
- CN
- China
- Prior art keywords
- crack
- fracture network
- fracture
- permeability
- dimensionless
- 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
- 230000035515 penetration Effects 0.000 title claims abstract description 27
- 238000000205 computational method Methods 0.000 title claims abstract description 12
- 230000035699 permeability Effects 0.000 claims abstract description 77
- 239000011435 rock Substances 0.000 claims abstract description 50
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000002591 computed tomography Methods 0.000 claims abstract description 8
- 238000012512 characterization method Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 230000006698 induction Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 3
- 238000012876 topography Methods 0.000 claims description 3
- 229910017435 S2 In Inorganic materials 0.000 claims 1
- 230000000877 morphologic effect Effects 0.000 abstract description 6
- 238000011161 development Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 206010017076 Fracture Diseases 0.000 description 103
- 208000010392 Bone Fractures Diseases 0.000 description 101
- 230000008859 change Effects 0.000 description 11
- 239000007789 gas Substances 0.000 description 10
- 238000010008 shearing Methods 0.000 description 10
- 238000012360 testing method Methods 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000001764 infiltration Methods 0.000 description 3
- 230000008595 infiltration Effects 0.000 description 3
- 208000002565 Open Fractures Diseases 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 238000009738 saturating Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000009918 complex formation Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000000007 visual effect Effects 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
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
本发明属于非常规油气开发技术领域,具体而言,涉及一种复杂裂缝网络渗透率计算方法,包括以下步骤:S1、从页岩储层段钻取岩心,并将钻取的岩心制成直径为2.5cm、长度为5.0cm的标准岩心;S2、对岩心端面进行CT扫描,扫描后图像进行灰度处理,识别出裂缝条数、位置及形态特征;S3、基于S2中的图像,运用盒子法来获取表征岩石复杂裂缝网络的分形维数;S4、基于S3中的分形维数,利用蒙特卡洛随机建模理论生成二维离散裂缝网络模型;S5、基于二维离散裂缝网络模型,通过改进的立方定律来计算复杂裂缝网络渗透率。本发明的有益效果:通过扫描获得分形维数,利用蒙特卡洛随机建模,考虑地层实际因素,计算复杂缝网渗透率,计算结果更加接近储层真实情况。
Description
技术领域
本发明涉及一种复杂裂缝网络渗透率计算方法,属于非常规油气开发技术领域。
背景技术
21世纪经济快速发展,能源需求不断攀升,煤、石油、天然气等传统能源的消耗也逐步增加,能源压力日益增大。在能源需求的新形势下,作为常规油气能源的有效补充,页岩气以非常规天然气的形式进入人们的视野。近年来,页岩气凭借其巨大的资源潜力与经济价值而备受关注,并且随着水平井钻井技术和体积压裂技术的发展,日益成为油气资源勘探开发的重要领域和目标。
页岩气储集层在成岩及烃类演化过程中,容易在有机质边缘或薄弱面处形成规则、不规则天然裂缝,大量发育的天然裂缝对页岩气质量传输影响显著。而在页岩气开采过程中,都会用到压裂施工,由于页岩储层本身具有较为发育的天然裂缝,并且页岩由石英等脆性矿物组成,在压裂后改造区域会形成人工裂缝与天然裂缝相互交叉的复杂裂缝网络。为了能够准确模拟出页岩气压裂施工后的产能,关键在于建立符合现场实际的页岩气复杂裂缝网络,并且准确计算出其渗透率。而如何准确获取复杂裂缝网络渗透率是一个亟待解决的难题。
目前国内外有大量学者开展了关于裂缝网络建模的相关研究,对于流体在裂缝中的流动主要有连续体模型、离散裂缝网络模型以及两者的混合模型。针对页岩压裂后形成复杂裂缝问题的特殊性,通常选择离散裂缝网络模型进行研究。
分形维数在分析岩石裂缝中应用较为广泛,通常情况下分形维数量化了曲线或者曲面在一定范围内填充空间的程度,可以更接近实际地来描述裂缝的分布及其形态特征。目前,在离散裂缝网络建模过程中对页岩压裂后裂缝形态的表征有所欠缺,大多数的裂缝形态主要是根据经验参数取得分形维数的大小。同样地,在复杂裂缝网络渗透率计算时,目前存在不能反映真实储层情况,不能同时考虑张、剪切裂缝以及围压等因素的缺点。
本发明模型通过分形维数生成的裂缝网络更加接近实际地层条件下裂缝的情况,且考虑裂缝密度、裂缝开度、围压,以及张开裂缝和剪切裂缝同时存在的多种地层条件下的因素,通过在在裂缝复杂网络生成的基础上,利用改进的立方定律计算出渗透率也因此更加接近实际地层渗透率。
发明内容
本发明提供了一种复杂裂缝网络渗透率计算方法,其目的在于,解决现有技术中存在的上述问题。
本发明的技术方案如下:
一种复杂裂缝网络渗透率计算方法,其特征在于,包括以下步骤:
S1、从页岩储层段钻取岩心,并将钻取的岩心制成直径为2.5cm、长度为5.0cm的标准岩心;
S2、对岩心端面进行CT扫描,扫描后图像进行灰度处理,识别出裂缝条数、位置及形态特征;
S3、基于S2中的图像,运用盒子法来获取表征岩石复杂裂缝网络的分形维数;
S4、基于S3中的分形维数,利用蒙特卡洛随机建模理论建立二维离散裂缝网络模型;
S5、基于二维离散裂缝网络模型,通过改进的立方定律求解复杂裂缝网络渗透率。
进一步地,所述步骤S2中对岩心进行CT扫描即利用形貌扫描仪拍下岩心端面的裂缝细节;利用编制的程序,将岩石裂缝面间距图进行灰度处理,并将灰度图处理成二值图,完成数字化处理。
进一步地,所述步骤S3的盒子法是用正方形格子去覆盖裂缝间距图;给定盒子尺码δ,可以计算出覆盖住裂缝间距所需的总盒子数目N;运用盒子法求得的分形维数表达式为:
式中:D为分形维数,无量纲;N为盒子数目,无量纲;δ为盒子尺码,无量纲;i为计数符号。
进一步的,所述步骤S4的二维离散裂缝网络模型的建模方法为蒙特卡洛随机建模方法,所述二维离散裂缝网络模型的生成由以下步骤组成:A、模拟区域生成;B、裂缝数目确定;C、裂缝位置表征;D、裂缝几何参数;E、裂缝网络模型生成。
进一步的,所述步骤S5中复杂裂缝网络渗透率包括正常张开型裂缝渗透率、剪切型裂缝渗透率,计算渗透率过程中,由张力形成的张开型和由剪应力形成的剪切型裂缝渗透率分别进行计算,最后将正应力和剪应力诱导产生裂缝所具有渗透率进行叠加,所述正应力和剪应力诱导产生裂缝所具有渗透率的计算公式如下:
式中:kx为复杂缝网在x方向上的渗透率,10-3μm2;ky为复杂缝网在y方向上的渗透率,10-3μm2;fx为在x方向上正应力产生的裂缝等效频率,1/m;fy为在y方向上正应力产生的裂缝的等效频率,1/m;br为裂缝的初始平均开度,μm;bmax为受压下裂缝最大开度,μm;αx为x方向上的压力系数,无量纲;αy为y方向上的压力系数,无量纲;σx为x方向上正应力大小,MPa;σy为y方向上正应力大小,MPa;βx为在xy下的函数斜率,无量纲;βy为在yx下的函数斜率,无量纲;fdx为在x方向上剪应力产生的裂缝等效频率,1/m;fdy为在y方向上剪应力产生的裂缝等效频率,1/m;dmax为裂缝发生剪切破坏后的最大开度,μm;k—x与y方向应力比,无量纲;γx——x方向剪切压力系数,无量纲;γy——y方向剪切压力系数,无量纲;kc——裂缝变形的临界压力比,无量纲;kc为裂缝变形的临界压力比,无量纲。
本发明的有益效果为:
本发明通过在钻取地层实际岩心的基础上,通过扫描获取分形维数,并根据分形维数利用蒙特卡洛随机建模理论生成裂缝网络模拟地层复杂缝网情况;在此基础上,通过改进的立方定律定量计算复杂缝网渗透率。该方法充分考虑了裂缝密度、裂缝开度、围压,以及张开裂缝和剪切裂缝同时存在等因素,通过该方法计算出的复杂缝网渗透率也更加接近储层真实情况。
附图说明
为了更清楚地说明本发明实施方式的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它相关的附图。
图1本发明提供的流程图;
图2本发明提供的1号岩心扫描图;
图3本发明提供的2号岩心扫描图;
图4本发明提供的3号岩心扫描图;
图5本发明提供的模型在x方向上渗透率与岩心实验测试渗透率的对比图;
图6本发明提供的模型在y方向上渗透率与岩心实验测试渗透率的对比图;
图7本发明提供的模型计算得到的渗透率与岩心实验测试渗透率的对比图;
图8本发明提供的模型计算得到的不同水平主应力下,x方向上渗透率的变化图;
图9本发明提供的模型计算得到的不同水平主应力下,y方向上渗透率的变化图;
图10本发明提供的模型计算得到的不同水平主应力下,复杂裂缝网络渗透率的变化图;
图11本发明提供的模型计算得到的不同分形维数下,复杂裂缝网络渗透率的变化图;
图12本发明提供的模型计算得到的不同面密度,复杂裂缝网络渗透率的变化图;
图13本发明提供的模型计算得到的不同裂缝开度,复杂裂缝网络渗透率的变化图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施方式,都属于本发明保护的范围。
如图1所示,本发明提供的一种复杂裂缝网络渗透率计算方法,包括以下步骤:
S1、从页岩储层段钻取岩心,并将钻取的岩心制成直径为2.5cm、长度为5.0cm的标准岩心;
S2、对岩心端面进行CT扫描,扫描后图像进行灰度处理,识别出裂缝条数、位置及形态特征;
S3、基于S2中的图像,运用盒子法来获取表征岩石复杂裂缝网络的分形维数;
S4、基于S3中的分形维数,利用蒙特卡洛随机建模理论建立二维离散裂缝网络模型;
S5、基于二维离散裂缝网络模型,通过改进的立方定律求解复杂裂缝网络渗透率。
在本发明提供的实施例中,所述步骤S2中对岩心进行CT扫描即利用形貌扫描仪拍下岩心端面的裂缝细节;利用编制的程序,将岩石裂缝面间距图进行灰度处理,并将灰度图处理成二值图,完成数字化处理。
在本发明提供的实施例中,所述步骤S3的分形维数在分析岩石裂缝中应用较为广泛,可以更接近实际来描述裂缝的分布及其形态特征。盒子法是用正方形的格子去覆盖岩石裂缝间距图,盒子大小是变化的;给定盒子尺码δ,可以数出覆盖住岩石裂缝间距所需的总盒子数目N;假设第i步覆盖使用δi×δi的格子,所需盒子数目为Ni(δi),在第i+1步需要的格子(δi+1)×(δi+1),则需要盒子数目Ni+1(δi+1)。得到任意两个尺度下所需盒子数之比与尺码之比的关系如下:
推广到:
N=aδ-D (2)
得到分形维数表达式为:
式中:D为分形维数,无量纲;N为盒子数目,无量纲;δ为盒子尺码,无量纲;i为计数符号。
在本发明提供的实施例中,所述步骤S4的二维离散裂缝网络模型的建模方法为随机建模方法,所述二维离散裂缝网络模型的生成由以下步骤组成:
A、模拟区域生成:模拟区域可以根据地层中裂缝系统分布范围或者根据研究需要来选择。裂缝网络模拟时,一般先要生成裂中心点。考虑到有裂缝只有部分在模拟区域内,中心点位置可能在模拟区域外。因此,为减少边界效应的影响,应扩大模拟区域。为了便于区别,扩大后的模拟区域称为生成域,而原来模拟区域称为分析区。为保证模拟结果可靠,生成域应该足够大,原则上每一个侧面距分析域的距离应大于裂缝可能出现的最大半径。
B、裂缝数目确定:模拟区域的大小确定后,就可以根据裂缝的密度来确定该区域内每组裂缝的数目。若迹长服从指数分布,裂缝的面密度可由下式确定。
C、裂缝空间位置表征:对于二维离散裂缝网络模型,裂缝的位置可以由中心点及迹长,或者裂缝的两个端点完全确定,即
式中:θ为裂缝与x轴的夹角,°;l——裂缝长度,m;x,y——生成裂缝坐标随机变量,m;x0,y0——生成裂缝中心点坐标随机变量,m。
D、裂缝几何参数:a、按均匀分布生成裂缝中心点坐标随机变量x0,y0,做二维模拟时只需生成其中的两个变量;b、生成裂缝的倾向及倾角,生成二维离散裂缝模型时,根据研究目的来选择其中一个变量;c、生成裂缝的迹长;d、生成裂缝宽度,通常把每组裂缝的宽度设为一常量。
E、裂缝网络模型生成及简化。当各组裂缝的数目、中心点坐标、裂缝产状、裂缝半径及缝宽确定后,生成域内裂缝也就完全确定,生成域裂缝生成后,根据分析几何边界来生成裂缝网络模型。
在本发明提供的实施例中,所述步骤S5中复杂裂缝网络渗透率包括正常张开型裂缝渗透率、剪切型裂缝渗透率,计算渗透率过程中,正应力和剪应力分别进行计算,最后将正应力和剪应力诱导产生裂缝所具有渗透率进行叠加,表达式如下:
kx=knx+kdx (5)
ky=kny+kdy (6)
式中:kx——在x方向上复杂缝网的渗透率,10-3μm2;ky——在y方向上复杂缝网的渗透率,10-3μm2;knx——在x方向上由于裂缝正常闭合后的渗透率,10-3μm2;kny——在y方向上由于裂缝正常闭合后的渗透率,10-3μm2;kdx——在x方向上由于裂缝发生剪切破坏产生的渗透率,10-3μm2;kdy——在y方向上由于裂缝发生剪切破坏产生的渗透率,10-3μm2。
在离散裂缝网络中,由于单条裂缝并不能反映出不同长度不同开度的裂缝,因此,在发明中使用等效裂缝频率(fx,fy)和等效开度(bx,by)的概念,使用这种形式来相对接近的表述离散裂缝网络中的复杂裂缝变形。通过等效频率和等效开度,正应力产生的张开或者闭合裂缝在x和y方向上(knx,kny)的渗透率可以通过下式计算所得:
等效开度是指用一个值来近似代表整个离散裂缝网络所有裂缝的宽度,然后根据式(7、8)的立方定律来表达离散裂缝网络的水力特性。其中裂缝的等效频率可以通过数值实验和初始开度反推得到:
式中:bi为裂缝的初始平均开度,μm;
裂缝的开度与地应力之间的关系,可以由如下公式表示:
b=br+bm=br+bmaxexp(-α'σ) (11)
式中:br为裂缝的初始平均开度,μm;bmax为受压下裂缝最大开度,μm;α'为正常闭合下正常剪切下的函数斜率,无量纲;σ为正应力大小,MPa。
根据式(11)可以从一维的公式推导到二维条件,在x方向上和在y方向上裂缝开度的大小:
bx=br+bmax=br+bmaxexp{-(αxσx+αyσy)} (12)
式中:αx为x方向上应力系数,无量纲;αy为y方向上应力系数,无量纲;σx为x方向上正应力大小,MPa;σy为y方向上正应力大小,MPa。
同理:
by=br+bmax=br+bmaxexp{-(βxσx+βyσy)} (13)
式中:βx为正常闭合下正常剪切下的函数斜率,无量纲;βy为正常闭合下正常剪切下的函数斜率,无量纲;σx为x方向上正应力大小,MPa;σy为y方向上正应力大小,MPa。
接下来计算剪切应力条件下裂缝渗透率,计算公式如下:
式中:fdx为在x方向上剪切裂缝的等效频率,1/m;fdy为在y方向上剪切裂缝的等效频率,1/m;dx为在x方向上剪切裂缝的等效开度,μm;dy为在y方向上剪切裂缝的等效开度,μm。
由于只有在临界压力附近的裂缝才会剪切破坏,因此只有部分裂缝会出现剪切变形,其它大部分裂缝都不会出现剪切变形。所以发生剪切破坏的裂缝产生的渗透率会小于之前计算所得的正常闭合裂缝产生的渗透率。由裂缝剪切变形产生的开度大小可以通过下式计算:
对于k<kc:
dx=0 (16)
dy=0 (17)
对于k≥kc:
dx=dmax[1-exp{-γx(k-kc)}] (18)
dy=dmax[1-exp{-γy(k-kc)}] (19)
式中:kc为裂缝变形的临界压力比,无量纲;dmax为裂缝发生剪切破坏后的最大开度,μm。
导致裂缝破坏的临界压力比和临界方向值可以通过库伦破坏准则来计算:
式中:φ为岩石内摩擦角,°;φf为岩石临界破坏角,°。
最后,x和y方向上的渗透率(kx,ky)大小即可通过叠加正常闭合和剪切变形的裂缝渗透率得到,如下:
式中:kx为裂缝在x方向上的渗透率,10-3μm2;ky为裂缝在y方向上的渗透率,10-3μm2;fx为在x方向上正应力产生的裂缝的等效频率,1/m;fy为在y方向上正应力产生的裂缝的等效频率,1/m;br为裂缝的初始平均开度,μm;bmax为受压下裂缝最大开度,μm;αx为x方向上的压力系数,无量纲;αy为y方向上的压力系数,无量纲;σx为x方向上正应力大小,MPa;σy为y方向上正应力大小,MPa;βx为在xy下的函数斜率,无量纲;βy为在yx下的函数斜率,无量纲;fdx为在x方向上剪应力产生的裂缝的等效频率,1/m;fdy为在y方向上剪应力产生的裂缝等效频率,1/m;dmax为裂缝发生剪切破坏后的最大开度,μm;kc为裂缝变形的临界压力比,无量纲。
实施例:
如图2~4所示,选取三块岩心,首先对岩心进行CT扫描,对扫描图像进行灰度处理,识别出裂缝在岩石表面的分布位置及形态特征,并根据本发明模型获得的裂缝长度、裂缝密度和分形维数等参数建立二维离散裂缝网络。根据扫描结果获取的分形维数等参数见表1。
表1计算基础参数Ⅰ
为了进一步说明本发明方法的适用性,进一步使用表2的数据进行了敏感性分析。
表2计算基础参数Ⅱ
参数 | 单位 | 数值 |
分形维数 | - | 1.41 |
最大裂缝长度rmax | m | 150 |
最小裂缝长度rmin | m | 1 |
裂缝区域宽度 | m | 100 |
裂缝区域长度 | m | 100 |
裂缝面密度 | 条/m2 | 5 |
bx上的压力系数(αx,αy) | - | 0.035,0.132 |
by上的压力系数(βx,βy) | - | 0.113,0.069 |
dx,dy的剪切应力系数γx,γy | - | 1.88,1.95 |
剩余开度(br)* | μm | 50 |
最大机械变形开度(bmax)* | μm | 25 |
最大剪切开度(dmax)* | μm | 50 |
裂缝初始开度(bi)* | μm | 20 |
摩擦角 | ° | 24.9 |
临界压力比 | - | 2.66 |
水平主应力σx | MPa | 1~40 |
水平主应力σy | MPa | 1~40 |
注:*实验测试数据
如图5所示,根据式(22)求得x方向上渗透率大小,得到模型在x方向上渗透率与岩心实验测试渗透率的对比。
如图6所示,根据式(23)求得y方向上渗透率大小,得到模型在y方向上渗透率与岩心实验测试渗透率的对比。
如图7所示,根据式(24)计算得到的渗透率与岩心实验测试渗透率的对比。
如图8所示,根据式(22)求得x方向上渗透率大小,得到模型在x方向上渗透率随x与y方向应力比值增加的变化情况。
如图9所示,根据式(23)求得y方向上渗透率大小,得到模型在y方向上渗透率随x与y方向应力比值增加的变化情况。
如图10所示,根据式(24)求得复杂裂缝网络渗透率大小,得到该区域的平均渗透率随压力的变化情况。
如图11所示,复杂裂缝网络渗透率计算模型,得到不同分形维数下,复杂裂缝网络渗透率的变化情况;
如图12所示,复杂裂缝网络渗透率计算模型,得到不同面密度下,复杂裂缝网络渗透率的变化情况;
如图13所示,复杂裂缝网络渗透率计算模型,得到不同裂缝开度下,复杂裂缝网络渗透率的变化情况;
首先通过对比x与y方向渗透率与实验测试所得岩心渗透率的值进行对比可以看出,岩心1计算得到的渗透率为30.5×10-3μm2,实验测得的渗透率为31.45×10-3μm2。对于岩心2和岩心3,计算得出的渗透率大小分别为55.7×10-3μm2和11.2×0-3μm2,实验测得的渗透率分别为50.06×10-3μm2和13.13×10-3μm2。通过对比可以发现,本发明方法计算所得渗透率与实验室测试得到的渗透率误差较小。
从图10可以看出,在水平主应力不等比例增加的条件下,保持x方向上的水平主应力不变,然后增大y方向上水平主应力,应力的比值k从0.5增加到5。岩石的渗透率在k值增大的情况下,出现了先减小后增大的趋势,裂缝渗透率先从2.3×10-3μm2减小到了0.7×10-3μm2。当k值大于2.5以后,渗透率又增加到了2.1×10-3μm2;从图9可以看出,在y轴上渗透率的变化没有x轴方向明显,原因在于x轴上的应力在不断增大,y方向上的渗透率增加受限。
通过对比不同分形维数可以发现,离散裂缝网络的分形维数从1.15增加到1.45的过程中,离散裂缝网络的渗透率在不断的增加(图11);通过对比面密度可以发现随着区域内裂缝密度的增大,渗透率的增加非常明显,说明在区域内裂缝密度的增加对渗透率的贡献非常大,能够显著的改变区域内渗透率的大小。当分形维数从1.15增加到1.45时,在裂缝面密度等于5条/m2时,缝网渗透率从0.15×10-3μm2增加到0.92×10-3μm2;在裂缝面密度等于10条/m2时,缝网渗透率从0.97×10-3μm2增加到1.73×10-3μm2;在裂缝面密度等于15条/m2时,缝网渗透率从1.75×10-3μm2增加到2.71×10-3μm2(图12);
通过对比裂缝开度可以发现,当裂缝开度为30μm时,随着裂缝密度从2条/m2增加到16条/m2,复杂缝网渗透率从0.1×10-3μm2增加到了0.81×10-3μm2;当裂缝开度为120μm时,随着裂缝密度从2条/m2增加到16条/m2,渗透率随着面密度从2×10-3μm2增加到了16×10-3μm2。
本发明提供的一种复杂裂缝网络渗透率测试方法的有益效果为:本发明模型利用对实际岩心进行扫描,通过分形维数生成的裂缝网络,并且考虑了裂缝密度、裂缝开度、围压多种地层实际因素,使得本发明模型更加接近实际地层条件下裂缝的形态分布,通过该模型计算出渗透率也因此更加接近实际地层渗透率。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
1.一种复杂裂缝网络渗透率计算方法,其特征在于,包括以下步骤:
S1、从页岩储层段钻取岩心,并将钻取的岩心制成直径为2.5cm、长度为5.0cm的标准岩心;
S2、对岩心端面进行CT扫描,扫描后图像进行灰度处理,识别出裂缝条数、位置及形态特征;
S3、基于S2中的图像,运用盒子法来获取表征岩石复杂裂缝网络的分形维数;
S4、基于S3中的分形维数,利用蒙特卡洛随机建模理论建立二维离散裂缝网络模型;
S5、基于二维离散裂缝网络模型,通过改进的立方定律求解复杂裂缝网络渗透率。
2.根据权利要求1所述的一种复杂裂缝网络渗透率计算方法,其特征在于,所述步骤S2中对岩心进行CT扫描即利用形貌扫描仪拍下岩心端面的裂缝细节;利用编制的程序,将岩石裂缝面间距图进行灰度处理,并将灰度图处理成二值图,完成数字化处理。
3.根据权利要求1所述的一种复杂裂缝网络渗透率计算方法,其特征在于,所述步骤S3的盒子法是用正方形格子去覆盖裂缝间距图;给定盒子尺码δ,可以计算出覆盖住裂缝间距所需的总盒子数目N;运用盒子法求得的分形维数表达式为:
式中:D为分形维数,无量纲;N为盒子数目,无量纲;δ为盒子尺码,m;i为计数符号,无量纲。
4.根据权利要求1所述的一种复杂裂缝网络渗透率计算方法,其特征在于,所述步骤S4的二维离散裂缝网络模型的建模方法为蒙特卡洛随机建模方法,所述二维离散裂缝网络模型的生成由以下步骤组成:A、模拟区域生成;B、裂缝数目确定;C、裂缝位置表征;D、裂缝几何参数;E、裂缝网络模型生成。
5.根据权利要求1所述的一种复杂裂缝网络渗透率计算方法,其特征在于,所述步骤S5中复杂裂缝网络渗透率包括正常张开型裂缝渗透率、剪切型裂缝渗透率,计算渗透率过程中,由张力形成的张开型和由剪应力形成的剪切型裂缝渗透率分别进行计算,最后将正应力和剪应力诱导产生裂缝所具有渗透率进行叠加,所述正应力和剪应力诱导产生裂缝所具有渗透率的计算公式如下:
式中:kx为复杂缝网在x方向上的渗透率,10-3μm2;ky为复杂缝网在y方向上的渗透率,10-3μm2;fx为在x方向上正应力产生的裂缝等效频率,1/m;fy为在y方向上正应力产生的裂缝的等效频率,1/m;br为裂缝的初始平均开度,μm;bmax为受压下裂缝最大开度,μm;αx为x方向上的压力系数,无量纲;αy为y方向上的压力系数,无量纲;σx为x方向上正应力大小,MPa;σy为y方向上正应力大小,MPa;βx为在xy下的函数斜率,无量纲;βy为在yx下的函数斜率,无量纲;fdx为在x方向上剪应力产生的裂缝等效频率,1/m;fdy为在y方向上剪应力产生的裂缝等效频率,1/m;dmax为裂缝发生剪切破坏后的最大开度,μm;k—x与y方向应力比,无量纲;γx——x方向剪切应力系数,无量纲;γy——y方向剪切应力系数,无量纲;kc——裂缝变形的临界压力比,无量纲;kc为裂缝变形的临界压力比,无量纲。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810312664.XA CN108590640B (zh) | 2018-04-09 | 2018-04-09 | 一种复杂裂缝网络渗透率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810312664.XA CN108590640B (zh) | 2018-04-09 | 2018-04-09 | 一种复杂裂缝网络渗透率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108590640A true CN108590640A (zh) | 2018-09-28 |
CN108590640B CN108590640B (zh) | 2019-10-22 |
Family
ID=63621351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810312664.XA Active CN108590640B (zh) | 2018-04-09 | 2018-04-09 | 一种复杂裂缝网络渗透率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108590640B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933844A (zh) * | 2019-01-28 | 2019-06-25 | 西南石油大学 | 一种基于分形维数的岩石破裂复杂程度表征方法 |
CN110018290A (zh) * | 2019-04-16 | 2019-07-16 | 中国矿业大学(北京) | 一种基于二维模拟试验台覆岩裂隙演化分形研究方法 |
CN110567858A (zh) * | 2019-10-17 | 2019-12-13 | 西南石油大学 | 一种基于分形理论的预测页岩纳米孔隙渗透率的方法 |
CN110738001A (zh) * | 2019-10-17 | 2020-01-31 | 重庆科技学院 | 一种非常规储层压裂增产改造区计算方法 |
CN111272630A (zh) * | 2020-02-28 | 2020-06-12 | 西南石油大学 | 致密岩心人工裂缝参数的计算方法 |
CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
CN112329975A (zh) * | 2020-09-04 | 2021-02-05 | 中国石油天然气股份有限公司 | 油气资源预测方法、装置、计算机设备及存储介质 |
CN112730198A (zh) * | 2021-01-15 | 2021-04-30 | 中国地质大学(武汉) | 一种页岩天然裂缝渗透率计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102682184A (zh) * | 2011-03-08 | 2012-09-19 | 中国科学院研究生院 | 随机分布三维裂隙网络中裂隙对相交判断方法 |
RU2476670C1 (ru) * | 2011-09-15 | 2013-02-27 | Общество с ограниченной ответственностью "Газпромнефть Научно-Технический Центр" (ООО "Газпромнефть НТЦ") | Способ определения фильтрационных свойств совместно работающих пластов (варианты) |
CN105160050A (zh) * | 2015-07-01 | 2015-12-16 | 中国矿业大学(北京) | 一种基于ct扫描的受载煤岩损伤本构方程构建方法 |
CN105279790A (zh) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | 裂缝网络三维数字岩心建模方法 |
CN105512366A (zh) * | 2015-11-26 | 2016-04-20 | 东北石油大学 | 含天然裂缝的致密储层体积压裂树状随机缝网描述方法 |
WO2018026485A1 (en) * | 2016-08-02 | 2018-02-08 | Schlumberger Technology Corporation | Wellbore sealant using nanoparticles |
-
2018
- 2018-04-09 CN CN201810312664.XA patent/CN108590640B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102682184A (zh) * | 2011-03-08 | 2012-09-19 | 中国科学院研究生院 | 随机分布三维裂隙网络中裂隙对相交判断方法 |
RU2476670C1 (ru) * | 2011-09-15 | 2013-02-27 | Общество с ограниченной ответственностью "Газпромнефть Научно-Технический Центр" (ООО "Газпромнефть НТЦ") | Способ определения фильтрационных свойств совместно работающих пластов (варианты) |
CN105279790A (zh) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | 裂缝网络三维数字岩心建模方法 |
CN105160050A (zh) * | 2015-07-01 | 2015-12-16 | 中国矿业大学(北京) | 一种基于ct扫描的受载煤岩损伤本构方程构建方法 |
CN105512366A (zh) * | 2015-11-26 | 2016-04-20 | 东北石油大学 | 含天然裂缝的致密储层体积压裂树状随机缝网描述方法 |
WO2018026485A1 (en) * | 2016-08-02 | 2018-02-08 | Schlumberger Technology Corporation | Wellbore sealant using nanoparticles |
Non-Patent Citations (2)
Title |
---|
丁志文 等: "岩体分形离散裂隙网络系统中流体流动模拟研究进展", 《水利水电科技进展》 * |
李大奇 等: "裂缝网络地层钻井液漏失模拟", 《钻井液与完井液》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109933844A (zh) * | 2019-01-28 | 2019-06-25 | 西南石油大学 | 一种基于分形维数的岩石破裂复杂程度表征方法 |
CN110018290A (zh) * | 2019-04-16 | 2019-07-16 | 中国矿业大学(北京) | 一种基于二维模拟试验台覆岩裂隙演化分形研究方法 |
CN110738001B (zh) * | 2019-10-17 | 2022-04-22 | 重庆科技学院 | 一种非常规储层压裂增产改造区计算方法 |
CN110567858A (zh) * | 2019-10-17 | 2019-12-13 | 西南石油大学 | 一种基于分形理论的预测页岩纳米孔隙渗透率的方法 |
CN110738001A (zh) * | 2019-10-17 | 2020-01-31 | 重庆科技学院 | 一种非常规储层压裂增产改造区计算方法 |
CN110567858B (zh) * | 2019-10-17 | 2020-11-03 | 西南石油大学 | 一种基于分形理论的预测页岩纳米孔隙渗透率的方法 |
CN111272630A (zh) * | 2020-02-28 | 2020-06-12 | 西南石油大学 | 致密岩心人工裂缝参数的计算方法 |
CN111963158A (zh) * | 2020-08-12 | 2020-11-20 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
CN111963158B (zh) * | 2020-08-12 | 2022-04-01 | 西南石油大学 | 一种碳酸盐岩酸压后基质渗透率计算方法 |
CN112329975A (zh) * | 2020-09-04 | 2021-02-05 | 中国石油天然气股份有限公司 | 油气资源预测方法、装置、计算机设备及存储介质 |
CN112329975B (zh) * | 2020-09-04 | 2024-04-30 | 中国石油天然气股份有限公司 | 油气资源预测方法、装置、计算机设备及存储介质 |
CN112730198A (zh) * | 2021-01-15 | 2021-04-30 | 中国地质大学(武汉) | 一种页岩天然裂缝渗透率计算方法 |
CN112730198B (zh) * | 2021-01-15 | 2022-04-01 | 中国地质大学(武汉) | 一种页岩天然裂缝渗透率计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108590640B (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108590640B (zh) | 一种复杂裂缝网络渗透率计算方法 | |
CN104992468B (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
CN108763650B (zh) | 一种覆岩采动裂隙网络模型构建方法 | |
Aimene et al. | Geomechanical modeling of hydraulic fractures interacting with natural fractures—Validation with microseismic and tracer data from the Marcellus and Eagle Ford | |
CN110147561B (zh) | 一种含天然裂缝致密油气储层体积压裂缝网预测方法 | |
CN106501856B (zh) | 亚地震断层定量预测方法 | |
CN107220493A (zh) | 基于微地震事件的页岩气水平井网络裂缝建模方法 | |
CN105740563A (zh) | 一种成熟油田二次开发之优势通道识别技术 | |
KR101591430B1 (ko) | 지하구조 파악방법 및 이를 이용한 셰일가스 시추방법 | |
CN106769463A (zh) | 一种岩心压后裂缝复杂程度定量表征方法 | |
Ding et al. | The analytical model of hydraulic fracture initiation for perforated borehole in fractured formation | |
Ouenes et al. | Using geomechanical modeling to quantify the impact of natural fractures on well performance and microseismicity: Application to the Wolfcamp, Permian Basin, Reagan County, Texas | |
CN108397130A (zh) | 钻井方法和装置 | |
CN104865610A (zh) | 一种深部储层岩心的数值表征方法 | |
Li et al. | A novel approach of two-dimensional representation of rock fracture network characterization and connectivity analysis | |
CN104751002A (zh) | 判别有效砂岩储层的方法 | |
CN112861890A (zh) | 一种储层评价模型构建方法及储层识别方法 | |
Jia et al. | Exploration on the characteristics of 3D crack network expansion induced by hydraulic fracturing: a hybrid approach combining experiments and algorithms | |
CN111077053A (zh) | 一种裂缝型多孔介质渗流分析可视化方法 | |
Youjun et al. | Analysis of pore scale fluid migration in a porous medium-application to coal rock seam | |
CN109190218A (zh) | 一种致密储层有效缝网的动态识别方法 | |
Peng et al. | A numerical simulation of seepage structure surface and its feasibility | |
CN117371161A (zh) | 一种模拟压裂缝在砾岩中扩展的数值模拟方法和系统 | |
Chernyshov et al. | Improving the efficiency of terrigenous oil-saturated reservoir development by the system of oriented selective slotted channels | |
Jiang et al. | A Numerical Simulation Approach for Shale Fracture Network Characterization Using Hybrid EDFM Method |
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 |