一种适用于测井约束反演的波阻抗曲线去压实处理方法技术

本发明专利技术提供一种适用于测井约束反演的测井纵波阻抗曲线去压实处理方法，包括分岩性计算并剔除各岩性的测井纵波阻抗压实趋势，以及利用测井曲线统计得到的各岩性测井纵波阻抗基线对各自岩性进行测井纵波阻抗基线赋值，从而得到去压实处理的测井纵波阻抗曲线，并应用于测井约束反演中，进行整个工区的储层纵波阻抗定量预测。该适用于测井约束反演的测井纵波阻抗曲线去压实处理方法能够合理地放大储层与围岩之间的纵波阻抗差异，减少了储层与围岩在原始测井纵波阻抗曲线上因纵波阻抗值叠合区间较大而造成反演结果的多解性，使得最终的纵波阻抗反演结果能够准确地预测储层的分布情况。

【技术实现步骤摘要】
一种适用于测井约束反演的波阻抗曲线去压实处理方法
本专利技术涉及地球物理勘探数据处理
，特别是涉及到一种适用于测井约束反演的测井纵波阻抗去压实处理方法。
技术介绍
测井约束反演是目前储层预测中的一种重要方法，应用范围较为广泛。在浊积岩、砂砾岩体和滩坝砂岩等岩性、构造-岩性油藏勘探开发中发挥了重要的作用。通常意义下的测井约束反演所得到的结果是地震的纵波阻抗反演体，具有明确的物理意义，在储层预测中，通过结合储层在测井纵波阻抗曲线上的分布特征来指导整个工区的储层预测。但大量的实际应用表明测井约束反演所得到的纵波阻抗反演体在储层预测中有时存在着较大多解性，主要表现为储层与围岩在纵波阻抗分布区间上存在较大的重叠区间，难以有效区分。造成这一现象的根本原因在于，在原始的测井纵波阻抗曲线上，储层与围岩之间本身就存在着较大的重叠区间，因而在测井约束下所得到的地震纵波阻抗反演体也必然存在这一问题。针对这一问题，现有的解决方法是采取声波曲线重构的方法来放大储层与围岩之间的速度差异，具体做法是在原始的声波曲线中引入自然伽马、电阻率、自然电位等此类对储层响应敏感的测井曲线，通过数学拟合公式，将这些原始测井曲线重新拟合为一条新的曲线，称为“伪声波”曲线，再与密度曲线相乘，得到一条“伪纵波阻抗”曲线，从而参与到最终的测井约束反演中。分析认为，这种声波曲线重构方法所得到的“伪声波”曲线虽然能放大储层与围岩之间的速度差异，但采取这种“伪声波”曲线最终得到的纵波阻抗反演体的预测结果并不准确，原因在于采取这种重构方式得到的“伪声波”曲线会严重的影响到由此而最终计算得到的反射系数曲线，这种曲线处理方式会明显改变测井的反射系数序列，致使其与子波褶积时所得到的合成记录道与原始地震道之间的匹配性变差而影响到最终的反演结果。
技术实现思路
本专利技术的目的是解决了以上技术问题，提供一种适用于测井约束反演的测井纵波阻抗去压实处理方法，通过对原始测井的纵波阻抗曲线进行合理的处理，从而有效地放大储层与围岩之间纵波阻抗差异，使得最终的纵波阻抗反演体能够准确地预测储层的分布情况。本专利技术的目的通过如下的技术方案来实现：本适用于测井约束反演的测井纵波阻抗去压实处理方法包括：步骤1，将测井纵波阻抗曲线与岩性曲线对应，分岩性提取出各岩性的测井纵波阻抗曲线；步骤2，计算各岩性的纵波阻抗非压实分量；步骤3，剔除步骤2中计算得到的各岩性纵波阻抗非压实分量，用各岩性的纵波阻抗值减去各自岩性的纵波阻抗非压实分量，得到剔除纵波阻抗非压实分量后的各岩性纵波阻抗；步骤4，分岩性分别拟合计算出剔除了纵波阻抗非压实分量的各岩性纵波阻抗的压实趋势，并将其剔除掉。即利用剔除了纵波阻抗非压实分量的各岩性纵波阻抗减去各自岩性的纵波阻抗压实趋势，得到去压实趋势后的各岩性纵波阻抗；步骤5，将步骤2中计算得到的各岩性纵波阻抗非压实分量回归到步骤4中得到的去压实趋势后的各岩性纵波阻抗中，从而得到各岩性的去压实纵波阻抗；步骤6，将步骤5中得到的各岩性去压实趋势后的纵波阻抗分别进行各自岩性的纵波阻抗基线赋值，得到基线赋值后的各岩性纵波阻抗；步骤7，将基线赋值后的各岩性纵波阻抗进行合并，得到一条去压实趋势处理后的测井纵波阻抗曲线，并参与到后续的纵波阻抗反演中。上述技术方案进一步说明：在步骤1中，通过岩性曲线与测井的纵波阻抗曲线之间的对应关系，将每种岩性的测井纵波阻抗单独提取出来，从而将连续型的测井组波阻抗曲线分解成多个离散型的单岩性测井纵波阻抗曲线。在步骤2中，各岩性的纵波阻抗值随着埋深的增加会显著增大，但不同岩性之间影响纵波阻抗变化的因素并不完全相同。以砂、泥岩为例，泥岩纵波阻抗的变化主要受压实的影响，非压实因素的影响较小，可忽略不计，而砂岩的纵波阻抗的变化除了受压实因素的影响外，还受到含流体性、孔隙性等非压实因素的影响，将这部分影响因素统称为纵波阻抗非压实分量，纵波阻抗非压实分量可由岩石物理实验室模拟计算得出。在步骤3中，用各岩性的纵波阻抗值减去各自岩性的纵波阻抗非压实分量后，各岩性的纵波阻抗值随着埋深的变化的影响因素只剩下压实因素。剔除纵波阻抗非压实分量的目的是为了后续能够更准确的计算出各岩性的纵波阻抗压实趋势。在纵波阻抗非压实分量的剔除过程中，因为泥岩的纵波阻抗非压实分量较小，可以忽略不计，但砂岩、灰岩等岩性的纵波阻抗非压实分量需要剔除。在步骤4中，对各岩性纵波阻抗的压实趋势进行拟合计算。泥岩的纵波阻抗压实趋势，可直接由泥岩的纵波阻抗曲线进行拟合计算，而砂岩、灰岩等岩性的纵波阻抗压实趋势，只有在剔除了纵波阻抗非压实分量之后才能进行拟合计算。在拟合计算出各岩性的纵波阻抗压实趋势并将其剔除后，就得到了去压实趋势后的各岩性纵波阻抗曲线，此时需要注意的是，泥岩仅剔除了泥岩的纵波阻抗压实趋势，而砂岩、灰岩等岩性不仅剔除了它们的纵波阻抗压实趋势，在此之前还剔除了它们的纵波阻抗非压实分量。在步骤5中，需要将各岩性的纵波阻抗非压实分量进行回归，其原因在于，泥岩的去压实趋势纵波阻抗仅仅是剔除了它的纵波阻抗压实趋势，与之相对应，砂岩、灰岩等岩性也只能是剔除了它们的纵波阻抗压实趋势，因而需要将步骤2中计算得到的砂岩、灰岩等岩性的纵波阻抗非压实分量进行回归，这样得到的各岩性的去压实纵波阻抗曲线相对于各自岩性的原始纵波阻抗曲线而言，均只剔除了各自岩性的纵波阻抗压实趋势。在步骤5中，处理得到的各岩性的去压实纵波阻抗具有以下特征：处理之后的各岩性纵波阻抗的值域分布范围较为集中，并且由于消除了纵波阻抗压实趋势，同一岩性的纵波阻抗在浅层、中层、深层都比较接近。在步骤6中，各岩性的测井纵波阻抗基线是从原始的测井纵波阻抗曲线中选取埋深相对较浅的层段中统计出来的各岩性纵波阻抗平均值，因为在埋深相对较浅的层段中，各岩性的纵波阻抗受到压实的影响相对较小，各岩性之间，如泥岩、砂岩、灰岩等岩性在浅层埋深段中叠合区间较小，纵波阻抗平均值存在较大的差异。在步骤6中，在选取埋深相对较浅的层段时，应遵循以下原则：优选目的层系的相对浅层部分，若沉积相带变化不大，可取得较好的效果，如果沉积相带差异较大，则可考虑取目的层系之上的相对浅层层系；在目的层系之上的相对浅层层系选择中，优先选择与目的层系具有相似沉积背景的层系。遵循以上原则有利于经过处理而最终得到的纵波阻抗曲线所计算的反射系数序列与子波褶积所得到的合成地震记录道到原始地震之间具有最佳的匹配性。在步骤7中，基线赋值后的各岩性纵波阻抗在深度域上都是单一岩性的离散纵波阻抗值，需要进行合并，而合并后得到的是一条消除压实影响的测井纵波阻抗曲线，各岩性之间的纵波阻抗具有较大的差异，因而在后续的测井约束纵波阻抗反演中能够取得较好的效果。本专利技术的适用于测井约束反演的测井纵波阻抗曲线去压实处理方法能够合理地放大储层与围岩之间的纵波阻抗差异，减少了储层与围岩在原始测井纵波阻抗曲线上因纵波阻抗值叠合区间较大而造成反演结果的多解性，使得最终的纵波阻抗反演结果能够准确地预测储层的分布情况。附图说明图1为本专利技术的具体实施例的流程图；图2是实测井的原始纵波阻抗(A)曲线图；图3是实测井原始纵波阻抗(A)的概率密度分布图；图4是从实测井原始纵波阻抗(A)分岩性剥离出来的砂岩纵波阻抗(I)曲线和泥岩纵波阻抗(a)曲线图；图5是砂本文档来自技高网...

【技术保护点】
一种适用于测井约束反演的测井纵波阻抗去压实处理方法，其特征在于包括：步骤1，将测井纵波阻抗曲线与岩性曲线对应，分岩性提取出各岩性的测井纵波阻抗曲线；步骤2，计算各岩性的纵波阻抗非压实分量；

【技术特征摘要】
1.一种适用于测井约束反演的测井纵波阻抗去压实处理方法，其特征在于包括：步骤1，将测井纵波阻抗曲线与岩性曲线对应，分岩性提取出各岩性的测井纵波阻抗曲线；步骤2，计算各岩性的纵波阻抗非压实分量；步骤3，剔除步骤2中计算得到的各岩性纵波阻抗非压实分量，用各岩性的纵波阻抗值减去各自岩性的纵波阻抗非压实分量，得到剔除纵波阻抗非压实分量后的各岩性纵波阻抗；步骤4，分岩性分别拟合计算出剔除了纵波阻抗非压实分量的各岩性纵波阻抗的压实趋势，并将其剔除掉；即利用剔除了纵波阻抗非压实分量的各岩性纵波阻抗减去各自岩性的纵波阻抗压实趋势，得到去压实后的各岩性纵波阻抗；步骤5，将步骤2中计算得到的各岩性纵波阻抗非压实分量回归到步骤4中得到的去压实趋势后的各岩性纵波阻抗中，从而得到各岩性的去压实纵波阻抗；步骤6，将步骤5中得到的各岩性去压实趋势后的纵波阻抗分别进行各自岩性的纵波阻抗基线赋值，得到基线赋值后的各岩性纵波阻抗；步骤7，将基线赋值后的各岩性纵波阻抗进行合并，得到一条去压实趋势处理后的纵波阻抗曲线，并参与到后续的测井约束反演中。2.根据权利要求1所述...

【专利技术属性】
技术研发人员：罗平平，崔世凌，谭明友，韩宏伟，潘兴祥，曲志鹏，于景强，杨旭，蕾蕾，徐元刚，李晓晨，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：北京,11