CN115434698A - 基于光电吸收截面指数的地层钙质胶结物含量计算方法 - Google Patents
基于光电吸收截面指数的地层钙质胶结物含量计算方法 Download PDFInfo
- Publication number
- CN115434698A CN115434698A CN202211076515.0A CN202211076515A CN115434698A CN 115434698 A CN115434698 A CN 115434698A CN 202211076515 A CN202211076515 A CN 202211076515A CN 115434698 A CN115434698 A CN 115434698A
- Authority
- CN
- China
- Prior art keywords
- quartz
- feldspar
- clay
- content
- 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.)
- Granted
Links
- 239000004568 cement Substances 0.000 title claims abstract description 111
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 title claims abstract description 99
- 239000011575 calcium Substances 0.000 title claims abstract description 99
- 229910052791 calcium Inorganic materials 0.000 title claims abstract description 99
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 20
- 239000010433 feldspar Substances 0.000 claims abstract description 93
- 238000012937 correction Methods 0.000 claims abstract description 71
- 239000010453 quartz Substances 0.000 claims abstract description 66
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims abstract description 66
- 239000004927 clay Substances 0.000 claims abstract description 61
- 239000012530 fluid Substances 0.000 claims abstract description 56
- 238000005553 drilling Methods 0.000 claims abstract description 55
- TZCXTZWJZNENPQ-UHFFFAOYSA-L barium sulfate Chemical compound [Ba+2].[O-]S([O-])(=O)=O TZCXTZWJZNENPQ-UHFFFAOYSA-L 0.000 claims abstract description 53
- 239000010428 baryte Substances 0.000 claims abstract description 53
- 229910052601 baryte Inorganic materials 0.000 claims abstract description 53
- 230000004044 response Effects 0.000 claims abstract description 52
- 238000004458 analytical method Methods 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000002360 preparation method Methods 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims abstract description 4
- 239000004698 Polyethylene Substances 0.000 claims description 266
- 229920000573 polyethylene Polymers 0.000 claims description 257
- -1 polyethylene Polymers 0.000 claims description 22
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 21
- 239000011707 mineral Substances 0.000 claims description 21
- 239000010459 dolomite Substances 0.000 claims description 12
- 229910000514 dolomite Inorganic materials 0.000 claims description 12
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 claims description 12
- AXCZMVOFGPJBDE-UHFFFAOYSA-L calcium dihydroxide Chemical compound [OH-].[OH-].[Ca+2] AXCZMVOFGPJBDE-UHFFFAOYSA-L 0.000 claims description 11
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 229910021532 Calcite Inorganic materials 0.000 claims description 7
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 claims description 7
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 claims description 7
- 229910052900 illite Inorganic materials 0.000 claims description 7
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 claims description 7
- 229910052622 kaolinite Inorganic materials 0.000 claims description 7
- 229910052901 montmorillonite Inorganic materials 0.000 claims description 7
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 claims description 7
- 150000003839 salts Chemical class 0.000 claims description 7
- 239000011148 porous material Substances 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 229910001919 chlorite Inorganic materials 0.000 claims description 5
- 229910052619 chlorite group Inorganic materials 0.000 claims description 5
- 229910052742 iron Inorganic materials 0.000 claims description 5
- 239000012267 brine Substances 0.000 claims description 4
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 3
- 239000011780 sodium chloride Substances 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 2
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 208000010378 Pulmonary Embolism Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009545 invasion 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
- 230000000704 physical effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- 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
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了基于光电吸收截面指数的地层钙质胶结物含量计算方法，具体包括以下步骤：S1、数据准备，S2、PE原始曲线响应分析，S3、分段式PE基线确定，S4、最小PE基线段的石英与长石骨架PE值确定，S5、重晶石钻井液泥饼校正，S6、其他PE基线段的石英与长石骨架PE值确定，S7、粘土骨架PE值和钙质胶结物骨架PE值确定，S8、扩径校正，S9、钙质胶结物含量计算，本发明涉及油气勘探开发技术领域。该基于光电吸收截面指数的地层钙质胶结物含量计算方法，通过建立分段式PE基线，确定岩石中石英与长石骨架、粘土骨架和钙质胶结物的骨架PE值，对PE曲线进行重晶石钻井液泥饼校正和扩径校正，最终实现全井段的钙质胶结物含量计算。
Description
技术领域
本发明涉及油气勘探开发技术领域，具体为基于光电吸收截面指数的地层钙质胶结物含量计算方法。
背景技术
砂岩油气储层一直都是各大油气田增储上产的主要目标，但地层中的钙质胶结物极大地影响了砂岩储层的物性和有效性。如何准确计算钙质胶结物含量成为测井评价的一个重要内容。目前常用的钙质胶结物含量计算方法主要依靠自然伽马、声波时差、密度、中子及电阻率等测井曲线，但这些测井曲线受到泥质含量、孔隙度、流体性质、钻井液侵入、井眼扩径等因素的综合影响，其评价效果受到很大的制约，其准确度不高。
光电吸收截面指数(简称PE)是油气田勘探阶段较为普遍的测井曲线之一，由于PE对岩性(矿物组分)非常敏感，因此利用PE曲线进行钙质胶结物含量计算具有较好的基础。然而，在实际钻井过程中，钻井液中经常添加重晶石，且重晶石对PE曲线有非常严重的影响，导致PE曲线失真。这使得PE曲线的作用被长期忽略，没有发挥其应有的作用。如果能解决钻井液中重晶石的影响，并研发一种基于PE曲线的钙质胶结物含量计算方法，将会更好的确定砂岩储层有效性并提高砂岩储层勘探开发效率。
发明内容
(一)解决的技术问题
针对现有技术的不足，本发明提供了基于光电吸收截面指数的地层钙质胶结物含量计算方法，通过针对砂岩储层油气勘探过程中，在测量有光电吸收截面指数(简称PE)曲线的情况下，提出了一套方法来校正钻井液中重晶石对PE曲线的影响，并利用校正后的PE曲线进行钙质胶结物含量的计算，以更好的确定砂岩储层有效性并提高砂岩储层勘探开发效率。
(二)技术方案
为实现以上目的，本发明通过以下技术方案予以实现：基于光电吸收截面指数的地层钙质胶结物含量计算方法，具体包括以下步骤：
S1、数据准备：准备好研究区包含PE原始曲线的测井数据，如果有研究区少量的岩心矿物含量分析数据更好；
S2、PE原始曲线响应分析：PE原始数据由6部分组成，包括重晶石钻井液泥饼的响应、扩径的响应、盐水的响应、粘土的响应、石英与长石矿物的响应以及钙质胶结物的响应，按照井的深度由浅至深分别命名为第①段PE基线、第②段PE基线、第③段PE基线……等，依此类推，一般不超过10个平台段；
S3、分段式PE基线确定：利用测井数据，以深度为横坐标并以PE原始数据为纵坐标绘制连线图，根据PE原始曲线的形态将其分为多个平台段，在每一段选择若干个PE值最低的纯砂岩层，取其PE平均值确定每段的PE基线值；
S4、最小PE基线段的石英与长石骨架PE值确定：PE基线由两部分构成：重晶石钻井液泥饼的PE响应和石英与长石的PE响应，PE基线即表示为以下公式：
PE基线＝PE重晶石钻井液泥饼+PE石英与长石
式中：PE基线为平台段内几个PE值最低纯砂岩层的PE平均值(B/e)，在所划分的多段PE基线中，选择数值最小的那段PE基线，由于石英骨架的PE值(1.81B/e)低于长石骨架的PE值(2.86B/e)(表1)，因此，数值最小的那段PE基线所对应的地层具有全井段最高的石英矿物含量；
表1主要矿物及流体的PE骨架值
S5、重晶石钻井液泥饼校正：针对数值最小的PE基线段，读取其PE基线值(PE基线min)，并令该段的PE石英与长石＝PE石英与长石骨架，将这两个参数代入步骤S4中的PE基线表示公式中，可得PE重晶石钻井液泥饼：
PE重晶石钻井液泥饼＝PE基线min-PE石英与长石骨架；
利用以下公式计算得到重晶石钻井液泥饼校正后的PE曲线，即PE校正1曲线：
PE校正1＝PE原始-PE重晶石钻井液泥饼
式中：PE校正1为重晶石钻井液泥饼校正后的PE曲线(B/e)；
S6、其他PE基线段的石英与长石骨架PE值确定：在不同平台段，PE重晶石钻井液泥饼的值保持不变，但在不同平台段，石英与长石的比例不同，则其PE石英与长石骨架也不同，在数值最小的PE基线段计算出PE重晶石钻井液泥饼的基础上，利用以下公式计算其他段的PE石英与长石骨架：
PE石英与长石骨架＝PE基线-PE重晶石钻井液泥饼；
S7、粘土骨架PE值和钙质胶结物骨架PE值确定；
S8、扩径校正：找到典型扩径段(一般CAL-BIT>0.2in)，典型扩径段的PE校正1明显高于正常井径段的PE校正1，选择该段数据，以CAL-BIT为横坐标，并以PE校正1为纵坐标，绘制交会图并获得线性拟合关系式：
PE校正1＝a×(CAL-BIT)+b
式中：CAL为井径(in)；BIT为钻头直径(in)；a、b为拟合系数；
S9、钙质胶结物含量计算：PE曲线经过重晶石钻井液泥饼校正和扩径校正后得到PE校正2曲线，PE校正2主要由4部分构成：盐水的响应、粘土的响应、石英与长石的响应、钙质胶结物的响应。PE校正2表示为下式：
PE校正2＝PE盐水+PE粘土+PE石英与长石+PE钙质胶结物
PE校正2可进一步表示为下式：
PE校正2＝POR×PE盐水骨架+V粘土×PE粘土骨架+V石英与长石×PE石英与长石骨架+V钙质胶结物×PE钙质胶结物骨架
式中：POR为孔隙度(小数)，可通过其他测井曲线计算得到；V粘土为粘土含量(小数)，通过其他测井曲线计算得到；V石英与长石为石英与长石含量(小数)；V钙质胶结物为钙质胶结物含量(小数)；PE盐水骨架为孔隙中盐水骨架的PE值(参考表1，PE盐水骨架＝1.64)。
在岩石中，孔隙、粘土、石英与长石、钙质胶结物的总体积为1，即：
POR+V粘土+V石英与长石+V钙质胶结物＝1
将上式代入PE校正2可进一步表示的公式，得到钙质胶结物含量(V钙质胶结物)的计算公式：
优选的，所述步骤S1中PE原始数据表示为以下公式：
PE原始＝PE重晶石钻井液泥饼+PE扩径+PE盐水+PE粘土+PE石英与长石+PE钙质胶结物
式中：PE原始为PE原始数据(B/e)、PE重晶石钻井液泥饼为重晶石钻井液泥饼对PE的响应(B/e)、PE扩径为扩径对PE的响应(B/e)、PE盐水为孔隙中盐水对PE的响应(B/e)、PE粘土为粘土对PE的响应(B/e)、PE石英与长石为石英与长石对PE的响应(B/e)以及PE钙质胶结物为钙质胶结物对PE的响应(B/e)。
优选的，所述步骤S2中在每一段选择3-5个PE值最低的纯砂岩层。
优选的，所述步骤S4中确定数值最小PE基线段的石英与长石骨架PE值(PE石英与长石骨架)的方法为：在该平台段有岩心矿物含量分析数据的情况下，优先利用以下公式计算该段的PE石英与长石骨架：；
式中：PE石英与长石骨架为平台段内石英与长石骨架的PE值(B/e)；V石英C为岩心分析的石英含量(小数)；V长石C为岩心分析的长石含量(小数)；V石英与长石C为岩心分析的石英与长石总含量(小数)。
优选的，所述步骤S4中确定数值最小PE基线段的石英与长石骨架PE值(PE石英与长石骨架)的方法为：在该平台段没有岩心矿物含量分析数据的情况下，默认该PE基线对应的骨架颗粒全为石英，则该段的PE石英与长石骨架＝1.81B/e。
优选的，所述步骤S7中确定粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)的方法为：当研究区有岩心矿物含量分析数据的情况下(少量岩心分析数据即可)，参考表1中的各种主要矿物骨架PE值，利用以下公式分别计算粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)：
式中：PE粘土骨架为粘土骨架的PE值(B/e)；V伊利石C为岩心分析的伊利石含量(小数)；V蒙脱石C为岩心分析的蒙脱石含量(小数)；V高岭石C为岩心分析的高岭石含量(小数)；V绿泥石C为岩心分析的绿泥石含量(小数)；V粘土C为岩心分析的粘土含量(即岩心分析伊利石、蒙脱石、高岭石、绿泥石含量之和)(小数)。
式中：PE钙质胶结物骨架为钙质胶结物骨架的PE值(B/e)；V方解石C为岩心分析的方解石含量(小数)；V白云石C为岩心分析的白云石含量(小数)；V铁白云石C为岩心分析的铁白云石含量(小数)；V钙质胶结物C为岩心分析的钙质胶结物含量(即岩心分析方解石、白云石、铁白云石含量之和)(小数)。
优选的，所述步骤S7中确定粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)的方法为：当研究区没有岩心矿物含量分析数据的情况下，粘土骨架PE值取伊利石、蒙脱石、高岭石、绿泥石骨架PE值的平均值(PE粘土骨架＝3.41B/e)，钙质胶结物骨架PE值取方解石、白云石、铁白云石骨架PE值的平均值(PE钙质胶结物骨架＝5.85B/e)，或者取研究区的经验值.
优选的，所述步骤S8中利用以下公式计算扩径校正后的PE曲线，即PE校正2曲线：
PE校正2＝PE校正1-a×(CAL-BIT)
式中：PE校正2为重晶石钻井液泥饼校正和扩径校正后的PE曲线(B/e)。
(三)有益效果
本发明提供了基于光电吸收截面指数的地层钙质胶结物含量计算方法。与现有技术相比具备以下有益效果：该基于光电吸收截面指数的地层钙质胶结物含量计算方法，通过基于光电吸收截面指数(PE)原始曲线及其他测井数据，通过建立分段式PE基线，确定岩石中石英与长石骨架、粘土骨架和钙质胶结物的骨架PE值，对PE曲线进行重晶石钻井液泥饼校正和扩径校正，最终实现全井段的钙质胶结物含量计算。本发明计算的钙质胶结物含量与岩心分析钙质胶结物含量计算结果吻合很好，证明了本发明的准确性(图2)。本发明对于确定砂岩储层的有效性、提高砂岩储层的勘探开发效率具有重要现实意义。
附图说明
图1为本发明基于光电吸收截面指数的地层钙质胶结物含量计算方法流程图；
图2为本发明PE曲线校正与钙质胶结物含量计算结果图；
图3为本发明深度-PE原始数据绘制的连线图；
图4为本发明深度-PE校正1数据绘制的连线图；
图5为本发明深度-PE校正1与CAL数据绘制的连线图；
图6为本发明PE校正1-(CAL-BIT)交会图；
图7为本发明深度-PE校正2与CAL数据绘制的连线图；
图8为本发明第①段PE曲线校正与钙质胶结物含量计算结果图；
图9为本发明第②段PE曲线校正与钙质胶结物含量计算结果图；
图10为本发明第③段PE曲线校正与钙质胶结物含量计算结果图；
图11为本发明第④段PE曲线校正与钙质胶结物含量计算结果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1-11，本发明实施例提供一种技术方案：基于光电吸收截面指数的地层钙质胶结物含量计算方法，以我国某油田的A井为例，说明本发明的具体实施方式，具体步骤如下：、
(1)数据准备
准备好A井包含PE原始曲线的测井数据，以及少量研究区的岩心矿物含量分析数据。
(2)PE原始曲线响应分析
PE原始数据由6部分组成：重晶石钻井液泥饼的响应、扩径的响应、盐水的响应、粘土的响应、石英与长石矿物的响应、钙质胶结物的响应，PE原始数据可表示为以下公式：
PE原始＝PE重晶石钻井液泥饼+PE扩径+PE盐水+PE粘土+PE石英与长石+PE钙质胶结物
(3)分段式PE基线确定
利用测井数据，以深度为横坐标并以PE原始数据为纵坐标绘制连线图(图1)。根据PE原始曲线形态将其分为4段(第①段，809～1453m；第②段，1453～2335m；第③段，2335～2670m；第④段，2670～3459m)，在每一段选择3-5个PE值最低的纯砂岩层，取其PE平均值确定每段的PE基线值(图3)。这4段的PE基线值确定如下：第①段，PE基线＝3.60B/e；第②段，PE基线＝4.00B/e；第③段，PE基线＝3.10B/e；第④段，PE基线＝3.40B/e。
(4)最小PE基线段的石英与长石骨架PE值确定
PE基线由2部分构成：重晶石钻井液泥饼的PE响应、石英与长石的PE响应，PE基线即可表示为以下公式：
PE基线＝PE重晶石钻井液泥饼+PE石英与长石
在所划分的多段PE基线中，第③段PE基线处的PE基线值最小，因此第③段PE基线所对应的地层具有全井段最高的石英矿物含量。本井在第③段PE基线处有纯砂岩岩心，岩心分析表明，
因此优先利用以下公式计算第③段PE基线处的PE石英与长石骨架：
(5)重晶石钻井液泥饼校正
针对数值最小的PE基线段(第③段PE基线)，读取其PE基线值(PE基线min＝3.10B/e)，利用以下公式计算PE重晶石钻井液泥饼：
PE重晶石钻井液泥饼＝PE基线min-PE石英与长石骨架＝3.10-1.86＝1.24B/e
进一步，利用以下公式计算得到重晶石钻井液泥饼校正后的PE曲线，即PE校正1曲线(图4)：
PE校正1＝PE原始-PE重晶石钻井液泥饼。
(6)其他PE基线段的石英与长石骨架PE值确定
在不同平台段，PE重晶石钻井液泥饼的值保持不变。但在不同平台段，石英与长石的比例不同，则其PE石英与长石骨架也不同。在数值最小的PE基线段(第③段PE基线)计算出PE重晶石钻井液泥饼的基础上(PE重晶石钻井液泥饼＝1.24B/e)，利用以下公式计算其他PE基线段(第①段、第②段、第④段PE基线)的PE石英与长石骨架：
PE石英与长石骨架＝PE基线-1.24
所有PE基线段的石英与长石骨架PE值(PE石英与长石骨架)计算结果为：第①段PE基线，PE石英与长石响应＝2.36B/e；第②段PE基线，PE石英与长石响应＝2.76B/e；第③段PE基线，PE石英与长石响应＝1.86B/e；第④段PE基线，PE石英与长石响应＝2.16B/e。
(7)粘土骨架PE值和钙质胶结物骨架PE值确定
研究区有少量岩心矿物含量分析数据
参考表1中的各种主要矿物骨架PE值，优先利用以下公式分别计算粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物骨架)：
粘土骨架PE值和钙质胶结物骨架PE值计算结果为：PE粘土骨架＝3.47B/e，PE钙质胶结物骨架＝7.23B/e。
(8)扩径校正
找到典型扩径段(例如3430.8米-3453.6米)，该典型扩径段的PE校正1明显高于正常井径段的PE校正1(图5)。
选择该段数据，以CAL-BIT为横坐标，并以PE校正1为纵坐标，绘制交会图(图6)并获得线性拟合关系式：
PE校正1＝2.0916×(CAL-BIT)+4.3278
进一步，利用以下公式计算扩径校正后的PE曲线，即PE校正2曲线(图7。
PE校正2＝PE校正1-2.0916×(CAL-BIT)。
(9)钙质胶结物含量计算
PE曲线经过重晶石钻井液泥饼校正和扩径校正后得到PE校正2曲线，以PE校正2曲线为基础，通过以下公式即可计算钙质胶结物含量(V钙质胶结物)：
上式中，POR、V粘土可通过其他测井曲线计算得到；PE盐水骨架＝1.64B/e(表1)；PE石英与长石骨架、PE粘土骨架、PE钙质胶结物骨架已在前面步骤中计算出来，因此钙质胶结物含量(V钙质胶结物)的计算公式可进一步表示为：
第①段PE基线(809～1453m)：
第②段PE基线(1453～2335m)：
第③段PE基线(2335～2670m)：
第④段PE基线(2670～3459m)：
根据上述公式，计算得到全井段的钙质胶结物含量(V钙质胶结物)，图8-图11分别展示了第①段PE基线-第④段PE基线中的部分深度段的钙质胶结物含量(V钙质胶结物)计算结果。本发明计算的钙质胶结物含量(V钙质胶结物)与岩心分析钙质胶结物含量计算结果吻合很好，证明了本发明的准确性。
同时本说明书中未作详细描述的内容均属于本领域技术人员公知的现有技术。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。
Claims (6)
1.基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：具体包括以下步骤：
S1、数据准备：准备好研究区包含PE原始曲线的测井数据；
S2、PE原始曲线响应分析：PE原始数据由6部分组成，包括重晶石钻井液泥饼的响应、扩径的响应、盐水的响应、粘土的响应、石英与长石矿物的响应以及钙质胶结物的响应；
S3、分段式PE基线确定：利用测井数据，以深度为横坐标并以PE原始数据为纵坐标绘制连线图，根据PE原始曲线的形态将其分为多个平台段，在每一段选择若干个PE值最低的纯砂岩层，取其PE平均值确定每段的PE基线值；
S4、最小PE基线段的石英与长石骨架PE值确定：PE基线由两部分构成：重晶石钻井液泥饼的PE响应和石英与长石的PE响应，PE基线即表示为以下公式：
PE基线＝PE重晶石钻井液泥饼+PE石英与长石
式中：PE基线为平台段内几个PE值最低纯砂岩层的PE平均值(B/e)，在所划分的多段PE基线中，选择数值最小的那段PE基线，由于石英骨架的PE值低于长石骨架的PE值，因此，数值最小的那段PE基线所对应的地层具有全井段最高的石英矿物含量；
确定数值最小PE基线段的石英与长石骨架PE值(PE石英与长石骨架)的方法为：在该平台段有岩心矿物含量分析数据的情况下，优先利用以下公式计算该段的PE石英与长石骨架：
式中：PE石英与长石骨架为平台段内石英与长石骨架的PE值(B/e)；V石英C为岩心分析的石英含量；V长石C为岩心分析的长石含量；V石英与长石C为岩心分析的石英与长石总含量；
若在该平台段没有岩心矿物含量分析数据的情况下，默认该PE基线对应的骨架颗粒全为石英，则该段的PE石英与长石骨架＝1.81B/e；
S5、重晶石钻井液泥饼校正：针对数值最小的PE基线段，读取其PE基线值(PE基线min)，并令该段的PE石英与长石＝PE石英与长石骨架，将这两个参数代入步骤S4中的PE基线表示公式中，可得PE重晶石钻井液泥饼：
PE重晶石钻井液泥饼＝PE基线min-PE石英与长石骨架；
利用以下公式计算得到重晶石钻井液泥饼校正后的PE曲线，即PE校正1曲线：
PE校正1＝PE原始-PE重晶石钻井液泥饼
式中：PE校正1为重晶石钻井液泥饼校正后的PE曲线(B/e)；
S6、其他PE基线段的石英与长石骨架PE值确定：在不同平台段，PE重晶石钻井液泥饼的值保持不变，但在不同平台段，石英与长石的比例不同，则其PE石英与长石骨架也不同，在数值最小的PE基线段计算出PE重晶石钻井液泥饼的基础上，利用以下公式计算其他段的PE石英与长石骨架：
PE石英与长石骨架＝PE基线-PE重晶石钻井液泥饼；
S7、粘土骨架PE值和钙质胶结物骨架PE值确定；
S8、扩径校正：找到典型扩径段，典型扩径段的PE校正1明显高于正常井径段的PE校正1，选择该段数据，以CAL-BIT为横坐标，并以PE校正1为纵坐标，绘制交会图并获得线性拟合关系式：
PE校正1＝a×(CAL-BIT)+b
式中：CAL为井径(in)；BIT为钻头直径(in)；a、b为拟合系数；
S9、钙质胶结物含量计算：PE曲线经过重晶石钻井液泥饼校正和扩径校正后得到PE校正2曲线，PE校正2主要由4部分构成：盐水的响应、粘土的响应、石英与长石的响应、钙质胶结物的响应。PE校正2表示为下式：
PE校正2＝PE盐水+PE粘土+PE石英与长石+PE钙质胶结物
PE校正2可进一步表示为下式：
PE校正2＝POR×PE盐水骨架+V粘土×PE粘土骨架+V石英与长石×PE石英与长石骨架+V钙质胶结物×PE钙质胶结物骨架
式中：POR为孔隙度，可通过其他测井曲线计算得到；V粘土为粘土含量，通过其他测井曲线计算得到；V石英与长石为石英与长石含量；V钙质胶结物为钙质胶结物含量；PE盐水骨架为孔隙中盐水骨架的PE值。
在岩石中，孔隙、粘土、石英与长石、钙质胶结物的总体积为1，即：
POR+V粘土+V石英与长石+V钙质胶结物＝1
将上式代入PE校正2可进一步表示的公式，得到钙质胶结物含量(V钙质胶结物)的计算公式：
2.根据权利要求1所述的基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：所述步骤S1中PE原始数据表示为以下公式：
PE原始＝PE重晶石钻井液泥饼+PE扩径+PE盐水+PE粘土+PE石英与长石+PE钙质胶结物
式中：PE原始为PE原始数据(B/e)、PE重晶石钻井液泥饼为重晶石钻井液泥饼对PE的响应(B/e)、PE扩径为扩径对PE的响应(B/e)、PE盐水为孔隙中盐水对PE的响应(B/e)、PE粘土为粘土对PE的响应(B/e)、PE石英与长石为石英与长石对PE的响应(B/e)以及PE钙质胶结物为钙质胶结物对PE的响应(B/e)。
3.根据权利要求1所述的基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：所述步骤S2中在每一段选择3-5个PE值最低的纯砂岩层。
4.根据权利要求1所述的基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：所述步骤S7中确定粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)的方法为：当研究区有岩心矿物含量分析数据的情况下，利用以下公式分别计算粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)：
式中：PE粘土骨架为粘土骨架的PE值(B/e)；V伊利石C为岩心分析的伊利石含量；V蒙脱石C为岩心分析的蒙脱石含量；V高岭石C为岩心分析的高岭石含量；V绿泥石C为岩心分析的绿泥石含量；V粘土C为岩心分析的粘土含量。
式中：PE钙质胶结物骨架为钙质胶结物骨架的PE值(B/e)；V方解石C为岩心分析的方解石含量；V白云石C为岩心分析的白云石含量；V铁白云石C为岩心分析的铁白云石含量；V钙质胶结物C为岩心分析的钙质胶结物含量。
5.根据权利要求1所述的基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：所述步骤S7中确定粘土骨架PE值(PE粘土骨架)和钙质胶结物骨架PE值(PE钙质胶结物)的方法为：当研究区没有岩心矿物含量分析数据的情况下，粘土骨架PE值取伊利石、蒙脱石、高岭石、绿泥石骨架PE值的平均值(PE粘土骨架＝3.41B/e)，钙质胶结物骨架PE值取方解石、白云石、铁白云石骨架PE值的平均值(PE钙质胶结物骨架＝5.85B/e)，或者取研究区的经验值。
6.根据权利要求1所述的基于光电吸收截面指数的地层钙质胶结物含量计算方法，其特征在于：所述步骤S8中利用以下公式计算扩径校正后的PE曲线，即PE校正2曲线：
PE校正2＝PE校正1-a×(CAL-BIT)
式中：PE校正2为重晶石钻井液泥饼校正和扩径校正后的PE曲线(B/e)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211076515.0A CN115434698B (zh) | 2022-09-01 | 2022-09-01 | 基于光电吸收截面指数的地层钙质胶结物含量计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211076515.0A CN115434698B (zh) | 2022-09-01 | 2022-09-01 | 基于光电吸收截面指数的地层钙质胶结物含量计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115434698A true CN115434698A (zh) | 2022-12-06 |
| CN115434698B CN115434698B (zh) | 2024-05-28 |
Family
ID=84246698
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211076515.0A Active CN115434698B (zh) | 2022-09-01 | 2022-09-01 | 基于光电吸收截面指数的地层钙质胶结物含量计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115434698B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116502557A (zh) * | 2023-05-04 | 2023-07-28 | 西南石油大学 | 基于网格分析法的储层压裂后可动流体饱和度计算方法 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6150655A (en) * | 1998-03-06 | 2000-11-21 | Computalog Research, Inc. | Inferential measurement of photoelectric absorption cross-section of geologic formations from neutron-induced, gamma-ray spectroscopy |
| CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
| CN104564037A (zh) * | 2013-10-10 | 2015-04-29 | 中国石油天然气股份有限公司 | 一种页岩气储层脆性矿物含量测井计算方法 |
| CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
| CN105425309A (zh) * | 2015-11-17 | 2016-03-23 | 山东科技大学 | 确定高伽马砂岩凝灰质含量的方法 |
| CN105443122A (zh) * | 2015-12-28 | 2016-03-30 | 中国石油天然气股份有限公司 | 一种测井解释模型的处理方法及装置 |
| CN109061752A (zh) * | 2018-06-26 | 2018-12-21 | 西南石油大学 | 一种含灰质地层的电阻率曲线校正方法 |
| CN109424360A (zh) * | 2017-08-28 | 2019-03-05 | 中国石油天然气股份有限公司 | 确定地层中矿物含量的方法及装置 |
| CN111797546A (zh) * | 2020-07-29 | 2020-10-20 | 重庆科技学院 | 一种页岩油气储层矿物组分模型最优化反演方法 |
| CN112836163A (zh) * | 2019-11-25 | 2021-05-25 | 中国石油大学(北京) | 成岩相测井定量划分方法、装置和设备 |
| CN113031102A (zh) * | 2019-12-09 | 2021-06-25 | 中国石油大学(北京) | 碳酸盐岩储层岩性岩相测井识别方法、装置及存储介质 |
-
2022
- 2022-09-01 CN CN202211076515.0A patent/CN115434698B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6150655A (en) * | 1998-03-06 | 2000-11-21 | Computalog Research, Inc. | Inferential measurement of photoelectric absorption cross-section of geologic formations from neutron-induced, gamma-ray spectroscopy |
| CN103375166A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种确定地层中多种矿物组分含量的方法 |
| CN104564037A (zh) * | 2013-10-10 | 2015-04-29 | 中国石油天然气股份有限公司 | 一种页岩气储层脆性矿物含量测井计算方法 |
| CN105298476A (zh) * | 2014-07-04 | 2016-02-03 | 中国石油化工股份有限公司 | 一种识别页岩地层矿物体积含量的方法 |
| CN105425309A (zh) * | 2015-11-17 | 2016-03-23 | 山东科技大学 | 确定高伽马砂岩凝灰质含量的方法 |
| CN105443122A (zh) * | 2015-12-28 | 2016-03-30 | 中国石油天然气股份有限公司 | 一种测井解释模型的处理方法及装置 |
| CN109424360A (zh) * | 2017-08-28 | 2019-03-05 | 中国石油天然气股份有限公司 | 确定地层中矿物含量的方法及装置 |
| CN109061752A (zh) * | 2018-06-26 | 2018-12-21 | 西南石油大学 | 一种含灰质地层的电阻率曲线校正方法 |
| CN112836163A (zh) * | 2019-11-25 | 2021-05-25 | 中国石油大学(北京) | 成岩相测井定量划分方法、装置和设备 |
| CN113031102A (zh) * | 2019-12-09 | 2021-06-25 | 中国石油大学(北京) | 碳酸盐岩储层岩性岩相测井识别方法、装置及存储介质 |
| CN111797546A (zh) * | 2020-07-29 | 2020-10-20 | 重庆科技学院 | 一种页岩油气储层矿物组分模型最优化反演方法 |
Non-Patent Citations (3)
| Title |
|---|
| FENG WU ET AL.: "Numerical Simulation of the Influence of Pore Structure on Resistivity, Formation Factor and Cementation Index in Tight Sandstone", 《ACTA GEOLOGICA SINICA》, 28 February 2020 (2020-02-28), pages 290 * |
| 韩浩东等: "《复杂碎屑岩矿物组分最优化测井定量分析》", 新疆石油地质, 30 June 2016 (2016-06-30), pages 726 - 731 * |
| 魏吉凯等: "华庆油田东部长6储层沉积微相及其测井评价研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》, 15 March 2013 (2013-03-15), pages 1 - 84 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116502557A (zh) * | 2023-05-04 | 2023-07-28 | 西南石油大学 | 基于网格分析法的储层压裂后可动流体饱和度计算方法 |
| CN116502557B (zh) * | 2023-05-04 | 2024-01-30 | 西南石油大学 | 基于网格分析法的储层压裂后可动流体饱和度计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115434698B (zh) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Morris et al. | Using log-derived values of water saturation and porosity | |
| CN108252709B (zh) | 一种致密砂岩油藏的油水性质识别方法及系统 | |
| CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
| CN107346455A (zh) | 一种用于识别页岩气产能的方法 | |
| CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
| CN115434698A (zh) | 基于光电吸收截面指数的地层钙质胶结物含量计算方法 | |
| Ehrenberg | Factors controlling porosity in Upper Carboniferous–Lower Permian carbonate strata of the Barents Sea | |
| CN112145165A (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
| CN110593855A (zh) | 一种测井电阻率的校正方法、油层的识别方法 | |
| CN110208874B (zh) | 一种致密砂岩油藏有效储层识别方法 | |
| CN110688781B (zh) | 一种低渗非均质气藏储层测井解释方法 | |
| CN110410068B (zh) | 一种确定地层破裂压力梯度的测井方法 | |
| CN112814667A (zh) | 基于热中子测井计数率比值评价致密气层含水性的方法 | |
| CN113720745A (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
| CN112228050A (zh) | 一种定量评价致密油储层宏观非均质性的方法及其应用 | |
| CN112346147A (zh) | 一种基于中子密度孔隙度差的储层评价方法 | |
| CN114692360A (zh) | 一种用于碳酸盐岩储层的测井评价方法及系统 | |
| CN113622904A (zh) | 用元素俘获谱测井氯离子含量识别测井流体类型的方法 | |
| Rathmell et al. | Application of Low Invasion Coring and Outcrop Studies to Reservoir Development Planning for the Villano Field | |
| CN116774279B (zh) | 基于气水平衡岩性组合单元的页岩含水饱和度预测方法 | |
| CN110469321B (zh) | 一种确定地层破裂压力梯度的录井方法 | |
| CN117784278B (zh) | 一种致密砂岩气甜点的预测方法和预测系统 | |
| Jia et al. | Study On Thin Differential Interpretation Method Of 0.2 m High Resolution Logging Series in Sabei Development Zone | |
| Khatimah et al. | Petrophysical Analysis Using Deterministic and Probabilistic Methods in Low Resistivity Zone of Gumai and Talang Akar Formations in Jambi Subbasin South Sumatra | |
| CN116291373A (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 |