CN106777649A - 一种裂缝型储集层孔隙结构定量评价方法 - Google Patents
一种裂缝型储集层孔隙结构定量评价方法 Download PDFInfo
- Publication number
- CN106777649A CN106777649A CN201611128885.9A CN201611128885A CN106777649A CN 106777649 A CN106777649 A CN 106777649A CN 201611128885 A CN201611128885 A CN 201611128885A CN 106777649 A CN106777649 A CN 106777649A
- Authority
- CN
- China
- Prior art keywords
- porosity
- reservoir
- throat radius
- capillary pressure
- pore
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种裂缝型储集层孔隙结构定量评价方法,步骤为:利用微电阻率扫描成像测井仪器采集数据,并对采集的数据进行处理,得到反映储集层基质和裂缝孔隙度发育情况的孔隙度频率谱;通过对常规的孔隙度测井资料进行处理,计算储集层孔隙度曲线,并利用计算得到的孔隙度曲线将储集层划分为三种不同类型;针对三种不同类型的储集层,在大孔隙部分和小孔隙部分,分别采用分段幂函数刻度方法,将孔隙度频率谱转化为储集层孔喉半径分布谱;根据毛管压力与孔喉半径之间的关系,将孔喉半径分布谱转化为伪毛管压力曲线;根据构造的伪毛管压力曲线和孔喉半径分布谱,计算储集层孔隙结构评价参数,实现定量评价裂缝型储集层孔隙结构。
Description
技术领域
本发明属于储集层评价技术领域,具体涉及一种裂缝型储集层孔隙结构定量评价方法。
背景技术
在自然界中,将具有一定储集空间并能使储集在其中的流体在一定压差下流动的岩石称为储集岩。由储集岩所构成的地层称为储集层,简称储层。
裂缝型储集层是按储集空间性质划分的一种储集层类型。储集油气的空间和渗滤通道是裂缝,或者是以孔隙为主要储集空间而裂缝为渗滤通道的储集层。其储集空间的形成多以次生作用为主。碳酸盐岩、结晶岩类、泥质岩类多具有此类储集空间,并往往形成不规则的储集体,少部分坚硬致密的碎屑岩也可能属于这一类。
对于裂缝型储集层而言,由于储集层致密,导致基质孔隙度和渗透率低、孔隙结构极差,裂缝对储集层孔隙结构的改造作用至关重要。裂缝不仅可以作为油气储集的空间,而且还能疏导基质孔隙,提高储集层的孔隙连通性。一般含有裂缝的地层,其孔隙连通性将大大提高,储集层的生产潜力也增大。因此,如果能够准确的评价裂缝型储集层的孔隙结构,在此基础上,划分储集层类型,对于在普遍发育较差的储集层中寻找孔隙连通性较好的优质储集层、提高储集层的勘探效率、降低开发风险将具有重要意义。
对储集层孔隙结构进行定量评价最有效的资料是毛管压力曲线。通过对毛管压力曲线进行处理,可以获取反映储集层孔隙大小及连通性的孔喉半径分布谱以及计算储集层孔隙结构评价参数,以利用其划分储集层类型,确定优质储集层。
然而,通过钻取岩心开展毛管压力实验获取的毛管压力曲线资料非常有限,无法实现连续定量评价储集层孔隙结构的目的。
对于裂缝型储集层评价而言,最有效的测井方法是微电阻率扫描成像测井。一般情况下,利用微电阻率扫描成像测井仪器测量并处理得到的是采用变密度形式反映储集层电导率高低的成像图,通过观察成像图的变密度图像颜色的深浅,可以反映储集层的裂缝发育状况。采用特殊的方法对该图像进行处理后,可以提取反映裂缝发育状况的裂缝孔隙度、裂缝密度、裂缝宽度等参数以及反映储集层基质和裂缝孔隙度分布情况的孔隙度频率谱。通过分析这些参数,可以了解储集层裂缝的发育状况和判断储集层的有效性。一般情况下,裂缝孔隙度越高、裂缝密度越大、裂缝宽度越宽的储集层,其孔隙度频率谱往往表现为双峰分布,且谱的分布较宽,该储集层往往是优质储集层。裂缝孔隙度越低、裂缝密度越小、裂缝宽度越窄的储集层,其孔隙度频率谱往往表现为单峰分布,且谱的分布较窄,反映储集层的有效性较差。
上述分析表明,通过分析储集层裂缝孔隙度等参数的大小及其孔隙度频率谱的形态,可以判断储集层的基质和裂缝发育状况,以识别优质储集层和研究储集层的有效性。然而,这只能算是一种定性或半定量的评价储集层有效性的方法。
对于如何利用孔隙度频率谱来定量评价裂缝型储集层孔隙结构,进而评价储集层的有效性,目前尚未有相关文献报道。
发明内容
为了克服现有技术只能定性或半定量地评价储集层有效性的缺陷,本发明提供一种裂缝型储集层孔隙结构定量评价方法,本发明利用微电阻率扫描成像测井孔隙度频率谱获取储集层孔喉半径分布谱、连续构造伪毛管压力曲线并计算储集层孔隙结构参数,实现根据实际测量的微电阻率扫描成像测井资料、连续定量评价裂缝型储集层孔隙结构的目的。
本发明解决其技术问题所采用的技术方案是:
本发明一种裂缝型储集层孔隙结构定量评价方法,所述定量评价方法包括以下步骤:
1)利用微电阻率扫描成像测井仪器采集数据,并对采集的数据进行处理,得到反映储集层基质和裂缝孔隙度发育情况的孔隙度频率谱;
2)通过对常规的孔隙度测井资料进行处理,计算储集层孔隙度曲线,并利用计算得到的孔隙度曲线将储集层划分为三种不同类型;
3)针对三种不同类型的储集层,在大孔隙部分和小孔隙部分,分别采用分段幂函数刻度方法,将孔隙度频率谱转化为储集层孔喉半径分布谱;
4)根据毛管压力与孔喉半径之间的关系,将孔喉半径分布谱转化为伪毛管压力曲线;
5)根据构造的伪毛管压力曲线和孔喉半径分布谱,计算储集层孔隙结构评价参数,实现定量评价裂缝型储集层孔隙结构。
步骤1)中按照如下公式对微电阻率扫描成像测井资料进行处理,得到储集层孔隙度频率谱:
式中:Φ为孔隙度,小数;m为胶结指数;n为饱和度指数;a和b为与岩性有关的常数;Rmf为泥浆电阻率,单位为Ω.m;Rxo为冲洗带电阻率,在微电阻率扫描成像测井资料中即为测量的电阻率结果,单位为Ω.m;Sxo为冲洗带含水饱和度,微电阻率扫描成像测井的探测深度较浅,主要反映泥浆侵入带的信息,Sxo等于1。
步骤2)中计算储集层孔隙度曲线的方法为利用岩石体积物理模型或地区的经验性公式,三种不同类型的储集层按照如下标准划分:I类储集层:孔隙度≥6.0%;II类储集层:2.0%≤孔隙度<6.0%;III类储集层:孔隙度<2.0%。
步骤3)中对于三种不同类型的储集层,采用如下方法将孔隙度频率谱转化为储集层孔喉半径分布谱:
a、对于I类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m11为待定的系数,由岩心分析得到;n11为待定的系数,由岩心分析得到;m12为待定的系数,由岩心分析得到;n12为待定的系数,由岩心分析得到;
b、对于II类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m21为待定的系数,由岩心分析得到;n21为待定的系数,由岩心分析得到;m22为待定的系数,由岩心分析得到;n22为待定的系数,由岩心分析得到;
c、对于III类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m31为待定的系数,由岩心分析得到;n31为待定的系数,由岩心分析得到;m32为待定的系数,由岩心分析得到;n32为待定的系数,由岩心分析得到。
步骤4)中按照如下式所述的关系将孔喉半径分布谱转化为毛管压力曲线:
式中:Pc(i)为第i个毛管压力,MPa;Rc(i)为第i个毛管压力所对应的孔喉半径,μm。
步骤4)中获取伪毛管压力曲线的步骤为:
a、将孔隙度频率谱的幅度按照储集层孔隙度频率从大到小的顺序进行反向累加,并进行归一化,得到一条在物理意义上与毛管压力曲线上非润湿相饱和度相似的孔隙度频率谱幅度反向累加曲线;
b、以归一化后的孔隙度频率谱幅度反向累加曲线为线性横坐标,以转化的毛管压力为对数纵坐标作图,得到储集层伪毛管压力曲线。
步骤5)中利用构造的伪毛管压力曲线和孔喉半径分布谱计算储集层孔隙结构评价参数方法,按照杨胜来等著作的普通高等学校“十五”规划教材《油层物理学》中209-233页所述的方法进行。
本发明的优点是:在根据孔隙度曲线对储集层进行分类的基础上,利用微电阻率扫描成像测井孔隙度频率谱转换得到储集层孔喉半径分布谱,依据储集层孔喉半径分布谱构造出储集层伪毛管压力曲线,并计算储集层孔隙结构评价参数,可以实现利用微电阻率扫描成像测井资料连续定量评价裂缝型储集层孔隙结构的目的。
附图说明
图1是本发明提供的一种裂缝型储集层孔隙结构定量评价方法流程图;
图2是本发明实施例提供的裂缝发育储集层的微电阻率扫描成像测井孔隙度频率谱示意图;
图3是本发明实施例提供的裂缝不发育储集层的微电阻率扫描成像测井孔隙度频率谱示意图;
图4是本发明实施例提供的裂缝发育储集层的岩心孔喉半径分布谱示意图;
图5是本发明实施例提供的裂缝不发育储集层的岩心孔喉半径分布谱示意图;
图6是本发明实施例提供的第I类储集层典型岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心实验的毛管压力曲线对比示意图;
图7是本发明实施例提供的第II类储集层典型岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心实验的毛管压力曲线对比示意图;
图8是本发明实施例提供的第III类储集层典型岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心实验的毛管压力曲线对比示意图;
图9所示为根据本发明所述方法,利用孔隙度频率谱构造的伪毛管压力曲线及计算的储集层孔隙结构评价参数与岩心压汞实验获取的毛管压力曲线和孔隙结构评价参数的对比图。
具体实施方式
下面将结合附图对本发明作进一步的说明:
通过对微电阻率扫描成像测井资料进行处理,获取储集层孔隙度频率谱。利用Alphalog软件上开发的微电阻率扫描成像测井处理模块,对微电阻率扫描成像测井资料进行处理,可以获取反映储集层基质和裂缝孔隙度发育情况的孔隙度频率谱。
通过对比图2和图3可以看到,在裂缝发育井段,其孔隙度频率谱的分布较宽,且表现为双峰分布,其中,左边的峰为基质孔隙度响应,而右边的峰则是孔喉半径较宽的裂缝孔隙度响应;在裂缝不发育井段,其孔隙度频率谱的分布较窄,且表现为单峰分布,主要反映的是储集层基质孔隙度特征。
图4和图5则分别为与图2和图3相同层段上,根据压汞毛管压力曲线获取的孔喉半径分布谱。从图2与图4以及图3与图5的对比可以看到,孔隙度频率谱的形态与相应深度上孔喉半径分布谱的形态非常的类似,裂缝较发育的层段,相应的孔喉半径分布谱也较宽,且也表现为双峰分布。而裂缝不发育的层段,其孔喉半径分布谱也较窄,为单峰分布。
因此,从上述对实验结果的分析可以看到,储集层孔喉半径分布谱与孔隙度频率谱在形态上具有很好地相似性,但二者所反映的物理意义存在差异,利用孔隙度频率谱只能定性的反映储集层的孔隙分布特征,而不能定量表征储集层的孔隙结构。为了实现利用孔隙度频率谱定量表征裂缝型储集层的孔隙结构,最直接有效的方法就是建立适当的模型,将孔隙度频率谱转化为孔喉半径分布谱,再利用转化得到的孔喉半径分布谱,得到储集层的伪毛管压力曲线,并计算出评价储集层孔隙结构的平均孔喉半径、最大孔喉半径、排驱压力等参数。
本发明在上述实验结果分析的基础上,利用实际测量的从岩心压汞毛管压力实验中得到的孔喉半径分布谱和通过对微电阻率扫描成像测井进行处理得到的孔隙度频率谱,提出一种利用孔隙度频率谱构造伪毛管压力曲线,以定量评价裂缝型储集层孔隙结构的方法。
如图1所示,一种裂缝型储集层孔隙结构定量评价方法,所述定量评价方法包括以下步骤:
1)利用微电阻率扫描成像测井仪器采集数据,并对采集的数据进行处理,得到反映储集层基质和裂缝孔隙度发育情况的孔隙度频率谱;
2)通过对常规的孔隙度测井资料进行处理,计算储集层孔隙度曲线,并利用计算得到的孔隙度曲线将储集层划分为三种不同类型;
3)针对三种不同类型的储集层,在大孔隙部分和小孔隙部分,分别采用分段幂函数刻度方法,将孔隙度频率谱转化为储集层孔喉半径分布谱;
4)根据毛管压力与孔喉半径之间的关系,将孔喉半径分布谱转化为伪毛管压力曲线;
5)根据构造的伪毛管压力曲线和孔喉半径分布谱,计算储集层孔隙结构评价参数,实现定量评价裂缝型储集层孔隙结构。
步骤1)中按照如下公式对微电阻率扫描成像测井资料进行处理,得到储集层孔隙度频率谱:
式中:Φ为孔隙度,小数;m为胶结指数;n为饱和度指数;a和b为与岩性有关的常数;Rmf为泥浆电阻率,单位为Ω.m;Rxo为冲洗带电阻率,在微电阻率扫描成像测井资料中即为测量的电阻率结果,单位为Ω.m;Sxo为冲洗带含水饱和度,微电阻率扫描成像测井的探测深度较浅,主要反映泥浆侵入带的信息,Sxo等于1。
步骤2)中计算储集层孔隙度曲线的方法为利用岩石体积物理模型或地区的经验性公式,三种不同类型的储集层按照如下标准划分:I类储集层:孔隙度≥6.0%;II类储集层:2.0%≤孔隙度<6.0%;III类储集层:孔隙度<2.0%。
步骤3)中对于三种不同类型的储集层,采用如下方法将孔隙度频率谱转化为储集层孔喉半径分布谱:
a、对于I类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m11为待定的系数,由岩心分析得到;n11为待定的系数,由岩心分析得到;m12为待定的系数,由岩心分析得到;n12为待定的系数,由岩心分析得到;
b、对于II类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m21为待定的系数,由岩心分析得到;n21为待定的系数,由岩心分析得到;m22为待定的系数,由岩心分析得到;n22为待定的系数,由岩心分析得到;
c、对于III类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m31为待定的系数,由岩心分析得到;n31为待定的系数,由岩心分析得到;m32为待定的系数,由岩心分析得到;n32为待定的系数,由岩心分析得到。
步骤4)中按照如下式所述的关系将孔喉半径分布谱转化为毛管压力曲线:
式中:Pc(i)为第i个毛管压力,MPa;Rc(i)为第i个毛管压力所对应的孔喉半径,μm。
步骤4)中获取伪毛管压力曲线的步骤为:
a、将孔隙度频率谱的幅度按照储集层孔隙度频率从大到小的顺序进行反向累加,并进行归一化,得到一条在物理意义上与毛管压力曲线上非润湿相饱和度相似的孔隙度频率谱幅度反向累加曲线;
b、以归一化后的孔隙度频率谱幅度反向累加曲线为线性横坐标,以转化的毛管压力为对数纵坐标作图,得到储集层伪毛管压力曲线。
步骤5)中利用构造的伪毛管压力曲线和孔喉半径分布谱计算储集层孔隙结构评价参数方法,按照杨胜来等著作的普通高等学校“十五”规划教材《油层物理学》中209-233页所述的方法进行。
根据本发明所述的一种裂缝型储集层孔隙结构定量评价方法,对国内东部地区某油田29块岩心毛管压力资料和实际测量的微电阻率扫描成像测井资料进行了处理,分别标定了三类储集层利用孔隙度频率谱转换孔喉半径分布谱的模型参数,获取了利用孔隙度频率谱连续构造孔喉半径分布谱的模型。通过孔喉半径与毛管压力之间的关系,得到了连续分布的伪毛管压力曲线并计算出平均孔喉半径、最大孔喉半径、排驱压力等孔隙结构评价参数。
为了定量表征本发明所述方法的可靠性,选取了3块代表性岩心样品,将利用本发明所述方法,从微电阻率扫描成像测井孔隙度频率谱中获取的伪毛管压力曲线与相应深度上的岩心压汞毛管压力曲线进行了对比。
图6至图8列举了3块代表性岩心压汞毛管压力曲线与构造的伪毛管压力曲线对比图。其中,图6所示为第I类储集层代表性岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心压汞毛管压力曲线的对比结果;图7所示为第II类储集层代表性岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心压汞毛管压力曲线的对比结果;图8所示为第III类储集层代表性岩心利用孔隙度频率谱构造的伪毛管压力曲线与岩心压汞毛管压力曲线的对比结果。从图中所示的对比结果可以明显看到,对于3类不同孔隙结构储集层的岩石而言,利用本发明所述方法,从微电阻率扫描成像测井中获取的伪毛管压力曲线与岩心实验得到的压汞毛管压力曲线之间吻合较好,充分证明了本发明所述的利用孔隙度频率谱构造伪毛管压力曲线,以评价裂缝型储集层孔隙结构方法的可靠性。
图9所示的效果图共分为十道,图中第一道包括自然伽马(GR)、自然电位(SP)和井径曲线(CAL),主要用于识别有效储集层;第二道包括密度测井(DEN)曲线、中子测井(CNL)曲线和声波时差测井(DT)曲线,主要用于计算储集层的孔隙度;第三道为深侧向电阻率曲线(RT)和浅侧向电阻率曲线(RXO);第四道为深度道,单位m;第五道为根据常规孔隙度测井曲线,利用体积物理模型计算的储集层孔隙度(PHIT);第六道列举了实际测量的微电阻率扫描成像测井变密度图像;第七道列举了通过对微电阻率扫描成像测井变密度图进行处理获取的储集层孔隙度频率谱图像;第八道为根据本发明所述方法利用孔隙度频率谱连续构造的伪毛管压力曲线与岩心压汞毛管压力曲线的对比,图9中黑色细线为利用孔隙度频率谱连续构造的伪毛管压力曲线,为了曲线显示的方便,本发明中将构造的伪毛管压力曲线累加到最大进汞压力29.70MPa的状态,而离散的黑色粗线为岩心压汞毛管压力曲线;从图9中所示的对比结果可以看出,利用本发明所述方法构造的伪毛管压力曲线与岩心压汞实验结果具有较好的一致性。这说明,利用本发明所述方法可以将孔隙度频率谱连续地转化成伪毛管压力曲线。图9中第九道至第十三道列举了根据构造的伪毛管压力曲线计算的储集层孔隙结构评价参数与岩心实验结果的对比。第九道中RMAX为利用本发明实施例提供的方法连续构造的伪毛管压力曲线计算的最大孔喉半径,CRMAX为岩心压汞实验得到的最大孔喉半径;第十道中R50为利用本发明实施例提供的方法连续构造的伪毛管压力曲线计算的中值半径,CR50为岩心压汞实验得到的中值半径;第十一道中RM为利用本发明实施例提供的方法连续构造的伪毛管压力曲线计算的平均孔喉半径,CRM为岩心压汞实验得到的平均孔喉半径;第十二道中PC50为利用本发明实施例提供的方法连续构造的伪毛管压力曲线计算的中值压力,CP50为岩心压汞实验得到的中值压力;第十三道中PD为利用本发明实施例提供的方法连续构造的伪毛管压力曲线计算的排驱压力,CPD为岩心压汞实验得到的排驱压力。
从图9上可以看出,利用本发明所述方法连续构造的伪毛管压力曲线计算的储集层孔隙结构评价参数与岩心压汞实验结果非常吻合,这说明,利用本发明所述方法可以将孔隙度频率谱连续地转化成伪毛管压力曲线,以计算出准确的储集层孔隙结构评价参数。
最后应说明的是:显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明的保护范围之中。
Claims (6)
1.一种裂缝型储集层孔隙结构定量评价方法,其特征在于,所述定量评价方法包括以下步骤:
1)利用微电阻率扫描成像测井仪器采集数据,并对采集的数据进行处理,得到反映储集层基质和裂缝孔隙度发育情况的孔隙度频率谱;
2)通过对常规的孔隙度测井资料进行处理,计算储集层孔隙度曲线,并利用计算得到的孔隙度曲线将储集层划分为三种不同类型;
3)针对三种不同类型的储集层,在大孔隙部分和小孔隙部分,分别采用分段幂函数刻度方法,将孔隙度频率谱转化为储集层孔喉半径分布谱;
4)根据毛管压力与孔喉半径之间的关系,将孔喉半径分布谱转化为伪毛管压力曲线;
5)根据构造的伪毛管压力曲线和孔喉半径分布谱,计算储集层孔隙结构评价参数,实现定量评价裂缝型储集层孔隙结构。
2.如权利要求1所述的一种裂缝型储集层孔隙结构定量评价方法,其特征在于,步骤1)中按照如下公式对微电阻率扫描成像测井资料进行处理,得到储集层孔隙度频率谱:
式中:Φ为孔隙度,小数;m为胶结指数;n为饱和度指数;a和b为与岩性有关的常数;Rmf为泥浆电阻率,单位为Ω.m;Rxo为冲洗带电阻率,在微电阻率扫描成像测井资料中即为测量的电阻率结果,单位为Ω.m;Sxo为冲洗带含水饱和度,微电阻率扫描成像测井的探测深度较浅,主要反映泥浆侵入带的信息,Sxo等于1。
3.如权利要求2所述的一种裂缝型储集层孔隙结构定量评价方法,其特征在于,步骤2)中计算储集层孔隙度曲线的方法为利用岩石体积物理模型或地区的经验性公式,三种不同类型的储集层按照如下标准划分:I类储集层:孔隙度≥6.0%;II类储集层:2.0%≤孔隙度<6.0%;III类储集层:孔隙度<2.0%。
4.如权利要求3所述的一种裂缝型储集层孔隙结构定量评价方法,其特征在于,步骤3)中对于三种不同类型的储集层,采用如下方法将孔隙度频率谱转化为储集层孔喉半径分布谱:
a、对于I类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m11为待定的系数,由岩心分析得到;n11为待定的系数,由岩心分析得到;m12为待定的系数,由岩心分析得到;n12为待定的系数,由岩心分析得到;
b、对于II类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m21为待定的系数,由岩心分析得到;n21为待定的系数,由岩心分析得到;m22为待定的系数,由岩心分析得到;n22为待定的系数,由岩心分析得到;
c、对于III类储集层,在大孔隙部分和小孔隙部分,分别采用如下式所述的幂函数将孔隙度频率谱转化为储集层孔喉半径分布谱:
大孔隙部分:
小孔隙部分:
式中:Rc为储集层孔喉半径分布谱;POR_dist为储集层孔隙度频率分布谱;m31为待定的系数,由岩心分析得到;n31为待定的系数,由岩心分析得到;m32为待定的系数,由岩心分析得到;n32为待定的系数,由岩心分析得到。
5.如权利要求4所述的一种裂缝型储集层孔隙结构定量评价方法,其特征在于,步骤4)中按照如下式所述的关系将孔喉半径分布谱转化为毛管压力曲线:
式中:Pc(i)为第i个毛管压力,MPa;Rc(i)为第i个毛管压力所对应的孔喉半径,μm。
6.如权利要求5所述的一种裂缝型储集层孔隙结构定量评价方法,其特征在于,步骤4)中获取伪毛管压力曲线的步骤为:
a、将孔隙度频率谱的幅度按照储集层孔隙度频率从大到小的顺序进行反向累加,并进行归一化,得到一条在物理意义上与毛管压力曲线上非润湿相饱和度相似的孔隙度频率谱幅度反向累加曲线;
b、以归一化后的孔隙度频率谱幅度反向累加曲线为线性横坐标,以转化的毛管压力为对数纵坐标作图,得到储集层伪毛管压力曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611128885.9A CN106777649B (zh) | 2016-12-09 | 2016-12-09 | 一种裂缝型储集层孔隙结构定量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611128885.9A CN106777649B (zh) | 2016-12-09 | 2016-12-09 | 一种裂缝型储集层孔隙结构定量评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777649A true CN106777649A (zh) | 2017-05-31 |
CN106777649B CN106777649B (zh) | 2020-06-12 |
Family
ID=58874964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611128885.9A Active CN106777649B (zh) | 2016-12-09 | 2016-12-09 | 一种裂缝型储集层孔隙结构定量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777649B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525753A (zh) * | 2017-07-25 | 2017-12-29 | 中国石油大学(北京) | 确定孔隙曲折度的方法及装置 |
CN109386285A (zh) * | 2018-10-12 | 2019-02-26 | 西安石油大学 | 一种特低渗透油层泥浆侵入程度及其影响的评价方法 |
CN112098291A (zh) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | 基于电成像孔隙度谱的渗透率计算方法及系统 |
CN112377175A (zh) * | 2020-11-03 | 2021-02-19 | 长江大学 | 一种优化钻井泥浆快速识别低阻油气层的方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN103279647A (zh) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法 |
CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
-
2016
- 2016-12-09 CN CN201611128885.9A patent/CN106777649B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN103279647A (zh) * | 2013-05-03 | 2013-09-04 | 中国石油天然气集团公司 | 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法 |
CN103592690A (zh) * | 2013-10-24 | 2014-02-19 | 长江大学 | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 |
CN105386753A (zh) * | 2015-10-28 | 2016-03-09 | 中国地质大学(北京) | 一种利用核磁共振测井构造伪毛管压力曲线方法 |
CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学(北京) | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
Non-Patent Citations (5)
Title |
---|
刘晓鹏, 肖亮, 张伟: "储层毛管压力曲线构造方法及其应用", 《 西南石油大学学报(自然科学版)》 * |
李宁,潘保芝: "岩心核磁几谱转换为岩心孔径分布谱及其应用", 《中国地球物理·2009》 * |
李爱芬,任晓霞,王桂娟,王永政,江凯亮: "核磁共振研究致密砂岩孔隙结构的方法及应用", 《 中国石油大学学报(自然科学版)》 * |
肖亮,毛志强,邹长春: "一种利用核磁共振测井构造毛管压力曲线的新方法", 《2015中国地球科学联合学术年会论文集(二十八)——专题64应用地球物理学前沿、专题65地球生物学》 * |
运华云,赵文杰,刘兵开,周灿灿,周凤鸣: "利用T2 分布进行岩石孔隙结构研究", 《测井技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525753A (zh) * | 2017-07-25 | 2017-12-29 | 中国石油大学(北京) | 确定孔隙曲折度的方法及装置 |
CN109386285A (zh) * | 2018-10-12 | 2019-02-26 | 西安石油大学 | 一种特低渗透油层泥浆侵入程度及其影响的评价方法 |
CN109386285B (zh) * | 2018-10-12 | 2022-04-01 | 西安石油大学 | 一种特低渗透油层泥浆侵入程度及其影响的评价方法 |
CN112098291A (zh) * | 2019-06-18 | 2020-12-18 | 中国石油化工股份有限公司 | 基于电成像孔隙度谱的渗透率计算方法及系统 |
CN112377175A (zh) * | 2020-11-03 | 2021-02-19 | 长江大学 | 一种优化钻井泥浆快速识别低阻油气层的方法及系统 |
CN112377175B (zh) * | 2020-11-03 | 2023-07-25 | 长江大学 | 一种优化钻井泥浆快速识别低阻油气层的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106777649B (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
CN105353419A (zh) | 一种基于储集层分类的核磁毛管压力曲线构造方法 | |
CN106777649A (zh) | 一种裂缝型储集层孔隙结构定量评价方法 | |
CN104374685B (zh) | 一种数字岩心三孔隙组分求取储层渗透率的方法 | |
CN103670384B (zh) | 一种储层分类方法及系统 | |
CN104775810A (zh) | 一种页岩气储层可压性评价方法 | |
CN103026202A (zh) | 获取多孔介质的相容和综合物理性质的方法 | |
CN105781539A (zh) | 一种致密油气储层饱和度测井计算方法 | |
CN103592690A (zh) | 基于电成像测井孔隙度谱信息自动识别储层裂缝的方法 | |
CN101881153B (zh) | 一种常规测井信息融合可视化方法及系统 | |
CN103628871A (zh) | 一种基于阿尔奇公式的电阻率侵入校正的新方法 | |
CN104712330A (zh) | 一种测井渗透率解释方法 | |
CN112363242B (zh) | 基于测录井融合的储层流体识别方法与装置 | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CN107346455A (zh) | 一种用于识别页岩气产能的方法 | |
CN103630939A (zh) | 一种气层识别评价方法 | |
CN107829731A (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN110320569B (zh) | 一种致密砂岩储层单井裂缝发育强度定量评价方法 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN105093333A (zh) | 基于地震-测井标定的缝洞型碳酸盐岩储层评价方法 | |
Kharrat et al. | Rock typeand permeability prediction of a heterogeneous carbonate reservoir using artificial neural networks based on flow zone index approach | |
CN106351652A (zh) | 一种含烃储集层核磁共振测井t2谱形态校正方法 | |
CN110344825A (zh) | 一种低孔低渗低电阻率砂岩油层综合判识方法 | |
CN105590018A (zh) | 砂泥岩薄互层油藏的油水层判别方法 | |
CN104991280A (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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20200518 Address after: 100086 Beijing City, Haidian District Zhongguancun South Street No. 6 Zhucheng building room 507 Applicant after: BOOMING (BEIJING) ENERGY TECHNOLOGY Co.,Ltd. Applicant after: China University of Geosciences (Beijing) Address before: 100086 Beijing City, Haidian District Zhongguancun South Street No. 6 Zhucheng building room 507 Applicant before: BOOMING (BEIJING) ENERGY TECHNOLOGY Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |