CN113945497B - 油气藏储层流体的可动性的评价方法 - Google Patents
油气藏储层流体的可动性的评价方法 Download PDFInfo
- Publication number
- CN113945497B CN113945497B CN202010681539.3A CN202010681539A CN113945497B CN 113945497 B CN113945497 B CN 113945497B CN 202010681539 A CN202010681539 A CN 202010681539A CN 113945497 B CN113945497 B CN 113945497B
- Authority
- CN
- China
- Prior art keywords
- core
- curve
- saturated
- fluid
- spectrum
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 51
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 103
- 238000001228 spectrum Methods 0.000 claims abstract description 74
- 238000009826 distribution Methods 0.000 claims abstract description 66
- 238000006243 chemical reaction Methods 0.000 claims abstract description 59
- 239000011148 porous material Substances 0.000 claims abstract description 59
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 48
- 238000006073 displacement reaction Methods 0.000 claims abstract description 45
- 238000009825 accumulation Methods 0.000 claims abstract description 38
- 238000012360 testing method Methods 0.000 claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 13
- 230000001186 cumulative effect Effects 0.000 claims description 34
- 238000004519 manufacturing process Methods 0.000 claims description 21
- 239000007789 gas Substances 0.000 claims description 20
- 239000011435 rock Substances 0.000 claims description 19
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 9
- 239000011521 glass Substances 0.000 claims description 8
- 238000011282 treatment Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 239000004215 Carbon black (E152) Substances 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 5
- 229930195733 hydrocarbon Natural products 0.000 claims description 5
- 150000002430 hydrocarbons Chemical class 0.000 claims description 5
- 229910052757 nitrogen Inorganic materials 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000002791 soaking Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims 1
- 238000009738 saturating Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 5
- 238000012512 characterization method Methods 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 abstract 1
- 230000035699 permeability Effects 0.000 description 7
- 230000018109 developmental process Effects 0.000 description 6
- 230000003068 static effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000008398 formation water Substances 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000005303 weighing Methods 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
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/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Dispersion Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Fluid Mechanics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请公开了一种油气藏储层流体的可动性的评价方法，属于油田开发技术领域。本申请提供的评价方法，首先采用CT孔径分布转换饱和水T2谱曲线，并结合误差分析法求取最佳转换系数；利用相似转换原则获取室内实验压差和流速，并测试得到驱替T2谱曲线；然后根据最佳转换系数C0、饱和水T2谱曲线和驱替T2谱曲线，进一步转换得到饱和流体累计曲线和束缚流体累计曲线，进而评价流体的可动性。本申请通过将CT孔径分布、核磁共振T2谱以及误差分析等结合，可充分模拟不同作业制度对流体的动用过程，不需要依赖大量实验样品，能够实现不同作业制度下流体孔喉半径动用分布的定量化表征，避免人工主观评价，为油气藏合理开发提供依据。
Description
技术领域
本申请涉及油气田开发技术领域，特别涉及一种油气藏储层流体的可动程度和可动孔径下限的评价方法。
背景技术
油气藏储量规模、储量可动性以及不同技术措施下的生产效果是影响油气藏能否投入开发生产直接因素。而影响油气藏开发的经济效益的关键因素在于油气藏的储量可动性，评价储量可动性包括对油气藏储层流体的可动孔径下限、可动用程度以及不同孔喉半径的流体可动用程度的评价。
目前，对储量可动性的评价方法包括静态法和动态法。静态法所用参数为实验结果和经验取值；而动态法所用参数来自实际生产过程中，可根据不同时期调整物性下限值。静态法主要包括压汞实验法、储层物性参数统计法、束缚水饱和度法等，动态法包括测井资料法、产能模拟实验法。
储量可动性评价是油气藏合理开发的关键基础，也是油气藏评价中很重要的一个环节。然而，采用目前的方法进行评价需要大量的岩心物性数据，评价结果主观性强，且不能反映不同压差、流速下不同孔径内的流体动用程度。
发明内容
鉴于此，本申请提供一种油气藏储层流体的可动性的评价方法。
具体而言，包括以下的技术方案：
提供一种油气藏储层流体的可动性的评价方法包括以下步骤：
获取油气藏储层的岩心；利用岩心CT扫描仪测试所述岩心，得到所述岩心的CT孔径分布曲线；
采用模拟地层水浸泡所述岩心，获取饱和岩心；利用核磁设备测试所述饱和岩心，得到所述岩心的饱和水T2谱曲线；
根据所述CT孔径分布曲线和所述饱和水T2谱曲线，得到CT孔径分布累计曲线和饱和水T2谱累计曲线；采用不同的转换系数Ci转换所述饱和水T2谱累计曲线；
根据所述CT孔径分布累计曲线和不同的转换系数Ci转换后的所述饱和水T2谱累计曲线，采用误差分析法计算最佳转换系数Co；
根据所述单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心；采用所述核磁设备测试所述驱替岩心，得到驱替T2谱曲线；
根据所述最佳转换系数C0、所述饱和水T2谱曲线和所述驱替T2谱曲线，得到饱和流体累计曲线和束缚流体累计曲线；
根据所述饱和流体累计曲线和所述束缚流体累计曲线，确定所述油气藏储层的可动孔径下限。
在一种可能的实施例中，所述利用岩心CT扫描仪测试所述岩心，得到所述岩心的CT孔径分布曲线，包括：
利用所述岩心CT扫描仪测试所述岩心，获取所述岩心的三维切片灰度图像；
利用数字岩心软件对所述灰度图像进行图像滤波处理、孔隙三维重构处理和孔隙分割处理，获取所述岩心的CT孔径分布曲线。
在一种可能的实施例中，所述利用核磁设备测试所述饱和岩心，得到所述岩心的饱和水T2谱曲线，包括：
将所述饱和岩心放入玻璃注射器内，使所述饱和岩心与注射器内壁、注射器活塞贴合；
设置所述核磁设备的实验参数，将装有所述饱和岩心的所述玻璃注射器放入核磁共振测试线圈中，保持所述饱和岩心位于所述测试线圈的中部位置；
对所述饱和岩心施加电磁射频信号，得到所述岩心的饱和水T2谱曲线。
在一种可能的实施例中，所述根据所述CT孔径分布累计曲线和不同的转换系数Ci转换后的所述饱和水T2谱累计曲线，计算最佳转换系数Co，包括：
对于Ci转换后的饱和水T2谱累计曲线，分别利用差值法求取T2谱的不同累计百分数对应的弛豫时间Tij；利用差值法求取T2谱百分比分布曲线在不同区间内的分布百分比ω(Tij)；对于所述CT孔径分布累计曲线，利用差值法分别求取不同累计百分数对应的孔喉半径Rj；
利用曲线拟合误差分析式，计算CT孔喉半径与不同Ci转换后的核磁共振T2谱之间的拟合误差δi，所述曲线拟合误差分析式为：
其中，N表示对T2谱百分比分布曲线所划分的区间的总个数，Rj表示对应不同区间的孔喉半径的取值，i表示所采用的不同转换系数的次数，i，j为正整数。
绘制所述拟合误差δi与所述转换系数Ci的关系曲线，并拟合所述拟合误差δi与所述转换系数Ci的关系式，关系式为：δi＝f(Ci)
求取δi′＝0时，对应的C值，该值为所述最佳转换系数Co。
在一种可能的实施例中，所述根据所述单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心，包括：
根据单井的生产资料，确定气井地层压力PR、井底压力Pwf、泄流半径re、井筒半径rw、油气层有效厚度h、产量Q，计算室内实验岩心等效流量QC，计算式为：
其中，D为所述岩心的直径；
根据岩心等效流量QC，利用氮气顺岩心标记方向驱替所述饱和岩心，得到所述驱替岩心。
在一种可能的实施例中，所述根据所述单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心，包括：
根据单井的生产资料，确定气井地层压力PR、井底压力Pwf、泄流半径re，计算室内实验岩心等效压差ΔP，计算式为：
其中L为所述岩心的长度；
根据室内实验岩心等效压差ΔP，利用氮气顺岩心标记方向驱替所述饱和岩心，得到所述驱替岩心。
在一种可能的实施例中，所述根据最佳转换系数C0、饱和水T2谱曲线和驱替T2谱曲线，得到饱和流体累计曲线和束缚流体累计曲线，包括：
利用所述最佳转换系数Co转换所述饱和水T2谱曲线和所述驱替T2谱曲线，得到饱和流体孔径分布曲线和束缚流体孔径分布曲线；将所述饱和流体孔径分布曲线和所述束缚流体孔径分布曲线转换为饱和流体累计曲线和束缚流体累计曲线。
在一种可能的实施例中，根据所述饱和流体累计曲线和所述束缚流体累计曲线，确定所述油气藏储层的可动孔径下限，包括：
对比所述饱和流体累计曲线和所述束缚流体累计曲线，将饱和流体累计饱和度大于束缚流体累计饱和度预设值作为可动界限，对应的孔喉半径作为所述可动孔径下限。
在一种可能的实施例中，所述预设值为5％。
在一种可能的实施例中，所述评价方法还包括：
根据所述饱和流体孔径分布曲线和所述束缚流体孔径分布曲线，计算不同孔喉半径范围内初始流体饱和度和束缚流体饱和度；
根据所述初始流体饱和度和所述束缚流体饱和度，计算不同孔喉半径范围内的流体动用程度。
本申请实施例提供的技术方案的有益效果至少包括：
本申请提供的油气藏储层流体的可动性的评价方法，首先采用CT孔径分布转换饱和水T2谱曲线，并结合误差分析法求取最佳转换系数；利用相似转换原则获取室内实验压差和流速，并测试得到驱替T2谱曲线；然后根据最佳转换系数C0、饱和水T2谱曲线和驱替T2谱曲线，进一步转换得到饱和流体累计曲线和束缚流体累计曲线，进而评价流体的可动性。本申请通过将CT孔径分布、核磁共振T2谱以及误差分析等结合，可充分模拟不同作业制度对流体的动用过程，不需要依赖大量实验样品，能够实现不同作业制度下流体孔喉半径动用分布的定量化表征，避免人工主观评价，为油气藏合理开发提供依据。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种油气藏储层流体的可动性的评价方法的流程图；
图2为本申请实施例提供的一种玻璃注射器的结构示意图；
图3为本申请实施例提供的述CT孔径分布累计曲线和不同的转换系数Ci转换后的饱和水T2谱累计曲线示意图；
图4为本申请实施例提供的不同转换系数Ci与对应的误差值δi的关系曲线示意图；
图5为本申请实施例提供的驱替前后流体分布曲线以及累计曲线示意图；
图6为本申请实施例提供的不同孔喉半径内流体分布及动用程度的示意图。
通过上述附图，已示出本申请明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本申请构思的范围，而是通过参考特定实施例为本领域技术人员说明本申请的概念。
具体实施方式
除非另有定义，本申请实施例所用的所有技术术语均具有与本领域技术人员通常理解的相同的含义。在对申请实施方式作进一步地详细描述之前，对理解本申请实施例一些术语进行说明。
压汞实验法是根据中值压力与渗透率的关系得到渗透率下限值，然后根据岩心化验得到孔渗关系，有此求出对应的孔隙度下限。储层物性参数统计法是以岩心所测的孔隙度、渗透率为基础，以渗透率或孔隙度累计数据5％丢失为界限。束缚水饱和度法是在孔隙度与束缚水饱和度相关系数大于70％时，选取束缚水饱和度80％，对应的孔隙度下限值即为下限值。测井资料法是在试油气资料指导下，得到其对应储层的声波时差等，依据测井解释，得到解释孔隙度，再依据孔渗关系求出渗透率下限。产能模拟实验法是通过取产层岩心，模拟不同物性条件下产能的变化，以此确定储层物性下限。
然而，在流体可动性及可动孔径下限评价过程中，这些方法存在不合理之处：静态法依赖于大量的实验岩心的物性数据，对于取芯困难的油气储层操作性不强；通过岩心物性数据统计分析下限值，选取具有一定的人工主观性，不同人员的评价结果一致性不强；评价结果为宏观物性参数，而储量动用主要受不同孔径内的储量分布控制，目前评价方法不能反映孔径内的动用差异；分析结果为孔隙度下限或渗透率下限，属于静态参数，而储量的动用下限受不同作业制度的影响，是动态的评价结果。
基于此，本申请实施例提供了一种油气藏储层流体的可动性的评价方法，该评价方法包括：
步骤101、获取油气藏储层的岩心；利用岩心CT扫描仪测试所述岩心，得到所述岩心的CT孔径分布曲线；
步骤102、采用模拟地层水浸泡所述岩心，获取饱和岩心；利用核磁设备测试所述饱和岩心，得到所述岩心的饱和水T2谱曲线；
步骤103、根据所述CT孔径分布曲线和所述饱和水T2谱曲线，得到CT孔径分布累计曲线和饱和水T2谱累计曲线；采用不同的转换系数Ci转换所述饱和水T2谱累计曲线；
步骤104、根据所述CT孔径分布累计曲线和不同的转换系数Ci转换后的所述饱和水T2谱累计曲线，采用误差分析法计算最佳转换系数Co；
步骤105、根据所述单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心；采用所述核磁设备测试所述驱替岩心，得到驱替T2谱曲线；
步骤106、根据所述最佳转换系数C0、所述饱和水T2谱曲线和所述驱替T2谱曲线，得到饱和流体累计曲线和束缚流体累计曲线；
步骤107、根据所述饱和流体累计曲线和所述束缚流体累计曲线，确定所述油气藏储层的可动孔径下限。
本申请提供的油气藏储层流体的可动性的评价方法，首先采用CT孔径分布转换饱和水T2谱曲线，并结合误差分析法求取最佳转换系数；利用相似转换原则获取室内实验压差和流速，并测试得到驱替T2谱曲线；然后根据最佳转换系数C0、饱和水T2谱曲线和驱替T2谱曲线，进一步转换得到饱和流体累计曲线和束缚流体累计曲线，进而评价流体的可动性。本申请通过将CT孔径分布、核磁共振T2谱以及误差分析等结合，可充分模拟不同作业制度对流体的动用过程，不需要依赖大量实验样品，能够实现不同作业制度下流体孔喉半径动用分布的定量化表征，避免人工主观评价，为油气藏合理开发提供依据。
为使本申请的技术方案和优点更加清楚，下面将结合附图对本申请具体实施方式作进一步地详细描述。
以西南油气田的SY1-1井为例，在对其进行开发前，对油气藏储层流体的可动性进行评价，评价步骤如下：
步骤1、沿储层渗流方向钻取柱塞岩心，柱塞岩心直径1英寸，长度约为2英寸；钻取岩心标记5-2，记录取芯井号ST3、层位栖霞组和深度7651.12m；将钻取的岩心切平，保持岩样端面与柱面平整；采用乙醇等试剂对岩心进行洗油洗盐处理；然后烘干岩心，称取岩心干重66.55g、长度2.515cm、直径5.094cm，标记岩心驱替方向，并测试岩心孔隙度5.80％、渗透率0.434mD；将制备完成的岩心样品放入干燥皿中储存备用。
步骤2、利用岩心CT扫描仪测试岩心，获取岩心三维不同切片灰度图像；利用数字岩心软件对灰度图像进行图像滤波处理、孔隙三维重构、孔隙分割等处理；对处理完成后的孔隙进行孔喉半径统计，获取不同孔喉半径对应的孔隙体积分布，即得到CT孔径分布曲线。
步骤3、将岩心放入真空-饱和装置，真空-饱和装置抽真空6h，保持装置内真空压力小于-0.1MPa，将脱气后的模拟地层水放入真空室内侵没岩心，并继续抽真空2h；再采用加压装置对饱和岩心加压6h，饱和压力12MPa；称取岩心的湿重为68.02g，将饱和岩心浸泡在模拟地层水中备用。
步骤4、利用玻璃材料制作活塞注射器，使注射器的内径与岩心直径一致，注射器内部长度10cm；将注射器完全侵入饱和流体，将饱和岩心放入注射器内，推入注射器活塞，排出多余流体，直至注射器内壁、活塞与岩心紧密贴合。活塞注射器可如图2所示。其中，21为饱和岩心，22为注射器活塞，23为注射器排液出口，24为饱和模拟地层水，25为烧杯。
步骤5、对核磁设备进行调试，设置核磁共振测试的回波时间、采样等待时间和叠加次数等各项参数；将装有岩心的玻璃注射器放入核磁共振测试线圈，保持岩心位于线圈中部位置；对岩心施加电磁射频信号，同一样品可重复采集多次(例如，3次)信号，得到岩心的饱和水T2谱曲线。
步骤6、将CT孔径分布曲线转化为CT孔径分布累计曲线；将T2谱的信号强度换算为百分比，并将饱和水T2谱曲线转化为饱和水T2谱累计曲线；选取不同的转换系数Ci，如0.01、0.1、1、10等，利用转换系数Ci乘以T2谱累计曲线的弛豫时间，得到不同的转换系数Ci转换后的饱和水T2谱累计曲线，如图3所示；
步骤7、对于不同Ci转换后的累计曲线，利用差值法求取T2谱累计百分数5％、10％、15％....90％、95％对应的弛豫时间Tij；利用差值法求取T2谱百分比分布曲线在2.5％～7.5％、7.5％～12.5％....87.5％～92.5％、92.5％～97.5％区间内分布百分比ω(Tij)。将CT孔径分布百分比曲线转化为累计曲线，利用差值法求取累计百分数5％、10％、15％....90％、95％对应的孔喉半径Rj。
步骤8、利用曲线拟合误差分析式，计算CT孔喉半径与不同转换系数Ci转换后的核磁共振T2谱之间的拟合误差δi，曲线拟合误差分析式为：
其中，N表示对T2谱百分比分布曲线所划分的区间的总个数，Rj表示对应不同区间的孔喉半径的取值，i表示所采用的不同转换系数的次数，i，j为正整数。
步骤9、绘制拟合误差δi与转换系数Ci的关系曲线，并拟合拟合误差δi与转换系数Ci的关系式，拟合符合率R2>0.99，关系式为：δi＝f(Ci)，如图4所示；对关系式求导，在定义域内，求取δi′＝0时，对应的C＝1.1736值，确定该值为最佳转换系数Co。
步骤10、根据SY1-1井生产动态等资料，气井地层压力PR＝95.32MPa、井底压力Pwf＝92.16MPa、泄流半径re＝67.11m；按照室内实验压力降梯度与矿场地层压力降梯度等效，根据计算室内实验岩心等效压差ΔP为2.18kPa。
步骤11、利用氮气沿岩心标记方向驱替岩心，驱替条件按照等效压差2.18kPa进行；驱替过程中岩心称重，若岩心3次称重质量不变，确定驱替后岩心重量为67.46g；通过称重计算流体动用程度为38.09％，如表1所示；岩心放入活塞注射器，排出多余的液体，直至注射器内壁、活塞与岩心紧密贴合。
表1岩心的流体动用程度
其中，流体动用程度的计算公式为：De＝100％-(Gw-Go)/(Gwi-Go)*100％。
步骤12、对核磁设备进行调试，设置核磁共振测试的回波时间、采样等待时间和叠加次数各项参数与步骤5相同；将装有岩心的玻璃注射器放入核磁共振测试线圈，保持岩心位于线圈中部位置；对岩心施加电磁射频信号，测试岩心T2谱，同一样品可重复采集多次(例如，3次)信号，得到驱替岩心的驱替T2谱曲线。
步骤13、利用最佳转换系数Co＝1.1736转换饱和水T2谱曲线和驱替T2谱曲线，得到饱和流体孔径分布曲线和束缚流体孔径分布曲线；将饱和流体孔径分布曲线和束缚流体孔径分布曲线转换为饱和流体累计曲线和束缚流体累计曲线；对比饱和流体累计曲线和束缚流体累计曲线，将饱和流体累计饱和度大于束缚流体累计饱和度5％作为可动界限，对应的孔喉半径作为可动孔径下限，可动孔径下限为21.99μm。如图5所示。
步骤14、将饱和流体孔径分布曲线和束缚流体孔径分布曲线按照不同孔径分布区间进行统计，分别计算不同孔喉半径范围内初始流体饱和度(饱和岩心的相应孔径的含水量)和束缚流体饱和度(驱替后岩心的含水量)，根据初始流体饱和度和束缚流体饱和度，计算不同孔喉半径范围内的流体动用程度，结果如图6所示和表2所示。
表2岩心在不同孔径范围内的流体动用程度
可见，本申请提供的油气藏储层流体的可动性的评价方法，至少具有以下
有益效果：
(1)利用CT孔径分布转换核磁共振T2谱，相比利用压汞转换T2谱，转换结果更加可靠；
(2)采用误差分析法求取不同转换系数的误差量，并采用数学分析方法求取最佳转换系数，转换精度较高；
(3)利用相似转换原则，将气井生产压差和产气速度转换为室内实验压差和流速，可对不同生产方式下储量的可动用程度进行评价；
(4)对不同孔径区间内的流体饱和度进行统计，获取不同孔径内的流体可动性，实现不同生产方式下动用下限的确定。
上述实施例中所使用的实验方法如无特殊说明，均为常规方法。
上述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
以上所述仅是为了便于本领域的技术人员理解本申请的技术方案，并不用以限制本申请。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (8)
1.一种油气藏储层流体的可动性的评价方法，其特征在于，包括以下步骤：
获取油气藏储层的岩心；利用岩心CT扫描仪测试所述岩心，得到所述岩心的CT孔径分布曲线；
采用模拟地层水浸泡所述岩心，获取饱和岩心；利用核磁设备测试所述饱和岩心，得到所述岩心的饱和水T2谱曲线；
根据所述CT孔径分布曲线和所述饱和水T2谱曲线，得到CT孔径分布累计曲线和饱和水T2谱累计曲线；采用不同的转换系数Ci转换所述饱和水T2谱累计曲线；
根据所述CT孔径分布累计曲线和不同的转换系数Ci转换后的所述饱和水T2谱累计曲线，采用误差分析法计算最佳转换系数Co；
根据单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心；采用所述核磁设备测试所述驱替岩心，得到驱替T2谱曲线；
利用所述最佳转换系数Co转换所述饱和水T2谱曲线和所述驱替T2谱曲线，得到饱和流体孔径分布曲线和束缚流体孔径分布曲线；将所述饱和流体孔径分布曲线和所述束缚流体孔径分布曲线转换为饱和流体累计曲线和束缚流体累计曲线；
对比所述饱和流体累计曲线和所述束缚流体累计曲线，将饱和流体累计饱和度大于束缚流体累计饱和度预设值作为可动界限，对应的孔喉半径作为可动孔径下限。
2.根据权利要求1所述的评价方法，其特征在于，所述利用岩心CT扫描仪测试所述岩心，得到所述岩心的CT孔径分布曲线，包括：
利用所述岩心CT扫描仪测试所述岩心，获取所述岩心的三维切片灰度图像；
利用数字岩心软件对所述灰度图像进行图像滤波处理、孔隙三维重构处理和孔隙分割处理，获取所述岩心的CT孔径分布曲线。
3.根据权利要求1所述的评价方法，其特征在于，所述利用核磁设备测试所述饱和岩心，得到所述岩心的饱和水T2谱曲线，包括：
将所述饱和岩心放入玻璃注射器内，使所述饱和岩心与注射器内壁、注射器活塞贴合；
设置所述核磁设备的实验参数，将装有所述饱和岩心的所述玻璃注射器放入核磁共振测试线圈中，保持所述饱和岩心位于所述测试线圈的中部位置；
对所述饱和岩心施加电磁射频信号，得到所述岩心的饱和水T2谱曲线。
4.根据权利要求1所述的评价方法，其特征在于，所述根据所述CT孔径分布累计曲线和不同的转换系数Ci转换后的所述饱和水T2谱累计曲线，计算最佳转换系数Co，包括：
对于Ci转换后的饱和水T2谱累计曲线，分别利用差值法求取T2谱的不同累计百分数对应的弛豫时间Tij；利用差值法求取T2谱百分比分布曲线在不同区间内的分布百分比ω(Tij)；对于所述CT孔径分布累计曲线，利用差值法分别求取不同累计百分数对应的孔喉半径Rj；
利用曲线拟合误差分析式，计算CT孔喉半径与不同Ci转换后的核磁共振T2谱之间的拟合误差δi，所述曲线拟合误差分析式为：
其中，N表示对T2谱百分比分布曲线所划分的区间的总个数，Rj表示对应不同区间的孔喉半径的取值，i表示所采用的不同转换系数的次数，i，j为正整数；
绘制所述拟合误差δi与所述转换系数Ci的关系曲线，并拟合所述拟合误差δi与所述转换系数Ci的关系式，关系式为：δi＝f(Ci)；
求取δi′＝0时，对应的C值，该值为所述最佳转换系数Co。
5.根据权利要求1所述的评价方法，其特征在于，所述根据单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心，包括：
根据单井的生产资料，确定气井地层压力PR、井底压力Pwf、泄流半径re、井筒半径rw、油气层有效厚度h、产量Q，计算室内实验岩心等效流量QC，计算式为：
其中，D为所述岩心的直径；
根据岩心等效流量QC，利用氮气顺岩心标记方向驱替所述饱和岩心，得到所述驱替岩心。
6.根据权利要求1所述的评价方法，其特征在于，所述根据单井的生产资料，基于相似转换原则确定实验驱替条件，并在所述驱替条件下驱替所述饱和岩心，获取驱替岩心，包括：
根据单井的生产资料，确定气井地层压力PR、井底压力Pwf、泄流半径re，计算室内实验岩心等效压差ΔP，计算式为：
其中L为所述岩心的长度；
根据室内实验岩心等效压差ΔP，利用氮气顺岩心标记方向驱替所述饱和岩心，得到所述驱替岩心。
7.根据权利要求1所述的评价方法，其特征在于，所述预设值为5％。
8.根据权利要求1所述的评价方法，还包括：
根据所述饱和流体孔径分布曲线和所述束缚流体孔径分布曲线，计算不同孔喉半径范围内初始流体饱和度和束缚流体饱和度；
根据所述初始流体饱和度和所述束缚流体饱和度，计算不同孔喉半径范围内的流体动用程度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010681539.3A CN113945497B (zh) | 2020-07-15 | 2020-07-15 | 油气藏储层流体的可动性的评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010681539.3A CN113945497B (zh) | 2020-07-15 | 2020-07-15 | 油气藏储层流体的可动性的评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113945497A CN113945497A (zh) | 2022-01-18 |
| CN113945497B true CN113945497B (zh) | 2023-09-26 |
Family
ID=79325827
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010681539.3A Active CN113945497B (zh) | 2020-07-15 | 2020-07-15 | 油气藏储层流体的可动性的评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113945497B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115127970B (zh) * | 2022-08-30 | 2022-11-15 | 西南石油大学 | 一种高温高压储层等效井下渗流特征及产能实验模拟方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760082A (zh) * | 2014-01-08 | 2014-04-30 | 中国石油天然气股份有限公司 | 页岩层系致密储层原油有效可动空间确定方法和装置 |
| CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
| CN104990851A (zh) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | 一种新的页岩敏感性实验研究方法 |
| CN111157424A (zh) * | 2020-01-07 | 2020-05-15 | 中南大学 | 一种岩石材料孔径分布测定方法 |
-
2020
- 2020-07-15 CN CN202010681539.3A patent/CN113945497B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760082A (zh) * | 2014-01-08 | 2014-04-30 | 中国石油天然气股份有限公司 | 页岩层系致密储层原油有效可动空间确定方法和装置 |
| CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
| CN104990851A (zh) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | 一种新的页岩敏感性实验研究方法 |
| CN111157424A (zh) * | 2020-01-07 | 2020-05-15 | 中南大学 | 一种岩石材料孔径分布测定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 致密砂岩气层毛细管自吸水锁损害及控制因素;王俊杰 等;断块油气田;第26卷(第5期);626-631 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113945497A (zh) | 2022-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107727679B (zh) | 一种表征深层碳酸盐岩岩石物理学特征方法 | |
| Walls | Tight gas sands-permeability, pore structure, and clay | |
| CN110231272B (zh) | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 | |
| RU2331057C2 (ru) | Способ и устройство для оценки физических параметров подземного месторождения полезных ископаемых на основе исследования отобранных из этого месторождения фрагментов горной породы | |
| CN110927035A (zh) | 一种低渗致密砂岩束缚水饱和度计算方法 | |
| US7092822B2 (en) | Method of evaluating the capillary pressure curve of an underground deposit rocks based on rock cuttings measurements | |
| CN110346258A (zh) | 一种致密岩石油相相对渗透率测定的方法 | |
| CN106501144A (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| CN113740515B (zh) | 综合表征深层海相碳酸盐岩油藏赋存下限临界条件的方法 | |
| CN108827999B (zh) | 低孔渗砂岩储层可动油比例以及可动油资源量的评价方法 | |
| CN109682850B (zh) | 一种在线渗吸实验用核磁共振测试装置及实验方法 | |
| CN106777515B (zh) | 一种基于岩心实验资料分析致密气井产能的方法 | |
| CN113945497B (zh) | 油气藏储层流体的可动性的评价方法 | |
| CN114370269B (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
| CN115266800A (zh) | 一种基于二维核磁共振的凝析油饱和度测试方法 | |
| CN116027453A (zh) | 一种水合物混合层饱和度定量评价方法及装置 | |
| CN112487620A (zh) | 一种页岩油可动资源量评价模型、评价方法、应用 | |
| CN111239372A (zh) | 一种基于覆压渗流实验的碳酸盐岩孔隙结构分类的方法 | |
| US11852576B2 (en) | Method for determining the pore size distribution in a reservoir | |
| Jacobsen | New oedometer and new triaxial apparatus for firm soils | |
| CN113138205B (zh) | 确定多孔介质内部气水渗吸情况的方法及系统 | |
| CN112012727B (zh) | 获得气相有效渗透率的方法及储层产能的预测方法 | |
| CN111693427B (zh) | 油气藏流体可动性的分析方法 | |
| CN111364982B (zh) | 针对于致密砂岩含油饱和度的取值方法 | |
| CN108717036B (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 |