CN113075105B - 一种利用综合分形维Dp定量评价砂岩储层的方法 - Google Patents
一种利用综合分形维Dp定量评价砂岩储层的方法 Download PDFInfo
- Publication number
- CN113075105B CN113075105B CN202110320805.4A CN202110320805A CN113075105B CN 113075105 B CN113075105 B CN 113075105B CN 202110320805 A CN202110320805 A CN 202110320805A CN 113075105 B CN113075105 B CN 113075105B
- Authority
- CN
- China
- Prior art keywords
- fractal dimension
- reservoir
- pressure
- sandstone
- comprehensive
- 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 35
- 230000000704 physical effect Effects 0.000 claims abstract description 54
- 239000011148 porous material Substances 0.000 claims abstract description 43
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 230000035699 permeability Effects 0.000 claims description 46
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 42
- 229910052753 mercury Inorganic materials 0.000 claims description 42
- 238000006073 displacement reaction Methods 0.000 claims description 29
- 230000008859 change Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000011548 physical evaluation Methods 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 description 28
- 239000007789 gas Substances 0.000 description 11
- 238000010219 correlation analysis Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提出一种利用综合分形维Dp定量评价砂岩储层的方法，其步骤包括分段求出砂岩中较大孔隙分形维D1和小孔吼分形维D2，然后结合大孔隙分形维D1和小孔吼分形维D2求出综合分形维DP，最后选出与储层物性相关系数最高的综合分形维作为最优储层物性评价指标，再根据综合分形维与储层物性参数之间的关系确定出不同类型储层对应的分形维指标阈值范围，来作为评价不同砂岩储层的标准。本发明综合考虑大孔、中小孔分形维值对储层物性的贡献度，能够定量描述孔隙复杂程度，并通过建立储层物性评价指标与储层类型对应的量化关系，确定出不同类型储层对应的分形维指标阈值范围，可为油气田砂岩储层类型划分提供参考。
Description
技术领域
本发明涉及油田勘探技术领域，具体涉及一种利用综合分形维Dp定量评价砂岩储层的方法。
背景技术
陆相含油气盆地中，砂岩是最常见的油气储层，而储层的物性好坏直接关系到储层储集油气的能力以及内部流体的渗流特征，是油气田勘探开发过程中重点评价的参数之一。砂岩储层的物性特征受控于岩石中孔隙、吼道的组成，实验室内利用压汞曲线评价砂岩孔吼，压汞进行过程中，随着压力增大，液态汞率先进入较大的孔隙，在压力增幅相对较小的条件下，含汞饱和度增大速度较快；当含汞饱和度达到一定程度之后，含汞饱和度增大的速度明显放缓，该阶段对应储层孔隙多对应小孔、微孔。非润湿相的汞驱替孔隙内流体对应的启动压力、含汞饱和度达到50％对应的中值压力、不同压力对应的孔隙半径等指标常常被应用在储层物性评价过程中，国内外学者将以上指标与压汞曲线形态相结合，综合评价储层孔喉特征，并划分成不同类型的储层。
然而，早期对压汞曲线的描述仅仅停留在定性分析方面，后期随着分形维的引入，对储层孔隙不均一性的描述逐渐进入定量化。目前利用分形处理的压汞曲线的方法主要存在如下问题：(1)仅描述中值压力之后的部分的分维值(如文献“徐守余,王淑萍.砂岩储层微观结构分形特征研究——以胜坨油田古近系沙河街组储层为例[J].天然气地球科学,2013,24(5):886-893”)，未考虑大孔隙部分对孔隙分维值的贡献，而该部分孔隙对砂岩孔隙度、渗透率贡献值较大，不考虑该部分孔隙的分维指标缺乏对砂岩孔隙的客观评价；(2)分段求取了孔隙的分维值，同时利用2段指标评价储层，缺乏整体评价储层孔隙的分维指标(如文献“王付勇,程辉.鄂尔多斯盆地延长组致密砂岩孔喉结构与油藏物性表征.吉林大学学报(地球科学版),2020,50(3):721–731；刘凯,石万忠,王任,等.鄂尔多斯盆地杭锦旗地区盒1段致密砂岩孔隙结构分形特征及其与储层物性的关系[J].地质科技通报,2021,40(1):57–68”)；(3)综合考虑了大孔部分对应的分维值D1，和小孔部分对应的分维值D2，但过分夸大了D1在综合指标的作用，求取到的综合分维指标与渗透率、孔隙度指标正相关，明显关系不符合地质规律(如文献“Su Penghui,Xia Zhaohui,Qu Liangchao,et al.Fractalcharacteristics of low-permeability gas sandstones based on anew model formercury intrusion porosimetry.Journal of Natural Gas ence and Engineering,2018,60：246-255”)。本发明提供一种利用综合分形维Dp评价砂岩储层物性的方法，综合考虑大孔、中小孔分维值对储层物性的贡献度，提出一种新的综合分形维Dp计算模型，并建立了该指标与储层类型对应的量化关系，为油气田砂岩储层定量评价提供参考。
发明内容
针对目前分形处理压汞曲线的描述中缺少中值压力与排驱压力之间的大孔部分的描述以及分维描述方法不清晰的技术问题，本发明提出一种利用综合分形维Dp定量评价砂岩储层的方法，综合考虑大孔、中小孔分维值对储层物性的贡献度，能够定量描述孔隙复杂程度，并通过建立储层物性评价指标与储层类型对应的量化关系，确定出不同类型储层对应的分形维指标阈值范围，为油气田砂岩储层类型划分提供参考。
为解决上述技术问题，本发明采用以下技术方案：一种利用综合分形维Dp定量评价砂岩储层的方法，包括以下步骤：
步骤S1：分段求取压汞曲线的分形维，所述的分形维包括砂岩中大孔隙分形维D1和小孔吼分形维D2；
步骤S2：根据所得的大孔隙分形维D1和小孔吼分形维D2，求取压汞曲线的综合分形维DP；
步骤S3：采用线性拟合的方法分别建立大孔隙分形维D1、小孔吼分形维D2和综合分形维DP与储层物性参数之间的关系，然后从大孔隙分形维D1、小孔吼分形维D2和综合分形维DP中选出与储层物性相关系数最高的分形维作为最优储层物性评价指标；
步骤S4：根据储层物性参数划分储层类型，然后根据最优储层物性评价指标与储层物性参数之间的关系，确定出不同类型储层对应的分形维指标阈值范围、并作为评价不同砂岩储层的标准。
所述步骤S1中求取大孔隙分形维D1和小孔吼分形维D2的步骤为：
步骤S1.1：分别对压汞获取的毛管压力Pc和进汞汞饱和度SHg求取相应的对数Log(Pc)、Log(SHg)，并以Log(Pc)作为横坐标、Log(SHg)作为纵坐标，然后将计算后的数据点投影到坐标系中进行拟合得到压汞曲线，并根据坐标系中拟合曲线斜率变化设置拐点；
步骤S1.2：取排驱压力至拐点压力之间的拟合线段的斜率K1，然后按照公式D1＝3+K1计算出砂岩中大孔隙分形维D1；
步骤S1.3：取大于拐点压力部分拟合曲线的斜率K2，然后按照公式D2＝3+K2计算出砂岩中小孔吼分形维D2。
所述步骤S2中求取压汞曲线的综合分形维DP的计算公式为：
DP＝D1*ΔPD1/P+D2*ΔPD2/P
式中：ΔPD1为拐点压力与排驱压力的差值；ΔPD2为大于拐点压力的高压区毛细管压力差；P为整体毛细管压力差；D1*ΔPD1/P表示砂岩中大孔隙对孔隙复杂程度的贡献值；D2*ΔPD2/P表示砂岩中小孔喉对孔隙复杂程度的贡献值。
步骤S3中所述的储层物性参数包括渗透率、孔隙度、饱和度中值压力和排驱压力。
所述步骤S3中选出最优储层物性评价指标的步骤为：
步骤S3.1：以大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值分别作为横坐标，以储层物性相关的渗透率、孔隙度、饱和度中值压力和排驱压力这四项参数分别作为纵坐标，然后采用线性拟合的方法分别拟合各分维值与四项储层物性参数之间的关系，得出各分维值与四项储层物性参数的相关系数；
步骤S3.2：分别对比大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值与四项储层物性参数之间相关系数的大小，最终得出相关系数最大的综合分形维DP作为最优储层物性评价指标。
所述步骤S4中根据渗透率和孔隙度的分布范围划分储层类型，包括Ⅰ类储层、Ⅱ类储层和Ⅲ类储层，其中Ⅰ类储层对应孔隙度≥15％、渗透率≥100md，Ⅱ类储层对应10％≤孔隙度<15％、10md≤渗透率<100md，Ⅲ类储层对应孔隙度<10％、渗透率<10md。
根据综合分形维DP分别与渗透率、孔隙度之间的关系，确定Ⅰ类储层、Ⅱ类储层和Ⅲ类储层所对应的综合分形维DP的分维值，最终得出综合分形维DP的分维值小于3.14的储层为Ⅰ类储层，综合分形维DP的分维值大于3.25的储层为Ⅲ类储层，综合分形维DP的分维值介于3.14～3.25之间的储层为Ⅱ类储层。
本发明首先通过分段分别求取出砂岩中较大孔隙的分形维数D1和小孔吼的分形维数D2，然后结合大孔隙分形维D1和小孔吼分形维D2求出综合分形维DP，采用线性拟合的方法选出与储层物性相关系数最高的综合分形维作为最优储层物性评价指标，然后再根据综合分形维与储层物性参数之间的关系，确定出不同类型储层对应的分形维指标阈值范围来作为评价不同砂岩储层的标准。本发明综合考虑大孔、中小孔分维值对储层物性的贡献度，能够定量描述孔隙复杂程度，并通过建立储层物性评价指标与储层类型对应的量化关系，确定出不同类型储层对应的分形维指标阈值范围，可为油气田砂岩储层类型划分提供参考。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明综合分形维DP评价砂岩物性的流程图；
图2为本发明分段求取大孔隙分形维D1和小孔吼分形维D2的方法示意图；
图3为本发明大孔隙分形维D1和小孔吼分形维D2的地质含义分析图；
图4为本发明大孔隙分形维D1、小孔吼分形维D2和综合分形维DP分别与储层物性参数相关性分析图，其中(a1)～(d1)为大孔隙分形维D1分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图；(a2)～(d2)为小孔吼分形维D2分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图；(a3)～(d3)为综合分形维DP分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图；可以看出，综合分形维DP作为物性评价指标，全面优于大孔隙分形维D1、小孔吼分形维D2。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明提供了一种利用综合分形维Dp定量评价砂岩储层的方法，包括以下步骤，步骤S1：分段求取压汞曲线的分形维，所述的分形维包括砂岩中大孔隙分形维D1和小孔吼分形维D2。步骤S2：根据所得的大孔隙分形维D1和小孔吼分形维D2，求取压汞曲线的综合分形维DP。步骤S3：采用线性拟合的方法分别建立大孔隙分形维D1、小孔吼分形维D2和综合分形维DP与储层物性参数之间的关系，然后从大孔隙分形维D1、小孔吼分形维D2和综合分形维DP中选出与储层物性相关系数最高的分形维作为最优储层物性评价指标。本实施例中所述的储层物性参数包括渗透率、孔隙度、饱和度中值压力和排驱压力。步骤S4：根据储层物性参数划分储层类型，然后根据最优储层物性评价指标与储层物性参数之间的关系，确定出不同类型储层对应的分形维指标阈值范围、并作为评价不同砂岩储层的标准，可为油气田砂岩储层类型划分提供参考。
具体地，所述步骤S1中求取大孔隙分形维D1和小孔吼分形维D2的具体步骤如下：
步骤S1.1：分别对压汞获取的毛管压力Pc和进汞汞饱和度SHg求出相应的对数Log(Pc)、Log(SHg)，本实施例中的数据为白云凹陷15组样品中的6号样品。毛管压力Pc(即进汞压力)和进汞汞饱和度SHg在实验室压汞实验中直接读取。然后以Log(Pc)作为横坐标、Log(SHg)作为纵坐标，将计算后的数据点(Log(Pc)，Log(SHg))投影到坐标系中进行拟合得到压汞曲线。根据坐标系中拟合曲线斜率变化设置拐点，拐点选择在斜率明显变化弯曲程度最大的点，保证拐点两侧的点相关系数之和R1 2+R2 2最大。其中R1 2代表拐点与排驱压力点之间进汞压力对数值Log(Pc)与进汞汞饱和度对数值Log(SHg)之间的相关程度，数值可以从Excel表格中直接获取；R2分别表示代表拐点与最大进汞压力之间进汞压力对数值Log(Pc)与进汞汞饱和度对数值Log(SHg)之间的相关程度，数值可以从Excel表格中直接获取。该样品压汞测试得到的排驱压力为0.097MPa，拐点压力为0.347MPa(两压力点对应的进汞汞饱和度SHg分别为5.840和34.124，可以从压汞实验数据中直接获取)，如图2所示，拟合后两侧的相关系数分别为0.969、0.8929，拟合后相关系数全部大于0.8，证明两侧点在双对数坐标系中符合线性关系，具有分形特征；拐点两侧的点相关指数R1 2+R2 2为1.8619，向两侧移动拐点位置，对应左侧压力0.2098MPa处或者右侧0.483MPa处，重新计算两侧相关系数小于1.86，故确定(Log(0.3470)，Log(34.124))为该样品的拐点位置。
步骤S1.2：取排驱压力至拐点压力之间的拟合线段的斜率K1，即排驱压力的对数值Log(0.097)至拐点压力Log(0.347)之间的拟合曲线斜率K1，该斜率数值为1.4238，小数点后两位为1.42，然后按照公式D1＝3+K1计算出砂岩中大孔隙分形维D1为4.42，该段孔隙对应的孔径范围可以由压汞测试中的毛细管压力Pc换算得到，孔径取值范围为2.12-7.56μm。
步骤S1.3：取大于拐点压力部分拟合曲线的斜率K2，即拐点压力Log(0.347)右侧部分数据点拟合后的斜率为K2，数值上等于0.1887，取小数点后两位为0.19，然后按照公式D2＝3+K2计算出砂岩中小孔吼分形维D2为3.19，该部分孔径范围0.02-2.12μm。由此可以看出，该部分孔径明显小于拐点左侧部分，即小孔吼分形维D2对应中小孔隙，大孔隙分形维D1对应较大孔隙。
按照图2中所采用的方法，分别计算出白云凹陷其余14组样品的分段分维数，计算得到的分形维D1、D2列于表1中，如表1所示：
表1分维数与储层物性对应关系
从表1中也可以看出，分段求取到压汞曲线的分形维D1、D2，并且已经明确分形维D1与砂岩的大孔有关，分形维D2与砂岩的中小孔隙有关，这两个参数代表同一样品中不同孔径的孔隙，未能实现样品孔隙的整体评价。因此，需要进行步骤S2中所述压汞曲线的综合分形维DP的定量求取，然而综合分形维DP的求取过程中必须综合考虑上述两部分孔隙对砂岩物性的影响。
所述的分形维D1、D2是否能与砂岩物性关联，需要结合图3进行分析，图3为本实施例中分形维D1、D2的地质含义分析图。图3中的四条曲线为本次实验的15组样品中的4组，样品分别为11号、12号、13号、14号，选取该部分样品分析的主要原因是排驱压力和拐点压力之间的线段几乎平行，计算得到的分形维D1数值分别为5.70、5.73、6.04、5.55，即压汞曲线的分形维D1对应段大致平行时，求取的分形维D1的数值几乎相等，该数值仅能代表该样品的大孔部分均一程度接近；而11号、12号、13号、14号样品对应的分形维D1段代表的压力段分别为0.207～0.686MPa、0.345～1.032MPa、0.693～1.724MPa、7.182～16.153MPa，代表的孔径段分别为1.07～3.54μm、0.54～2.13μm、0.43～1.06μm、0.046～0.10μm，即分形维D1只能代表大孔的均匀程度，不能反映孔径变化，故无法直接表征储层的物性。
所述的步骤S2中求取压汞曲线的综合分形维DP的计算公式为：
DP＝D1*ΔPD1/P+D2*ΔPD2/P
式中：D1为样品中大孔的分维值，代表大孔的均一程度，与孔径的大小无关；ΔPD1为拐点压力与排驱压力的差值，代表孔径较大的低压区毛细管压力差；ΔPD2为大于拐点压力的高压区毛细管压力差；P为整体毛细管压力差，数值上等于为ΔPD1与ΔPD2之和，同时等于最大毛细管压力减去排驱压力；D1*ΔPD1/P表示砂岩中大孔隙对孔隙复杂程度的贡献值；D2*ΔPD2/P表示砂岩中小孔喉对孔隙复杂程度的贡献值。
在本实施例的步骤S3中分维值与储层孔隙度、渗透率的关系确立采用如下方法：
步骤S3.1：以大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值分别作为横坐标，以储层物性相关的渗透率、孔隙度、饱和度中值压力和排驱压力这四项参数分别作为纵坐标，拟合分维值与物性参数之间的关系，排驱压力和饱和度中值压力由压汞实验中直接提供，砂岩的孔隙度、渗透率分别由气体孔隙度仪器、气体渗透率仪器直接测试获得，得出各分维值与四项储层物性参数的相关系数。
步骤S3.2：分别对比大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值与四项储层物性参数之间相关系数的大小，最终得出相关系数最大的综合分形维DP作为最优储层物性评价指标。具体以图3中的四组样品(11号、12号、13号、14号)为例，计算得到的DP数值分别为3.11、3.18、3.29、3.98，从表1中可知对应的孔隙度依次为13.7％、10.8％、8.5％、3.3％，渗透率分别为3.59md、0.96md、0.32md、0.02md，可以表明综合分形维DP可以作为评价储层物性(孔隙度、渗透率)的有效指标，DP数值增大，孔隙的复杂程度增加，相应的储层孔隙度、渗透率降低。
图4为本实施例中大孔隙分形维D1、小孔吼分形维D2和综合分形维DP分别与储层物性参数相关性分析图，其中(a1)～(d1)为大孔隙分形维D1分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图；(a2)～(d2)为小孔吼分形维D2分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图；(a3)～(d3)为综合分形维DP分别与排驱压力、饱和度中值压力、渗透率和孔隙度相对应的相关性分析图。证实综合分形维Dp作为评价指标的有效性，不能仅靠15组数据中的4组，需要拟合分形维D1、D2、Dp与常用储层物性评价参数排驱压力、饱和度中值压力、渗透率、孔隙度等物性参数，相关系数较高的分形维可以作为量化评价储层好坏的定量指标。从图4可以看出，D1与排驱压力、饱和度中值压力、渗透率、孔隙度的相关系数(保留小数点后2位)分别为0.09、0.05、0.46、0.11，即分形维D1与储层物性基本不具备相关性。D2与排驱压力、饱和度中值压力、渗透率、孔隙度的相关系数(保留小数点后2位)分别为0.75、0.77、0.51、0.73，综合分形维Dp与排驱压力、饱和度中值压力、渗透率、孔隙度的相关系数(保留小数点后2位)分别为0.79、0.85、0.77、0.64、0.74，证实综合分形维DP作为物性评价指标，全面优于分形维D1、D2。
步骤S4中根据渗透率和孔隙度的分布范围划分储层类型，包括Ⅰ类储层、Ⅱ类储层和Ⅲ类储层，其中Ⅰ类储层对应孔隙度≥15％、渗透率≥100md，Ⅱ类储层对应10％≤孔隙度<15％、10md≤渗透率<100md，Ⅲ类储层对应孔隙度<10％、渗透率<10md。根据实验室测试得到孔隙度与渗透率的数据以及综合分形维DP分别与渗透率、孔隙度之间的关系，确定Ⅰ类储层、Ⅱ类储层和Ⅲ类储层所对应的综合分形维DP的分维值，最终得出综合分形维DP的分维值小于3.14的储层为优质的Ⅰ类储层，综合分形维DP的分维值大于3.25的储层为较差的Ⅲ类储层，综合分形维DP的分维值介于3.14～3.25之间的储层为Ⅱ类储层。按照上述指标分类方案，处理白云凹陷15组样品之后(见表2)，完全符合DP分类标准的样品共7组(孔隙度、渗透率全部符合)，分别为1号、3号、4号、5号、13号、14号、15号；基本符合DP标准(孔隙度、渗透率单项符合)的样品共6组，分别为2号、7号、8号、10号、11号、12号，符合率达到86.67％；对于不符合DP分类标准(孔隙度、渗透率2项指标全不符合)的6号，按照DP标准划分为Ⅱ类储层，然而按照孔隙度、渗透率标准划分为Ⅰ类储层，该样品的孔隙度15.08％与储层分类时Ⅰ类储层、Ⅱ类储层之间的界限15％极为接近，属于两种类型的过渡型储层。样品9号按照孔隙度、渗透率分类属于Ⅰ类储层，而按照DP大于3.14小于3.25，则被划分为Ⅱ类储层，该样品的孔隙度15.9％与分类界限15％接近，DP数值的3.17与DP分类界限的3.14同样比较接近。表2为15组样品最终的评价结果，具体如下：
表2评价结果
总之，利用DP指标划分储层类型符合度达到86.6％，剩余13.4％样品属于过度型储层，在未能直接测试获取孔隙度、渗透率的条件下，DP可以作为一个较好的储层评价指标。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种利用综合分形维Dp定量评价砂岩储层的方法，其特征在于：包括以下步骤：
步骤S1：分段求取压汞曲线的分形维，所述的分形维包括砂岩中大孔隙分形维D1和小孔吼分形维D2；
所述步骤S1中求取大孔隙分形维D1和小孔吼分形维D2的步骤为：
步骤S1.1：分别对压汞获取的毛细管压力Pc和进汞饱和度SHg求取相应的对数Log（Pc）、Log（SHg），并以Log（Pc）作为横坐标、Log（SHg）作为纵坐标，然后将计算后的数据点投影到坐标系中进行拟合得到压汞曲线，并根据坐标系中拟合曲线斜率变化设置拐点；
步骤S1.2：取排驱压力至拐点压力之间的拟合线段的斜率K1，然后按照公式D1=3+ K1计算出砂岩中大孔隙分形维D1；
步骤S1.3：取大于拐点压力部分拟合曲线的斜率K2，然后按照公式D2=3+ K2计算出砂岩中小孔吼分形维D2；
步骤S2：根据所得的大孔隙分形维D1和小孔吼分形维D2，求取压汞曲线的综合分形维DP；
所述步骤S2中求取压汞曲线的综合分形维DP的计算公式为：
DP=D1*ΔPD1/P+ D2*ΔPD2/P；
式中：ΔPD1为拐点压力与排驱压力的差值；ΔPD2为大于拐点压力的高压区毛细管压力差；P为整体毛细管压力差；D1*ΔPD1/P表示砂岩中大孔隙对孔隙复杂程度的贡献值；D2*ΔPD2/P表示砂岩中小孔喉对孔隙复杂程度的贡献值；
步骤S3：采用线性拟合的方法分别建立大孔隙分形维D1、小孔吼分形维D2和综合分形维DP与储层物性参数之间的关系，然后从大孔隙分形维D1、小孔吼分形维D2和综合分形维DP中选出与储层物性相关系数最高的分形维作为最优储层物性评价指标；
步骤S3中所述的储层物性参数包括渗透率、孔隙度、饱和度中值压力和排驱压力；
所述步骤S3中选出最优储层物性评价指标的步骤为：
步骤S3.1：以大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值分别作为横坐标，以储层物性相关的渗透率、孔隙度、饱和度中值压力和排驱压力这四项参数分别作为纵坐标，然后采用线性拟合的方法分别拟合各分维值与四项储层物性参数之间的关系，得出各分维值与四项储层物性参数的相关系数；
步骤S3.2：分别对比大孔隙分形维D1、小孔吼分形维D2和综合分形维DP的分维值与四项储层物性参数之间相关系数的大小，最终得出相关系数最大的综合分形维DP作为最优储层物性评价指标；
步骤S4：根据储层物性参数划分储层类型，然后根据最优储层物性评价指标与储层物性参数之间的关系，确定出不同类型储层对应的分形维指标阈值范围、并作为评价不同砂岩储层的标准。
2.根据权利要求1所述的利用综合分形维Dp定量评价砂岩储层的方法，其特征在于：所述步骤S4中根据渗透率和孔隙度的分布范围划分储层类型，包括Ⅰ类储层、Ⅱ类储层和Ⅲ类储层，其中Ⅰ类储层对应孔隙度≥15%、渗透率≥100md，Ⅱ类储层对应10%≤孔隙度<15%、10md≤渗透率<100md，Ⅲ类储层对应孔隙度<10%、渗透率<10md。
3.根据权利要求2所述的利用综合分形维Dp定量评价砂岩储层的方法，其特征在于：根据综合分形维DP分别与渗透率、孔隙度之间的关系，确定Ⅰ类储层、Ⅱ类储层和Ⅲ类储层所对应的综合分形维DP的分维值，最终得出综合分形维DP的分维值小于3.14的储层为Ⅰ类储层，综合分形维DP的分维值大于3.25的储层为Ⅲ类储层，综合分形维DP的分维值介于3.14～3.25之间的储层为Ⅱ类储层。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110320805.4A CN113075105B (zh) | 2021-03-25 | 2021-03-25 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110320805.4A CN113075105B (zh) | 2021-03-25 | 2021-03-25 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113075105A CN113075105A (zh) | 2021-07-06 |
| CN113075105B true CN113075105B (zh) | 2023-06-23 |
Family
ID=76611616
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110320805.4A Active CN113075105B (zh) | 2021-03-25 | 2021-03-25 | 一种利用综合分形维Dp定量评价砂岩储层的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113075105B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113343502B (zh) * | 2021-07-09 | 2022-11-15 | 成都地铁运营有限公司 | 一种轨道打磨确定方法 |
| CN113642656B (zh) * | 2021-08-18 | 2023-09-05 | 中国石油大学(北京) | 一种低渗砂岩储层开采方式确定方法及相关装置 |
| CN114386881B (zh) * | 2022-03-22 | 2022-06-10 | 中国石油大学(华东) | 一种中高渗油气藏出砂模式划分及固体控制与管理方法 |
| CN115711120B (zh) * | 2022-12-02 | 2024-04-12 | 西南石油大学 | 一种油气井水泥环完整性及损伤演化定量评价方法 |
| CN116297111B (zh) * | 2023-05-24 | 2023-09-05 | 中国石油大学(华东) | 一种基于压汞及分形理论的致密砂岩渗透率预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106124378A (zh) * | 2016-06-16 | 2016-11-16 | 中国石油大学(华东) | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 |
| CN111398122A (zh) * | 2020-04-03 | 2020-07-10 | 中国矿业大学 | 一种页岩全尺度孔隙结构非均质性特征的综合表征方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2811760B1 (fr) * | 2000-07-17 | 2002-09-13 | Inst Francais Du Petrole | Methode pour modeliser des deplacements de fluides dans un milieu poreux tenant compte d'effets d'hysteresis |
-
2021
- 2021-03-25 CN CN202110320805.4A patent/CN113075105B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106124378A (zh) * | 2016-06-16 | 2016-11-16 | 中国石油大学(华东) | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 |
| CN111398122A (zh) * | 2020-04-03 | 2020-07-10 | 中国矿业大学 | 一种页岩全尺度孔隙结构非均质性特征的综合表征方法 |
Non-Patent Citations (2)
| Title |
|---|
| 低渗透碎屑岩储层孔隙结构分形维数计算方法――以川中地区须家河组储层41块岩样为例;赖锦 等;东北石油大学学报;37(01);全文 * |
| 准噶尔盆地玛湖凹陷百口泉组致密砂砾岩孔隙分形特征及影响因素探讨;司马立强 等;测井技术;40(05);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113075105A (zh) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113075105B (zh) | 一种利用综合分形维Dp定量评价砂岩储层的方法 | |
| Liu et al. | Nanoscale pore structure characterization of the Bakken shale in the USA | |
| CN107132171B (zh) | 一种基于压汞-氮吸附联测数据确定致密储层孔径分布的方法 | |
| Xiao et al. | A full-scale characterization method and application for pore-throat radius distribution in tight oil reservoirs | |
| CN110032744B (zh) | 一种非常规油气储量计算方法及*** | |
| Zhao et al. | Fractal study on pore structure of tight sandstone based on full-scale map | |
| CN111827996B (zh) | 基于力学性质的多参数综合定性致密砂岩储层分类方法 | |
| CN115078210B (zh) | 页岩孔隙结构测试方法 | |
| CN106932312B (zh) | 多孔介质中微球状态动态分级表征方法及装置 | |
| CN111537544A (zh) | 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法 | |
| CN110909466B (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
| CN110715879B (zh) | 基于气水分布的高演化页岩储层微孔隙吸附气量评价方法 | |
| CN110029989B (zh) | 一种非常规油气采出程度计算方法及*** | |
| CN106370577B (zh) | 用于缝洞油藏的模拟试验的方法 | |
| CN106285658B (zh) | 一种毛管力曲线校正的方法及装置 | |
| CN112014906B (zh) | 一种致密储层评价方法 | |
| CN112031753B (zh) | 一种随钻识别低渗气水层的方法 | |
| Song et al. | Pore fractal characteristics of lignite at different temperatures based on mercury intrusion test | |
| Bennion et al. | Determination of true effective in situ gas permeability in subnormally water-saturated tight gas reservoirs | |
| CN113670791A (zh) | 一种对储层孔隙及颗粒表面积的定量分析方法 | |
| CN113567319A (zh) | 一种页岩油气储层微观孔隙形态识别方法及装置 | |
| CN113239955A (zh) | 一种碳酸盐岩储层岩石分类方法 | |
| CN113176186A (zh) | 碳酸盐岩原油注入压力和孔喉下限的判断方法 | |
| CN117607388A (zh) | 一种定量表征煤体裂隙率和孔隙率的方法 | |
| CN115753866B (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 |