CN109779611A - 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 - Google Patents
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 Download PDFInfo
- Publication number
- CN109779611A CN109779611A CN201910129076.7A CN201910129076A CN109779611A CN 109779611 A CN109779611 A CN 109779611A CN 201910129076 A CN201910129076 A CN 201910129076A CN 109779611 A CN109779611 A CN 109779611A
- Authority
- CN
- China
- Prior art keywords
- cement
- cement sheath
- thickness
- amendment
- counting rate
- 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.)
- Withdrawn
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 54
- 239000003208 petroleum Substances 0.000 title claims abstract description 31
- 239000004568 cement Substances 0.000 claims abstract description 157
- 238000012937 correction Methods 0.000 claims abstract description 22
- 238000010521 absorption reaction Methods 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims abstract description 11
- 238000013459 approach Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims description 18
- 230000000694 effects Effects 0.000 claims description 18
- 230000002285 radioactive effect Effects 0.000 claims description 13
- 230000002159 abnormal effect Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000000149 penetrating effect Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 8
- 238000004519 manufacturing process Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 3
- 229910052770 Uranium Inorganic materials 0.000 description 10
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 10
- 230000004044 response Effects 0.000 description 9
- 230000005251 gamma ray Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000012417 linear regression Methods 0.000 description 3
- 239000003129 oil well Substances 0.000 description 3
- 239000002002 slurry Substances 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000001568 sexual effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000726221 Gemma Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，本发明通过制作多个不同厚度不同密度的水泥环模型孔，模拟井下情况，建立模型，利用不同模型水泥环吸收系数值进行拟合，最终找出固井水泥环对伽玛测井计数率的影响系数，并进行修正论证实验，使修正后的数据无限接近地层真实值。文中结合内蒙某油田区块实际资料，通过水泥环实验测井修正系数，对资料进行修正，探讨了水泥环密度和厚度与计数率的关系，进一步验证了方法技术的可行性，从而为指导下一步石油老井复测结果提供修正依据。
Description
技术领域
本发明涉及地球物理伽玛测井，尤其涉及石油老井复测的伽马测井结果修正方法。
背景技术
近年来，寻找地浸砂岩型铀矿床己成为铀矿工作的主攻方向，并取得重大进展。地球物理伽玛测井是铀矿地质勘查的一个重要组成部分，其目的是定量确定铀矿(层)体的空间位置、含量及其厚度，为地质研究、钻探施工、储量计算等工作提供科学依据，同时为提交铀矿藏的储量报告及后期的开发利用提供有力证据。
运用铀矿勘查技术在内蒙某油田区块进行选区选井研究，为了进一步查明工作区成矿潜力，充分收集目的层放射性异常信息，同时为了弥补石油井浅部资料的缺失，对工作区多个凹陷内石油钻孔进行老井复测工作，并进行解释和总结。在已完成的伽马测井中发现了新的问题：即自然伽马测井曲线的幅度不仅与地层的放射性有关，而且还受井眼条件(井径、泥浆比重、套管、水泥环等参数)的影响。泥浆、套管、水泥环均吸收伽马射线会使自然伽马测井值降低。
在做定性解释时，如果井内泥浆稳定，则整个曲线的相对趋势反映地层性质，可不做修正。在大井眼和套管井中，定量解释自然伽马资料时，要进行必要的修正。由中国核工业地质局《测井资料自动化处理解释系统》软件中γ测井解释系统软件模块可直接对井径、泥浆比重、套管等参数进行校正。老井复测时所面对的石油井，除了考虑铁套管吸收因素，还需要考虑固井后，铁套管外壁与岩层之间固井水泥环的影响因素，这种影响因素目前在工作中没法修正，无法最大程度的接近真实数据。针对此种情况，经过多次研究论证，制定了详细的修正实验方案，考虑利用模型响应拟合修正系数，从而为指导下一步石油老井复测的伽马测井结果提供修正依据。
本着“循序渐进”的原则，本次工作通过制作多个不同厚度不同密度的水泥环，模拟井下情况，建立模型，研究水泥环影响因素的测井响应，最终找出固井水泥环对测井的影响系数，并进行修正论证实验，使修正后的数据无限接近地层真实值。
发明内容
本发明的目的是提供了石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，本发明通过制作多个不同厚度不同密度的水泥环模型孔，模拟井下情况，建立模型，利用不同模型水泥环吸收系数值进行拟合，最终找出固井水泥环对伽玛测井计数率的影响系数，并进行修正论证实验，使修正后的数据无限接近地层真实值。
本发明的技术方案是：一种石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：基于放射性测井的康普顿效应，建立多个不同厚度不同密度的水泥环模型，模拟井下情况，利用不同模型水泥环吸收系数值进行拟合，通过模拟计算得出放射性计数率和水泥环密度以及水泥环厚度的函数关系，分析可知放射性计数率和水泥环密度是线性关系，和水泥环厚度是三阶多项式非线性关系，并推导出了计数率与水泥环属性参数的综合函数关系，最终找出固井水泥环对伽玛测井计数率的影响系数，得出具体修正系数，并进行修正论证实验，使修正后的数据无限接近地层真实值，对过去的和即将测量的石油老井进行重新解释修正，还原真实地层放射性效应。
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：分析了不同厚度、不同密度水泥环对放射性计数率的影响，通过对实际资料处理解释，得出水泥环对计数率的影响最高为30％。
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：建立的水泥环模型，从井眼开始依次是套管，水泥环，地层；γ源为137Cs点状源，仪器推靠井轴测量。
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：对于任意点，只要已知水泥环密度、水泥环厚度、水泥返深，就可以求出各个参数对计数率修正值，其和就是总体修正值；
最终水泥环综合修正公式为：
Y＝100*Y1/(100-0.0007h(s) 3+0.0519h(s) 2-0.9652h(s))+100*Y1/(190.461-64.344*ρ(s))-Y1
(1-13)
Y1为任意一点实际测量计数率，h(s)为固井水泥环厚度，ρ(s)为固井水泥环密度，Y为总体修正值。
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：对于某一口井，只要已知水泥环密度值和实测测量值，就可以对其异常层段计数率进行修正。
Y2＝100*Y1/(190.461-64.344*ρ(s)) (1-5)
公式(1-5)是水泥环密度修正公式，Y1为任意一点实际测量计数率，Y2为希望得到的修正值，ρ(s)为固井水泥环密度。
石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：
计数率和水泥环厚度表达式如下：
Y＝4087.716-39.457*h(s)+2.120*h(s) 2-0.027*h(s) 3 (1-7)
其中Y1为任意一点实际测量计数率，h(s)为固井水泥环厚度。
附图说明
图1是系数修正实验流程框架图。
图2是水泥环计算模型示意图。
图3是相同水泥环厚度不同水泥环密度计数率关系图。
图4是水泥环密度与计数率关系散点图。
图5是相同水泥环密度不同水泥环厚度计数率关系图。
图6是水泥环厚度与计数率关系散点图。
图7是水泥环厚度定点测量曲线估算图。
图8是伽玛复测系数修正程序流程图。
图9是ZKM修正前后散点对比图(340米-370米)。
图10是ZKM修正前后散点对比图(438米-447米)。
图11是ZKM修正前后和验证孔曲线对比图。
具体实施方式
以下结合附图和实施例对本发明内容进行详细描述。
一、系数修正实验理论基础
如图1所示，首先建立水泥环模型，方法是基于放射性测井的康普顿效应。固井水泥环厚度及密度是如何影响放射性强度需要研究各个参数的测井响应。通过各个测井响应数值模拟方法可以全面认识水泥环与放射性关系特征，把握放射性衰变的规律，从而为找出固井水泥环对测井的影响系数，并进行修正论证实验，使修正后的数据无限接近地层真实值。
放射性元素核衰变放出的伽马射线，能量一般在0.5MeV到5.3MeV之间，当伽马射线能量介于0.25MeV到2.0MeV之间时，能量较高伽马射线与物质的原子核外电子碰撞时，一部分能量转交给电子，使之脱离原子电子壳层而飞出，同时伽马射线改变自己运动方向，继续与其它电子相撞。每碰撞一次，能量损失一部分，并改变其运动方向，形成康普顿效应。伽马射线经多次碰撞之后，能量不断降低，最后以光电效应结束。本次实验使用的伽马源与岩石作用时，主要产生康普顿效应，并散射伽马射线。
一个物体的质量和长度是具有固定数值的，只要量具和测量仪器质量可靠、环境条件不变、工作仔细、多次重复测量，总可以测得相同的数值，放射性测量却完全不同。不管测量仪器多么稳定精确，环境条件控制极好，工作非常细心，每次测量仍然得不到完全相同的结果，甚至相差很大。因为一个放射性源含有很多不稳定的原子核，每个原子核只有衰变时才放出射线被探测器接收，每个核衰变是完全独立的，完全不与其他核衰变相关。哪个核先衰变，哪个核后衰变，没有任何规定的顺序，纯属随机。每个单位时间的核衰变数，不可能完全相同。然而虽然有差别，但他们是有规律的，放射性核素衰变符合高斯分布，对放射性进行无限多次测量得到的平均值，才是最准确的期望值。基于上述规律，本次水泥环实验定点测量数据不少于300个，去掉最大值、最小值以及野值求取平均值作为期望值。
二、水泥环模型的建立
在实际生产中，井下水泥环的状态，是处于套管与外壁围岩之间，水泥浆凝固后，套管、水泥环及井壁围岩为一个组合弹性体。考虑到实际情况，建立的模型，从井眼开始依次是套管，水泥环，地层。γ源为137Cs点状源，仪器推靠井轴测量。图2为水泥环计算模型示意图。
(1)水泥环密度模型
水泥环密度模型目的是研究相同套管类型、相同水泥环厚度、不同水泥环密度对计数率的影响，采用定点测量和连续测量两种测量方式，为了研究计数率与水泥环密度的关系，本次实验设计六种模型参数进行实验。采用密度定向源，源距15厘米，统一放置井深为1米的位置，模型一是裸眼孔，模型一和模型二均没有制作水泥环，模型三～模型六铁套管规格一致，水泥环厚度均为26.15毫米，密度有明显差异，具体见下表1。
表1水泥环密度模型参数表
(2)水泥环厚度模型
水泥环厚度模型目的是研究相同套管类型、相同水泥环密度、不同水泥环厚度对计数率的影响，采用定点测量和连续测量两种测量方式，放射源贴水泥环外壁和源距为15cm分别测量，水泥环厚度设计为6.15mm、20.00mm、26.15mm、50.15mm等四组数据，厚度利用PVC管进行控制。计算得到放射性计数率与水泥环厚度测井响应关系，具体见下表2。
表2水泥环厚度模型参数表
| 序号 | 模型井名 | 套管规格 | 水泥环密度(g/cm<sup>3</sup>) | 水泥环厚度(mm) |
| 模型一 | ZK-PVC | 无 | 无 | 无 |
| 模型二 | ZK-PVC-TG | N80-139.7-7.72 | 无 | 无 |
| 模型三 | ZKSH1 | N80-139.7-7.72 | 1.82 | 6.15 |
| 模型四 | ZKSH2 | N80-139.7-7.72 | 1.82 | 50.15 |
| 模型五 | ZKSH3 | N80-139.7-7.72 | 1.82 | 20.00 |
| 模型六 | ZKSH4 | N80-139.7-7.72 | 1.82 | 26.15 |
(3)结合前面各种模型关系确定计数率与变量关系
分析对于油田所有井，忽略井内流体的影响因素，需要确定水泥环密度ρ(s)、水泥环厚度h(s)；其中，主要研究计数率与水泥环密度、水泥环厚度的关系函数，这里，将计数率与变量的关系表示为：N∝f(ρ(s),h(s))；结合前面各种响应分析，利用SPSS软件绘制散点图，研究自变量和因变量之间的相关系数。限定另一个自变量，在此基础上进行相关性分析，看是否呈线性正相关，还是别的关系。进一步对不同模型数据进行单自变量线性(非线性)回归分析，建立回归方程。
三、实验计算实例分析
本实验投入北京中地英捷公司生产的FD-3019闪烁伽玛测井探管，仪器性能良好，实验开展前均在国防科技工业1313二级计量站进行了标定。选择一个放射性相对平稳而不受人文干扰场影响的地区，面积大约30平方米，一共建立模型孔12个，深度统一为1.5米，孔径统一为296mm，孔与孔之间距离1米以上。在测量之前，将所有孔的本底均加以测量，以防有的孔自身有放射性引起误差。
3.1水泥环密度模型分析
对于连续测量，对表1水泥环密度模型中各个模型进行测井响应，以测井深度为横坐标，各组计数率为纵坐标，绘制出了不同模型计数率随与密度关系散点图，见图3。
由图3可见，没有铁套管和水泥环时，计数率峰值最高，加入铁套管后计数率峰值次之，水泥环厚度统一时，随着水泥环密度的增加，计数率呈现出衰减的趋势，且峰值衰减表现出一定的规律。
对于定点测量，利用各个模型峰值大小差异做出散点图，分析随水泥环密度的增加，计数率的变化情况，得出的散点图如下图4所示：
由图4可见，当水泥环厚度一定时，3019探测器的计数率随着水泥环密度的增加呈现逐渐衰减的近线性关系，在这个比较明显的近线性趋势里，没有偏离很远的点。由于自变量水泥环密度与因变量计数率均为连续变量，考虑建立单自变量线性回归方程。假设计数率Y与水泥环密度ρ(s)有如下线性关系：Y＝a+b*ρ(s)；式中：a为常数项，b为Y对应于ρ(s)的回归系数，ρ(s)代表水泥环密度，单位为g/cm3，Y代表计数率值，单位为cps，利用SPSS软件进行单自变量线性回归模型，结果如下表3。
表3模型摘要表
由于调整的判定系数R2较接近1，德宾-沃森模型结果值为2.175，说明数据之间是互相独立的，因此认为该模型拟合优度较高。建立线性模型如下表4所示。
表4模型系数表
通过采用逐步法回归分析得到表4中自变量水泥环密度ρ(s)系数的概率P值远远小于0.05，差异显著，即对方程有用，所以最后得到的模型为：
Y＝42249.030-18718.301*ρ(s) (1-1)
其中，ρ(s)代表水泥环密度，单位为g/cm3。
实际上，假定水泥环厚度不发生变化，在遇到不同放射性强度的目的层时，水泥环对放射源的吸收值是随密度值变化而变化的，而实际测量值是原始值(希望得到的修正值)经过衰减后的值，假设水泥环吸收系数值为Q％，任意一点实际测量值为Y1，希望得到的修正值为Y2，则有：
(Y2-Y1)/Y2＝Q％ (1-2)
Y2＝Y1/(1-Q％) (1-3)
由上述模型可知，当仅有铁套管而没有水泥环的时候，按照理论上应该为29091cps；当加上一层厚度21.65mm密度1.63g/cm3水泥环时，计数率变成了24517cps；衰减了15.72％；当加上同一厚度，密度1.72g/cm3水泥环时，计数率变成了23659cps，衰减了18.67％；当加上同一厚度，密度1.82g/cm3水泥环时，计数率变成了21499cps，衰减了26.09％；当加上同一厚度，密度1.94g/cm3水泥环时，计数率变成了18866cps，衰减了35.15％。
尝试建立水泥环吸收系数值和水泥环密度的响应函数关系，则归一化后的函数关系为：
Q＝-90.461+64.344*ρ(s) (1-4)
将公式(1-4)代入(1-3)可知：
Y2＝100*Y1/(190.461-64.344*ρ(s)) (1-5)
公式(1-5)是水泥环密度修正公式，对于某一口井，只要已知水泥环密度值和实测测量值，就可以对其异常层段计数率进行修正。
3.2水泥环厚度模型分析
对于连续测量，对表2水泥环厚度模型中各个模型进行测井响应，采样间隔0.05米，以测井深度为横坐标，各组计数率为纵坐标，绘制出了不同模型计数率随水泥环厚度的变化图，见图5。
由图5可见，相同水泥环密度时，水泥环厚度为6.15mm时计数率峰值最小，水泥环厚度为20mm时次之，水泥环厚度为50.15mm时最大。随水泥环厚度增加，计数率没有衰减反而增加了。利用各个模型定点测量结果差异做出散点图，分析随水泥环厚度的增加，计数率的变化情况，得出的散点图如下图6所示。
由图6可见，当水泥环密度一定时，3019探测器的计数率随着水泥环厚度的增加呈现增加的非线性关系，利用SPSS回归模块里的曲线估算进行计算，结果如下图7和表5所示：
由图7知，拟合曲线为三次多项式时最接近，假设计数率Y与水泥环厚度h(s)有如下关系：
Y＝a+b1*h(s)+b2*h(s) 2+b3*h(s) 3 (1-6)
由表5可见，拟合曲线为三次非线性曲线时最接近，且R2为1，显著性为0，所以计数率和水泥环厚度表达式如下：
Y＝4087.716-39.457*h(s)+2.120*h(s) 2-0.027*h(s) 3 (1-7)
表5模型摘要和参数估算表
同样假设水泥环厚度吸收系数值为Qh％，任意一点实际测量值为Y1，希望得到的修正值为Y2，则：
(Y2-Y1)/Y2＝Qh％ (1-8)
所以Y2＝Y1/(1-Qh％) (1-9)
尝试建立水泥环厚度吸收系数值和水泥环厚度的响应函数关系，由上述模型可知，当没有水泥环的时候，按照理论上应该为4087.716cps；铁套管规格均相同前提下，当加上密度相同，厚度为6.15mm水泥环时，归一化处理后计数率变成了3919cps；衰减了4.127％；当加上密度相同，厚度为20.00mm水泥环时，归一化处理后计数率变成了3932cps；衰减了3.809％；当加上密度相同，厚度为26.15mm水泥环时，归一化处理后计数率变成了4026cps；衰减了1.509％；当加上密度相同，厚度为50.15mm水泥环时，归一化处理后计数率变成了4058cps；衰减了0.726％；同样，利用SPSS软件对水泥环厚度吸收系数和水泥环厚度进行拟合，建立水泥环厚度吸收系数值和铁套管厚度的响应函数关系：
Qh＝0.0007h(s) 3-0.0519h(s) 2+0.9652h(s)+2E-6 (1-10)
所以，对于任意点，只要已知套管厚度和得到实测测量值，就可以求出水泥环厚度吸收系数值。即Y2＝Y1/(1-Qh％)；Y2为希望得到的修正值，Y1为任意一点实际工作测量计数率。将公式(1-10)代入(1-9)即可得到最终修正公式为：
Y2＝100*Y1/(100-0.0007h(s) 3+0.0519h(s) 2-0.9652h(s)) (1-11)
3.3综合确定计数率与参数的关系
经过前面的研究发现，计数率与水泥环密度成线性关系，与水泥环厚度成三次非线性关系，即表达式应该为如下公式：
Y＝a+b*ρ(s)+c*h(s)+d*h(s) 2+e*h(s) 3 (1-12)
对于任意点，只要已知水泥环密度、水泥环厚度、水泥返深，对于我们实测计数率值，就可以求出各个参数对计数率修正值，其和就是总体修正值。Y1为任意一点实际测量计数率，h(s)为固井水泥环厚度，ρ(s)为固井水泥环密度。即可得到最终水泥环综合修正公式为：
Y＝100*Y1/(100-0.0007h(s) 3+0.0519h(h) 2-0.9652h(s))+100*Y1/(190.461-64.344*ρ(x))-Y1
(1-13)
四、实际资料应用效果分析
4.1计数率对比分析
以内蒙某油田区块某口已测矿化孔ZKM石油老井为研究对象，在油田作业队检泵时已获取该孔井深结构以及固井材料参数。石油老井ZKM套管及水泥浆数据如下表6所示：
表6套管及水泥浆数据统计表
利用模块对ZKM其中显示最好的两层矿化层计数率进行修正，取其中一部分数据作对比如表7所示：
表7实际资料修正结果对比表
修正前后异常层数据散点图对比图如下图9和图10所示.
由上述修正结果表以及修正前后散点对比图可见，由于水泥自身也具有放射性，所以厚度对放射性计数率的吸收是很小的，而密度占主要影响因素。修正之后计数率大约为之前的1.44倍，且在高异常背景区越明显，经计算，对于水泥环厚度38.15mm，密度1.88g/cm3时，计数率大约被衰减了30.68％。分析可见，水泥环对计数率的吸收影响较明显，其变化符合模型规律。
4.2解释结果对比分析
当验证孔距离石油老井非常接近时，可以将验证孔测井数据与石油老井复测数据进行对比，由于铀矿具有特殊性，即使在最理想状况下，也不会是同一位置，因此，数据上会有一定差距，下面将老井复测3019与修正后3019和验证孔3019数据进行对比分析。
在该石油复测钻孔ZKM(未修正之前为矿化孔)布设的验证ZKY(工业孔),主要为了验证系数修正结果，本次经过系数修正后为工业孔，且数据具有一定说服力，对比曲线与数据如下图11所示.
重点对矿化层进行研究分析，发现对于ZKM复测孔，修正前仅有4层矿化层和11层异常，修正后变成了4层矿化，1层工业层和14层异常层，异常层值整体得到了抬高，使得一些层位得到了真实合并，未修正前表外矿层厚度很大，修正之后表外矿层转化成了表内矿层，对比工业层修正前后解释结果如下表8所示：
表8修正前后与验证孔解释成果对比表
由上表可见：由于验证孔与原老井复测孔相距15m，井下异常存在差异，因此有一定的差异，但在修正之后，二者已经很接近了，经计算，修正后的工业层品位相差0.002％，平米铀量相差0.0094kg/m2；假设以验证孔为基准值，则品位误差约为5.433％，平米铀量误差为0.87％。这个误差结果在铀矿勘查中符合规范要求。上述分析可见，水泥环对计数率的吸收影响较明显，其变化符合模型规律，可以使修正后的数据无限接近地层真实值，为指导下一步石油老井伽马复测结果提供修正依据。
五、结论和意义
(1)通过实验制作多个水泥环模型孔，建立不同厚度不同密度水泥环的测井响应分析，并且建立了合理的解释模型，得出具体修正系数，利用验证孔和修正后解释结果进行对比证明了这种方法的准确性，同时可以对过去的和即将测的石油老井进行重新解释修正，还原真实地层放射性效应。
(2)本文基于大量实验数据的基础上，通过模拟计算得出放射性计数率和水泥环密度以及水泥环厚度的函数关系，分析可知和水泥环密度是线性关系，和水泥环厚度是三阶多项式非线性关系。在此基础上，推导出了计数率与水泥环属性参数的综合函数关系，有助于研究放射性及其衰变规律。
(3)对于油田任一采油井，钻井水泥固井参数均是严格按照设计执行的，水泥环厚度及密度也是确定的参数，本文分析了不同厚度、不同密度水泥环对放射性计数率的影响，通过对实际资料处理解释，得出水泥环对计数率的影响最高约为30％。
(4)得出修正系数能将定性解释最大限度的向定量-半定量解释靠近，工作更加科学实用。利用系数修正结果进行修正之后，能使得一些矿化孔经过修正解释后达到工业品位(岩性具体分析)，一些异常孔经过修正之后达到矿化标准，一些无异常孔经过修正之后达到异常孔标准，有利于我们准确地搜集工作区放射性资料。
(5)由于在全套管的石油井，别的测井曲线无意义。石油老井伽马复测在此背景下意义重大，系数修正有助于最大程度接近真实数据，为后期查明工作区成矿潜力提供有力证据。利用实验结果能将石油老井复测这种经济、可行、可靠、高效的方法在其他油田加以推广。
Claims (6)
1.石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：基于放射性测井的康普顿效应，建立多个不同厚度不同密度的水泥环模型，模拟井下情况，利用不同模型水泥环吸收系数值进行拟合，通过模拟计算得出放射性计数率和水泥环密度以及水泥环厚度的函数关系，分析可知放射性计数率和水泥环密度是线性关系，和水泥环厚度是三阶多项式非线性关系，并推导出了计数率与水泥环属性参数的综合函数关系，最终找出固井水泥环对伽玛测井计数率的影响系数，得出具体修正系数，并进行修正论证实验，使修正后的数据无限接近地层真实值，对过去的和即将测量的石油老井进行重新解释修正，还原真实地层放射性效应。
2.如权利要求1所述的石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：分析了不同厚度、不同密度水泥环对放射性计数率的影响，通过对实际资料处理解释，得出水泥环对计数率的影响最高为30％。
3.如权利要求1所述的石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：建立的水泥环模型，从井眼开始依次是套管，水泥环，地层；γ源为137Cs点状源，仪器推靠井轴测量。
4.如权利要求1所述的石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：对于任意点，只要已知水泥环密度、水泥环厚度、水泥返深，就可以求出各个参数对计数率修正值，其和就是总体修正值；
最终水泥环综合修正公式为：
Y＝100*Y1/(100-0.0007h(s) 3+0.0519h(s) 2-0.9652h(s))+100*/Y1(190.461-64.344*ρ(s))Y1 (1-13)
Y1为任意一点实际测量计数率，h(s)为固井水泥环厚度，ρ(s)为固井水泥环密度，Y为总体修正值。
5.如权利要求1所述的石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：对于某一口井，只要已知水泥环密度值和实测测量值，就可以对其异常层段计数率进行修正；
Y2＝100*Y1/(190.461-64.344*ρ(s)) (1-5)
公式(1-5)是水泥环密度修正公式，Y1为任意一点实际测量计数率，Y2为希望得到的修正值，ρ(s)为固井水泥环密度。
6.如权利要求1所述的石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用，其特征是：
计数率和水泥环厚度表达式如下：
Y＝4087.716-39.457*h(s)+2.120*h(s) 2-0.027*h(s) 3 (1-7)
其中Y1为任意一点实际测量计数率，h(s)为固井水泥环厚度。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910129076.7A CN109779611A (zh) | 2019-02-21 | 2019-02-21 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
| CN201910645297.XA CN110344817B (zh) | 2019-02-21 | 2019-07-17 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910129076.7A CN109779611A (zh) | 2019-02-21 | 2019-02-21 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109779611A true CN109779611A (zh) | 2019-05-21 |
Family
ID=66503539
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910129076.7A Withdrawn CN109779611A (zh) | 2019-02-21 | 2019-02-21 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
| CN201910645297.XA Active CN110344817B (zh) | 2019-02-21 | 2019-07-17 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910645297.XA Active CN110344817B (zh) | 2019-02-21 | 2019-07-17 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN109779611A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534283A (zh) * | 2021-06-10 | 2021-10-22 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿要素特征的定量评价方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1044286C (zh) * | 1994-07-15 | 1999-07-21 | 西安石油勘探仪器总厂 | 放射性能谱示踪吸水剖面测井法 |
| CN100516458C (zh) * | 2006-10-25 | 2009-07-22 | 大庆油田有限责任公司 | 低渗透薄互层压裂层段分层方法 |
| CN102094639B (zh) * | 2009-12-14 | 2013-12-04 | 西安威尔罗根能源科技有限公司 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
| US9383473B2 (en) * | 2012-06-26 | 2016-07-05 | Exxonmobil Upstream Research Company | Method for cement evaluation with neutron logs |
| CN105370273B (zh) * | 2015-12-10 | 2018-01-26 | 中国海洋石油集团有限公司 | 一种碳氧比测井环境影响因素校正的方法 |
| CN106761666B (zh) * | 2016-11-09 | 2019-12-17 | 中国石油大学(北京) | 一种四探头散射伽马测井及非线性数据反演的方法及装置 |
-
2019
- 2019-02-21 CN CN201910129076.7A patent/CN109779611A/zh not_active Withdrawn
- 2019-07-17 CN CN201910645297.XA patent/CN110344817B/zh active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534283A (zh) * | 2021-06-10 | 2021-10-22 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿要素特征的定量评价方法 |
| CN113534283B (zh) * | 2021-06-10 | 2024-03-19 | 核工业北京地质研究院 | 一种砂岩型铀矿成矿要素特征的定量评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110344817B (zh) | 2020-10-16 |
| CN110344817A (zh) | 2019-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108457640B (zh) | 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法 | |
| CN104360415B (zh) | 一种致密砂岩储层裂缝识别的方法 | |
| US8694261B1 (en) | 3D-well log invention | |
| CN108825220B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井仪及铀定量方程 | |
| US11099289B2 (en) | Multivariate analysis of seismic data, microseismic data, and petrophysical properties in fracture modeling | |
| Mendoza et al. | Linear iterative refinement method for the rapid simulation of borehole nuclear measurements: Part I—Vertical wells | |
| US10690802B2 (en) | Cement evaluation using neutron tool | |
| CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
| CN108825219B (zh) | 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法 | |
| CN105301639A (zh) | 基于vsp旅行时双加权层析反演速度场的方法及其装置 | |
| CN103748486A (zh) | 使用多个标准和空间区域的核频谱分析校正 | |
| CN118265927A (zh) | 用于获得真实井眼西格玛和真实地层西格玛的神经网络 | |
| CN109779611A (zh) | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 | |
| Zhang et al. | A new gamma density measurement method for cased-hole formation evaluation | |
| CN110469324A (zh) | 一种基于脉冲中子测井的计算地层密度方法 | |
| CA3040439C (en) | Petrophysical field evaluation using self-organized map | |
| CN107489418B (zh) | 密度测井正演模拟方法及装置 | |
| Ge et al. | Porosity evaluation of igneous rocks based on deep learning | |
| Zhou et al. | Experimental study of the influence of borehole parameters on prompt fission neutron uranium logging and its corrections | |
| Ma et al. | Searching Method of Critical Slip Surface of Slope Based on Improved Wolf Swarm Algorithm | |
| CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
| CN116500694B (zh) | 套后物理量反演方法、装置、计算设备及存储介质 | |
| Geier et al. | Prediction of Inflow into the D-Holes at the Stripa Mine | |
| Cai et al. | Theoretical Method and Experimental Research on High Precision Four Detector Density Logging | |
| Baxter et al. | Exploring conditioned simulations of discrete fracture networks in support of hydraulic acceptance of deposition holes |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20190521 |
|
| WW01 | Invention patent application withdrawn after publication |