CN113969777A - 一种储层含油饱和度确定方法、存储介质及电子设备 - Google Patents
一种储层含油饱和度确定方法、存储介质及电子设备 Download PDFInfo
- Publication number
- CN113969777A CN113969777A CN202111364543.8A CN202111364543A CN113969777A CN 113969777 A CN113969777 A CN 113969777A CN 202111364543 A CN202111364543 A CN 202111364543A CN 113969777 A CN113969777 A CN 113969777A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- determining
- effective
- porosity
- water
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 136
- 239000011435 rock Substances 0.000 claims abstract description 92
- 239000011148 porous material Substances 0.000 claims description 53
- 238000005481 NMR spectroscopy Methods 0.000 claims description 38
- 238000002474 experimental method Methods 0.000 claims description 25
- 238000001228 spectrum Methods 0.000 claims description 21
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 229920006395 saturated elastomer Polymers 0.000 claims description 10
- 230000008569 process Effects 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 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
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 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
- E21B47/00—Survey of boreholes 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
- 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
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
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)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种储层含油饱和度确定方法及相关设备，该方法包括：确定储层的第一束缚水孔隙度和第一有效孔隙度；利用岩心实验数据，确定储层的岩电参数；根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度；根据储层的第一有效含水饱和度，确定储层的有效含油饱和度。本发明提供的储层含油饱和度确定方法在确定储层的有效含油饱和度的过程中，通过引入储层的第一束缚水孔隙度和第一有效孔隙度，并利用岩心实验数据对储层的岩电参数进行确定，可以大幅提升对储层含油饱和度确定时的精度。
Description
技术领域
本发明涉及油气勘探储层评价技术领域，尤其涉及一种储层含油饱和度确定方法、存储介质及电子设备。
背景技术
随着油田勘探开发工作日趋深入，以泥质砂岩储层为主的岩性油气藏成为油田勘探开发的重点。然而，适合常规分选性好的纯砂岩储层的含油饱和度解释模型已经不能满足复杂储层评价的需要，如何提高对泥质砂岩储层含油饱和度的计算精度成为储层评价中面临的重要问题。因此，需要建立一种适合泥质砂岩储层的含油饱和度计算方法。
发明内容
本发明旨在至少解决现有技术或相关技术中存在的技术问题之一。
为此，本发明的第一方面提供了一种储层含油饱和度确定方法。
本发明的第二方面提供了一种存储介质。
本发明的第三方面提供了一种电子设备。
有鉴于此，根据本申请实施例的第一方面提出了一种储层含油饱和度确定方法，包括：
确定储层的第一束缚水孔隙度和第一有效孔隙度；
利用岩心实验数据，确定储层的岩电参数；
根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度；
根据第一有效含水饱和度，确定储层的有效含油饱和度。
在一种可行的实施方式中，确定储层的第一束缚水孔隙度和第一有效孔隙度的步骤，包括：
利用核磁共振测井数据，获取储层的第一核磁共振标准T2谱；
利用核磁共振岩心实验数据，确定第一T2截止值；
根据第一核磁共振标准T2谱和第一T2截止值，确定第一束缚水孔隙度，并确定第一有效孔隙度。
在一种可行的实施方式中，利用核磁共振岩心实验数据，确定第一T2截止值的步骤，包括：
利用核磁共振岩心实验数据，获取多个岩心样品的第二T2截止值；
确定多个第二T2截止值的平均值作为第一T2截止值。
在一种可行的实施方式中，根据第一核磁共振标准T2谱和第一T2截止值，确定第一束缚水孔隙度，并确定第一有效孔隙度的步骤，包括：
将第一核磁共振标准T2谱中，横向弛豫时间小于或等于第一T2截止值的各个T2区间对应的孔隙度之和确定为第一束缚水孔隙度；
将第一核磁共振标准T2谱中，横向弛豫时间大于第一T2截止值的各个T2区间对应的孔隙度之和确定为第一有效孔隙度。
在一种可行的实施方式中，储层的岩电参数包括：
有效孔隙的胶结指数，微孔隙的胶结指数，第一有效孔隙的岩性系数，微孔隙的岩性系数，第二有效孔隙的岩性系数和饱和度指数中的至少一种系数。
在一种可行的实施方式中，利用岩心实验数据，确定储层的岩电参数的步骤，包括：
获取岩心实验数据中多个岩心样品的饱含水岩石电阻率；
获取岩心实验数据中实验用水电阻率；
根据饱含水岩石电阻率和实验用水电阻率，确定岩心样品的地层因数；
获取岩心实验数据中多个岩心样品的第二有效孔隙度；
获取岩心实验数据中多个岩心样品的第二束缚水孔隙度；
建立岩心样品的地层因数、第二有效孔隙度与第二束缚水孔隙度之间的第一关系；
根据第一关系，分析得到有效孔隙的胶结指数、微孔隙的胶结指数、第一有效孔隙的岩性系数和微孔隙的岩性系数。
在一种可行的实施方式中，利用岩心实验数据，确定储层的岩电参数的步骤，还包括：
获取岩心实验数据中，每个岩心样品在多个不同的总含水饱和度下的相应的第一岩石电阻率；
根据第一岩石电阻率，确定相应的岩心样品的电阻增大率；
根据岩心样品的总含水饱和度，确定相应的岩心样品的第二有效含水饱和度；
建立岩心样品的电阻增大率与第二有效含水饱和度之间的第二关系；
根据第二关系，分析得到第二有效孔隙的岩性系数和饱和度指数。
在一种可行的实施方式中，储层含油饱和度确定方法还包括：
获取储层的第二岩石电阻率和储层水电阻率；
根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度的步骤可以为：
根据第一束缚水孔隙度、第一有效孔隙度、岩电参数、第二岩石电阻率和储层水电阻率，确定储层的第一有效含水饱和度。
在一种可行的实施方式中，第二岩石电阻率可采用通过储层的电阻率测井得到的深探测电阻率；
储层水电阻率可采用储层的区块经验参数值。
为了实现上述目的，根据本发明的第二方面，提供了一种存储介质，存储介质包括存储的程序，其中，在程序运行时控制存储介质所在设备执行如上述第一方面中任一项提出的储层含油饱和度确定方法。
为了实现上述目的，根据本发明的第三方面，提供了一种电子设备，电子设备包括至少一个处理器、以及与处理器连接的至少一个存储器；其中，处理器用于调用存储器中的程序指令，执行如上述第一方面中任一项提出的储层含油饱和度确定方法。
相比现有技术，本发明至少包括以下有益效果：本发明提供的储层含油饱和度确定方法包括：确定储层的第一束缚水孔隙度和第一有效孔隙度；利用岩心实验数据，确定储层的岩电参数；根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度；根据储层的第一有效含水饱和度，确定储层的有效含油饱和度。该储层含油饱和度确定方法在确定储层的有效含油饱和度的过程中，通过引入储层的第一束缚水孔隙度和第一有效孔隙度，从而综合考虑了储层中微孔隙和有效孔隙形成的不同导电路径对储层整体导电特性的影响；并且，通过利用岩心实验数据对储层的岩电参数进行确定，可以进一步使得确定过程中涉及到的变量和参数能够与岩心实验数据进行对比。进而，本发明提供的储层含油饱和度确定方法可以大幅提升对储层含油饱和度确定时的精度，尤其适用于对如泥质砂岩储层等复杂储层的含油饱和度的确定。
附图说明
通过阅读下文示例性实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出示例性实施方式的目的，而并不认为是对本申请的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1为本申请提供的一种实施例的储层含油饱和度确定方法的示意性流程图；
图2为本申请提供的一种电子设备的示意性结构框图。
具体实施方式
下面将参照附图更详细地描述本申请的示例性实施例。虽然附图中显示了本申请的示例性实施例，然而应当理解，可以以各种形式实现本申请而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本申请，并且能够将本申请的范围完整的传达给本领域的技术人员。
如图1所示，根据本申请实施例的第一方面提出了一种储层含油饱和度确定方法，包括：
步骤101：确定储层的第一束缚水孔隙度和第一有效孔隙度；
具体地，第一束缚水孔隙度Φic和第一有效孔隙度Φfc可通过储层的测井数据计算确定。测井种类包括但不限于电测井、声波测井和核测井。
步骤102：利用岩心实验数据，确定储层的岩电参数；
具体地，根据对储层进行岩心实验得到的实验数据，对储层的岩电参数进行确定，从而可以保证得到的岩电参数能够更加准确地反应储层中的岩石特性，进而确保确定储层有效含油饱和度的过程中所涉及到的变量和参数能够与通过岩心实验得到的实测数据进行对比，提升确定的储层有效含油饱和度的精度。
步骤103：根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度；
具体地，根据前述确定的第一束缚水孔隙度Φic、第一有效孔隙度Φfc和岩电参数，对储层的第一有效含水饱和度Swfc进行计算确定。其中，由于引入了第一束缚水孔隙度Φic和第一有效孔隙度Φfc，从而综合考虑了储层中微孔隙和有效孔隙对储层导电特性的影响，提升了对储层含水情况的确定精度。
步骤104：根据第一有效含水饱和度，确定储层的有效含油饱和度。
具体地，有效含水饱和度Swf和有效含油饱和度Sof之间具有如下关系：
Sof＝1-Swf (1)
式中，有效含水饱和度Swf和有效含油饱和度Sof的单位均为v/v，也即体积比。
将前述确定的第一有效含水饱和度Swfc代入式(1)，可计算得到储层的有效含油饱和度Sofc，即：
Sofc＝1-Swfc (2)
需要说明的是，一些含油饱和度计算方法的理论基础主要包括经典的Archie公式和考虑多重孔隙影响的饱和度解释模型。经典的Archie公式没有考虑束缚水的影响，考虑多重孔隙影响的饱和度解释模型中的参数和变量无法与岩心实测数据进行对比，进而降低了含油饱和度的计算精度。
本发明提供的储层含油饱和度确定方法是基于以下基本思想建立的：
其一是，岩石电性的主要影响因素有孔隙空间中可以流动的有效孔隙和完全由束缚水占据的微孔隙；
其二是，有效孔隙和微孔隙形成的不同导电路径以并联的方式影响储层整体的导电特性；
其三是，有效孔隙和微孔隙具有不同的地层因数和胶结指数；
其四是，油气只能取代孔隙空间中的有效孔隙。
综上，在确定第一有效含水饱和度Swfc的过程中，引入了第一束缚水孔隙度Φic、第一有效孔隙度Φfc以及通过岩心实验数据确定的岩电参数，从而综合考虑了储层中微孔隙和有效孔隙对储层导电特性的影响，并使得确定过程中涉及到的参数与变量能够和岩心实验的实测数据进行对比，进而大幅提升了确定得到的储层的有效含油饱和度Sofc的精确度。
在一些示例中，获取储层的第一束缚水孔隙度和第一有效孔隙度的步骤，包括：
利用核磁共振测井数据，获取储层的第一核磁共振标准T2谱；
利用核磁共振岩心实验数据，确定第一T2截止值；
根据第一核磁共振标准T2谱和第一T2截止值，确定第一束缚水孔隙度，并确定第一有效孔隙度。
具体地，以储层的核磁共振测井数据为依据，可获取储层的第一核磁共振标准T2谱；以核磁共振岩心实验数据为依据，可确定第一T2截止值。进而根据第一核磁共振标准T2谱和第一T2截止值，对第一束缚水孔隙度Φic和第一有效孔隙度Φfc进行确定。
需要说明的是，通过核磁共振测井确定储层孔隙度Φ的依据来自观测信号强度与孔隙流体中氢核数的对应关系，观测信号在零时刻的数值大小与储层孔隙中的含氢总量成正比，经过刻度，可把零时刻的信号强度标定为储层孔隙度Φ，从而通过核磁共振标准T2谱可以得到T2(横向弛豫时间)与孔隙度P之间的关系，并通过如下公式计算储层孔隙度Φ：
式中，Φ为储层孔隙度，Pi为第一核磁共振标准T2谱中第i个T2区间对应的孔隙度，i为正整数。其中，Φ和Pi的单位均为v/v。
进而，根据第一T2截止值，认为横向弛豫时间T2小于或等于该值的所有孔隙中的流体均为束缚状态的，并在储层压力条件下是不能流动的；认为横向弛豫时间T2大于该值的所有孔隙中的流体在储层压力条件下是可以流动的，可以分别得到储层的第一束缚水孔隙度Φic和第一有效孔隙度Φfc。
在一些示例中，利用核磁共振岩心实验数据，确定第一T2截止值的步骤，包括：
利用核磁共振岩心实验数据，获取多个岩心样品的第二T2截止值；
确定多个第二T2截止值的平均值作为第一T2截止值。
具体的，储层的第一束缚水孔隙度Φic和第一有效孔隙度Φfc的确定结果的准确度与储层的第一T2截止值存在密切联系。根据核磁共振岩心实验数据，可以获得储层多个岩心样品的第二T2截止值，并将多个第二T2截止值的平均值作为第一T2截止值，从而可以使储层的第一T2截止值的确定结果更加准确，有利于提高第一束缚水孔隙度Φic和第一有效孔隙度Φfc的确定结果的准确性。
在一些可行的示例中，多个岩心样品可以来自于同一核磁共振测井的同一深度范围。
在一些可行的示例中，根据第一核磁共振标准T2谱和第一T2截止值，确定第一束缚水孔隙度，并确定第一有效孔隙度的步骤，包括：
将第一核磁共振标准T2谱中，横向弛豫时间小于或等于第一T2截止值的各个T2区间对应的孔隙度之和确定为第一束缚水孔隙度；
将第一核磁共振标准T2谱中，横向弛豫时间大于第一T2截止值的各个T2区间对应的孔隙度之和确定为第一有效孔隙度。
具体地，第一束缚水孔隙度Φic可通过下式计算：
式中，T2cutoff为第一T2截止值。其中，横向弛豫时间T2和第一T2截止值T2cutoff的单位均为ms，第一束缚水孔隙度Φic的单位为v/v。
也即对第一核磁共振标准T2谱中，横向弛豫时间T2小于或等于第一T2截止值T2cutoff的各T2区间对应的孔隙度Pi进行求和，得到第一束缚水孔隙度Φic。
第一有效孔隙度Φfc可通过下式计算：
其中，第一有效孔隙度Φfc的单位为v/v。
也即对第一核磁共振标准T2谱中，横向弛豫时间T2大于第一T2截止值T2cutoff的各T2区间对应的孔隙度Pi进行求和，得到第一有效孔隙度Φfc。
在一些示例中，储层的岩电参数包括：
有效孔隙的胶结指数，微孔隙的胶结指数，第一有效孔隙的岩性系数，微孔隙的岩性系数，第二有效孔隙的岩性系数和饱和度指数中的至少一种系数。
具体地，如上面提到的本方法所基于的基本思想，在根据岩心实验数据确定储层的岩电参数时，岩电参数可以包括有效孔隙的胶结指数mf，微孔隙的胶结指数mi，第一有效孔隙的岩性系数af，微孔隙的岩性系数ai，第二有效孔隙的岩性系数b和饱和度指数n，从而对储层中的有效孔隙和微孔隙带来的影响，进行针对性的考量，有助于提升对储层含油饱和度的确定结果的准确性。
其中，有效孔隙的胶结指数mf与储层中有效孔隙部分岩石的胶结情况和孔隙结构有关；微孔隙的胶结指数mi与储层中微孔隙部分岩石的胶结情况和孔隙结构有关；第一有效孔隙的岩性系数af和第二有效孔隙的岩性系数b与储层中有效孔隙部分岩石的岩性有关；微孔隙的岩性系数ai与储层中微孔隙部分岩石的岩性有关；饱和度指数n与储层中油、气、水在有效孔隙中的分布状况有关。
在一些示例中，利用岩心实验数据，确定储层的岩电参数的步骤，包括：
获取岩心实验数据中多个岩心样品的饱含水岩石电阻率；
获取岩心实验数据中实验用水电阻率；
根据饱含水岩石电阻率和实验用水电阻率，确定岩心样品的地层因数；
获取岩心实验数据中多个岩心样品的第二有效孔隙度；
获取岩心实验数据中多个岩心样品的第二束缚水孔隙度；
建立岩心样品的地层因数、第二有效孔隙度与第二束缚水孔隙度之间的第一关系；
根据第一关系，分析得到有效孔隙的胶结指数、微孔隙的胶结指数、第一有效孔隙的岩性系数和微孔隙的岩性系数。
具体地，本发明基于前述基本思想，提出了一种基于双孔隙的导电体积解释模型，该模型如下式所示：
式中，Rt为含油气层岩石电阻率，单位为Ω·m；Φf为有效孔隙度，单位为v/v；Φi为束缚水孔隙度，单位为v/v；Rw为水电阻率，单位为Ω·m。
当岩石饱含水时，式(6)可写成：
式中，Ro为饱含水岩石电阻率，单位为Ω·m。
令：
式中，F为地层因数，为无因次项。
通过式(7)和式(8)可以得到：
将第二有效孔隙度Φfy和第二束缚水孔隙度Φiy代入式(9)，可得到：
式中，Fy为岩心样品的地层因数；Φfy为第二有效孔隙度；Φiy为第二束缚水孔隙度。
式(10)即岩心样品的地层因数Fy、第二有效孔隙度Φfy与第二束缚水孔隙度Φiy之间的第一关系。
从而，通过岩心实验数据中的岩石电阻率实验数据，获取其中的多个岩心样品的饱含水岩石电阻率Roy以及实验用水电阻率Rwy，根据式(8)计算出各岩心样品的地层因数Fy，即：
并通过岩心实验数据中的核磁共振实验数据，确定多个岩心样品的第二有效孔隙度Φfy和第二束缚水孔隙度Φiy；根据式(10)所示出的第一关系，结合多个岩心样品的第二有效孔隙度Φfy和第二束缚水孔隙度Φiy以及多个岩心样品的地层因数Fy，可以分析得到有效孔隙的胶结指数mf，微孔隙的胶结指数mi，第一有效孔隙的岩性系数af和微孔隙的岩性系数ai。
同时，根据岩石电阻率实验的测量原理，当岩心样品的第二有效孔隙度Φfy为1.0000v/v，并且第二束缚水孔隙度Φiy为0.0000v/v时，岩石电阻率实验测量岩心样品的饱含水岩石电阻率Roy为实验用水电阻率Rwy，根据公式(11)计算地层因数Fy为1.00。在这一情形下的Fy、Φfy和Φiy值亦可引入对mf、mi、af和ai这4种系数的分析之中。
在一些可行的示例中，根据第一关系，分析得到有效孔隙的胶结指数、微孔隙的胶结指数、第一有效孔隙的岩性系数和微孔隙的岩性系数的方法可以为多元非线性回归法。
在一些示例中，利用岩心实验数据，确定储层的岩电参数的步骤，还包括：
获取岩心实验数据中，每个岩心样品在多个不同的总含水饱和度下的相应的第一岩石电阻率；
根据第一岩石电阻率，确定相应的岩心样品的电阻增大率；
根据岩心样品的总含水饱和度，确定相应的岩心样品的第二有效含水饱和度；
建立岩心样品的电阻增大率与第二有效含水饱和度之间的第二关系；
根据第二关系，分析得到第二有效孔隙的岩性系数和饱和度指数。
具体地，当岩石含油气时，式(6)可写成：
令：
式中，I为电阻增大率，为无因次项。
通过式(12)和式(13)可以得到：
将第二有效含水饱和度代入式(14)，可得到：
式中，Iy为岩心样品的电阻增大率；Swfy为第二有效含水饱和度。
式(15)即岩心样品的电阻增大率Iy与第二有效含水饱和度Swfy之间的第二关系。
从而，通过岩心实验数据中的岩石电阻率实验数据，获取其中的多个岩心样品在多个不同总含水饱和度Swy下的第一岩石电阻率Rty；在有效孔隙的胶结指数mf，微孔隙的胶结指数mi，第一有效孔隙的岩性系数af和微孔隙的岩性系数ai，以及每个岩心样品的第二有效孔隙度Φfy、第二束缚水孔隙度Φiy和实验用水电阻率Rwy通过前述步骤确定了的情况下，根据式(13)计算出每个岩心样品在各个总含水饱和度Swy下，相应的岩心样品的电阻增大率Iy，即：
同时，总含水饱和度Sw与有效含水饱和度Swf之间具有如下关系：
进而，根据岩心样品的总含水饱和度Swy，结合式(17)，可以确定相应的岩心样品的第二有效含水饱和度Swfy，即：
随即根据如式(15)所示出的第二关系，结合每个岩心样品的多个电阻增大率Iy和相应的多个第二有效含水饱和度Swfy，可以分析得到第二有效孔隙的岩性系数b和饱和度指数n。
在一些示例中，储层含油饱和度确定方法还包括：
获取储层的第二岩石电阻率和储层水电阻率；
根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度的步骤可以为：
根据第一束缚水孔隙度、第一有效孔隙度、岩电参数、第二岩石电阻率和储层水电阻率，确定储层的第一有效含水饱和度。
具体地，在获取了储层的第二岩石电阻率Rtc和储层水电阻率Rwc的情况下，可以根据第一束缚水孔隙度Φic、第一有效孔隙度Φfc、岩电参数、第二岩石电阻率Rtc和储层水电阻率Rwc，确定储层的第一有效含水饱和度Swfc。
确定储层的第一有效含水饱和度Swfc的计算式可以通过式(6)得到，即：
将第一束缚水孔隙度Φic、第一有效孔隙度Φfc、岩电参数、第二岩石电阻率Rtc和储层水电阻率Rwc代入式(19)，得到：
在一些示例中，第二岩石电阻率可采用通过储层的电阻率测井得到的深探测电阻率；储层水电阻率可采用储层的区块经验参数值。
具体地，将储层的电阻率测井得到的深探测电阻率作为第二岩石电阻率Rtc，可以保证第二岩石电阻率Rtc的数值准确度；将储层的区块经验参数值作为储层水电阻率Rwc，可以在保证储层水电阻率Rwc具有一定准确度的同时，更加易于获得，节省了实验测定的过程。
根据本申请实施例的第二方面还提出了一种存储介质，其上存储有程序，该程序运行时控制该存储介质所在设备执行如上述第一方面中任一项提出的储层含油饱和度确定方法的步骤。
根据本申请实施例的第三方面还提出了一种电子设备200，如图2所示，设备包括至少一个处理器201、以及与处理器201连接的至少一个存储器202；其中，处理器201用于调用存储器202中的程序指令，以执行实现如上述第一方面中任一项提出的储层含油饱和度确定方法的步骤。
借由上述技术方案，本申请实施例提供了一种储层含油饱和度确定方法、存储介质及电子设备。储层含油饱和度确定方法包括：确定储层的第一束缚水孔隙度和第一有效孔隙度；利用岩心实验数据，确定储层的岩电参数；根据第一束缚水孔隙度、第一有效孔隙度和岩电参数，确定储层的第一有效含水饱和度；根据储层的第一有效含水饱和度，确定储层的有效含油饱和度。该储层含油饱和度确定方法在确定储层的有效含油饱和度的过程中，通过引入储层的第一束缚水孔隙度和第一有效孔隙度，从而综合考虑了储层中微孔隙和有效孔隙形成的不同导电路径对储层整体导电特性的影响；并且，通过利用岩心实验数据对储层的岩电参数进行确定，可以进一步使得确定过程中涉及到的变量和参数能够与岩心实验数据进行对比。进而，本发明提供的储层含油饱和度确定方法可以大幅提升对储层含油饱和度确定时的精度，尤其适用于对如泥质砂岩储层等复杂储层的含油饱和度的确定。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的；应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合；可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程流程管理设备的处理器以产生一个机器，使得通过计算机或其他可编程流程管理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
在一个典型的配置中，电子设备可以包括一个或多个处理器(CPU)、存储器和总线；电子设备还可以包括输入/输出接口、网络接口等。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片；存储器是存储介质的示例。
存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储；信息可以是计算机可读指令、数据结构、程序的模块或其他数据；计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息；按照本文中的界定，存储介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本申请的说明书和权利要求书及附图中的术语“第一”、“第二”、“第三”、“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序；应该理解这样使用的数据在适当情况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素；在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品；因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式；而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
可以由一种或多种程序设计语言或其组合来编写用于执行本申请实施例操作的计算机程序代码，程序设计语言包括面向对象的程序设计语言——诸如Common Lisp、Python、C++、Objective-C、Smalltalk、Delphi、Java、Swift、C#、Perl、Ruby、JavaScript和PHP等，还包括常规的过程式程序设计语言——诸如Fortran、ALGOL、COBOL、PL/I、BASIC、Pascal和C等，还包括其他任意一种编程语言——诸如Lisp、Tcl、Prolog、VisualBasic.NET、SQL和R等；程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行；在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)——连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
以上仅为本申请的实施例而已，并不用于限制本申请；对于本领域技术人员来说，本申请可以有各种更改和变化；凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (11)
1.一种储层含油饱和度确定方法，其特征在于，包括：
确定所述储层的第一束缚水孔隙度和第一有效孔隙度；
利用岩心实验数据，确定所述储层的岩电参数；
根据所述第一束缚水孔隙度、所述第一有效孔隙度和所述岩电参数，确定所述储层的第一有效含水饱和度；
根据所述第一有效含水饱和度，确定所述储层的有效含油饱和度。
2.根据权利要求1所述的储层含油饱和度确定方法，其特征在于，所述确定所述储层的第一束缚水孔隙度和第一有效孔隙度的步骤，包括：
利用核磁共振测井数据，获取所述储层的第一核磁共振标准T2谱；
利用核磁共振岩心实验数据，确定第一T2截止值；
根据所述第一核磁共振标准T2谱和第一T2截止值，确定所述第一束缚水孔隙度，并确定所述第一有效孔隙度。
3.根据权利要求2所述的储层含油饱和度确定方法，其特征在于，所述利用核磁共振岩心实验数据，确定第一T2截止值的步骤，包括：
利用所述核磁共振岩心实验数据，获取多个岩心样品的第二T2截止值；
确定多个所述第二T2截止值的平均值作为所述第一T2截止值。
4.根据权利要求2所述的储层含油饱和度确定方法，其特征在于，所述根据所述第一核磁共振标准T2谱和第一T2截止值，确定所述第一束缚水孔隙度，并确定所述第一有效孔隙度的步骤，包括：
将所述第一核磁共振标准T2谱中，横向弛豫时间小于或等于所述第一T2截止值的各个T2区间对应的孔隙度之和确定为所述第一束缚水孔隙度；
将所述第一核磁共振标准T2谱中，所述横向弛豫时间大于所述第一T2截止值的各个T2区间对应的孔隙度之和确定为所述第一有效孔隙度。
5.根据权利要求1所述的储层含油饱和度确定方法，其特征在于，所述储层的岩电参数包括：
有效孔隙的胶结指数，微孔隙的胶结指数，第一有效孔隙的岩性系数，微孔隙的岩性系数，第二有效孔隙的岩性系数和饱和度指数中的至少一种系数。
6.根据权利要求5所述的储层含油饱和度确定方法，其特征在于，所述利用岩心实验数据，确定所述储层的岩电参数的步骤，包括：
获取所述岩心实验数据中多个岩心样品的饱含水岩石电阻率；
获取所述岩心实验数据中实验用水电阻率；
根据所述饱含水岩石电阻率和所述实验用水电阻率，确定所述岩心样品的地层因数；
获取所述岩心实验数据中多个所述岩心样品的第二有效孔隙度；
获取所述岩心实验数据中多个所述岩心样品的第二束缚水孔隙度；
建立所述岩心样品的地层因数、所述第二有效孔隙度与所述第二束缚水孔隙度之间的第一关系；
根据所述第一关系，分析得到所述有效孔隙的胶结指数、所述微孔隙的胶结指数、所述第一有效孔隙的岩性系数和所述微孔隙的岩性系数。
7.根据权利要求6所述的储层含油饱和度确定方法，其特征在于，所述利用岩心实验数据，确定所述储层的岩电参数的步骤，还包括：
获取所述岩心实验数据中，每个所述岩心样品在多个不同的总含水饱和度下的相应的第一岩石电阻率；
根据所述第一岩石电阻率，确定相应的所述岩心样品的电阻增大率；
根据所述岩心样品的总含水饱和度，确定相应的所述岩心样品的第二有效含水饱和度；
建立所述岩心样品的电阻增大率与所述第二有效含水饱和度之间的第二关系；
根据所述第二关系，分析得到所述第二有效孔隙的岩性系数和所述饱和度指数。
8.根据权利要求7所述的储层含油饱和度确定方法，其特征在于，还包括：
获取所述储层的第二岩石电阻率和储层水电阻率；
所述根据所述第一束缚水孔隙度、所述第一有效孔隙度和所述岩电参数，确定所述储层的第一有效含水饱和度的步骤可以为：
根据所述第一束缚水孔隙度、所述第一有效孔隙度、所述岩电参数、所述第二岩石电阻率和所述储层水电阻率，确定所述储层的第一有效含水饱和度。
9.根据权利要求8所述的储层含油饱和度确定方法，其特征在于，
所述第二岩石电阻率为通过所述储层的电阻率测井得到的深探测电阻率；
所述储层水电阻率为所述储层的区块经验参数值。
10.一种存储介质，其特征在于，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行如权利要求1至权利要求9中任一项所述的储层含油饱和度确定方法。
11.一种电子设备，其特征在于，所述电子设备包括至少一个处理器、以及与所述处理器连接的至少一个存储器，其中，所述处理器用于调用所述存储器中的程序指令，执行如权利要求1至权利要求9中任一项所述的储层含油饱和度确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111364543.8A CN113969777A (zh) | 2021-11-17 | 2021-11-17 | 一种储层含油饱和度确定方法、存储介质及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111364543.8A CN113969777A (zh) | 2021-11-17 | 2021-11-17 | 一种储层含油饱和度确定方法、存储介质及电子设备 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113969777A true CN113969777A (zh) | 2022-01-25 |
Family
ID=79589908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111364543.8A Pending CN113969777A (zh) | 2021-11-17 | 2021-11-17 | 一种储层含油饱和度确定方法、存储介质及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113969777A (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN104278989A (zh) * | 2013-07-02 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种获取低孔低渗储层饱和度指数的方法 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| US20190137648A1 (en) * | 2016-08-10 | 2019-05-09 | Petrochina Company Limited | Method and system for determining heterogeneous carbonate reservoir saturation exponent |
-
2021
- 2021-11-17 CN CN202111364543.8A patent/CN113969777A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN104278989A (zh) * | 2013-07-02 | 2015-01-14 | 中国石油天然气股份有限公司 | 一种获取低孔低渗储层饱和度指数的方法 |
| CN105114064A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油天然气股份有限公司 | 确定致密砂岩储层饱和度的方法 |
| US20190137648A1 (en) * | 2016-08-10 | 2019-05-09 | Petrochina Company Limited | Method and system for determining heterogeneous carbonate reservoir saturation exponent |
Non-Patent Citations (1)
| Title |
|---|
| 司马立强;王超;吴丰;王亮;马力;王紫娟;: "川西马井气田蓬莱镇组致密砂岩储层可动水饱和度计算方法", 测井技术, no. 02, pages 77 - 82 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110717249B (zh) | 页岩气储层测井孔隙度快速预测方法及系统 | |
| Guarracino et al. | A physically based analytical model to describe effective excess charge for streaming potential generation in water saturated porous media | |
| US8138763B2 (en) | Method for quantitative separation of electromagnetic induction and induced polarization effects | |
| CN103760082A (zh) | 页岩层系致密储层原油有效可动空间确定方法和装置 | |
| CN103485770B (zh) | 一种基于人工神经网络求取含油饱和度的方法及系统 | |
| Quintana-Murillo et al. | An explicit numerical method for the fractional cable equation | |
| Walker et al. | Permeability models of porous media: Characteristic length scales, scaling constants and time-dependent electrokinetic coupling | |
| CN111665560A (zh) | 油气层识别方法、装置、计算机设备及可读存储介质 | |
| Soleymanzadeh et al. | A new technique for determining water saturation based on conventional logs using dynamic electrical rock typing | |
| Wang et al. | The evaluation of rock permeability with streaming current measurements | |
| CN108681793B (zh) | 深层油藏采油指数预测方法及装置 | |
| CN113969777A (zh) | 一种储层含油饱和度确定方法、存储介质及电子设备 | |
| CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| CN111999331B (zh) | 基于核磁测试的毛管压力确定方法及装置 | |
| Biggin et al. | How many paleointensity determinations are required from a single lava flow to constitute a reliable average? | |
| CN111562629B (zh) | 一种基于等效孔隙截面指数的饱和度确定方法及装置 | |
| Cassiani et al. | A note on in situ estimates of sorption using push‐pull tests | |
| CN115901563A (zh) | 用核磁共振和电阻率联合计算低孔低渗储层渗透率的方法 | |
| Liu et al. | Permeability estimation using relaxation time spectra derived from differential evolution inversion | |
| CN108647417B (zh) | 一种确定页岩气储层含气饱和度的简易方法 | |
| McPhee | The core analysis elephant in the formation evaluation room | |
| CN112147172A (zh) | 基于核磁共振t2谱评估含水饱和度的方法和装置 | |
| CN112861061A (zh) | 地质分层方法及系统、处理设备 | |
| Yeganeh et al. | Using a LDG method for solving an inverse source problem of the time-fractional diffusion equation | |
| CN114165226B (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 |