CN114075965B - 一种河道砂岩粒度的确定方法 - Google Patents
一种河道砂岩粒度的确定方法 Download PDFInfo
- Publication number
- CN114075965B CN114075965B CN202010777823.0A CN202010777823A CN114075965B CN 114075965 B CN114075965 B CN 114075965B CN 202010777823 A CN202010777823 A CN 202010777823A CN 114075965 B CN114075965 B CN 114075965B
- Authority
- CN
- China
- Prior art keywords
- granularity
- sandstone
- logging curve
- core
- core sample
- 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 26
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 239000011229 interlayer Substances 0.000 claims description 8
- 238000000611 regression analysis Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 description 14
- 239000011435 rock Substances 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003921 particle size analysis Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Computer Hardware Design (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Agronomy & Crop Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种河道砂岩粒度的确定方法，属于砂岩油气藏地质勘探技术领域。本发明利用河道砂岩取心段的岩心样品，通过实验分析，确定岩心样品的对应的粒度，基于该取心段对应的测井数据，建立岩心粒度与测井数据之间的关系，利用该关系，根据没有取心段的测井数据确定出没有取心段的砂岩的粒度。因此，本发明能够根据测井曲线实现对砂岩粒度的确定，无需获取大量的岩心数据，能够快速、准确地实现对砂岩粒度的确定。
Description
技术领域
本发明涉及一种河道砂岩粒度的确定方法，属于砂岩油气藏地质勘探技术领域。
背景技术
储层评价是油气藏地质勘探的一个重要工作，在评价储集层的时候，往往从储层的岩性、物性等方面入手。其中，储层岩性、物性主要受岩石粒度的影响。有效的分析河道砂岩粒度、评价储层是地质勘探和油气开发共同需求。
目前常用的砂岩粒度确定大都采用岩心数据，即通过获取大量岩心样品，对岩心样品进行试验测量得到对应的砂岩粒度，但是由于取心数量有限，所以粒度分析并不能有效的覆盖所有储层段，同时因需要在实验室操作所以粒度分析效率低、不方便在现场使用。例如申请公布号为CN107688037A的专利文件，该文件公开了一种利用核磁测井T2分布确定井下岩石粒度曲线的方法，该方法首先对岩心样品进行实验分析获取每块岩心的核磁T2分布及岩石粒度分布曲线，然后利用孔隙度分析仪测量多块岩心的孔隙度参数，并通过核磁共振测井装置采集和记录地层核磁T2分布数据以及地层孔隙度参数；最后选取区域T2分布转粒度分布曲线的转换系数，确定地层的粒度分布。本发明虽能够实现对岩石粒度的确定，但是该方法需要孔隙度分析仪、核磁共振实验分析仪等设备，实现起来比较复杂。
发明内容
本发明的目的是提供一种河道砂岩粒度的确定方法，以解决目前河道砂岩粒度分析存在复杂度高的问题。
本发明为解决上述技术问题而提供一种河道砂岩粒度的确定方法，该方法包括以下步骤：
1)获取目标工区测井曲线，所述测井曲线包括自然伽马测井曲线、声波时差测井曲线、补偿中子曲线和密度测井曲线；
2)将获取的测井曲线代入到砂岩粒度与测井曲线关系式中，确定出对应的砂岩粒度，所述砂岩粒度与测井曲线关系式是通过对目标工区中的岩心样品粒度与对应的测井曲线进行多元回归分析得到。
本发明利用河道砂岩取心段的岩心样品，通过实验分析，确定岩心样品的对应的粒度，基于该取心段对应的测井数据，建立岩心粒度与测井数据之间的关系，利用该关系，根据没有取心段的测井数据确定出没有取心段的砂岩的粒度。因此，本发明能够根据测井曲线实现对砂岩粒度的确定，无需获取大量的岩心数据，能够快速、准确地实现对砂岩粒度的确定。
进一步地，所述砂岩粒度与测井曲线关系式为：
式中，为砂岩粒度，GR为实测的自然伽马值，单位为API；DEN为实测的补偿密度，单位为g/cm3；AC为实测的声波时差，单位为μm/s；CNL为实测的补偿中子，单位为％。
进一步地，为了准确获取砂岩粒度与测井曲线的关系，所述砂岩粒度与测井曲线关系的建立过程如下：
a.获取目标工区岩心样品的粒度；
b.在单井上进行岩心样品深度归位，并获取归位后岩心样品深度对应的测井曲线；
c.将岩心样品的粒度和对应深度的测井曲线进行多元回归分析，建立砂岩粒度与测井曲线之间的关系。
进一步地，所述步骤a中的岩心样品的粒度是通过激光粒度分析仪对岩心样品进行试验分析得到。
进一步地，为了提高归位的准确性，所述步骤b是通过选取取心段泥岩夹层，利用泥岩夹层的测井曲线相应特征进行深度归位，其中，泥岩夹层与砂岩段相比较测井曲线特征表现为高自然伽马、高声波时差、高补偿中子、低补偿密度。
附图说明
图1是本发明河道砂岩粒度的确定方法的流程图；
图2是本发明河道砂岩粒度的确定方法的应用成果图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步地说明。
本发明利用河道砂岩取心段的岩心样品，通过实验分析，确定岩心样品的对应的粒度，基于该取心段对应的测井数据，建立岩心粒度与测井数据之间的关系，利用该关系，根据没有取心段的测井数据确定出没有取心段的砂岩的粒度，该方法的具体实现流程如图1所示，具体实现过程如下。
1.获取河道砂岩取心段的岩心样品，确定每个岩心样品的平均粒径。
对本实施例而言，对河道砂岩取心段进行采样，每米4～5个采样点，将获取的各采样点的岩心样品拿到实验室进行粒度分析，确定出各岩心样品的粒度，本实施例采用激光粒度分析仪测量出砂岩颗粒的直径，砂岩粒度划分的方案采用值标度，/>值的计算公式如下：
其中D为颗粒直径，通过对岩心样品进行实验测试得到。
2.获取河道砂岩取心段的测井数据，在单井上对岩心样品进行深度归位，获取每个岩心样品所在位置的测井曲线。
由于泥岩夹层与砂岩段相比较测井曲线特征表现为高自然伽马、高声波时差、高补偿中子、低补偿密度。因此本发明选取取心段泥岩夹层，利用泥岩测井曲线响应特征进行深度归位。由于自然伽马测井曲线(GR)、声波时差测井曲线(AC)、补偿中子曲线(CNL)和密度测井曲线(DEN)能够反映岩性及物性的变化，对本实施例而言，选取的测井曲线包括GR、AC、CNL、DEN曲线。
3.建立岩性粒径与所选测井曲线之间的关系。
根据实验测得的取心段的粒径，与对应深度的GR、AC、CN、DEN曲线进行多元回归，得到砂岩粒径与GR、AC、CN、DEN曲线的多元回归关系，对本实施例而言，得到的多元回归关系为：
式中，GR为实测的自然伽马值，单位为API；DEN为实测的补偿密度，单位为g/cm3；AC为实测的声波时差，单位为μm/s；CNL为实测的补偿中子，单位为％。
4.获取没有取心段的测井数据，利用岩性粒径与所选测井曲线之间的关系确定对应的岩性粒径。
考虑到成本和时间问题，取心的数量是有限的，因此，对于没有取心数据的砂岩，砂岩的粒度需要基于步骤3中所建立的关系来确定。本实施例中根据上述关系确定的岩性粒径结果如图2所示，可以看出，本发明能够快速准确地得到砂岩粒径。
Claims (4)
1.一种河道砂岩粒度的确定方法，其特征在于，该方法包括以下步骤：
1）获取目标工区测井曲线，所述测井曲线包括自然伽马测井曲线、声波时差测井曲线、补偿中子曲线和密度测井曲线；
2）将获取的测井曲线代入到砂岩粒度与测井曲线关系式中，确定出对应的砂岩粒度，所述砂岩粒度与测井曲线关系式是通过对目标工区中的岩心样品粒度与对应的测井曲线进行多元回归分析得到；
所述砂岩粒度与测井曲线关系式为：
φ=0.01*GR+16.84*DEN-0.93CNL-0.1AC
式中，φ为砂岩粒度，GR为实测的自然伽马值，单位为API；DEN为实测的补偿密度，单位为g/cm3；AC为实测的声波时差，单位为μm/s；CNL为实测的补偿中子，单位为%。
2.根据权利要求1所述的河道砂岩粒度的确定方法，其特征在于，所述砂岩粒度与测井曲线关系式的建立过程如下：
a.获取目标工区岩心样品的粒度；
b.在单井上进行岩心样品深度归位，并获取归位后岩心样品深度对应的测井曲线；
c.将岩心样品的粒度和对应深度的测井曲线进行多元回归分析，建立砂岩粒度与测井曲线之间的关系式。
3.根据权利要求2所述的河道砂岩粒度的确定方法，其特征在于，步骤a中的岩心样品的粒度是通过激光粒度分析仪对岩心样品进行试验分析得到。
4.根据权利要求2所述的河道砂岩粒度的确定方法，其特征在于，步骤b是通过选取取心段泥岩夹层，利用泥岩夹层的测井曲线相应特征进行深度归位，其中，泥岩夹层与砂岩段相比较测井曲线特征表现为高自然伽马、高声波时差、高补偿中子、低补偿密度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010777823.0A CN114075965B (zh) | 2020-08-05 | 2020-08-05 | 一种河道砂岩粒度的确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010777823.0A CN114075965B (zh) | 2020-08-05 | 2020-08-05 | 一种河道砂岩粒度的确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114075965A CN114075965A (zh) | 2022-02-22 |
| CN114075965B true CN114075965B (zh) | 2024-05-07 |
Family
ID=80279568
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010777823.0A Active CN114075965B (zh) | 2020-08-05 | 2020-08-05 | 一种河道砂岩粒度的确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114075965B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606150A (zh) * | 2012-03-30 | 2012-07-25 | 中国海洋石油总公司 | 一种碎屑岩岩性识别方法及*** |
| US8812236B1 (en) * | 2014-04-11 | 2014-08-19 | Particle Size Engineering, LLC | Method for using particle size analysis in near time or real time to create a proper particle size distribution within a drilling fluid management system for improved well drilling efficiency |
| CN107688037A (zh) * | 2017-08-17 | 2018-02-13 | 中国海洋石油总公司 | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 |
| CN108303752A (zh) * | 2018-02-11 | 2018-07-20 | 中国石油化工股份有限公司 | 砂砾岩有效储层常规测井定量识别方法 |
| CN110318745A (zh) * | 2019-06-10 | 2019-10-11 | 中国石油大学(华东) | 一种沉积微相约束下粒径岩性测井评价方法 |
| CN111323844A (zh) * | 2020-03-14 | 2020-06-23 | 长江大学 | 一种基于曲线重构的复杂砂砾岩体的岩性识别方法及*** |
| CN111425190A (zh) * | 2020-03-19 | 2020-07-17 | 中国石油大学(华东) | 一种页岩气地层岩性识别方法、***、存储介质、终端 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190266501A1 (en) * | 2018-02-27 | 2019-08-29 | Cgg Services Sas | System and method for predicting mineralogical, textural, petrophysical and elastic properties at locations without rock samples |
-
2020
- 2020-08-05 CN CN202010777823.0A patent/CN114075965B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606150A (zh) * | 2012-03-30 | 2012-07-25 | 中国海洋石油总公司 | 一种碎屑岩岩性识别方法及*** |
| US8812236B1 (en) * | 2014-04-11 | 2014-08-19 | Particle Size Engineering, LLC | Method for using particle size analysis in near time or real time to create a proper particle size distribution within a drilling fluid management system for improved well drilling efficiency |
| CN107688037A (zh) * | 2017-08-17 | 2018-02-13 | 中国海洋石油总公司 | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 |
| CN108303752A (zh) * | 2018-02-11 | 2018-07-20 | 中国石油化工股份有限公司 | 砂砾岩有效储层常规测井定量识别方法 |
| CN110318745A (zh) * | 2019-06-10 | 2019-10-11 | 中国石油大学(华东) | 一种沉积微相约束下粒径岩性测井评价方法 |
| CN111323844A (zh) * | 2020-03-14 | 2020-06-23 | 长江大学 | 一种基于曲线重构的复杂砂砾岩体的岩性识别方法及*** |
| CN111425190A (zh) * | 2020-03-19 | 2020-07-17 | 中国石油大学(华东) | 一种页岩气地层岩性识别方法、***、存储介质、终端 |
Non-Patent Citations (3)
| Title |
|---|
| Bayes判别法在麻黄山西探区油层识别中的应用;陈利雯等;勘探地球物理进展;31(04);第301-304页 * |
| 储层粒度神经网络预测模型研究;王利华等;西南石油大学学报(自然科学版);38(01);第53-59页 * |
| 吕国祥等.《计算机地质制图及应用》.电子科技拟大学出版社,2009,第13页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114075965A (zh) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106468172B (zh) | 一种超低渗砂岩油藏低阻储层测井解释方法 | |
| CA2903451C (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| WO2017084454A1 (zh) | 一种地层组分最优化确定方法及装置 | |
| CN1693893A (zh) | 确定地层流体性质的方法 | |
| CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及*** | |
| CN104075974A (zh) | 一种利用低场核磁共振精确测定页岩孔隙度的方法 | |
| CN111007230B (zh) | 定量评价陆相湖盆低孔隙度致密油储层含油量的方法 | |
| CN107688037B (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
| CN108561126B (zh) | 一种确定页岩气储层有机孔隙度的简易方法 | |
| CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
| CN111812736B (zh) | 一种致密砂岩无水气藏含气性评价方法 | |
| Liu et al. | Pore structure and physical properties of sandy conglomerate reservoirs in the Xujiaweizi depression, northern Songliao Basin, China | |
| CN110905493B (zh) | 一种测量井下地层流体污染率的方法 | |
| CN106290103A (zh) | 一种页岩气储层中粘土微孔孔隙度的测定方法 | |
| CN113216945B (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
| Zhou et al. | Pore structure analysis and classification of pyroclastic reservoirs in the Dehui fault depression based on experimental and well-logging data | |
| CN114075965B (zh) | 一种河道砂岩粒度的确定方法 | |
| CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
| CN108412488B (zh) | 快速确定页岩气储层有机孔隙度的测井方法 | |
| CN114086938B (zh) | 一种非均质性砂岩储层的含气饱和度预测方法 | |
| Spain et al. | Petrophysical rock typing of tight gas sands-beyond porosity and saturation: example from the cotton valley formation, East Texas | |
| CN103197348A (zh) | 利用各层内部样品进行加权编制测井交会图的方法 | |
| WO2022204297A1 (en) | The combined discrete gaussian analysis of micp and nmr t2 distributions of multi-modal carbonate rocks | |
| CN108647417B (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 |