CN111946336A - 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 - Google Patents
基于矿物分析数据获取页岩油储层核磁孔隙度的方法 Download PDFInfo
- Publication number
- CN111946336A CN111946336A CN202010925211.1A CN202010925211A CN111946336A CN 111946336 A CN111946336 A CN 111946336A CN 202010925211 A CN202010925211 A CN 202010925211A CN 111946336 A CN111946336 A CN 111946336A
- Authority
- CN
- China
- Prior art keywords
- nuclear magnetic
- shale oil
- content
- porosity
- magnetic porosity
- 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
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims abstract description 48
- 239000011707 mineral Substances 0.000 title claims abstract description 48
- 239000003079 shale oil Substances 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011435 rock Substances 0.000 claims abstract description 19
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 229910021532 Calcite Inorganic materials 0.000 claims description 21
- 229910000514 dolomite Inorganic materials 0.000 claims description 21
- 239000010459 dolomite Substances 0.000 claims description 21
- 229910052655 plagioclase feldspar Inorganic materials 0.000 claims description 11
- 229910052661 anorthite Inorganic materials 0.000 claims description 10
- GWWPLLOVYSCJIO-UHFFFAOYSA-N dialuminum;calcium;disilicate Chemical compound [Al+3].[Al+3].[Ca+2].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] GWWPLLOVYSCJIO-UHFFFAOYSA-N 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- 230000005311 nuclear magnetism Effects 0.000 abstract description 5
- 239000000463 material Substances 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000012937 correction Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 5
- 239000011148 porous material Substances 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 3
- 230000007935 neutral effect Effects 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 2
- 239000002366 mineral element Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000004058 oil shale Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012795 verification 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
- 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
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- Mathematical Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Data Mining & Analysis (AREA)
- Geophysics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及储层核磁孔隙度计算方法技术领域，是一种基于矿物分析数据获取页岩油储层核磁孔隙度的方法，其通过公式Ip＝100*(Kp*Pla+Kc*Cal+Kd*Dol+K)计算页岩油储层的核磁孔隙度。本发明基于矿物分析数据获取页岩油储层核磁孔隙度的方法，只需要利用成本较低的随钻录井采集的岩屑进行矿物分析，进而应用矿物含量拟合计算核磁孔隙度；由于是对岩屑实物进行矿物含量测量，测量值不受井筒环境影响，矿物含量测量精度较高，进而使得计算得到的核磁孔隙度数值精度也较高。
Description
技术领域
本发明涉及储层核磁孔隙度计算方法技术领域，是一种基于矿物分析数据获取页岩油储层核磁孔隙度的方法。
背景技术
在新疆油田页岩油区块(例如吉木萨尔页岩油区块)面临的问题是：通常是每口水平井完钻后进行核磁测井，得到核磁孔隙度参数用于储层评价；但是由于核磁测井费用较高，同时因井筒条件所限，部分水平井无法进行核磁测井，导致无法获得完整的核磁孔隙度参数用于储层评价。
授权公开号为106443800的中国专利文献公开了一种中基性火山岩的核磁孔隙度校正方法，该方法包括：结合ECS测井对实验数据进行统计分析，优选影响中基性火山岩核磁孔隙度的敏感矿物元素，建立核磁孔隙度和ECS测井矿物元素含量的校正模型，模型中的参数采用最优化拟合法确定；最后利用校正模型实现对核磁测井孔隙度的校正，得到反映地层真实情况的孔隙度曲线。该发明在火山岩岩石物理实验的基础上，根据ECS测井建立核磁孔隙度的校正模型，实现了对中基性火山岩核磁测井孔隙度的连续校正，提高了核磁测井孔隙度的计算精度。
而页岩油储层特征明显不同于火成岩(火山岩)；火成岩储层属于基质孔隙与缝洞双重储集空间，而页岩油储层以基质孔隙储集空间为主，基质孔隙中又以溶蚀孔隙类型占主导地位，溶蚀孔隙发育程度与特定矿物含量密切相关；火成岩核磁孔隙度校正方法是根据ECS元素测井测得的元素含量拟合计算得到校正核磁孔隙度模型，进而修正核磁孔隙度数值，该种方法首先必须在单井实施成本较高的ECS特殊测井项目才能得到元素含量测量值，其次，ECS元素测井测量的元素含量精度易受井筒环境影响，进而影响最后拟合计算的核磁孔隙度数值。
发明内容
本发明提供了一种基于矿物分析数据获取页岩油储层核磁孔隙度的方法，克服了上述现有技术之不足，其能有效解决现有核磁孔隙度测量方法成本较高，且火山岩储层核磁孔隙度计算方法不适用于页岩油储层核磁孔隙度计算的问题；同时能解决因井筒条件所限导致部分水平井无法进行核磁测井，从而难以获得完整的核磁孔隙度数据的问题。
本发明的技术方案是通过以下措施来实现的：一种基于矿物分析数据获取页岩油储层核磁孔隙度的方法，按下述方法进行：根据公式一计算页岩油储层的核磁孔隙度，
公式一：Ip＝100*(Kp*Pla+Kc*Cal+Kd*Dol+K)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol为白云石含量，％；
Kp、Kc、Kd、K为与斜长石、方解石和白云石矿物相关的系数；
将页岩油区块试验井的储层特征矿物中的斜长石、方解石和白云石的含量数据以及相应的测井核磁孔隙度数据代入公式一，求得相应的Kp、Kc、Kd、K数值，Kp、Kc、Kd、K数值确定后,得出确定系数的公式一，对页岩油区块正钻井随钻录井采集的岩屑进行矿物分析，得到斜长石、方解石和白云石矿物的质量百分含量，将斜长石、方解石和白云石矿物的质量百分含量代入Kp、Kc、Kd、K数值确定后的公式一中，就能计算得相应井的核磁孔隙度数据。
下面是对上述发明技术方案的进一步优化或/和改进：
上述对于新疆油田吉木萨尔页岩油区块，Kp＝0.002，Kc＝0.004，Kd＝0.001，K＝-0.016，其核磁孔隙度按下述公式计算，
Ip＝100*(0.002*Pla+0.004*Cal-0.001*Dol-0.016)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol：白云石含量，％。
本发明所述基于矿物分析数据获取页岩油储层核磁孔隙度的方法，只需要利用成本较低的随钻录井采集的岩屑进行矿物分析，进而应用矿物含量拟合计算核磁孔隙度；由于是对岩屑实物进行矿物含量测量，测量值不受井筒环境影响，矿物含量测量精度较高，进而使得计算得到的核磁孔隙度数值精度也较高。
具体实施方式
本发明不受下述实施例的限制，可根据本发明的技术方案与实际情况来确定具体的实施方式。本发明中的百分数如没有特殊说明，均为质量百分数。
下面结合实施例对本发明作进一步描述：
实施例1：该基于矿物分析数据获取页岩油储层核磁孔隙度的方法，按下述方法进行：根据公式一计算页岩油储层的核磁孔隙度，
公式一：Ip＝100*(Kp*Pla+Kc*Cal+Kd*Dol+K)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol为白云石含量，％；
Kp、Kc、Kd、K分别为与斜长石、方解石和白云石矿物相关的系数；
将页岩油区块试验井的储层特征矿物中的斜长石、方解石和白云石的含量数据以及相应的测井核磁孔隙度数据代入公式一，求得相应的Kp、Kc、Kd、K数值，Kp、Kc、Kd、K数值确定后,得出确定系数的公式一，对页岩油区块正钻井随钻录井采集的岩屑进行矿物分析，得到斜长石、方解石和白云石矿物的质量百分含量，将斜长石、方解石和白云石矿物的质量百分含量代入Kp、Kc、Kd、K数值确定后的公式一中，就能计算得相应井的核磁孔隙度数据。
通过对大量页岩油水平井储层的矿物分析数据和核磁测井数据进行相关性分析研究，找到了影响页岩油储层核磁孔隙度的敏感矿物组合，即斜长石、方解石和白云石的矿物组合。
不同区块的页岩油储层特性不同，那么对应的Kp、Kc、Kd、K也不相同。
对于某个区块的页岩油储层，按照实施例1所述方法，先将该区块个别井筒的矿物组合含量和核磁孔隙度数据代入上述公式后，得到相应的Kp、Kc、Kd、K数值，即可得到该区块的核磁孔隙度计算公式，对于该区块的其它井筒，无需核磁测井，只需将随钻录井采集的岩屑进行矿物分析，得到斜长石、方解石和白云石矿物的质量百分含量，即可计算得到该井筒处的核磁孔隙度，以此方法获得该区块完整的核磁孔隙度数据，以满足页岩油储层的评价需求。
通过多口水平井的测井核磁孔隙度与根据本发明方法计算的核磁孔隙度进行验证对比，根据本发明方法计算得到的核磁孔隙度精度达到90％以上，满足页岩油储层评价需求。
本发明基于矿物分析数据获取页岩油储层核磁孔隙度的方法，只需要利用成本较低的随钻录井采集的岩屑进行矿物分析，进而应用矿物含量拟合计算核磁孔隙度；由于是对岩屑实物进行矿物含量测量，测量值不受井筒环境影响，矿物含量测量精度较高，进而使得计算得到的核磁孔隙度数值精度也较高。
实施例2：对于新疆油田吉木萨尔页岩油区块，Kp＝0.002，Kc＝0.004，Kd＝0.001，K＝-0.016，其核磁孔隙度按下述公式计算，
Ip＝100*(0.002*Pla+0.004*Cal-0.001*Dol-0.016)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol：白云石含量，％。
根据上述公式，只需将新疆油田吉木萨尔页岩油区块各个井筒的随钻录井采集的岩屑进行矿物分析，得到斜长石、方解石和白云石矿物的质量百分含量，即可计算得到该井筒处的核磁孔隙度，以此方法获得该区块完整的核磁孔隙度数据；相比于以往的核磁孔隙度获取方法(核磁测井等)，采用本方法能够直接计算出新疆油田吉木萨尔页岩油区块的核磁孔隙度，计算精度能达到90％以上，计算出的核磁孔隙度参数精度不仅能满足储层评价的需求，同时获取孔隙度参数的成本相比核磁测井降低50％以上。
以上技术特征构成了本发明的实施例，其具有较强的适应性和实施效果，可根据实际需要增减非必要的技术特征，来满足不同情况的需求。
Claims (2)
1.一种基于矿物分析数据获取页岩油储层核磁孔隙度的方法，其特征在于按下述方法进行：根据公式一计算页岩油储层的核磁孔隙度，
公式一：Ip＝100*(Kp*Pla+Kc*Cal+Kd*Dol+K)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol为白云石含量，％；
Kp、Kc、Kd、K为与斜长石、方解石和白云石矿物以及区块地质特征相关的系数；
将页岩油区块试验井的储层特征矿物中的斜长石、方解石和白云石的含量数据以及相应的测井核磁孔隙度数据代入公式一，求得相应的Kp、Kc、Kd、K数值，Kp、Kc、Kd、K数值确定后,得出确定系数的公式一，对页岩油区块正钻井随钻录井采集的岩屑进行矿物分析，得到斜长石、方解石和白云石矿物的质量百分含量，将斜长石、方解石和白云石矿物的质量百分含量代入Kp、Kc、Kd、K数值确定后的公式一中，就能计算得相应井的核磁孔隙度数据。
2.根据权利要求1所述的基于矿物分析数据获取页岩油储层核磁孔隙度的方法，其特征在于对于新疆油田吉木萨尔页岩油区块，Kp＝0.002，Kc＝0.004，Kd＝0.001，K＝-0.016，其核磁孔隙度按下述公式计算，
Ip＝100*(0.002*Pla+0.004*Cal-0.001*Dol-0.016)
式中，Ip为核磁孔隙度，％；
Pla为斜长石含量，％；
Cal为方解石含量，％；
Dol：白云石含量，％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010925211.1A CN111946336A (zh) | 2020-09-06 | 2020-09-06 | 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010925211.1A CN111946336A (zh) | 2020-09-06 | 2020-09-06 | 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111946336A true CN111946336A (zh) | 2020-11-17 |
Family
ID=73356236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010925211.1A Pending CN111946336A (zh) | 2020-09-06 | 2020-09-06 | 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111946336A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114035227A (zh) * | 2021-11-11 | 2022-02-11 | 中国海洋石油集团有限公司 | 基于随钻xrd全岩录井的变质岩潜山储层孔隙度预测方法 |
| CN114791482A (zh) * | 2021-01-25 | 2022-07-26 | 中国石油化工股份有限公司 | 一种确定岩石矿物含量系数的方法及装置、存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2360108C1 (ru) * | 2007-12-14 | 2009-06-27 | Назия Мингалиевна Данилова | Способ определения пористости и проницаемости пласта месторождений нефти и газа |
| CN106290103A (zh) * | 2015-06-12 | 2017-01-04 | 中国石油化工股份有限公司 | 一种页岩气储层中粘土微孔孔隙度的测定方法 |
| CN106600436A (zh) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | 计算页岩气地层矿物组分含量及孔隙度的方法 |
| CN106841001A (zh) * | 2017-01-17 | 2017-06-13 | 西南石油大学 | 一种基于储层质量主控因素分析的致密砂岩孔隙度、渗透率预测方法 |
| CN106950347A (zh) * | 2017-03-03 | 2017-07-14 | 中国石油大学(华东) | 一种评价泥页岩各组分体积的方法 |
| CN109100815A (zh) * | 2018-09-13 | 2018-12-28 | 中国石油天然气股份有限公司 | 一种确定地层孔隙度的方法、装置及系统 |
-
2020
- 2020-09-06 CN CN202010925211.1A patent/CN111946336A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2360108C1 (ru) * | 2007-12-14 | 2009-06-27 | Назия Мингалиевна Данилова | Способ определения пористости и проницаемости пласта месторождений нефти и газа |
| CN106290103A (zh) * | 2015-06-12 | 2017-01-04 | 中国石油化工股份有限公司 | 一种页岩气储层中粘土微孔孔隙度的测定方法 |
| CN106600436A (zh) * | 2016-11-28 | 2017-04-26 | 中国石油集团川庆钻探工程有限公司 | 计算页岩气地层矿物组分含量及孔隙度的方法 |
| CN106841001A (zh) * | 2017-01-17 | 2017-06-13 | 西南石油大学 | 一种基于储层质量主控因素分析的致密砂岩孔隙度、渗透率预测方法 |
| CN106950347A (zh) * | 2017-03-03 | 2017-07-14 | 中国石油大学(华东) | 一种评价泥页岩各组分体积的方法 |
| CN109100815A (zh) * | 2018-09-13 | 2018-12-28 | 中国石油天然气股份有限公司 | 一种确定地层孔隙度的方法、装置及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 王道富;王玉满;董大忠;王世谦;黄金亮;黄勇斌;王淑芳;李新景;: "川南下寒武统筇竹寺组页岩储集空间定量表征", 天然气工业, no. 07, pages 1 - 10 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114791482A (zh) * | 2021-01-25 | 2022-07-26 | 中国石油化工股份有限公司 | 一种确定岩石矿物含量系数的方法及装置、存储介质 |
| CN114791482B (zh) * | 2021-01-25 | 2024-03-08 | 中国石油化工股份有限公司 | 一种确定岩石矿物含量系数的方法及装置、存储介质 |
| CN114035227A (zh) * | 2021-11-11 | 2022-02-11 | 中国海洋石油集团有限公司 | 基于随钻xrd全岩录井的变质岩潜山储层孔隙度预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110644980B (zh) | 一种超低渗透油藏储层综合分类评价方法 | |
| CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
| CN111946336A (zh) | 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 | |
| CN104750884B (zh) | 基于多因素非线性回归的页岩油气富集指数定量评价方法 | |
| CN111946398A (zh) | 一种复合地层盾构掘进效率现场预测计算方法 | |
| CN113419284B (zh) | 一种基于聚类分析测井岩石物理相双甜点识别方法 | |
| CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
| CN102900432B (zh) | 利用微钻时数据随钻计算录井孔隙度进行储层评价的方法 | |
| CN110471127B (zh) | 一种测井渗透率解释方法 | |
| CN115030707A (zh) | 一种油页岩“甜点”快速评价方法 | |
| CN110593855B (zh) | 一种测井电阻率的校正方法、油层的识别方法 | |
| CN107025358A (zh) | 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 | |
| CN104675391B (zh) | 计算地层含油饱和度的方法 | |
| CN113777655A (zh) | 一种盐湖盆地小层烃源岩平面分布预测方法 | |
| CN114677460A (zh) | 一种合成gr曲线的方法及应用 | |
| CN107808055A (zh) | 一种基于双重影响的页岩气饱和度测井计算方法 | |
| CN109343133B (zh) | 基于二维核磁共振测井技术的扩径段孔隙度校正方法 | |
| Brindle et al. | An integrated approach to unconventional resource play reservoir characterization, Thistleton-1 case study, NW England | |
| CN114429015A (zh) | 页岩气藏地层模型建立方法、电子设备及介质 | |
| CN110570102A (zh) | 一种储层评价方法 | |
| CN108150158A (zh) | 一种深层裂缝性致密砂岩气藏早期水体分析与预测方法 | |
| CN113763194B (zh) | 页岩单井含气丰度预测方法及系统 | |
| CN111487157B (zh) | 一种基于现场解吸与核磁联测的页岩损失气计算方法 | |
| CN112184033B (zh) | 一种碳酸盐岩岩层断层封闭性评价方法 | |
| CN110513098B (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 |