CN106226215A - 一种碳酸盐岩储层孔洞定量评价的方法 - Google Patents
一种碳酸盐岩储层孔洞定量评价的方法 Download PDFInfo
- Publication number
- CN106226215A CN106226215A CN201610121749.0A CN201610121749A CN106226215A CN 106226215 A CN106226215 A CN 106226215A CN 201610121749 A CN201610121749 A CN 201610121749A CN 106226215 A CN106226215 A CN 106226215A
- Authority
- CN
- China
- Prior art keywords
- rock core
- hole
- corresponding relation
- radius
- reservoir
- 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 56
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 22
- 238000011158 quantitative evaluation Methods 0.000 title 1
- 239000011435 rock Substances 0.000 claims abstract description 172
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 30
- 239000011148 porous material Substances 0.000 claims abstract description 27
- 239000011800 void material Substances 0.000 claims description 37
- 230000001186 cumulative effect Effects 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 8
- 238000002591 computed tomography Methods 0.000 claims description 6
- 238000002791 soaking Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 abstract description 8
- 238000011161 development Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000003860 storage Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 241000272525 Anas platyrhynchos Species 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000002247 constant time method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- DHNCFAWJNPJGHS-UHFFFAOYSA-J [C+4].[O-]C([O-])=O.[O-]C([O-])=O Chemical compound [C+4].[O-]C([O-])=O.[O-]C([O-])=O DHNCFAWJNPJGHS-UHFFFAOYSA-J 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 235000002639 sodium chloride Nutrition 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000006467 substitution reaction 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
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- High Energy & Nuclear Physics (AREA)
- Dispersion Chemistry (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和核磁共振技术对岩心进行实验，建立所述岩心的孔洞半径与横向弛豫时间之间的对应关系；获取储层的横向弛豫时间与孔隙度的对应关系；根据所述岩心的孔洞半径与横向弛豫时间之间的对应关系，由所述储层的横向弛豫时间与孔隙度的对应关系，得到储层的孔洞半径区间与孔洞体积百分比的对应关系。本发明通过采用CT和核磁共振技术对岩心进行实验，可以最终得到储层的孔洞半径区间与孔洞体积百分比的对应关系，进而直观了解到储层的储集性，满足了油气田勘探开发的需求。
Description
技术领域
本发明涉及油气田勘探开发技术领域，特别涉及一种碳酸盐岩储层孔洞定量评价的方法。
背景技术
储层是具有连通孔洞、允许油气在其中存储和渗透的岩层，储层的孔洞是影响油气分布和流动的重要因素，因此对储层的孔洞进行评价十分重要。碳酸盐岩储层在世界油气分布中占有重要地位，其油气储量约占全世界油气总储量的50％，油气产量达全世界油气总产量的60％以上。但碳酸盐岩储层的孔洞类型多、次生变化大、结构复杂，对碳酸盐岩储层的孔洞进行评价，特别是定量评价(采用数学的方法，收集和处理数据资料，对评价对象做出定量结果的价值判断)存在巨大困难。
目前研究碳酸盐岩储层的孔洞结构的方法包括：先对若干岩心进行核磁共振实验，得到若干岩心的核磁共振T2(横向弛豫时间)分布；再对若干岩心进行压泵实验，得到若干岩心的毛管压力曲线；根据若干岩心的核磁共振T2分布和若干岩心的毛管压力曲线，建立核磁共振T2分布与毛管压力曲线之间的对应关系；获取油田现场的T2分布；根据该对应关系，将油田现场的T2分布转换为油田现场的毛管压力曲线；由油田现场的毛管压力曲线得到油田现场的喉道半径。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题：
喉道是岩石中连通各个孔洞的狭窄通道，喉道半径侧重于反映各个孔洞的连通性，与储层(尤其是碳酸盐岩孔洞型储层)的储集性关联性不强，由于油气田勘探开发主要基于储层的储集性，因此目前研究碳酸盐岩储层的孔结构的方法无法全面满足油气田勘探开发的需求。
发明内容
为了解决现有技术无法满足油气田勘探开发的需求的问题，本发明实施例提供了一种碳酸盐岩储层孔洞定量评价的方法。所述技术方案如下：
本发明实施例提供了一种碳酸盐岩储层孔洞定量评价的方法，所述方法包括：
采用计算机断层成像技术CT和核磁共振技术对岩心进行实验，建立所述岩心的孔洞半径与横向弛豫时间之间的对应关系；
获取储层的横向弛豫时间与孔隙度的对应关系；
根据所述岩心的孔洞半径与横向弛豫时间之间的对应关系，由所述储层的横向弛豫时间与孔隙度的对应关系，得到储层的孔洞半径区间与孔洞体积百分比的对应关系。
具体地，所述采用计算机断层成像技术CT和核磁共振技术对岩心进行实验，建立所述岩心的孔洞半径与横向弛豫时间之间的对应关系，包括：
采用CT对所述岩心进行实验，得到所述岩心的所有孔洞的体积和对应的半径；
根据所述岩心的所有孔洞的体积和对应的半径，计算所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系；
采用核磁共振技术对所述岩心进行实验，得到所述岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系；
根据孔洞体积百分比的累加值与归一化的孔隙度的累加值之间的对应关系，得到所述岩心的孔洞半径与横向弛豫时间之间的对应关系。
可选地，所述根据所述岩心的所有孔洞半径，计算所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系，包括：
由设定的孔洞分级标准得到多个孔洞半径区间；
按照所述孔洞半径区间，计算各个所述孔洞半径区间的孔洞总体积；
计算各个所述孔洞半径区间的孔洞总体积与所述岩心的所有孔洞总体积的比值，得到各个所述孔洞半径区间的孔洞体积百分比；
将半径不超过所述孔洞半径区间的最大半径的所有所述孔洞半径区间的孔洞体积百分比进行累加，得到所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系。
可选地，所述采用核磁共振技术对所述岩心进行实验，得到所述岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系，包括：
将核磁共振设备参数按如下数值进行设置：
温度为35℃，回波间隔为0.2ms，等待时间为20s，扫描次数为512次，回波个数为16348个。
可选地，所述岩心的孔洞半径区间r与横向弛豫时间T2之间的对应关系为：
r＝0.0002×T2 1.4379，T2≤195.74ms；
r＝0.0391×T2 0.4455，T2＞195.74ms。
优选地，所述岩心的数量为三个。
可选地，所述方法还包括：
获取岩心原样；
对所述岩心原样进行加工，制成所述岩心。
优选地，所述岩心呈柱塞样，直径为2.5cm，长度为3-5cm。
优选地，所述获取岩心原样，包括：
选出具有半径大于5mm的孔洞的岩心原样。
优选地，所述方法还包括：
采用热缩膜包裹所述岩心的圆周面，滤纸包裹所述岩心的两端；
采用100000ppmNaCl溶液浸泡所述岩心5-6天。
本发明实施例提供的技术方案带来的有益效果是：
通过采用CT和核磁共振技术对岩心进行实验，建立岩心的孔洞半径与横向弛豫时间之间的对应关系，进而由储层的横向弛豫时间与孔隙度的对应关系得到储层的孔洞半径区间与孔洞体积百分比的对应关系，由于通过储层的孔洞半径区间与孔洞体积百分比的对应关系可以直观了解到储层的储集性，因此本发明实施例提供的碳酸盐岩储层孔洞定量评价的方法可以满足油气田勘探开发的需求。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的一种碳酸盐岩储层孔洞定量评价的方法的流程图；
图2是本发明实施例二提供的一种碳酸盐岩储层孔洞定量评价的方法的流程图；
图3是本发明实施例二提供的岩心的孔洞半径与孔洞体积百分比、孔洞体积百分比的累加值的对应关系的示意图；
图4是本发明实施例二提供的岩心的横向弛豫时间与归一化的孔隙度、孔隙度的累加值的对应关系的示意图；
图5是本发明实施例二提供的岩心的孔洞半径与横向弛豫时间之间的对应关系的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一
本发明实施例提供了一种碳酸盐岩储层孔洞定量评价的方法，参见图1，该方法包括：
步骤101：采用计算机断层成像技术(Computed Tomography，简称CT)和核磁共振技术对岩心进行实验，建立岩心的孔洞半径与横向弛豫时间(T2)之间的对应关系。
步骤102：获取储层的横向弛豫时间与孔隙度的对应关系。
步骤103：根据岩心的孔洞半径与横向弛豫时间之间的对应关系，由储层的横向弛豫时间与孔隙度的对应关系，得到储层的孔洞半径区间与孔洞体积百分比的对应关系。
本发明实施例通过采用CT和核磁共振技术对岩心进行实验，建立岩心的孔洞半径与横向弛豫时间之间的对应关系，进而由储层的横向弛豫时间与孔隙度的对应关系得到储层的孔洞半径区间与孔洞体积百分比的对应关系，由于通过储层的孔洞半径区间与孔洞体积百分比的对应关系可以直观了解到储层的储集性，因此本发明实施例提供的碳酸盐岩储层孔洞定量评价的方法可以满足油气田勘探开发的需求。
实施例二
本发明实施例提供了一种碳酸盐岩储层孔洞定量评价的方法，参见图2，该方法包括：
步骤201：采用CT对岩心进行实验，得到岩心的所有孔洞的体积和对应的半径。
在实际应用中，岩心是根据地质勘查工作或工程的需要，使用环状岩心钻头及其他取心工具，从孔内取出的圆柱状岩石样品。在本实施例中，岩心是来自于(碳酸盐岩)储层的样品。
在具体实现中，采用CT实验设备对岩心进行实验，CT实验设备可以列出岩心的所有孔洞体积和半径，即得到岩心的所有孔洞体积和半径。
可选地，岩心的数量可以为三个。
可以理解地，采用三个岩心进行实验，后续处理时(详见步骤202)可以取三个岩心的实验结果的平均值，提高实验结果的准确性，同时也不会造成较大的开销。
可以理解地，选取的岩心均具有孔洞结构，以实现对孔洞的测量。
步骤202：根据岩心的所有孔洞的体积和对应的半径，计算岩心的孔洞半径与孔洞体积百分比的累加值的对应关系。
在本实施例中，孔洞体积百分比为设定半径的所有孔洞与岩心的所有孔洞的体积比，孔洞体积百分比的累加值为不大于设定半径的所有孔洞与岩心的所有孔洞的体积比。设定半径为孔洞半径区间的最大半径。
在具体实现中，可根据岩心的孔洞体积计算属于同一个孔洞半径区间的所有孔洞总体积、以及岩心的所有孔洞的总体积，再计算出孔洞半径区间的所有孔洞与岩心的所有孔洞的体积比(实现归一化处理)，即孔洞体积百分比孔孔洞体积百分比，最后将半径不超过孔洞半径区间的最大半径的所有孔洞半径区间的体积百分比进行累加，得到孔洞体积百分比的累加值，进而由孔洞的体积与半径之间的对应关系得到孔洞半径区间与孔洞体积百分比的累加值的对应关系。具体地，可以利用Microsoft Office Excel软件对岩心的所有孔洞的体积和对应的半径进行统计和运算，得到岩心的孔洞半径区间与孔洞体积百分比的累加值的对应关系。
图3中的曲线为步骤202得到的岩心的孔洞半径区间与孔洞体积百分比、孔洞体积百分比的累加值的对应关系的示意图。图3中的横坐标表示孔洞半径，纵坐标表示孔洞体积百分比。横坐标按照孔洞半径区间进行了划分，划分后的方格内高度不等的柱形代表相应孔洞半径区间的孔洞体积百分比，曲线上的点表示小于相应孔洞半径区间的最大半径的所有孔洞半径区间的孔洞体积百分比的累加值。
在本实施例中，该步骤202可以包括：
由设定的孔洞分级标准得到多个孔洞半径区间；
按照孔洞半径区间，计算各个孔洞半径区间的孔洞总体积；
计算各个孔洞半径区间的孔洞总体积与岩心的所有孔洞总体积的比值，得到各个孔洞半径区间的孔洞体积百分比；
将半径不超过孔洞半径区间的最大半径的所有孔洞半径区间的孔洞体积百分比进行累加，得到岩心的孔洞半径与孔洞体积百分比的累加值的对应关系。
其中，设定的孔洞分级标准可以参考《综合录井岩心整理和描述规范中油测井技术服务有限责任公司企业标准ST/CNLC SOP 3005-2007》、《碳酸盐岩孔隙分类方法综述》等资料得到，如下表一所示：
表一
| 序号 | 级别 | 半径区间(mm) |
| 1 | 隐孔隙 | <0.005 |
| 2 | 微孔隙 | 0.05～0.005 |
| 3 | 细孔隙 | 0.125～0.05 |
| 4 | 中孔隙 | 0.25～0.125 |
| 5 | 粗孔隙 | 0.5～0.25 |
| 6 | 巨孔隙 | 1～0.5 |
| 7 | 小洞 | 2.5～1 |
| 8 | 中洞 | 5～2.5 |
| 9 | 大洞 | 250～5 |
| 10 | 洞穴 | ≥250 |
可以理解地，按照设定的孔洞分级标准进行统计，可以使统计结果更好地与实际应用相匹配。
步骤203：采用核磁共振技术对岩心进行实验，得到岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系。该步骤203可以在步骤201之后执行，也可以在步骤201之前执行，也可以同步执行。
在具体实现中，采用核磁共振实验设备对岩心进行实验，得到岩心的横向弛豫时间与孔隙度的对应关系，如依据SY/T 6490-2007《岩样核磁共振参数实验室测量规范》，采用MARAN DRX2核磁共振测量仪对岩心进行横向弛豫时间的测量。由于孔洞体积百分比的累加值的范围为0～100％，孔隙度的累加值的范围为0～孔洞体积/岩心体积(孔洞体积＜岩心体积，孔洞体积/岩心体积＜100％)，为了建立孔洞体积百分比的累加值和孔隙度的累加值之间的对应关系，会将孔隙度的累加值进行归一化，得到岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系。
横向弛豫时间和孔隙大小及形态有关，孔隙几何形态及孔隙表面弛豫强度一定情况下，孔隙越小，表面弛豫衰减越快，横向弛豫时间越小；孔隙越大，表面弛豫衰减越慢，横向弛豫时间越大，即横向弛豫时间与孔隙大小(即孔洞体积)之间具有对应关系。由于岩心体积是不变的，因此横向弛豫时间与孔隙度(孔洞体积与岩心体积之比，用百分比表示)之间也存在对应关系。可以按照横向弛豫时间的分布，将对应的横向弛豫时间不超过某个横向弛豫时间的孔隙度相加，得到孔隙度的累积值，进而得到各个横向弛豫时间与孔隙度的累积值的对应关系，最终得到岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系。
图4中的虚线为步骤203得到的岩心的横向弛豫时间与归一化的孔隙度、孔隙度的累加值的对应关系的示意图。图4中实线上的点表示横向弛豫时间对应的归一化的孔隙度，将不超过横向弛豫时间的孔隙度相加，即可得到虚线，虚线上的点表示横向弛豫时间对应的归一化的孔隙度的累加值。
可选地，该步骤203可以包括：
将核磁共振设备参数按如下数值进行设置：
温度为35℃，回波间隔为0.2ms，等待时间为20s，扫描次数为512次，回波个数为16348个。
需要说明的是，由于核磁共振设备的实验范围有限，因此需要设定参数对实验范围进行选择。考虑到采用CT对岩心进行实验时，受到CT分辨率的限制，对于直径为2.5cm的柱塞样岩心，无法识别微孔隙和隐孔隙，同时对于碳酸盐岩孔洞型储层，微孔隙和隐孔隙所占体积相对较小，其影响可以忽略。另外，进行实验的岩心体积有限(如直径为2.5cm、长度为3～5cm的柱塞样岩心)，上面不会有洞穴，因此按照上述参数进行设置，能够兼顾测量细孔隙、中孔隙、粗孔隙、巨孔隙、小洞、中洞、大洞的横向弛豫时间，与CT分析结果进行对比。
步骤204：根据孔洞体积百分比的累加值与归一化的孔隙度的累加值之间的对应关系，得到岩心的孔洞半径与横向弛豫时间之间的对应关系。
在具体实现中，孔洞体积百分比的累加值为不大于设定半径的所有孔洞体积与岩心的所有孔洞体积之比，即多个半径的所有孔洞体积与岩心的所有孔洞体积之比相加的结果。孔隙度的累加值为孔隙度(孔洞体积与岩心体积之比，用百分比表示)相加的结果，由于岩心的所有孔洞总体积、岩心体积均为固定不变的，因此孔洞体积的百分比与孔隙度之间是一一对应的，进而孔洞体积百分比的累加值与孔隙度的累加值之间具有对应关系。利用该对应关系、以及步骤202得到的岩心的孔洞半径与孔洞体积百分比的累加值的对应关系、步骤203得到的岩心的横向弛豫时间与孔隙度的累加值的对应关系，即可得到岩心的孔洞半径与横向弛豫时间之间的对应关系。
步骤204的执行结果可以如下表二所示：
表二
| 序号 | 级别 | 半径区间(mm) | 横向弛豫时间区间(ms) |
| 1 | 隐孔隙 | <0.005 | / |
| 2 | 微孔隙 | 0.05～0.005 | <2.1 |
| 3 | 细孔隙 | 0.125～0.05 | 2.1～14.6 |
| 4 | 中孔隙 | 0.25～0.125 | 14.6～64.2 |
| 5 | 粗孔隙 | 0.5～0.25 | 64.2～227.3 |
| 6 | 巨孔隙 | 1～0.5 | 227.3～365.6 |
| 7 | 小洞 | 2.5～1 | 365.6～685.4 |
| 8 | 中洞 | 5～2.5 | 685.4～1102.5 |
| 9 | 大洞 | 250～5 | ≥1102.5 |
| 10 | 洞穴 | ≥250 | / |
图5中的曲线为步骤204得到的岩心的孔洞半径与横向弛豫时间之间的对应关系的示意图。图5中的点表示各个孔洞半径区间的最大半径对应的横向弛豫时间，拟合出通过这些点的曲线即为岩心的孔洞半径与横向弛豫时间之间的对应关系。
在实际应用中，岩心的孔洞半径与横向弛豫时间之间的对应关系可以采用函数表示。由于岩心的孔洞半径与横向弛豫时间之间的对应关系中的数据点分布呈分段指数趋势，横向弛豫时间较大的数据点趋势较陡，横向弛豫时间较小的数据点趋势较缓，为了与得到的数据尽可能匹配，采用分段函数表示岩心的孔洞半径与横向弛豫时间之间的对应关系，如岩心的孔洞半径r和岩心的T2满足如下公式：
r＝0.0002×T2 1.4379，T2≤195.74ms；
r＝0.0391×T2 0.4455，T2＞195.74ms。
实验证明，岩心的孔洞半径r和岩心的T2之间的关系采用分段函数表示的准确度较高。
可以理解地，通过步骤201-步骤204实现了采用CT和核磁共振技术对岩心进行实验，建立岩心的孔洞半径与横向弛豫时间分布区间之间的对应关系。
步骤205：获取储层的横向弛豫时间与孔隙度的对应关系。
在实际应用中，可以通过对储层采用核磁共振技术进行测量，得到储层的横向弛豫时间与孔隙度的对应关系；也可以直接采用数据库中存储的数据，得到储层的横向弛豫时间与孔隙度的对应关系。
步骤206：根据岩心的孔洞半径与横向弛豫时间之间的对应关系，由储层的横向弛豫时间与孔隙度的对应关系，得到储层的孔洞半径区间与孔洞体积百分比的对应关系。
可以理解地，由于油气都是存储在储层的孔洞中的，得到储层的孔洞半径与孔洞体积百分比的对应关系，即可直观了解到储层的储集性，进而根据储层的储集性进行油气田的勘探开发，满足了油气田勘探开发的需求。
在实际应用中，还可以进一步得到储层的平均孔洞半径、最大孔洞半径、各个孔洞半径区间的孔洞体积等参数，实现碳酸盐岩储层孔洞的定量评价。
在本实施例的一种实现方式中，该方法还可以包括：
获取岩心原样；
对岩心原样进行加工，制成岩心。
可以理解地，上述实现方式在步骤201和步骤203之前执行。
可选地，岩心可以呈柱塞样，便于实验测量。
在实际应用中，可以从储层中钻取岩心样品，制成柱塞样，磨平端面，并进行清洗。
具体地，岩心的直径可以为2.5cm，岩心的长度可以为3-5cm，以保证CT实验的分辨率(直径越大，CT分辨率越差)。例如，对于直径为2.5cm的柱塞样岩心，CT理论分辨率为11μm，实际分辨率为40-55μm；对于直径为10cm的全直径岩心，CT理论分辨率为38μm，实际分辨率为60-100μm。
优选地，获取岩心原样，可以包括：
选出具有半径大于5mm的孔洞的岩心原样。
可以理解地，本发明主要是对储层的储集性进行评价，由于储层的储集能力主要由孔洞实现，孔洞越大，储层的储集能力越强，选择具有较大孔洞半径(大于5mm)的岩心原样，可以提高实验结果的准确性。
进一步地，选出具有半径大于5mm的孔洞的岩心原样，可以包括：
选出具有半径为5mm-12.5mm的孔洞的岩心原样。
如前所述，实验测量时针对的是细孔隙、中孔隙、粗孔隙、巨孔隙、小洞、中洞、大洞，因此直接选择具有半径为5mm-12.5mm的孔洞的岩心原样，可以提高实验的准确度。
在本实施例的另一种实现方式中，岩心可以呈全直径岩心(用取心技术从油气层中取出的岩心，不经过切割和劈分)，以节省操作。
可选地，该方法还可以包括：
采用热缩膜包裹岩心的圆周面，滤纸包裹岩心的两端；
采用100000ppmNaCl溶液浸泡岩心5-6天。
需要说明的是，采用热缩膜包裹岩心的圆周面，滤纸包裹岩心的两端，可以便于CT实验时测量岩心外部的孔洞。采用100000ppmNaCl溶液浸泡岩心5-6天，可以使岩心充分饱和，模拟真实油气位于岩心内的情形，提高CT实验的准确性。
需要说明的是，本发明实施例中的数据是由四川盆地磨溪地区龙王庙组碳酸盐岩储层得到的。
本发明实施例通过采用CT和核磁共振技术对岩心进行实验，建立岩心的孔洞半径与横向弛豫时间之间的对应关系，进而由储层的横向弛豫时间与孔隙度的对应关系得到储层的孔洞半径区间与孔洞体积百分比的对应关系，由于通过储层的孔洞半径区间与孔洞体积百分比的对应关系可以直观了解到储层的储集性，因此本发明实施例提供的碳酸盐岩储层孔洞定量评价的方法可以满足油气田勘探开发的需求。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种碳酸盐岩储层孔洞定量评价的方法，其特征在于，所述方法包括：
采用计算机断层成像技术CT和核磁共振技术对岩心进行实验，建立所述岩心的孔洞半径与横向弛豫时间之间的对应关系；
获取储层的横向弛豫时间与孔隙度的对应关系；
根据所述岩心的孔洞半径与横向弛豫时间之间的对应关系，由所述储层的横向弛豫时间与孔隙度的对应关系，得到储层的孔洞半径区间与孔洞体积百分比的对应关系。
2.根据权利要求1所述的方法，其特征在于，所述采用计算机断层成像技术CT和核磁共振技术对岩心进行实验，建立所述岩心的孔洞半径与横向弛豫时间之间的对应关系，包括：
采用CT对所述岩心进行实验，得到所述岩心的所有孔洞的体积和对应的半径；
根据所述岩心的所有孔洞的体积和对应的半径，计算所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系；
采用核磁共振技术对所述岩心进行实验，得到所述岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系；
根据孔洞体积百分比的累加值与归一化的孔隙度的累加值之间的对应关系，得到所述岩心的孔洞半径与横向弛豫时间之间的对应关系。
3.根据权利要求2所述的方法，其特征在于，所述根据所述岩心的所有孔洞半径，计算所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系，包括：
由设定的孔洞分级标准得到多个孔洞半径区间；
按照所述孔洞半径区间，计算各个所述孔洞半径区间的孔洞总体积；
计算各个所述孔洞半径区间的孔洞总体积与所述岩心的所有孔洞总体积的比值，得到各个所述孔洞半径区间的孔洞体积百分比；
将半径不超过所述孔洞半径区间的最大半径的所有所述孔洞半径区间的孔洞体积百分比进行累加，得到所述岩心的孔洞半径与孔洞体积百分比的累加值的对应关系。
4.根据权利要求2所述的方法，其特征在于，所述采用核磁共振技术对所述岩心进行实验，得到所述岩心的横向弛豫时间与归一化的孔隙度的累加值的对应关系，包括：
将核磁共振设备参数按如下数值进行设置：
温度为35℃，回波间隔为0.2ms，等待时间为20s，扫描次数为512次，回波个数为16348个。
5.根据权利要求2所述的方法，其特征在于，所述岩心的孔洞半径区间r与横向弛豫时间T2之间的对应关系为：
r＝0.0002×T2 1.4379，T2≤195.74ms；
r＝0.0391×T2 0.4455，T2＞195.74ms。
6.根据权利要求1-5任一项所述的方法，其特征在于，所述岩心的数量为三个。
7.根据权利要求1-5任一项所述的方法，其特征在于，所述方法还包括：
获取岩心原样；
对所述岩心原样进行加工，制成所述岩心。
8.根据权利要求7所述的方法，其特征在于，所述岩心呈柱塞样，直径为2.5cm，长度为3-5cm。
9.根据权利要求8所述的方法，其特征在于，所述获取岩心原样，包括：
选出具有半径大于5mm的孔洞的岩心原样。
10.根据权利要求7所述的方法，其特征在于，所述方法还包括：
采用热缩膜包裹所述岩心的圆周面，滤纸包裹所述岩心的两端；
采用100000ppmNaCl溶液浸泡所述岩心5-6天。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610121749.0A CN106226215A (zh) | 2016-03-03 | 2016-03-03 | 一种碳酸盐岩储层孔洞定量评价的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610121749.0A CN106226215A (zh) | 2016-03-03 | 2016-03-03 | 一种碳酸盐岩储层孔洞定量评价的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106226215A true CN106226215A (zh) | 2016-12-14 |
Family
ID=57518938
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610121749.0A Pending CN106226215A (zh) | 2016-03-03 | 2016-03-03 | 一种碳酸盐岩储层孔洞定量评价的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106226215A (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107727679A (zh) * | 2017-11-03 | 2018-02-23 | 中国科学院地质与地球物理研究所 | 一种表征深层碳酸盐岩岩石物理学特征方法 |
| CN109374489A (zh) * | 2018-08-21 | 2019-02-22 | 中国地质大学（武汉） | 联合x-ct技术的水合物沉积物nmr弛豫信号量标定装置和方法 |
| CN110133035A (zh) * | 2019-06-12 | 2019-08-16 | 陕西延长石油(集团)有限责任公司研究院 | 核磁共振耦合恒速压汞定量表征砂岩储层孔喉结构的方法 |
| CN111352155A (zh) * | 2018-12-21 | 2020-06-30 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
| CN112305000A (zh) * | 2020-10-22 | 2021-02-02 | 西南石油大学 | 一种基于页岩有机质孔形态表征页岩储集能力的方法 |
| CN113405962A (zh) * | 2020-03-16 | 2021-09-17 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的储集空间评价方法和装置 |
| CN114280686A (zh) * | 2020-09-27 | 2022-04-05 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的岩心物性分析方法及设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
| CN102262041A (zh) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
| CN103513285A (zh) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种确定横向表面弛豫速率的方法及装置 |
| CN104463362A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种预测储层流体产量的方法及装置 |
-
2016
- 2016-03-03 CN CN201610121749.0A patent/CN106226215A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
| CN102262041A (zh) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
| CN103513285A (zh) * | 2013-09-27 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种确定横向表面弛豫速率的方法及装置 |
| CN104463362A (zh) * | 2014-12-03 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种预测储层流体产量的方法及装置 |
Non-Patent Citations (5)
| Title |
|---|
| YANBIN YAO ET AL.: "《Petrophysical characterization of coals by low-field nuclear magnetic resonance (NMR)》", 《FUEL》 * |
| 何雨丹 等: "《利用核磁共振T2分布构造毛管压力曲线的新方法》", 《吉林大学学报（地球科学版）》 * |
| 吴运强 等: "《数字岩心重建与孔隙参数分析》", 《太赫兹科学与电子信息学报》 * |
| 李宁 等: "《应用CT分析及核磁测井预测碳酸盐岩产气量》", 《石油勘探与开发》 * |
| 胡荣泽: "《粉末颗粒和孔隙的测量》", 28 February 1982, 冶金工业出版社 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107727679A (zh) * | 2017-11-03 | 2018-02-23 | 中国科学院地质与地球物理研究所 | 一种表征深层碳酸盐岩岩石物理学特征方法 |
| CN107727679B (zh) * | 2017-11-03 | 2018-08-14 | 中国科学院地质与地球物理研究所 | 一种表征深层碳酸盐岩岩石物理学特征方法 |
| CN109374489A (zh) * | 2018-08-21 | 2019-02-22 | 中国地质大学（武汉） | 联合x-ct技术的水合物沉积物nmr弛豫信号量标定装置和方法 |
| CN109374489B (zh) * | 2018-08-21 | 2020-07-07 | 中国地质大学（武汉） | 联合x-ct技术的水合物沉积物nmr弛豫信号量标定装置和方法 |
| CN111352155A (zh) * | 2018-12-21 | 2020-06-30 | 中国石油天然气股份有限公司 | 流体识别方法及装置 |
| CN110133035A (zh) * | 2019-06-12 | 2019-08-16 | 陕西延长石油(集团)有限责任公司研究院 | 核磁共振耦合恒速压汞定量表征砂岩储层孔喉结构的方法 |
| CN113405962A (zh) * | 2020-03-16 | 2021-09-17 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的储集空间评价方法和装置 |
| CN113405962B (zh) * | 2020-03-16 | 2023-02-10 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的储集空间评价方法和装置 |
| CN114280686A (zh) * | 2020-09-27 | 2022-04-05 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的岩心物性分析方法及设备 |
| CN114280686B (zh) * | 2020-09-27 | 2024-05-28 | 中国石油天然气股份有限公司 | 碳酸盐岩储层的岩心物性分析方法及设备 |
| CN112305000A (zh) * | 2020-10-22 | 2021-02-02 | 西南石油大学 | 一种基于页岩有机质孔形态表征页岩储集能力的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106226215A (zh) | 一种碳酸盐岩储层孔洞定量评价的方法 | |
| EP2965126B1 (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| CN104533400B (zh) | 一种重构测井曲线的方法 | |
| EP2765409B1 (en) | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field | |
| CN105240001B (zh) | 核磁共振测井孔隙度校正方法及装置 | |
| US20060092766A1 (en) | Method and system for predicting production of a well | |
| CN111665560B (zh) | 油气层识别方法、装置、计算机设备及可读存储介质 | |
| CN110231272A (zh) | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 | |
| CN106644736B (zh) | 一种测量储层岩石的脆性指数的方法 | |
| CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
| CN103306671A (zh) | 一种四象限储层类型识别方法及系统 | |
| CN111827968B (zh) | 一种基于核磁共振测井的储层非均质性评价方法及装置 | |
| Wang et al. | Study on the pore structure and fractal dimension of tight sandstone in coal measures | |
| CN104932027A (zh) | 基于核磁共振测井的储层分类方法 | |
| CN110593857B (zh) | 油润湿致密砂岩核磁共振测井流体识别方法及其图版 | |
| CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
| CN112082918B (zh) | 孔隙度的确定方法、装置及设备 | |
| CN108518220A (zh) | 基于综合流体识别指数识别储层流体的方法与装置 | |
| CN111271054B (zh) | 裂缝型碳酸盐岩储层压裂效果的评估分析方法及装置 | |
| CN107843531A (zh) | 低渗透低粘油油层孔隙结构的核磁共振参数表征方法 | |
| CN111220639B (zh) | 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 | |
| CN104463362A (zh) | 一种预测储层流体产量的方法及装置 | |
| CN107605468B (zh) | 一种确定核磁共振测井的核磁有效孔隙度的方法及装置 | |
| CN117765286B (zh) | 一种水泥基材料的水分相态表征方法和装置 | |
| CN112710688B (zh) | 核磁共振纵向弛豫获取方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161214 |
|
| RJ01 | Rejection of invention patent application after publication |