CN107688037B - 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 - Google Patents
一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 Download PDFInfo
- Publication number
- CN107688037B CN107688037B CN201710705113.5A CN201710705113A CN107688037B CN 107688037 B CN107688037 B CN 107688037B CN 201710705113 A CN201710705113 A CN 201710705113A CN 107688037 B CN107688037 B CN 107688037B
- Authority
- CN
- China
- Prior art keywords
- particle size
- nuclear magnetic
- distribution
- size distribution
- value
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000002245 particle Substances 0.000 claims abstract description 45
- 238000006243 chemical reaction Methods 0.000 claims abstract description 21
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 238000001035 drying Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 235000011203 Origanum Nutrition 0.000 description 1
- 241001529744 Origanum Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000007771 core particle Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003921 particle size analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/34—Purifying; Cleaning
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/44—Sample treatment involving radiation, e.g. heat
-
- 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/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Dispersion Chemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，包括以下步骤：S1、实验岩心制备；S2、利用核磁共振实验分析仪及激光粒度分析仪对多块岩心样品进行试验分析，获取每块岩心的核磁T2分布及岩石粒度分布曲线；S3、利用孔隙度分析仪测量多块岩心的孔隙度参数；S4、通过核磁共振测井装置采集和记录地层核磁T2分布数据及地层孔隙度参数，记录步长为1点/0.1m；S5、选取区域T2分布转粒度分布曲线的转换系数C值，S6、求取地层的粒度分布曲线；S7、输出计算结果。本发明的有益效果为：解决了测井资料计算岩石粒度分布曲线的问题，成本低，较常规技术手段确定粒度分布更准确、适用范围更广。
Description
技术领域
本发明涉及地质储层岩石物理分析技术领域，尤其涉及一种利用核磁测井T2分布确定井下岩石粒度曲线的方法。
背景技术
在地层岩石物理研究中，传统的粒度分布曲线获取方法是岩心试验分析及基于常规测井资料建模计算。
岩心粒度试验分析技术主要包括直接测量法、筛析法及激光粒度分析仪法，井下钻井取心不连续且取心成本昂贵，这些粒度分析都是在实验室条件下获得的，时效性差。
国内现有的基于常规测井资料的粒度参数计算主要以自然伽玛、冲洗带电阻率等测井曲线与粒度参数建模，但上述方法获取的粒度参数较少，不能真实反映地层垂向组合特征。
发明内容
本发明提供了一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，为储层岩石物理分析提供连续的井下岩石粒度分布曲线，提高地质评价精度。
为解决上述技术问题，本发明提供了一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，包括以下步骤：
S1、实验岩心制备：选取不同孔隙度、渗透率的多块岩心样品，对岩心进行洗油、烘干预处理；
S2、利用核磁共振实验分析仪及激光粒度分析仪对多块岩心样品进行试验分析，获取每块岩心的核磁T2分布及岩石粒度分布曲线；
S3、利用孔隙度分析仪测量多块岩心的孔隙度参数；
S4、通过核磁共振测井装置采集和记录地层核磁T2分布数据及地层孔隙度参数，记录步长为1点/0.1m；
S5、选取区域T2分布转粒度分布曲线的转换系数C值，根据不同类型地层分类选取，范围为0.0003～0.0006；
S6、根据公式
求取地层的粒度分布曲线rgrain,核磁测井T2分布反演岩石粒度分布曲线公式：
式中：T2为核磁测井横向弛豫时间谱分布，用ms表示；
C为转换系数，无量纲；
rgrain为岩石粒度，用um表示；
Φ为地层孔隙度，用小数表示；
S7、输出计算结果。
作为本发明的一个优选的技术方案，步骤S5所述的选取区域T2分布转粒度分布曲线的转换系数C值，根据不同类型地层分类进行选取，其确定方法，包括：
对于多块岩心，获取核磁共振T2分布和岩石粒度分布曲线数据；确定转换系数C值取值区间；按照指定步长改变C值，得到伪粒度分布曲线；计算核磁模拟伪粒度分布曲线与实验粒度分布之间的误差，若误差小于给定的截止值0.0001，则输出C值，否则修改C值继续计算，直到计算误差小于截止值0.0001；获得每块岩心的T2分布转粒度分布曲线的转换系数C值。
本申请实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
解决了测井资料计算岩石粒度分布曲线的问题，成本低，较常规技术手段确定粒度分布更准确、适用范围更广。本发明在南海西部海域莺歌海盆地、琼东南盆地计算岩石粒度分布曲线30余井次，与岩心分析粒度结果相比较误差一般不超过10％，能较好满足地质研究需要。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明工作流程框图；
图2为本发明实际工作流程详图；
图3为发明应用成果图。
具体实施方式
为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
如图1、图2所示，本实施例所述的一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，包括以下步骤：
S1、实验岩心制备：选取不同孔隙度、渗透率的多块岩心样品，对岩心进行洗油、烘干预处理；
S2、利用核磁共振实验分析仪及激光粒度分析仪对多块岩心样品进行试验分析，获取每块岩心的核磁T2分布及岩石粒度分布曲线；
S3、利用孔隙度分析仪测量多块岩心的孔隙度参数；
S4、通过核磁共振测井装置采集和记录地层核磁T2分布数据及地层孔隙度参数，记录步长为1点/0.1m；
S5、选取区域T2分布转粒度分布曲线的转换系数C值，根据不同类型地层分类选取，范围为0.0003～0.0006；
S6、根据公式
求取地层的粒度分布曲线rgrain,核磁测井T2分布反演岩石粒度分布曲线公式：
式中：T2为核磁测井横向弛豫时间谱分布，用ms表示；
C为转换系数，无量纲；
rgrain为岩石粒度，用um表示；
Φ为地层孔隙度，用小数表示；
S7、输出计算结果。
2、根据权利要求1所述的一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，其特征在于：步骤S5所述的选取区域T2分布转粒度分布曲线的转换系数C值，根据不同类型地层分类进行选取，其确定方法，包括：
对于多块岩心，获取核磁共振T2分布和岩石粒度分布曲线数据；确定转换系数C值取值区间；按照指定步长改变C值，得到伪粒度分布曲线；计算核磁模拟伪粒度分布曲线与实验粒度分布之间的误差，若误差小于给定的截止值0.0001，则输出C值，否则修改C值继续计算，直到计算误差小于截止值0.0001；获得每块岩心的T2分布转粒度分布曲线的转换系数C值。
渗透率Perm<1mD，C＝0.00025；1mD≤Perm≤30mD，C＝0.0004；Perm>30mD，C＝0.0005。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
Claims (1)
1.一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，其特征在于：包括以下步骤：
S1、实验岩心制备：选取不同孔隙度、渗透率的多块岩心样品，对岩心进行洗油、烘干预处理；
S2、利用核磁共振实验分析仪及激光粒度分析仪对多块岩心样品进行试验分析，获取每块岩心的核磁T2分布及岩石粒度分布曲线；
S3、利用孔隙度分析仪测量多块岩心的孔隙度参数；
S4、通过核磁共振测井装置采集和记录地层核磁T2分布数据及地层孔隙度参数，记录步长为1点/0.1m；
S5、选取区域T2分布转粒度分布曲线的转换系数C值，根据不同类型地层分类选取，范围为0.0003～0.0006；
其确定方法为：对于多块岩心，获取核磁共振T2分布和岩石粒度分布曲线数据；确定转换系数C值取值区间；按照指定步长改变C值，得到伪粒度分布曲线；计算核磁模拟伪粒度分布曲线与实验粒度分布之间的误差，若误差小于给定的截止值0.0001，则输出C值，否则修改C值继续计算，直到计算误差小于截止值0.0001；获得每块岩心的T2分布转粒度分布曲线的转换系数C值；
S6、根据公式
求取地层的粒度分布曲线rgrain,核磁测井T2分布反演岩石粒度分布曲线公式：
式中：T2为核磁测井横向弛豫时间谱分布，用ms表示；
C为转换系数，无量纲；
rgrain为岩石粒度，用um表示；
Φ为地层孔隙度，用小数表示；
S7、输出计算结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710705113.5A CN107688037B (zh) | 2017-08-17 | 2017-08-17 | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710705113.5A CN107688037B (zh) | 2017-08-17 | 2017-08-17 | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107688037A CN107688037A (zh) | 2018-02-13 |
| CN107688037B true CN107688037B (zh) | 2020-04-07 |
Family
ID=61153427
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710705113.5A Active CN107688037B (zh) | 2017-08-17 | 2017-08-17 | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107688037B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117505A (zh) * | 2018-07-10 | 2019-01-01 | 中国海洋石油集团有限公司 | 一种基于介电实验的孔隙结构储层冲洗带含水饱和度计算方法 |
| CN108982320B (zh) * | 2018-07-10 | 2021-03-02 | 中国海洋石油集团有限公司 | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 |
| CN111581761B (zh) * | 2019-01-30 | 2022-11-04 | 中国石油天然气股份有限公司 | 基于岩石粒度的数据处理方法、装置及存储介质 |
| CN114075965B (zh) * | 2020-08-05 | 2024-05-07 | 中国石油化工股份有限公司 | 一种河道砂岩粒度的确定方法 |
| CN113125676A (zh) * | 2021-03-25 | 2021-07-16 | 中国石油天然气股份有限公司 | 一种油气常规储层样品联合测试方法 |
| CN112986309A (zh) * | 2021-04-01 | 2021-06-18 | 中海石油(中国)有限公司 | 一种利用岩屑破碎煤样测量煤层孔隙度的方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1896876B1 (en) * | 2005-06-03 | 2013-04-17 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| CN102253069B (zh) * | 2011-04-20 | 2014-01-15 | 中国石油天然气股份有限公司 | 一种根据核磁共振t2谱确定渗透率的方法和装置 |
| US20130257424A1 (en) * | 2012-03-27 | 2013-10-03 | Schlumberger Technology Corporation | Magnetic resonance rock analysis |
| CN104634804B (zh) * | 2013-11-08 | 2016-10-26 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
| CN103926267B (zh) * | 2014-04-28 | 2016-05-04 | 西安石油大学 | 一种定量评价应力敏感过程中孔喉变化程度的方法 |
| CN104634718B (zh) * | 2015-03-05 | 2017-03-22 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
-
2017
- 2017-08-17 CN CN201710705113.5A patent/CN107688037B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107688037A (zh) | 2018-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107688037B (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
| WO2017084454A1 (zh) | 一种地层组分最优化确定方法及装置 | |
| EP2834453B1 (en) | Methods for determining wettability from nmr measurements | |
| RU2315339C2 (ru) | Система петрофизической оценки в реальном времени | |
| EP2447469B1 (en) | Quantitative calculation method for hydrocarbon saturation of fractured reservoir | |
| CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN110318745B (zh) | 一种沉积微相约束下粒径岩性测井评价方法 | |
| CN108694264B (zh) | 一种确定页岩气储层渗透率的方法 | |
| CN105891089B (zh) | 储层渗透率确定方法及装置 | |
| CN106370814B (zh) | 基于成分-结构分类的湖相混积岩类储层测井识别方法 | |
| CA2903451A1 (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| CN108982320A (zh) | 一种利用粒度参数进行复杂孔隙结构储层渗透率计算方法 | |
| CN105064986A (zh) | 利用常规测录井资料构建储层四性关系谱的方法 | |
| CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
| CN102052074A (zh) | 一种岩心数字化确定岩电参数的方法 | |
| Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
| CN105205296A (zh) | 一种求取页岩气储层孔隙度的方法 | |
| Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
| CN106468169B (zh) | 计算油藏含油饱和度的方法 | |
| Newsham et al. | Tutorial: Organic-Mudstone Petrophysics: Workflow to Estimate Storage Capacity | |
| Feng et al. | A new empirical method for constructing capillary pressure curves from conventional logs in low-permeability sandstones | |
| CN103197348B (zh) | 利用各层内部样品进行加权编制测井交会图的方法 | |
| CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
| CN106837318B (zh) | 岩层稠油含量获取方法和装置 | |
| CN108412488B (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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: Zhanjiang Branch of CNOOC (China) Co., Ltd. Address before: Chinese CNOOC Building No. 25 Beijing 100010 Dongcheng District Chaoyangmen North Street Applicant before: CHINA NATIONAL OFFSHORE OIL Corp. Applicant before: Zhanjiang Branch of CNOOC (China) Co., Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |