CN107339099B - 一种确定储层岩性的方法和装置 - Google Patents
一种确定储层岩性的方法和装置 Download PDFInfo
- Publication number
- CN107339099B CN107339099B CN201710589979.4A CN201710589979A CN107339099B CN 107339099 B CN107339099 B CN 107339099B CN 201710589979 A CN201710589979 A CN 201710589979A CN 107339099 B CN107339099 B CN 107339099B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- detected
- time difference
- curve
- stratum
- 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 65
- 238000012937 correction Methods 0.000 claims abstract description 43
- 239000004576 sand Substances 0.000 claims abstract description 41
- 230000015572 biosynthetic process Effects 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000004048 modification Effects 0.000 claims description 6
- 238000012986 modification Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 239000011435 rock Substances 0.000 abstract description 3
- 238000005755 formation reaction Methods 0.000 description 21
- 238000001914 filtration Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 235000019738 Limestone Nutrition 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 239000010459 dolomite Substances 0.000 description 4
- 229910000514 dolomite Inorganic materials 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000006028 limestone Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 239000003245 coal Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000007704 transition Effects 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
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 And Detection Of Objects (AREA)
Abstract
本发明提供了一种确定储层岩性的方法和装置，其中，该方法包括：确定待测储层是否为砂泥岩地层；在确定待测储层为砂泥岩地层的情况下，获取待测储层的地层修正常量；获取待测储层的泥质含量和声波时差曲线的低频分量；根据待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到待测储层的储层特征曲线；根据储层特征曲线，确定待测储层的储层岩性。在本发明实施例中，能够大大提高了砂泥岩的区分精度，从而达到提高储层预测效果的目的。
Description
技术领域
本发明涉及地质勘探技术领域，特别涉及一种确定储层岩性的方法和装置。
背景技术
目前，一般可以采用波阻抗反演的方法来确定储层岩性。由于储层中砂岩与泥岩的纵波阻抗差异不明显，含气砂岩与泥岩的纵波阻抗叠置严重，利用纵波阻抗信息预测储层岩性难度较大。对于叠前波阻抗反演而言，可以根据工区中振幅随角度(即，偏移距)的变化情况来确定储层岩性。然而，对于振幅随角度或者偏移距的增加变化较小或者不变的砂泥岩工区而言，利用叠前波阻抗反演的方法预测储层岩性时会受到明显限制，反演结果也较难反映地下实际储层的含气情况。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本发明提供了一种确定储层岩性的方法和装置，以达到准确预测待测储层中岩性的目的。
本发明实施例提供了一种确定储层岩性的方法，可以包括：确定待测储层是否为砂泥岩地层；在确定所述待测储层为砂泥岩地层的情况下，获取所述待测储层的地层修正常量，其中，所述地层修正常量是根据参考储层的岩性特征确定的；获取所述待测储层的泥质含量和声波时差曲线的低频分量；根据所述待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到所述待测储层的储层特征曲线；根据所述储层特征曲线，确定所述待测储层的储层岩性。
在一个实施例中，获取所述待测储层的声波时差曲线的低频分量，可以包括：获取所述待测储层的声波时差曲线；将所述声波时差曲线校正至预设坐标系；从校正至预设坐标系的声波时差曲线中提取出低频分量。
在一个实施例中，确定待测储层是否为砂泥岩地层，可以包括：获取所述待测储层的声波时差曲线；将所述声波时差曲线校正至预设坐标系；获取校正至预设坐标系的声波时差曲线的取值范围；在所述取值范围在预设区间内的情况下，则确定所述待测储层为砂泥岩地层。
在一个实施例中，所述预设区间的取值范围可以为大于等于230us/m且小于等于300us/m。
在一个实施例中，可以按照以下公式计算得到所述待测储层的储层特征曲线：
AC_tezheng＝nAC_low+kSH_deal-P
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，n表示第二常量，k表示第三常量，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量，P表示所述地层修正常量。
在一个实施例中，所述地层修正常量的取值范围可以是大于等于25us/m且小于等于35us/m。
本发明实施例还提供了一种确定储层岩性的装置，可以包括：地层确定模块，用于确定待测储层是否为砂泥岩地层；常量获取模块，用于在确定所述待测储层为砂泥岩地层的情况下，获取所述待测储层的地层修正常量，其中，所述地层修正常量是根据参考储层的岩性特征确定的；储层参量获取模块，用于获取所述待测储层的泥质含量和声波时差曲线的低频分量；特征曲线计算模块，用于根据所述待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到所述待测储层的储层特征曲线；储层岩性确定模块，用于根据所述储层特征曲线，确定所述待测储层的储层岩性。
在一个实施例中，所述储层参量获取模块可以包括：时差曲线获取单元，可以用于获取所述待测储层的声波时差曲线；曲线校正单元，可以用于将所述声波时差曲线校正至预设坐标系；低频分量提取单元，可以用于从校正至预设坐标系的声波时差曲线中提取出低频分量。
在一个实施例中，所述地层确定模块可以包括：时差曲线获取单元，可以用于获取所述待测储层的声波时差曲线；曲线校正单元，可以用于将所述声波时差曲线校正至预设坐标系；取值范围获取单元，可以用于获取校正至预设坐标系的声波时差曲线的取值范围；砂泥岩地层确定单元，可以用于在所述取值范围在预设区间内的情况下，则确定所述待测储层为砂泥岩地层。
在一个实施例中，所述预设区间的取值范围可以为大于等于230us/m且小于等于300us/m。
在本发明实施例中，根据待测储层是否为砂泥岩地层，对区分后的地层进行储层特征曲线的计算。具体的，可以根据待测储层中声波时差曲线的低频分量、泥质含量以及地层修正常量，计算得到砂泥岩地层的特征曲线；可以根据声波时差曲线确定不是砂泥岩地层所对应的特征曲线。通过针对砂泥岩地层采用相应的特征曲线计算方法，解决了现有技术中由于对所有地层均采用相同的波阻抗反演方式而导致的砂泥岩反演结果差异不明显、重叠现象严重的缺陷，同时利用地层修正常量对不易区分的砂泥岩地层的储层特征曲线进行进一步的修正，大大提高了砂泥岩的区分精度，最终达到了提高储层预测效果的目的。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种确定储层岩性的方法流程图；
图2是本申请提供的三维地震资料区块示意图；
图3是本申请提供的AC曲线、特征曲线示意图，其中，图3(a)是AC曲线示意图，图3(b)是特征曲线示意图；
图4是本申请提供的一种确定储层岩性的装置的一种结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
考虑到利用叠前反演或者波阻抗反演的方法确定储层岩性时，由于对所有的待测储层采用相同的储层特征曲线而导致的砂泥岩反演结果差异不明显、重叠现象严重，难以反映地层真实情况的缺陷，鉴于泥质含量曲线在区分砂泥岩方面的优越性能，以及声波时差曲线在连接地震与测井曲线(时深关系)时的稳定性，发明人提出了先判断待测储层是否为砂泥岩地层，并针对待测储层为砂泥岩地层的情况下，根据砂泥岩地层中的声波时差曲线的低频分量、泥质含量以及地层修正常量来确定砂泥岩的岩性。基于此，提出了一种确定储层岩性的方法，如图1所示，可以包括以下步骤：
S101：确定待测储层是否为砂泥岩地层。
在本申请的一个实施例中，在得到待测储层之后，可以先确定待测储层是否为砂泥岩地层。其中，上述砂泥岩地层可以是砂岩地层和泥岩地层。
具体的，可以采用下述方式来确定待测储层是否为砂泥岩地层：
S1-1：获取待测储层的声波时差曲线；
S1-2：将声波时差曲线校正至预设坐标系；特别的，当只存在一条声波时差曲线时，可以不用校正，直接观察测得的这条声波时差曲线。
S1-3：获取校正至预设坐标系的声波时差曲线的取值范围；
S1-4：在取值范围在预设区间内的情况下，则确定待测储层为砂泥岩地层。
声波时差可以指的是接收声波的时间差值，同时，声波时差一般指的是纵波声波时差。利用这个时间差值可以进行相关运算，并求解各种量值。在地质勘探领域可以缩写为AC。常见岩石声波时差如下：砂岩为182us/m，灰岩为155us/m，白云岩为141us/m等。
受井径不规则、垮塌等问题的影响，声波、密度等测井曲线会发生畸变，需要进行环境校正。此外，受测井仪器、环境(如泥浆)、人为因素等问题的影响，不同井之间的测井曲线刻度可能发生变化，需要对其进行归一化校正处理。经过校正，各井测井曲线平均值与方差趋于一致。
本文针对测井曲线环境校正采用的方法主要为：运用电阻率拟合法校正井径垮塌井段的声波曲线。
在采用上述方式进行环境校正之后，可以通过归一化处理的方式，将不同井之间的测井曲线校正至预设坐标系。例如：可以假定不同井之间的测井曲线均满足正态分布，并计算多个测井曲线的平均值以及方差。根据计算得到的平均值以及方差调整每口井的测井曲线值，使得得到的每口井的测井曲线值的平均值以及方差等于上述计算得到的平均值以及方差。采用上述方法可以将不同井之间的多个测井曲线校正至相同的坐标系，从而保证对不同井之间的值可以在相同的刻度下进行对比。在本申请的一个实施例中，预设区间的取值范围可以为大于等于230us/m且小于等于300us/m。
即，当声波时差曲线取值范围大于等于230us/m且小于等于300us/m时，可以确定待测储层为砂泥岩储层；当声波时差曲线取值范围小于230us/m或者大于300us/m时，可以确定待测储层为其他岩性储层。具体的，运用声波时差曲线可以直接识别煤层(声波时差大于300us/m)，但砂泥岩相互叠置(声波时差在230～300us/m之间)，灰岩与白云岩亦相互叠置(声波时差均小于230us/m)。
S102：在确定待测储层为砂泥岩地层的情况下，获取待测储层的地层修正常量，其中，地层修正常量是根据参考储层的岩性特征确定的；
S103：获取待测储层的声波时差曲线的低频分量和泥质含量。
在本申请的一个实施例中，在确定待测储层为砂泥岩地层的情况下，可以获取待测储层的地层修正常量、获取待测储层的声波时差曲线的低频分量和泥质含量。
在本申请的一个实施例中，可以根据参考储层的岩性特征确定地层修正常量。具体的，地层修正常量的取值范围可以是大于等于25us/m、且小于等于35us/m。
在本申请的一个实施例中，可以按照以下方式获取待测储层的声波时差曲线的低频分量：
S2-1-1：获取待测储层的声波时差曲线；
S2-1-2：将声波时差曲线校正至预设坐标系；
S2-1-3：从校正至预设坐标系的声波时差曲线中提取出低频分量。
在采用S1-1和S1-2的方法得到校正至预设坐标系的多井的测井参数(声波时差曲线)之后，由于原始的AC曲线变化较快，可以对AC曲线进行低频滤波处理，从而得到声波时差曲线的低频分量。AC曲线不能有效区分砂泥岩，但可以正确标定时深关系，低频滤波后的AC曲线不仅保留了地层的变化趋势，并且有利于之后和泥质含量曲线结合起来，对地层岩性进行判断。
其中，低频滤波处理的方法包括但不限于以下处理方法中的一种：限幅滤波法、中位值滤波法、算术平均滤波法、递推平均滤波法、中位值平均滤波法、限幅平均滤波法、一阶滞后滤波法、加权递推平均滤波法、消抖滤波法以及限幅消抖滤波法。具体的，可以采用以下公式进行低频滤波运算：
AC_Low＝LowFilter(AC，overlap)
其中，AC表示AC曲线，AC_Low表示低频滤波后的AC曲线，overlap表示深度域滤波因子。在本申请的一个实施例中，overlap＝0.1。
指示砂泥岩的曲线有自然伽马、泥质含量等，而以SH(泥质含量)曲线最为直接。如果用自然伽马来区分岩性，不同井之间会存在一定的数据漂移，此外砂质泥岩、泥质砂岩与砂岩、泥岩数据之间均存在一些重叠。反之，泥质含量曲线定量指示砂岩、泥岩的百分比，不同井之间不存在数据漂移的问题。因此，特征曲线构建的关键是将泥质含量进行一定处理，使其在某常数值附近变化。在本申请的一个实施例中，可以按照以下方式获取待测储层的泥质含量：
S2-2-1：获取待测储层中的初始泥质含量以及预设泥质含量；
S2-2-2：对比初始泥质含量与预设泥质含量的大小；
S2-2-3：当初始泥质含量不大于预设泥质含量时，泥质含量等于初始泥质含量；
S2-2-4：当初始泥质含量大于预设泥质含量时，根据初始泥质含量与预设参量的差值计算得到泥质含量。
由于泥质含量曲线(SH)在0～100％之间变化，而砂泥岩的判识标准通常以25％为界。因此，以预设泥质含量为25％时为基线，对于岩性大于25％的预设泥质含量(泥质砂岩、砂质泥岩、泥岩)时，采取(SH-25)/3+25的计算方法，使得处理后的SH曲线在0～50％之间变化；当初始泥质含量不大于25％时，直接等于初始测得的泥质含量。即，可以采用以下公式进行计算得到泥质含量：
其中，SH_deal表示为泥质含量，SH表示初始泥质含量。
S104：根据待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到待测储层的储层特征曲线。
在本申请的一个实施例中，在确定待测储层为砂泥岩地层的情况下，可以按照以下公式计算得到待测储层的储层特征曲线：
AC_tezheng＝nAC_low+kSH_deal-P
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，n表示第二常量，k表示第三常量，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量，P表示所述地层修正常量。其中，所述地层修正常量的取值范围是大于等于25us/m、且小于等于35us/m。
在本申请的一个实施例中，n＝k＝1，P＝32us/m。可以按照以下公式计算得到待测储层的储层特征曲线：
AC_tezheng＝AC_low+SH_deal-32
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量。
在本申请的另一个实施例中，在确定待测储层为除砂泥岩地层以外的其他地层(例如：灰岩、白云岩及煤层等)的情况下，可以按照以下公式计算得到待测储层的储层特征曲线：
AC_tezheng＝AC
S105：根据储层特征曲线，确定待测储层的储层岩性。
波阻抗反演技术在鄂尔多斯盆地及其周边地区得到广泛的运用。然而，由于目的层段砂、泥岩纵波阻抗重叠现象非常普遍，利用该技术存在严重的多解性问题。在本申请的一个实施例中，反演确定待测储层的储层岩性时，可以将AC曲线替换成既能识别岩性又能保持AC曲线基本骨架的储层特征曲线。在得到待测储层的储层特征曲线之后，可以利用得到的储层特征曲线反演得到待测储层的储层岩性。即，当特征曲线的取值大于300us/m时，所对应位置处的待测储层为煤层；当特征曲线的取值范围为小于等于300us/m且大于等于280us/m时，所对应位置处的待测储层为泥岩；当特征曲线的取值范围为小于等于280us/m且大于等于250us/m时，所对应位置处的待测储层为砂质泥岩；当特征曲线的取值范围为小于等于250us/m且大于等于225us/m时，所对应位置处的待测处为砂岩；当特征曲线的取值范围为小于225us/m时，所对应位置处的待测储层为灰岩或者白云岩。
下面结合一个具体的实施例对上述确定储层岩性的方法进行具体说明，然而值得注意的是，该具体实施例仅是为了更好地说明本发明，并不构成对本发明的不当限定。
苏里格气田上古生界盒8段为主要产气层，由于声波在砂泥岩段传播速度十分接近，故而利用声波阻抗区分岩性的反演方法不适用于该地区。
为了检验特征曲线反演方法是否能够适应该区块储层预测的需要，选择一小块三维地震资料，运用该曲线进行波阻抗反演方法试验。如图2所示为所选取的三维地震资料区块示意图，其中，横坐标表示主测线号，纵坐标表示联络测线号。该区块位于苏里格南国际合作区1区，地表附近低降速厚度在20～50米左右，原始单炮记录视主频在15～20Hz之间，已有直井探井9口，测井曲线种类齐全。
为了构建既能反映地下储层岩性特征，又能与AC曲线变化趋势相一致的储层特征曲线。首先，选择SN108-05、106-05、143-05、145-05等4口井进行储层特征曲线构建，具体步骤为：
根据上述4口井的测井参数，确定上述4口井是否为砂泥岩地层；
1)在确定为砂泥岩地层的情况下，可以采用以下步骤构建储层特征曲线：S1：对AC曲线进行简单的环境校正与归一化处理，并对AC曲线进行低频滤波，得到特征曲线的低频分量；S2：获取待测储层中的泥质含量曲线；S3：基于现有的测井数据，确定待测储层的地层修正常量；S4：将泥质含量曲线以及低频分量进行叠加，得到储层特征曲线；
2)在确定为砂泥岩地层以外的地层时，将原始的AC曲线作为储层特征曲线。
通过上述特征曲线的构建，可以得到既能区分岩性，又能用来制作合成记录的特征曲线。
为了检验特征曲线在井震标定方面与AC曲线存在的差异，首先，可以运用30Hz零相位雷克子波分别对二者所制作的合成记录进行了比较，如图3所示的AC曲线、特征曲线示意图，图3中虚线表示井轨迹，对比图3中的波形形态、振幅强弱、以及转换时间，可知：苏里格气田所得到的AC曲线、特征曲线所制作的合成记录基本一致，这说明特征曲线完全遵从AC曲线所反映的时深关系。其中，图3(a)表示AC曲线，图3(b)表示特征曲线。
为了进一步测试特征曲线的稳定性，也可以运用两种曲线对工区内多口井分别提取子波，并进行综合，结果表明，运用特征曲线提取的子波在波形、相位方面要比AC曲线所提的子波更稳定。
综上所述，如果用特征曲线进行反演，那么所取得的效果自然更加理想。
在4口井特征曲线构建的基础上，分别进行基于模型和基于地震道的波阻抗反演，具体步骤如下：
1)井震标定
首先确立Tc2反射层为标志层，运用30Hz零相位雷克子波与特征曲线的褶积建立合成记录，并放大振幅系数，使之与地震振幅相匹配，然后通过同相轴的移动使得合成记录与井旁地震道的相关系数达到峰值。
2)子波提取
在进行子波提取时，需要考虑子波的长度，并对振幅、相位进行估算，理想的子波应单峰值，相位在0°左右，有效频带范围内保持稳定，如果相位出现一定的偏差，可以通过对时深关系进行上下微调来完成(若波峰出现在0°右侧则可以向下滑动，反之则可以向上滑动)。井与井之间所提取的子波应大体类似，这样随后得到的综合子波才能最大程度地兼顾各口井的合成记录，反演效果就越稳定，也越可信。
3)模型建立
通过特征曲线的内插外推，建立符合地质规律的波阻抗模型。
由于苏里格区块上古构造极其简单，没有明显的断层，且地层厚度横向变化不大，故而在建模时层与层之间选择平行顶底的方法是可以的。如果在构造复杂地区建模时，必须综合考虑如断层、超覆、削截、尖灭、盐丘等情况，这样才可以得到更加合理的模型。
4)反演
1、可以采用模型反演的方式确定待测储层的岩性：
模型反演是油气田进入开发阶段后运用最广泛的储层预测技术，多口井参与基础上的模型反演能细致反映储层的横向变化，具有纵向分辨率高的优势，然而受井控制的程度高，如果参与反演的井过少，会出现失真的问题。运用上述4口井进行模型反演，并利用剩余的5口井进行盲井检验，从而评估该算法的合理性。采用相同参数与流程，分别运用特征曲线与AC曲线进行模型反演结果表明：运用特征曲线进行的反演结果与盲井吻合更好。
2、可以采用地震道反演的方式确定待测储层的岩性：
地震道反演是油气田在勘探初期，钻井数量少、测井资料缺乏时所采用的储层预测技术，该反演方法忠实于地震信息，能识别大套的砂岩储层，但纵向分辨率低，不能有效分辨薄层。
通过对特征曲线与AC曲线地震道反演效果对比，可知：二者在纵向上对砂岩储层识别能力都要比模型反演弱很多，大约能识别20米左右的单砂层，而特征曲线与井旁地震道反演结果吻合性更好。
苏里格南区块的特征曲线反演试验证明，运用该方法，可以较快速地获得直观反映地下岩性变化的波阻抗反演数据体。与传统波阻抗反演方法相比，该方法不存在多解性，更容易解释；与叠前反演相比，该方法运算方便快捷，所需数据量小(不需要对叠前地震数据进行分角度或分偏移距叠加，不需要提取多个地震子波，不需要横波测井数据)，因此，可以预见：特征曲线的构建及相应的反演技术能在今后的储层预测中大显身手，发挥独特且不可替代的作用。
基于同一发明构思，本发明实施例中还提供了一种确定储层岩性的装置，如下面的实施例所述。由于确定储层岩性的装置解决问题的原理与确定储层岩性的方法相似，因此确定储层岩性的装置的实施可以参见确定储层岩性的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图4是本发明实施例的确定储层岩性的装置的一种结构框图，如图4所示，可以包括：地层确定模块401、常量获取模块402、储层参量获取模块403、特征曲线计算模块404、储层岩性确定模块405，下面对该结构进行说明。
地层确定模块401，可以用于确定待测储层是否为砂泥岩地层；
常量获取模块402，可以用于在确定所述待测储层为砂泥岩地层的情况下，获取所述待测储层的地层修正常量，其中，所述地层修正常量是根据参考储层的岩性特征确定的；
储层参量获取模块403，可以用于获取所述待测储层的泥质含量和声波时差曲线的低频分量；
特征曲线计算模块404，可以用于根据所述待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到所述待测储层的储层特征曲线；
储层岩性确定模块405，可以用于根据所述储层特征曲线，确定所述待测储层的储层岩性。
在一个实施例中，所述储层参量获取模块可以包括：时差曲线获取单元，可以用于获取所述待测储层的声波时差曲线；曲线校正单元，可以用于将所述声波时差曲线校正至预设坐标系；低频分量提取单元，可以用于从校正至预设坐标系的声波时差曲线中提取出低频分量。
在一个实施例中，所述地层确定模块可以包括：时差曲线获取单元，可以用于获取所述待测储层的声波时差曲线；曲线校正单元，可以用于将所述声波时差曲线校正至预设坐标系；取值范围获取单元，可以用于获取校正至预设坐标系的声波时差曲线的取值范围；砂泥岩地层确定单元，可以用于在所述取值范围在预设区间内的情况下，则确定所述待测储层为砂泥岩地层。
在一个实施例中，所述预设区间可以为大于等于230us/m且小于等于300us/m。
在一个实施例中，所述特征曲线计算模块具体可以用于按照以下公式计算得到所述待测储层的储层特征曲线：
AC_tezheng＝nAC_low+kSH_deal-P
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，n表示第二常量，k表示第三常量，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量，P表示所述地层修正常量。
在一个实施例中，所述地层修正常量的取值范围可以是大于等于25us/m、且小于等于35us/m。
利用上述各实施例所提供的确定储层岩性的装置的实施方式，可以自动实施计算得到储层特征曲线的方法，对储层进行岩性确定，可以不需要实施人员的具体参与，可以直接输出储层岩性确定结果，操作简单快捷，有效提高了用户体验。
确定储层岩性的装置中，确定待测储层是否为砂泥岩地层、声波时差曲线的低频分量的获取、储层特征曲线的计算的实施方式的扩展可以参照前述方法的相关描述。
尽管本申请内容中提到砂泥岩地层的确定方式、低频分量的获取方式、储层特征曲线的计算方式等描述，但是，本申请并不局限于必须是本申请实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取/计算存储/判断等获取的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的单元、装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机***环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器***、基于微处理器的***、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何***或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (9)
1.一种确定储层岩性的方法，其特征在于，包括：
确定待测储层是否为砂泥岩地层；
在确定所述待测储层为砂泥岩地层的情况下，获取所述待测储层的地层修正常量，其中，所述地层修正常量是根据参考储层的岩性特征确定的；
获取所述待测储层的泥质含量和声波时差曲线的低频分量；
根据所述待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到所述待测储层的储层特征曲线；
根据所述储层特征曲线，确定所述待测储层的储层岩性；
其中，按照以下公式计算得到所述待测储层的储层特征曲线：
AC_tezheng＝nAC_low+kSH_deal-P
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，n表示第二常量，k表示第三常量，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量，P表示所述地层修正常量。
2.如权利要求1所述的方法，其特征在于，获取所述待测储层的声波时差曲线的低频分量，包括：
获取所述待测储层的声波时差曲线；
将所述声波时差曲线校正至预设坐标系；
从校正至预设坐标系的声波时差曲线中提取出低频分量。
3.如权利要求1所述的方法，其特征在于，确定待测储层是否为砂泥岩地层，包括：
获取所述待测储层的声波时差曲线；
将所述声波时差曲线校正至预设坐标系；
获取校正至预设坐标系的声波时差曲线的取值范围；
在所述取值范围在预设区间内的情况下，则确定所述待测储层为砂泥岩地层。
4.如权利要求3所述的方法，其特征在于，所述预设区间的取值范围为大于等于230us/m且小于等于300us/m。
5.如权利要求1所述的方法，其特征在于，所述地层修正常量的取值范围是大于等于25us/m且小于等于35us/m。
6.一种确定储层岩性的装置，其特征在于，包括：
地层确定模块，用于确定待测储层是否为砂泥岩地层；
常量获取模块，用于在确定所述待测储层为砂泥岩地层的情况下，获取所述待测储层的地层修正常量，其中，所述地层修正常量是根据参考储层的岩性特征确定的；
储层参量获取模块，用于获取所述待测储层的泥质含量和声波时差曲线的低频分量；
特征曲线计算模块，用于根据所述待测储层的地层修正常量、泥质含量和声波时差曲线的低频分量，计算得到所述待测储层的储层特征曲线；
储层岩性确定模块，用于根据所述储层特征曲线，确定所述待测储层的储层岩性；
其中，所述特征曲线计算模块具体用于按照以下公式计算得到所述待测储层的储层特征曲线：
AC_tezheng＝nAC_low+kSH_deal-P
其中，AC_tezheng表示所述储层特征曲线，单位为us/m，n表示第二常量，k表示第三常量，AC_low表示所述声波时差曲线的低频分量，SH_deal表示所述待测储层中的泥质含量，P表示所述地层修正常量。
7.如权利要求6所述的装置，其特征在于，所述储层参量获取模块包括：
时差曲线获取单元，用于获取所述待测储层的声波时差曲线；
曲线校正单元，用于将所述声波时差曲线校正至预设坐标系；
低频分量提取单元，用于从校正至预设坐标系的声波时差曲线中提取出低频分量。
8.如权利要求6所述的装置，其特征在于，所述地层确定模块包括：
时差曲线获取单元，用于获取所述待测储层的声波时差曲线；
曲线校正单元，用于将所述声波时差曲线校正至预设坐标系；
取值范围获取单元，用于获取校正至预设坐标系的声波时差曲线的取值范围；
砂泥岩地层确定单元，用于在所述取值范围在预设区间内的情况下，则确定所述待测储层为砂泥岩地层。
9.如权利要求8所述的装置，其特征在于，所述预设区间的取值范围为大于等于230us/m且小于等于300us/m。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710589979.4A CN107339099B (zh) | 2017-07-19 | 2017-07-19 | 一种确定储层岩性的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710589979.4A CN107339099B (zh) | 2017-07-19 | 2017-07-19 | 一种确定储层岩性的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107339099A CN107339099A (zh) | 2017-11-10 |
| CN107339099B true CN107339099B (zh) | 2020-06-09 |
Family
ID=60215937
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710589979.4A Active CN107339099B (zh) | 2017-07-19 | 2017-07-19 | 一种确定储层岩性的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107339099B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111045108B (zh) * | 2018-10-11 | 2022-12-13 | 中国石油化工股份有限公司 | 随钻横波计算方法及*** |
| CN109343120B (zh) * | 2018-10-17 | 2019-10-01 | 吉林大学 | 融入约束稀疏脉冲反演低频补偿的声波曲线重构方法 |
| CN109386284B (zh) * | 2018-11-23 | 2021-09-28 | 中国石油化工股份有限公司 | 基于多曲线联合校正的储层物性分析方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3724589A (en) * | 1969-05-26 | 1973-04-03 | Schlumberger Technology Corp | Well logging methods and apparatus |
| US5349528A (en) * | 1990-03-12 | 1994-09-20 | Halliburton Logging Services, Inc. | Method apparatus for determination of porosity lithological composition |
| CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
| CN105301647A (zh) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | 区分灰质泥岩和砂岩的方法 |
-
2017
- 2017-07-19 CN CN201710589979.4A patent/CN107339099B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3724589A (en) * | 1969-05-26 | 1973-04-03 | Schlumberger Technology Corp | Well logging methods and apparatus |
| US5349528A (en) * | 1990-03-12 | 1994-09-20 | Halliburton Logging Services, Inc. | Method apparatus for determination of porosity lithological composition |
| CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
| CN105301647A (zh) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | 区分灰质泥岩和砂岩的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 储层特征曲线重构技术在储层预测中的应用研究;张静 等;《天然气地球科学》;20080630;第19卷(第3期);第396-401页 * |
| 灰质背景下基于敏感曲线的储层预测方法——以孤北洼陷沙三中亚段为例;乌洪翠;《中国石油勘探》;20170531;第22卷(第3期);第56-62页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107339099A (zh) | 2017-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10324211B2 (en) | Seismic spectral balancing | |
| US20180321404A1 (en) | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function | |
| US20220291418A1 (en) | An integrated geomechanics model for predicting hydrocarbon and migration pathways | |
| CN104675392B (zh) | 一种基于叠前多参数降维的储层岩性识别方法 | |
| Aarre et al. | Seismic detection of subtle faults and fractures | |
| CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
| US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| EP2420863A2 (en) | Determining a Position of a Geological Layer Relative to a Wavelet Response in Seismic Data | |
| Koesoemadinata et al. | Seismic reservoir characterization in Marcellus shale | |
| CN104678434A (zh) | 一种预测储层裂缝发育参数的方法 | |
| CN107339099B (zh) | 一种确定储层岩性的方法和装置 | |
| CN107045143A (zh) | 一种预测裂缝发育的方法及装置 | |
| CN105607120A (zh) | 一种基于时移测井的地震相约束构建初始模型的方法 | |
| Lefranc et al. | Fractured basement reservoir characterization for fracture distribution, porosity and permeability prediction | |
| US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
| CN111090117B (zh) | 一种相控正演约束下的有效储层预测方法及*** | |
| CN104820237B (zh) | 计算地层的方差体的方法 | |
| CN113376692B (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
| CN112305594B (zh) | 非均质储层的油气分布确定方法及*** | |
| CN107589457A (zh) | 测井资料尺度粗化的方法和装置 | |
| Chen et al. | Application of Digital Twin Technology to 3D Seismic Iterative Processing and Dynamic Geological Modeling | |
| CN117930380A (zh) | 基于地震沉积学的碳酸盐岩储层预测方法和*** | |
| Noghretab et al. | Simulation of VSP data based on DSI data and estimation of shear wave velocity and elastic Moduli for a well case study | |
| Ausburn et al. | Modern Seismic Methods-An Aid for the Petroleum Engineer | |
| CN105372708B (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 |