CN106597542B - 一种储层特征参数的预测方法及装置 - Google Patents
一种储层特征参数的预测方法及装置 Download PDFInfo
- Publication number
- CN106597542B CN106597542B CN201611046845.XA CN201611046845A CN106597542B CN 106597542 B CN106597542 B CN 106597542B CN 201611046845 A CN201611046845 A CN 201611046845A CN 106597542 B CN106597542 B CN 106597542B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- indicate
- node
- density
- modulus
- 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 57
- 239000011435 rock Substances 0.000 claims abstract description 255
- 230000004044 response Effects 0.000 claims abstract description 32
- 239000012530 fluid Substances 0.000 claims description 151
- 229920006395 saturated elastomer Polymers 0.000 claims description 95
- 239000011159 matrix material Substances 0.000 claims description 87
- 239000011148 porous material Substances 0.000 claims description 85
- 238000010008 shearing Methods 0.000 claims description 69
- 238000001615 p wave Methods 0.000 claims description 52
- 230000006835 compression Effects 0.000 claims description 43
- 238000007906 compression Methods 0.000 claims description 43
- 230000000153 supplemental effect Effects 0.000 claims description 39
- 238000009826 distribution Methods 0.000 claims description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 239000004575 stone Substances 0.000 claims description 10
- 230000010287 polarization Effects 0.000 claims description 6
- 238000012916 structural analysis Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 22
- 238000003860 storage Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 239000000047 product Substances 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 4
- 238000011017 operating method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 241001274660 Modulus Species 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种储层特征参数的预测方法及装置。所述方法包括：获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据通过地震叠前反演得到；将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测；输出预测得到的待预测位置的储层特征参数数据。利用本申请中各个实施例，可以通过构建的三维岩石物理模版，准确预测出储层的孔隙扁度、饱和度和孔隙度，可以更全面地反映出储层的结构特性，进而有效提高储层的结构分析结果的真实性和准确度。
Description
技术领域
本申请涉及地球物理勘探技术领域，特别涉及一种储层特征参数的预测方法及装置。
背景技术
储层特征参数的预测，是对储层进行预测的基础。准确而全面的储层特征参数数据，是能够精确预测储层的保证。因此，在储层预测过程中，储层特征参数的预测是一个重要环节。
现有技术中，对储层特征参数的预测，主要是采用交会图技术。利用所述交会图技术，可以较为准确地预测出储层的孔隙度和饱和度，进而可以较为准确地预测出储层的结构等特性。但是，对于一些致密储层来说，由于储层的孔隙结构比较复杂，可能具有两种甚至两种以上的孔隙结构。这个时候，仅仅依靠储层的孔隙度和饱和度两种参数，就无法更准确地预测储层的真实结构了。而现有技术中，没有同时预测储层的孔隙度、饱和度和储层孔隙结构特征参数的方法。
现有技术中至少存在如下问题：
现有技术只预测得到了储层的孔隙度和饱和度两种储层特征参数，没有得到反映储层的孔隙结构的特征参数，而且现有技术中也没有同时预测储层的孔隙度、饱和度和储层孔隙结构特征参数的方法。这样就导致了在利用得到的储层特征参数进行储层预测时，不能准确真实地预测出储层的真实结构。
发明内容
本申请实施例的目的是提供一种储层特征参数的预测方法及装置，以构建一种三维岩石物理模板，提供更全面更准确的用于储层预测的特征参数数据。
本申请实施例提供一种储层特征参数的预测方法及装置是这样实现的：
一种储层特征参数的预测方法，所述方法包括：
获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据通过地震叠前反演得到；
将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测，所述三维岩石物理模版包括采用下述构建方式生成：确定储层特征参数的分布范围，对所述分布范围进行网格剖分，生成三维网格，计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度，利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版；
输出预测得到的待预测位置的储层特征参数数据。
优选实施例中，所述对所述分布范围进行网格剖分，生成三维网格的方式包括：
分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系；
分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分，得到所述三维网格。
优选实施例中，所述利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版的方式，包括：
建立分别以纵波阻抗、横波阻抗和密度为三个坐标轴的第二个三维坐标系；
以所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度为坐标，将所述三维网格的各节点定位于所述第二个三维坐标系中，得到定位后的节点；
以所述定位后的节点确定出所述三维岩石物理模版。
优选实施例中，所述储层特征参数包括：
储层的孔隙扁度，储层的孔隙度，储层的饱和度。
优选实施例中，所述计算所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗的方式包括：
获取储层的岩石基质的体积模量、剪切模量和密度；
获取储层的孔隙中流体的体积模量和密度；
利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量；
利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量，进而计算得到储层孔隙中流体的密度；
利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量；
利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度；
利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
优选实施例中，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：
和表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
优选实施例中，采用包括下述计算公式计算得到储层孔隙中的体积模量：
Kf＝SWKW+(1-SW)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
优选实施例中，采用包括下述计算公式计算得到储层孔隙中流体的密度：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
优选实施例中，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
优选实施例中，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度；
φ表示各节点对应的孔隙度。
优选实施例中，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat
式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
优选实施例中，所述将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置进行储层特征参数预测的方式，包括：
利用所述待预测位置的三种地震响应特征参数数据，在所述三维岩石物理模版中，确定出所述待预测位置的储层特征参数数据。
优选实施例中，所述装置包括：
数据获取模块，用于获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据；
模版构建模块，用于构建预先构建好的三维岩石物理模版；
参数预测模块，用于将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测；
数据输出模块，用于输出预测得到的待预测位置的储层特征参数数据。
优选实施例中，所述模版构建模块包括：
分布范围确定模块，用于确定储层特征参数的分布范围；
网格生成模块，用于对所述储层特征参数的分布范围进行网格剖分，生成三维网格；
数据计算模块，用于计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度；
模版生成模块，用于利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，构建所述三维岩石物理模版。
优选实施例中，所述数据计算模块包括：
模量计算模块，用于利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量，还用于利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量，还用于利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量；
地震响应特征参数计算模块，用于利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度，还用于计算得到储层孔隙中流体的密度；
波阻抗计算模块，用于利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
优选实施例中，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：
和表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
优选实施例中，采用包括下述计算公式计算得到储层孔隙中的体积模量：
Kf＝SwKw+(1-Sw)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
优选实施例中，采用包括下述计算公式计算得到储层孔隙中流体的密度：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
优选实施例中，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
优选实施例中，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度；
φ表示各节点对应的孔隙度。
优选实施例中，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat
式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
利用本申请实施例提供的一种储层特征参数的预测方法，可以通过构建的三维岩石物理模版，准确预测出储层的孔隙扁度、饱和度和孔隙度。其中，所述储层的孔隙扁度可以反映储层的孔隙结构。利用所述三维岩石物理模版预测得到的三种储层特征参数，不仅可以准确地反映储层的孔隙大小以及储层对流体的储存特性，还可以准确预测出储层的孔隙结构，这样就可以更全面地反映出储层的结构特性。进而有效提高储层的结构分析结果的真实性和准确度。利用本申请实施例提供的一种储层特征参数的预测装置，可以自动实施所述储层特征参数的预测方法，并自动输出预测结果，操作简单快捷，不需要实施人员的过多参与，有效提高了用户体验。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一个实施例中提供的一种储层特征参数的预测方法的方法流程图；
图2是本申请一个实施例中提供的一种三维岩石物理模版的构建方式的方法流程图；
图3是本申请一个实施例中提供的一种确定出所述三维岩石物理模版的方法流程图；
图4是本申请一个实施例中提供的一种储层特征参数的预测装置的模块结构示意图；
图5是本申请一个实施例中提供的模版构建模块的模块结构示意图；
图6(a)是本申请一个实施例中得到的预先构建好的三维岩石物理模版；
图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比；
图7(a)是本申请另一个实施例中得到的预先构建好的三维岩石物理模版；
图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比；
图8(a)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙度连井剖面；
图8(b)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的饱和度连井剖面；
图8(c)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙扁度连井剖面；
图9是本申请一个实施例中提供的计算处理模块的模块结构示意图。
具体实施方式
本申请实施例提供一种储层特征参数的预测方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
图1是本申请所述一种储层特征参数的预测方法一种实施例的方法流程图。虽然本申请提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本申请实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理的实施环境)。
具体的，如图1所示，本申请提供的一种储层特征参数的预测方法的一种实施例可以包括：
S1：获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据包括通过地震叠前反演得到。
所述储层中待预测位置可以包括储层中的任意位置。
所述待预测位置可以只是一个位置，也可以是N个位置，N≥2。
在本申请一个实施例中，所述待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，可以通过地震叠前反演得到。
当然，具体的获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据的方式不必限定，在本申请其他实施例中，也可以通过其他方式获取，比如在本申请另一个实施例中，可以通过直接测量的方式得到所述三种地震响应特征数据。
S2：将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测。
所述三维岩石物理模版主要由三维网格中的各节点确定出来。
所述储层特征参数至少包括储层的孔隙扁度、储层的孔隙度、储层的饱和度。
图2是本申请一个实施例中提供的一种三维岩石物理模版的构建方式的方法流程图，具体的，如图2所示，所述预先构建好的三维岩石物理模版的构建方式可以包括：
S201：确定储层特征参数的分布范围，对所述分布范围进行网格剖分，生成三维网格。
所述储层特征参数的分布范围，可以认为是所述储层特征参数的变化区间，每个储层特征参数的变化区间都是可以确定的，比如，储层的孔隙扁度的分布范围可以是大于等于0且小于等于1。所述确定储层特征参数的分布范围属于本领域技术人员可以掌握的常规技术，具体的确定方式不必限定。
所述对所述分布范围进行网格剖分，可以包括对每种储层特征参数的分布范围进行划分，所述划分的方式，可以是等间隔划分，也可以是不等间隔划分，划分的间隔的宽度也不必限定。比如，对于上述的储层的孔隙扁度，划分的间隔的宽度可以均取为0.1，比如所述储层的孔隙扁度的分布范围为[0,1]，划分之后，每一个小间隔可以依次为[0,0.1)、(0.1,0.2)……(0.9,1.0]。当然，划分的间隔的宽度越小，最终得到的三维岩石物理模版的预测精度就越高。具体的，在本申请其他实施例中，可以根据所需要的预测精度以及实际实施情形，确定所述对每种储层特征参数的分布范围进行划分的方式以及划分的间隔的宽度。
最终由每种储层特征参数的分布范围划分出的间隔组成三维网格，所述三维网格的每个节点对应一组储层特征参数数据，所述一组储层特征参数可以包括储层的孔隙扁度、孔隙度和饱和度。
S202：计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。
所述三维网格中各节点的孔隙扁度、孔隙度和饱和度是已经确定的。
利用所述各节点的孔隙扁度、孔隙度和饱和度三种储层特征参数，结合储层的岩石基质的体积模量、剪切模量和密度，储层的孔隙中流体的体积模量和密度，就可以计算出所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。
S203：利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的各节点确定出所述三维岩石物理模版。
所述对所述三维网格中各节点进行定位，可以包括：将计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度分别作为所述各节点的三个坐标，利用所述坐标，将所述各节点在同一个分别以纵波阻抗、横波阻抗和密度为坐标轴的三维坐标系中定位。
利用定位后的各节点，就可以确定出所述三维岩石物理模版。
所述确定出所述三维岩石物理模版的方式属于常规的数据处理手段，可以采用拟合等方式确定出所述三维岩石物理模版，具体的确定方式不必限定。
利用上述实施例中得到的所述三维岩石物理模版，可以预测出储层的三种储层特征参数，可以包括：储层的孔隙扁度，储层的孔隙度，储层的饱和度。
S3：输出预测得到的待预测位置的储层特征参数数据。
所述待预测位置的储层特征参数可以包括：储层的孔隙扁度，储层的孔隙度和储层的饱和度。
具体的，输出所述预测得到的待预测位置的储层特征参数数据，可以根据实施人员的需要选择性输出，比如可以同时输出待预测位置的孔隙扁度、孔隙度和饱和度的数据，也可以只输出孔隙度和饱和度的数据，当然，也可以只输出孔隙扁度一种参数的数据。
同样的，在数据输出阶段，所述待预测位置可以只是一个位置，也可以是N个位置，N≥2。也就是说，在本申请其他实施例中，可以输出储层中一个位置的储层特征参数数据，也可以同时输出储层中N个位置的储层特征参数数据，N≥2。
利用上述各实施例提供的一种储层特征参数的预测方法可以利用所述三维岩石物理模版预测得到的储层的孔隙扁度、储层的孔隙度和储层的饱和度三种储层特征参数，不仅可以准确地反映储层的孔隙大小以及储层对流体的储存特性，还可以准确预测出储层的孔隙结构，这样就可以更全面地反映出储层的结构特性。进而有效提高储层的结构分析结果的真实性和准确度。
在本申请另一个实施例中，对所述分布范围进行网格剖分，生成三维网格的方式可以包括：
分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系。
本申请一个实施例中，所述三种预定储层特征参数可以包括储层的孔隙扁度、储层的孔隙度、储层的饱和度，选择上述三种储层特征参数，可以更全面地反映储层的结构特征，尤其是储层的孔隙特征。
当然，在本申请其他实施例中，实施人员也可以根据实际作业情况和作业需要，选择其他的储层特征参数作为三种预定储层特征参数，比如三种预定储层特征参数可以是储层的渗透率、储层的孔隙度和储层的饱和度。
分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分，得到所述三维网格。
所述分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分的方式，可以包括对在所述每种储层特征参数对应的坐标轴上，对所述每种储层特征参数的分布范围进行划分，所述划分的方式，可以是等间隔划分，也可以是不等间隔划分，划分的间隔的宽度也不必限定。
利用上述实施例，可以得到一个三维网格，并且三维网格中的每个节点对应着一组三种储层特征参数数据。
在本申请又一个实施例中，利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的各节点确定出所述三维岩石物理模版的方式如图3所示，具体的，可以包括：
S2031：建立分别以纵波阻抗、横波阻抗和密度为三个坐标轴的第二个三维坐标系。
S2032：以所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度为坐标，将所述三维网格的各节点定位于所述第二个三维坐标系中，得到定位后的各节点。
S2033：以所述定位后的各节点确定出所述三维岩石物理模版。
利用所述三维岩石物理模版，准确预测出储层的孔隙度、饱和度和孔隙扁度。
在本申请再一个实施例中，将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置进行储层特征参数预测的方式，可以包括：
利用所述待预测位置的三种地震响应特征参数数据，在所述三维岩石物理模版中，确定出所述待预测位置的储层特征参数数据。
具体的，确定出所述待预测位置的储层特征参数数据的方式，属于常规数学方法，可以以三种地震响应特征参数数据为坐标，在所述三维岩石物理模版中，确定出所述坐标对应的点，所述坐标对应的点对应的储层特征参数数据就是所述待预测位置的储层特征参数数据。
在本申请另一个实施例中，计算所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗的方式，可以包括：
获取储层的岩石基质的体积模量、剪切模量和密度。
获取储层的孔隙中流体的体积模量和密度。
利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量。
在本申请一个实施例中，可以包括采用下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：
和表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量，进而计算得到储层孔隙中流体的密度。
在本申请一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中流体的体积模量：
Kf＝SwKw+(1-Sw)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
在本申请一个实施例中，所述计算得到储层孔隙中流体的密度的计算公式可以包括：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量。
在本申请一个实施例中，可以包括采用下述计算公式计算得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度。
在本申请一个实施例中，可以包括采用下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度。
利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
在本申请一个实施例中，可以包括采用下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat
式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
利用上述实施例提供的计算方式，可以计算出所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗。
利用上述各实施例提供的一种储层特征参数的预测方法，可以得到的所述三维岩石物理模版，并利用所述三维岩石物理模版，准确预测出待预测位置的储层特征参数数据。
图6(a)是本申请一个实施例中得到的预先构建好的三维岩石物理模版，图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。
图7(a)是本申请另一个实施例中得到的预先构建好的三维岩石物理模版，图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。从图6和图7中，可以看出，利用构建的三维岩石物理模版，可以准确预测出储层的孔隙扁度、饱和度和孔隙度。并且预测得到的结果的真实性和准确度得到了有效提高。
图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。图中，纵坐标表示所述待预测位置的深度，单位是米(m)。图中黑色曲线表示的是测量得到的纵横波速度和密度曲线以及孔隙度和饱和度曲线，灰色曲线表示的是根据所述三维岩石物理模版预测得到的孔隙扁度及孔隙度和饱和度。可以看出，无论是储层的孔隙度还是储层的饱和度，测井解释和利用所述三维岩石物理模版预测得到的结果吻合度较高，这也证明了利用所述三维岩石物理模版进行储层特征参数预测具有较高的可靠性和准确度。从反演的结果看，储层段的孔隙扁度变化范围大，从0.1到1.0都存在分布，从图6(a)的三维岩石物理模版上也可以发现，这说明碳酸岩盐储层的纵横波速度不仅受孔隙度和饱和度的影响，也受孔隙扁度的影响，利用纵横波密度数据反演孔隙度和饱和度需要考虑孔隙扁度这一参数。
图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。图中，纵坐标表示所述待预测位置的深度，单位是米(m)。可以看到，利用所述三维岩石物理模版预测得到的孔隙度和饱和度与测井解释结果具有较高的吻合度。证明了利用所述三维岩石物理模版进行储层特征参数预测具有较高的可靠性和准确度。
图8是本申请一个实例中根据预测得到的储层特征参数数据反演得出的连井剖面。
其中，图8(a)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙度连井剖面。
图8(b)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的饱和度连井剖面。
图8(c)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙扁度连井剖面。
利用储层的孔隙扁度、饱和度和孔隙度，可以得到相应的连井剖面。
基于本申请所述的储层特征参数的预测方法，本申请提供一种储层特征参数的预测装置，所述装置可以集成在储层预测的功能组件中，进行储层预测。图4是本申请一个实施例中提供的储层特征参数的预测装置的模块结构示意图。如图4所示，所述装置可以包括：
数据获取模块101，可以用于获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据。
模版构建模块102，可以用于构建预先构建好的三维岩石物理模版。
参数预测模块103，可以用于将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测。
数据输出模块104，可以用于输出预测得到的待预测位置的储层特征参数数据。利用上述实施例提供的一种储层特征参数的预测装置，可以自动实施所述储层特征参数的预测方法，并自动输出预测结果，操作简单快捷，不需要实施人员的过多参与，有效提高了用户体验。
在本申请另一个实施例中，所述参数预测模块的模块结构示意图如图5所示，具体的，可以包括：
分布范围确定模块1031，可以用于确定储层特征参数的分布范围。
网格生成模块1032，可以用于对所述储层特征参数的分布范围进行网格剖分，生成三维网格。
数据计算模块1033，可以用于计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。
模版生成模块1034，可以用于利用计算得到的三维网格节点对应的纵波阻抗、横波阻抗和密度，构建所述三维岩石物理模版。
利用上述实施例提供的参数预测模块，可以自动构建所述三维岩石物理模版，操作简单快捷，提高了用户体验。
在本申请又一个实施例中，所述数据计算模块的模块结构示意图如图9所示，具体的，可以包括：
模量计算模块10331，可以用于利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量。
还可以用于利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量。
还可以用于利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量。
地震响应特征参数计算模块10332，可以用于利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度。还可以用于计算得到储层孔隙中流体的密度。
波阻抗计算模块10333，可以用于利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
利用上述实施例提供的一种数据计算模块，可以直接计算出所述各节点对应的密度、纵波阻抗和横波阻抗。
本申请另一个实施例中，所述计算得到储层的干岩石骨架的体积模量和剪切模量，可以包括采用下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：
和表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
本申请另一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中的体积模量：
Kf＝SwKw+(1-Sw)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
本申请另一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中流体的密度：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
本申请另一个实施例中，可以包括采用下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
本申请另一个实施例中，可以包括采用下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度。
本申请另一个实施例中，可以包括采用下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat
式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
所述储层特征参数的预测装置中，所述构建三维岩石物理模版、将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中对所述待预测位置的储层特征参数进行预测、输出预测得到的待预测位置的储层特征参数数据的实施方式的扩展可以参照前述方法的相关描述。
尽管本申请内容中提到不同的储层特征参数的预测方式，从获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据、构建三维岩石物理模版、对所述待预测位置的储层特征参数进行预测到输出预测得到的待预测位置的储层特征参数数据的各种时序方式、数据获取/处理/输出方式等的描述，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (19)
1.一种储层特征参数的预测方法，其特征在于，所述方法包括：
获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据通过地震叠前反演得到；
将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测，所述三维岩石物理模版包括采用下述构建方式生成：确定储层特征参数的分布范围，对所述分布范围进行网格剖分，生成三维网格，计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度，利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版；
输出预测得到的待预测位置的储层特征参数数据；
其中，所述对所述分布范围进行网格剖分，生成三维网格的方式包括：
分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系；
分别在所述三个坐标轴上对三种预定储层特征参数的分布范围进行划分，得到所述三维网格。
2.如权利要求1所述的一种储层特征参数的预测方法，其特征在于，所述利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版的方式，包括：
建立分别以纵波阻抗、横波阻抗和密度为三个坐标轴的第二个三维坐标系；
以所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度为坐标，将所述三维网格的各节点定位于所述第二个三维坐标系中，得到定位后的节点；
以所述定位后的节点确定出所述三维岩石物理模版。
3.如权利要求1所述的一种储层特征参数的预测方法，其特征在于，所述储层特征参数包括：
储层的孔隙扁度，储层的孔隙度，储层的饱和度。
4.如权利要求1所述的一种储层特征参数的预测方法，其特征在于，所述计算所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗的方式包括：
获取储层的岩石基质的体积模量、剪切模量和密度；
获取储层的孔隙中流体的体积模量和密度；
利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量；
利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量，进而计算得到储层孔隙中流体的密度；
利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量；
利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度；
利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
5.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：P*i-Q*i＝a+bK*(y)/G*(y)；
P*i和Q*i表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
6.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述计算公式计算得到储层孔隙中的体积模量：
Kf＝SWKW+(1-SW)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
7.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述计算公式计算得到储层孔隙中流体的密度：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
8.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
9.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度；
φ表示各节点对应的孔隙度。
10.如权利要求4所述的一种储层特征参数的预测方法，其特征在于，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat
式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
11.如权利要求1所述的一种储层特征参数的预测方法，其特征在于，所述将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置进行储层特征参数预测的方式，包括：
利用所述待预测位置的三种地震响应特征参数数据，在所述三维岩石物理模版中，确定出所述待预测位置的储层特征参数数据。
12.一种储层特征参数的预测装置，其特征在于，所述装置包括：
数据获取模块，用于获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据；
模版构建模块，用于构建预先构建好的三维岩石物理模版；
参数预测模块，用于将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测；
数据输出模块，用于输出预测得到的待预测位置的储层特征参数数据；
所述模版构建模块包括：
分布范围确定模块，用于确定储层特征参数的分布范围；
网格生成模块，用于分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系，分别在所述三个坐标轴上对三种预定储层特征参数的分布范围进行划分，得到所述三维网格；
数据计算模块，用于计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度；
模版生成模块，用于利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，构建所述三维岩石物理模版。
13.如权利要求12所述的一种储层特征参数的预测装置，其特征在于，所述数据计算模块包括：
模量计算模块，用于利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量，还用于利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量，还用于利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量；
地震响应特征参数计算模块，用于利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度，还用于计算得到储层孔隙中流体的密度；
波阻抗计算模块，用于利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。
14.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：
式中，Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
Gm表示储层的岩石基质的剪切模量；
φ表示各节点对应的孔隙度；
α表示各节点的孔隙扁度；
a和b满足：P*i-Q*i＝a+bK*(y)/G*(y)；
P*i和Q*i表示干岩石骨架的极化因子；
K*(y)和G*(y)表示岩石等效体积模量和剪切模量。
15.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述计算公式计算得到储层孔隙中的体积模量：
Kf＝SWKW+(1-SW)Kg
式中，Kf表示所述储层孔隙中流体的体积模量；
Sw表示所述各节点对应的饱和度；
Kw表示储层孔隙中水的体积模量；
Kg表示储层孔隙中气体的体积模量。
16.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述计算公式计算得到储层孔隙中流体的密度：
ρf＝Swρw+(1-Sw)ρg
式中，ρf表示所述储层孔隙中流体的密度；
Sw表示所述各节点对应的饱和度；
ρw表示储层孔隙中水的密度；
ρg表示储层孔隙中气体的密度。
17.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：
Gsat＝Gdry
式中，Ksat表示所述储层的饱和流体岩石的体积模量；
Gsat表示所述储层的饱和流体岩石的体积模量；
Kdry表示所述储层的干岩石骨架的体积模量；
Gdry表示所述储层的干岩石骨架的剪切模量；
Km表示储层的岩石基质的体积模量；
φ表示各节点对应的孔隙度；
Kf表示所述储层孔隙中流体的体积模量。
18.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：
ρsat＝(1-φ)ρm+φρf；
式中，ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度；
Gsat表示所述储层的饱和流体岩石的体积模量；
Ksat表示所述储层的饱和流体岩石的体积模量；
ρm表示储层的岩石基质的密度；
ρf表示所述储层孔隙中流体的密度；
φ表示各节点对应的孔隙度。
19.如权利要求13所述的一种储层特征参数的预测装置，其特征在于，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：
Ip＝ρVp；
Is＝ρVs；
ρ＝ρsat式中，Ip表示所述各节点对应的纵波阻抗；
Is表示所述各节点对应的横波阻抗；
ρ表示所述各节点对应的密度；
ρsat表示所述饱和流体岩石的密度；
Vp表示所述饱和流体岩石的纵波速度；
Vs表示所述饱和流体岩石的横波速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611046845.XA CN106597542B (zh) | 2016-11-23 | 2016-11-23 | 一种储层特征参数的预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611046845.XA CN106597542B (zh) | 2016-11-23 | 2016-11-23 | 一种储层特征参数的预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106597542A CN106597542A (zh) | 2017-04-26 |
| CN106597542B true CN106597542B (zh) | 2019-01-18 |
Family
ID=58592064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611046845.XA Active CN106597542B (zh) | 2016-11-23 | 2016-11-23 | 一种储层特征参数的预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106597542B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107367762B (zh) * | 2017-06-27 | 2019-06-11 | 中国石油天然气股份有限公司 | 一种确定储层参数的方法及装置 |
| CN108761529B (zh) * | 2018-05-03 | 2019-10-11 | 中国石油天然气股份有限公司 | 基于三维弹电岩石物理量版确定储层参数的方法及装置 |
| CN110426752B (zh) * | 2019-08-20 | 2021-05-25 | 赛哲尔能源科技(北京)有限公司 | 一种基于岩石物理模型的储层参数反演方法与系统 |
| CN113960664B (zh) * | 2021-10-25 | 2023-07-28 | 中国石油天然气股份有限公司 | 一种纵横波反射率模版估算水合物和游离气饱和度的方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9110190B2 (en) * | 2009-06-03 | 2015-08-18 | Geoscale, Inc. | Methods and systems for multicomponent time-lapse seismic measurement to calculate time strains and a system for verifying and calibrating a geomechanical reservoir simulator response |
| CN103278850B (zh) * | 2013-05-24 | 2014-08-20 | 西安石油大学 | 基于煤岩工业组分物理体积模型的横波时差曲线构建方法 |
| CN103984009B (zh) * | 2014-04-16 | 2016-04-13 | 刘致水 | 基于纵、横波测井资料的三种孔隙类型定量反演方法 |
| CN105182424B (zh) * | 2015-08-03 | 2017-09-01 | 中国石油天然气股份有限公司 | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 |
| CN105425280A (zh) * | 2015-11-21 | 2016-03-23 | 西南石油大学 | 一种矿物模量和孔隙结构预测方法 |
-
2016
- 2016-11-23 CN CN201611046845.XA patent/CN106597542B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106597542A (zh) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106597542B (zh) | 一种储层特征参数的预测方法及装置 | |
| Kim et al. | Estimation of fracture porosity of naturally fractured reservoirs with no matrix porosity using fractal discrete fracture networks | |
| CN105487122B (zh) | 用于确定储层流体识别因子的敏感性的方法及系统 | |
| Keogh et al. | The development of fluvial stochastic modelling in the Norwegian oil industry: A historical review, subsurface implementation and future directions | |
| Qu et al. | Incorporation of deformation band fault damage zones in reservoir models | |
| US10620339B2 (en) | Static earth model calibration methods and systems using tortuosity evaluations | |
| CN107219159B (zh) | 一种确定储层非均质性的方法和装置 | |
| CN109655903B (zh) | 页岩层横波速度预测方法及系统 | |
| CN105701319B (zh) | 一种水平井规则开发井网下的沉积微相建模方法 | |
| US10656295B2 (en) | Systems and methods for downscaling stress for seismic-driven stochastic geomechanical models | |
| CN106886046B (zh) | 确定缝洞型气藏未投产区块可动用储量的方法 | |
| CN106227957A (zh) | 等效裂缝建模的方法 | |
| CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
| CN107450108B (zh) | 甜点区的确定方法和装置 | |
| CN103946896A (zh) | 利用基于过程的模型和动态异质性估定感兴趣的地质体积的异质性的系统及方法 | |
| Liu et al. | Petrel2ANSYS: Accessible software for simulation of crustal stress fields using constraints provided by multiple 3D models employing different types of grids | |
| CN106772680A (zh) | 基于二维叠后地震资料微观裂缝表征确定方法和装置 | |
| CN106125133B (zh) | 一种基于气云区约束下的精细速度建模方法 | |
| CA3045880A1 (en) | Coupled reservoir-geomechanical models using compaction tables | |
| Suman et al. | Sensitivity study of rock-physics parameters for modeling time-lapse seismic response of Norne field | |
| CN106842316B (zh) | 裂缝确定方法和装置 | |
| CN107991716B (zh) | 薄层的沉积相图和沉积体厚度的确定方法和装置 | |
| CN106338768B (zh) | 一种生成储层预测属性数据的处理方法、装置及系统 | |
| CN111323814B (zh) | 基于岩石物理模板定量确定砂岩储层含水饱和度的方法 | |
| CN109001804B (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 |