CN111239810A - 基于结构模型的基岩储层预测方法及装置 - Google Patents
基于结构模型的基岩储层预测方法及装置 Download PDFInfo
- Publication number
- CN111239810A CN111239810A CN201811469385.0A CN201811469385A CN111239810A CN 111239810 A CN111239810 A CN 111239810A CN 201811469385 A CN201811469385 A CN 201811469385A CN 111239810 A CN111239810 A CN 111239810A
- Authority
- CN
- China
- Prior art keywords
- bedrock
- longitudinal wave
- wave impedance
- reservoir
- curve
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 35
- 239000011707 mineral Substances 0.000 claims abstract description 35
- 239000011148 porous material Substances 0.000 claims abstract description 33
- 239000007788 liquid Substances 0.000 claims abstract description 24
- 230000004044 response Effects 0.000 claims abstract description 16
- 239000011435 rock Substances 0.000 claims description 51
- 238000004458 analytical method Methods 0.000 claims description 41
- 239000000203 mixture Substances 0.000 claims description 32
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000004088 simulation Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 12
- 230000008859 change Effects 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 12
- 238000011161 development Methods 0.000 description 7
- 238000005553 drilling Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 239000002352 surface water Substances 0.000 description 4
- 238000002386 leaching Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- DLHONNLASJQAHX-UHFFFAOYSA-N aluminum;potassium;oxygen(2-);silicon(4+) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Si+4].[Si+4].[Si+4].[K+] DLHONNLASJQAHX-UHFFFAOYSA-N 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229910052652 orthoclase Inorganic materials 0.000 description 1
- 229910052655 plagioclase feldspar Inorganic materials 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 229910001868 water Inorganic materials 0.000 description 1
Images
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
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
技术领域
本发明涉及地球物理勘探技术领域，特别涉及一种基于结构模型的基岩储层预测方法及装置。
背景技术
基岩是一种盆地形成前的古老结晶岩，国内外已发现多个大型基岩油气藏，是非常有潜力的勘探目标层系。通常基岩由上至下可以划分为风化残积层(厚度1—10m)、强风化带(厚度30—150m)、弱风化带(100—500m)、母岩四部分。基岩储层主要发育在强风化带及弱风化带，其中强风化带储层受风化淋滤作用比较强，储层大范围连片分布，弱风化带储层主要受顺断层、裂缝向下渗流的地表水溶蚀作用形成，因此其发育位置受断层、裂缝发育位置的控制。目前勘探目标主要以强风化带为主，依据储层发育程度及勘探程度将风化残积层、强风化带定义为基岩浅层，弱风化带及母岩定义为基岩深层。
基岩储层预测通常是通过叠后稀疏脉冲反演预测基岩纵波阻抗来完成，反演过程中所需要的低频模型数据是利用测井纵波阻抗数据的内插来建立，然而通常以基岩为勘探目标的钻井都不会揭示深层基岩内幕，因此基岩深层的低频模型难以准确建立。目前，商业软件建立基岩深层低频模型的方法是利用纵波阻抗曲线最后一个样点值外推至基岩深层，这种方法具有很大的随机性，且不能够真实反映基岩深层纵波阻抗的变化趋势，因此难以准确建立基岩深层低频模型。
发明内容
本发明实施例提供了一种基于结构模型的基岩储层预测方法及装置，通过基于结构模型对基岩储层进行预测，能够更加准确的反映基岩浅层、深层储层的变化趋势，对储层进行精确预测。
该基于结构模型的基岩储层预测方法包括：
建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
根据测井数据确定储层测井响应特征和定量解释模板；
根据地震数据确定层位断层精细解释结果；
根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
该基于结构模型的基岩储层预测装置包括：
总孔隙度预测模块，用于建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
岩石骨架、流体成分及其弹性参数分析模块，用于确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
基岩深层纵波阻抗曲线模拟模块，用于根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
纵波阻抗曲线合并模块，用于将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
储层特征分析模块，用于根据测井数据确定储层测井响应特征和定量解释模板；
层位断层解释模块，用于根据地震数据确定层位断层精细解释结果；
稀疏脉冲反演模块，用于根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
储层精细预测模块，用于根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述基于结构模型的基岩储层预测方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述基于结构模型的基岩储层预测方法的计算机程序。
在本发明实施例中，基于结构模型的基岩储层预测方法是根据实测曲线分析基岩矿物成分组成，结合基岩深层的孔隙度预测结果，建立稳定的基岩结构模型，利用岩石物理模拟方法预测基岩深层纵波阻抗曲线，并与浅层实测曲线进行拼接，从而建立更加准确的基岩内部低频模型来约束稀疏脉冲反演运算。在保证浅层预测准确的前提下，更加准确的预测基岩深层的纵波阻抗变化，以取得更加准确的储层预测结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于结构模型的基岩储层预测方法流程图；
图2是本发明实施例提供的一种基于结构模型的基岩储层预测方法简化流程图；
图3是本发明实施例提供的一种基岩储层发育模式图；
图4是本发明实施例提供的一种叠前时间偏移成果地震剖面图；
图5是本发明实施例提供的一种基岩结构模型及原始纵波阻抗曲线及合成纵波阻抗曲线图；
图6是本发明实施例提供的一种应用基于结构模型方法前后的低频模型对比结果图；
图7是本发明实施例提供的一种应用基于结构模型方法前后的预测结果剖面对比结果图；
图8是本发明实施例提供的一种应用基于结构模型方法后的预测结果剖面图；
图9是本发明实施例提供的一种强风化带、弱风化带储层预测平面结果图；
图10是本发明实施例提供的一种基于结构模型的基岩储层预测装置结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明实施例中，提供了一种基于结构模型的基岩储层预测方法方法，如图1和图2所示，该方法包括：
步骤101：建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
具体的，基岩深层的纵波阻抗曲线模拟需要比较准确的确定基岩矿物成分及孔隙度。矿物成分可以通过实测曲线的分析结果来确定，孔隙度主要通过目标位置距基岩顶面距离(即基岩厚度)来进行计算。在暴露期，基岩由浅至深受风化淋滤作用的强度越来越弱，孔隙度越来越小，纵波阻抗越来越大，当基岩完全不受地表水风化淋滤作用时，孔隙度接近于零，此时纵波阻抗值最大且只与组成岩石的矿物成分相关。
可以按照如下方式建立基岩总孔隙度与基岩厚度的定量关系：
获取基岩浅层实测地质数据(即地质资料)；
根据所述基岩浅层实测地质数据计算孔隙度曲线；
将所述孔隙度曲线与基岩厚度进行交汇分析，获得分析结果；
根据所述分析结果和获取的区域地质数据(即区域地质资料)，建立基岩孔隙度同基岩厚度之间的定量关系，从而可以利用基岩厚度推算深层基岩总孔隙度。
具体的，上面基岩浅层实测数据只局限于局部地区(即要开展研究的目标区)，区域地质数据指的是在更广泛的范围内取得的地质数据(包括露头分析数据、测井数据、化验数据等)。通常受研究目标区范围局限的控制，利用基岩浅层实测数据很难建立孔隙度同基岩厚度之间的定量关系，需要参考区域地质数据，才可能获得更准确的基岩孔隙度同基岩厚度之间的定量关系。
步骤102：确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
具体的，按照如下方式确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数：
根据测井曲线(即测井资料)，确定基岩的矿物成分和孔隙液体成分；
根据完钻井的分析结果确定基岩矿物成分的体积含量；其中，完钻井分析结果包括测井曲线数据、录井数据、取芯数据、分析化验数据等，测井曲线只是其中一种；
根据区域分析结果确定孔隙流体的体积模量及密度；其中，流体的体积模量单独用测井曲线很难确定，且多数地区不做这方面的工作，这里的区域分析结果是指的相同盆地、相同沉积环境的邻区的流体体积模量分析结果；
通过纵波时差曲线、横波时差曲线和密度曲线，计算基岩不同组分的弹性参数，其中弹性参数包括体积模量、剪切模量和密度值。
步骤103：根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
其中，基岩结构模型整体框架是存在的，即基岩是由骨架、孔隙、孔隙液体组成，但骨架、孔隙、孔隙液体的参数是变化的(比如骨架的组成矿物、孔隙度的大小、孔隙液体的成分等)，需根据上述信息建立实际模型。
具体的，根据所述基岩结构模型预测基岩深层纵波阻抗曲线，包括：
根据所述基岩结构模型进行基岩深层纵波速度和密度的模拟，获得模拟结果；
根据所述模拟结果计算基岩深层纵波阻抗；
根据所述基岩深层纵波阻抗确定基岩深层纵波阻抗曲线。
步骤104：将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
步骤105：根据测井数据确定储层测井响应特征和定量解释模板；
具体的，根据测井解释结果，通过直方图、交会图等方法，确定储层测井响应特征，并建立定量解释模板。
定量解释模板指的是储层定量解释模板，即储层解释的标准。储层是根据测井数据解释出来的，解释储层的过程中，需要参考多种测井数据，设定不同数据的门槛值，门槛值以内是储层，以外不是储层(比如根据声波时差和伽马曲线解释储层，根据测井响应分析认为：伽马小于40gAPI，声波时差小于220us/m为储层，那么储层的伽马、声波时差的门槛值分别为40gAPI、220us/m)。通常储层解释是利用两种数据的交会分析来进行的(有时也用一种数据或两种以上的数据)，在交会图上确定储层门槛值即为储层定量解释模板。
步骤106：根据地震数据确定层位断层精细解释结果；
具体的，井震结合，确定主要目的层的地震反射特征，并据此进行层位、断层对比追踪。
步骤107：根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
步骤108：根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
具体的，①利用纵波阻抗反演结果，根据储层测井响应特征分析结果，提取属性平面图，对储层进行定性预测。
②利用纵波阻抗反演结果，根据储层定量解释模板，提取反演结果中的储层样点数并累积，对储层进行定量预测。
实施实例：
具体实施情况如下：
1)根据钻井分析结果及理论研究成果，基岩风化壳由上至下可以划分为风化残积层、强风化带、弱风化带、母岩四部分，基岩储层主要发育在强风化带及弱风化带，其中强风化带储层受风化淋滤作用比较强，储层大范围连片分布，弱风化带储层主要受顺断层、裂缝向下渗流的地表水溶蚀作用形成，因此其发育位置受断层、裂缝发育位置的控制(图3)。
2)根据钻井标定结果，基岩顶界解释为一强波谷，由于上覆地层与基岩之间的阻抗差异较明显，因此基岩顶界为全区连续稳定的波谷反射特征，容易对比追踪(图 4)。
3)根据基岩矿物成分分析结果，基岩骨架主要由石英、斜长石、正长石、暗色矿物、绿泥石组成，孔隙内充填水(图5a)。根据完钻井的分析结果，确定矿物成分的体积含量；利用纵波曲线、横波曲线、密度推算不同矿物的体积模量、剪切模量、密度值，根据区域分析结果确定孔隙流体的体积模量及密度；根据浅层实测曲线及区域地质资料，建立总孔隙度同基岩厚度之间的定量关系。利用上述数据进行岩石物理模拟，预测基岩深层纵波阻抗曲线变化趋势，并同实测曲线进行拼接。图5b为高切滤波后的原始纵波阻抗曲线，图5c为高切滤波后的合成纵波阻抗曲线，基岩深层上部随着深度的增加，纵波阻抗逐渐变大，下部孔隙不再发育，纵波阻抗保持不变。利用合成曲线可以更好的控制基岩深层的储层预测结果。
4)图6为低频模型的对比，其中图6的左面的图为利用合成纵波阻抗曲线建立的低频模型(0—10HZ)，图6的右面的图为利用原始纵波阻抗曲线建立的低频模型 (0—10Hz)。与原始模型相比，图6的左面的图由基岩浅层至深层，纵波阻抗为逐渐增大的变化规律，至基岩深层下部，纵波阻抗保持不变，符合实际地质认识。
5)图7为反演结果的对比，其中图7a为利用合成纵波阻抗曲线建立的低频模型作为约束计算而得的反演结果，图7b为利用原始纵波阻抗曲线建立的低频模型作为约束计算而得的反演结果。与原始结果(图7b)相比，新成果(图7a)在纵向更加能够准确反映储层由基岩浅层至深层的变化规律：基岩浅层受风化淋滤作用比较强，储层连片发育，深层受顺断层而下的地表水溶蚀作用的控制，储层呈带状发育，更深层的基岩储层基本不发育。同时在新成果上，容易划分出与地质结果相对应的强风化带、弱风化带和母岩(图8)，方便对不同的带分别开展储层预测工作。
6)分别提取强风化带与弱风化带的纵波阻抗属性平面图来表征其储层发育情况。图9a为强风化带储层预测结果，储层连片发育，实际工区内钻遇强风化带的钻井有 21口，预测结果与实钻结果有18口吻合，符合率85.71％；图9b为弱风化带储层预测结果，储层带状发育，实际工区内钻遇弱风化带的钻井有2口，预测结果与实钻结果均吻合，符合率100％。
基于同一发明构思，本发明实施例中还提供了一种基于结构模型的基岩储层预测装置，如下面的实施例所述。由于基于结构模型的基岩储层预测装置解决问题的原理与基于结构模型的基岩储层预测方法相似，因此基于结构模型的基岩储层预测装置的实施可以参见基于结构模型的基岩储层预测方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图10是本发明实施例的基于结构模型的基岩储层预测装置的结构框图，如图10所示，包括：
总孔隙度预测模块1001，用于建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
岩石骨架、流体成分及其弹性参数分析模块1002，用于确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
基岩深层纵波阻抗曲线模拟模块1003，用于根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
纵波阻抗曲线合并模块1004，用于将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
储层特征分析模块1005，用于根据测井数据确定储层测井响应特征和定量解释模板；
层位断层解释模块1006，用于根据地震数据确定层位断层精细解释结果；
稀疏脉冲反演模块1007，用于根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
储层精细预测模块1008，用于根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
下面对该结构进行说明。
在本发明实施例中，所述总孔隙度预测模块1001具体用于：
按照如下方式建立基岩总孔隙度与基岩厚度的定量关系：
获取基岩浅层实测地质数据；
根据所述基岩浅层实测地质数据计算孔隙度曲线；
将所述孔隙度曲线与基岩厚度进行交汇分析，获得分析结果；
根据所述分析结果和获取的区域地质数据，建立基岩孔隙度同基岩厚度之间的定量关系。
在本发明实施例中，所述岩石骨架、流体成分及其弹性参数分析模块1002具体用于：
按照如下方式确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数：
根据测井曲线，确定基岩的矿物成分和孔隙液体成分；
根据完钻井的分析结果确定基岩矿物成分的体积含量；
根据区域分析结果确定孔隙流体的体积模量及密度；
通过纵波时差曲线、横波时差曲线和密度曲线，计算基岩不同组分的弹性参数，其中弹性参数包括体积模量、剪切模量和密度值。
在本发明实施例中，所述基岩深层纵波阻抗曲线模拟模块1003具体用于：
按照如下方式根据所述基岩结构模型预测基岩深层纵波阻抗曲线：
根据所述基岩结构模型进行基岩深层纵波速度和密度的模拟，获得模拟结果；
根据所述模拟结果计算基岩深层纵波阻抗；
根据所述基岩深层纵波阻抗确定基岩深层纵波阻抗曲线。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述基于结构模型的基岩储层预测方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述基于结构模型的基岩储层预测方法的计算机程序。
综上所述，基于结构模型的基岩储层预测方法是根据实测曲线分析基岩矿物成分组成，结合基岩深层的孔隙度预测结果，建立稳定的基岩结构模型，利用岩石物理模拟方法预测基岩深层纵波阻抗曲线，并与浅层实测曲线进行拼接，从而建立更加准确的基岩内部低频模型来约束稀疏脉冲反演运算。在保证浅层预测准确的前提下，更加准确的预测基岩深层的纵波阻抗变化，以取得更加准确的储层预测结果。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等) 上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于结构模型的基岩储层预测方法，其特征在于，包括：
建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
根据测井数据确定储层测井响应特征和定量解释模板；
根据地震数据确定层位断层精细解释结果；
根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
2.如权利要求1所述的基于结构模型的基岩储层预测方法，其特征在于，按照如下方式建立基岩总孔隙度与基岩厚度的定量关系：
获取基岩浅层实测地质数据；
根据所述基岩浅层实测地质数据计算孔隙度曲线；
将所述孔隙度曲线与基岩厚度进行交汇分析，获得分析结果；
根据所述分析结果和获取的区域地质数据，建立基岩孔隙度同基岩厚度之间的定量关系。
3.如权利要求1所述的基于结构模型的基岩储层预测方法，其特征在于，按照如下方式确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数：
根据测井曲线，确定基岩的矿物成分和孔隙液体成分；
根据完钻井的分析结果确定基岩矿物成分的体积含量；
根据区域分析结果确定孔隙流体的体积模量及密度；
通过纵波时差曲线、横波时差曲线和密度曲线，计算基岩不同组分的弹性参数，其中弹性参数包括体积模量、剪切模量和密度值。
4.如权利要求1所述的基于结构模型的基岩储层预测方法，其特征在于，根据所述基岩结构模型预测基岩深层纵波阻抗曲线，包括：
根据所述基岩结构模型进行基岩深层纵波速度和密度的模拟，获得模拟结果；
根据所述模拟结果计算基岩深层纵波阻抗；
根据所述基岩深层纵波阻抗确定基岩深层纵波阻抗曲线。
5.一种基于结构模型的基岩储层预测装置，其特征在于，包括：
总孔隙度预测模块，用于建立基岩总孔隙度与基岩厚度的定量关系，根据所述基岩总孔隙度与基岩厚度的定量关系预测基岩总孔隙度；
岩石骨架、流体成分及其弹性参数分析模块，用于确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数；
基岩深层纵波阻抗曲线模拟模块，用于根据所述基岩总孔隙度、所述基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数，确定基岩结构模型，根据所述基岩结构模型预测基岩深层纵波阻抗曲线；
纵波阻抗曲线合并模块，用于将所述基岩深层纵波阻抗曲线与基岩浅层纵波阻抗曲线拼接，形成基岩内部完整纵波阻抗曲线；
储层特征分析模块，用于根据测井数据确定储层测井响应特征和定量解释模板；
层位断层解释模块，用于根据地震数据确定层位断层精细解释结果；
稀疏脉冲反演模块，用于根据所述测井数据、所述地震数据、基岩内部完整纵波阻抗曲线和层位断层精细解释结果，建立低频模型，根据所述低频模型和根据所述地震数据和所述测井数据提取的子波进行稀疏脉冲反演，获得纵波阻抗反演结果；
储层精细预测模块，用于根据所述纵波阻抗反演结果和所述储层测井响应特征对储层进行定性预测，根据所述纵波阻抗反演结果和定量解释模板对储层进行定量预测。
6.如权利要求5所述的基于结构模型的基岩储层预测装置，其特征在于，所述总孔隙度预测模块具体用于：
按照如下方式建立基岩总孔隙度与基岩厚度的定量关系：
获取基岩浅层实测地质数据；
根据所述基岩浅层实测地质数据计算孔隙度曲线；
将所述孔隙度曲线与基岩厚度进行交汇分析，获得分析结果；
根据所述分析结果和获取的区域地质数据，建立基岩孔隙度同基岩厚度之间的定量关系。
7.如权利要求5所述的基于结构模型的基岩储层预测装置，其特征在于，所述岩石骨架、流体成分及其弹性参数分析模块具体用于：
按照如下方式确定基岩的矿物成分、孔隙液体成分和基岩不同组分的弹性参数：
根据测井曲线，确定基岩的矿物成分和孔隙液体成分；
根据完钻井的分析结果确定基岩矿物成分的体积含量；
根据区域分析结果确定孔隙流体的体积模量及密度；
通过纵波时差曲线、横波时差曲线和密度曲线，计算基岩不同组分的弹性参数，其中弹性参数包括体积模量、剪切模量和密度值。
8.如权利要求5所述的基于结构模型的基岩储层预测装置，其特征在于，所述基岩深层纵波阻抗曲线模拟模块具体用于：
按照如下方式根据所述基岩结构模型预测基岩深层纵波阻抗曲线：
根据所述基岩结构模型进行基岩深层纵波速度和密度的模拟，获得模拟结果；
根据所述模拟结果计算基岩深层纵波阻抗；
根据所述基岩深层纵波阻抗确定基岩深层纵波阻抗曲线。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述基于结构模型的基岩储层预测方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述基于结构模型的基岩储层预测方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811469385.0A CN111239810B (zh) | 2018-11-28 | 2018-11-28 | 基于结构模型的基岩储层预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811469385.0A CN111239810B (zh) | 2018-11-28 | 2018-11-28 | 基于结构模型的基岩储层预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111239810A true CN111239810A (zh) | 2020-06-05 |
| CN111239810B CN111239810B (zh) | 2022-07-05 |
Family
ID=70879337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811469385.0A Active CN111239810B (zh) | 2018-11-28 | 2018-11-28 | 基于结构模型的基岩储层预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111239810B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114415238A (zh) * | 2022-01-25 | 2022-04-29 | 大庆油田有限责任公司 | 同期异相海相碳酸岩地震解释方法及装置、设备和介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020087272A1 (en) * | 2000-12-15 | 2002-07-04 | Dwight Mackie | Method for optimizing migration fields using time slice analysis |
| CN102043167A (zh) * | 2009-10-23 | 2011-05-04 | 中国石油化工股份有限公司 | 一种曲面地表地震资料速度分析的方法 |
| CN105298388A (zh) * | 2015-10-30 | 2016-02-03 | 中国石油天然气集团公司 | 一种水平井钻井轨迹导向方法及装置 |
| CN105301638A (zh) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | 一种提取风化层底界面的方法和装置 |
| CN105756591A (zh) * | 2014-12-17 | 2016-07-13 | 中国石油天然气股份有限公司 | 煤层气井完井方法 |
| CN106569266A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 页岩气储层脆性矿物含量预测方法 |
-
2018
- 2018-11-28 CN CN201811469385.0A patent/CN111239810B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020087272A1 (en) * | 2000-12-15 | 2002-07-04 | Dwight Mackie | Method for optimizing migration fields using time slice analysis |
| CN102043167A (zh) * | 2009-10-23 | 2011-05-04 | 中国石油化工股份有限公司 | 一种曲面地表地震资料速度分析的方法 |
| CN105756591A (zh) * | 2014-12-17 | 2016-07-13 | 中国石油天然气股份有限公司 | 煤层气井完井方法 |
| CN105301638A (zh) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | 一种提取风化层底界面的方法和装置 |
| CN106569266A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 页岩气储层脆性矿物含量预测方法 |
| CN105298388A (zh) * | 2015-10-30 | 2016-02-03 | 中国石油天然气集团公司 | 一种水平井钻井轨迹导向方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 刘家瑾，等: "《煤田测井资料数字处理》", 31 March 1991 * |
| 程平，等: "《现代爆破工程》", 31 October 2018 * |
| 陈国文，等: "基于风化壳结构特征的基岩孔缝储层预测技术", 《CPS/SEG北京2018国际地球物理会议暨展览电子论文集》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114415238A (zh) * | 2022-01-25 | 2022-04-29 | 大庆油田有限责任公司 | 同期异相海相碳酸岩地震解释方法及装置、设备和介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111239810B (zh) | 2022-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10572611B2 (en) | Method and system for characterizing fractures in a subsurface region | |
| US10712472B2 (en) | Method and system for forming and using a subsurface model in hydrocarbon operations | |
| RU2669948C2 (ru) | Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности | |
| EP2497900B1 (en) | Modeling hydraulic fractures | |
| US8712746B2 (en) | Method of modified facies proportions upon history matching of a geological model | |
| US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| US8868390B2 (en) | Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling | |
| WO2017041074A1 (en) | Method of integrating fracture, production, and reservoir operations into geomechanical operations of a wellsite | |
| WO2017027340A1 (en) | Method integrating fracture and reservoir operations into geomechanical operations of a wellsite | |
| US20100185424A1 (en) | Method, Program and Computer System for Conciliating Hydrocarbon Reservoir Model Data | |
| CN107558999B (zh) | 一种孤立缝洞体剩余油定量预测方法 | |
| EP3394647B1 (en) | A methodology for building realistic numerical forward stratigraphic models in data sparse environment | |
| CN110632657B (zh) | 泥岩涂抹型断层封闭性分析方法及装置 | |
| GB2595809A (en) | Method, apparatus, and system for quantitative analysis of sealing strength of trap fault | |
| US20190204464A1 (en) | Method and System for Modeling in a Subsurface Region | |
| US20190203593A1 (en) | Method and System for Modeling in a Subsurface Region | |
| US10605955B2 (en) | Multi-step subsidence inversion for modeling lithospheric layer thickness through geological time | |
| US10337294B2 (en) | Reservoir permeability upscaling | |
| US11555937B2 (en) | Method and system for generating simulation grids by mapping a grid from the design space | |
| CN111239810B (zh) | 基于结构模型的基岩储层预测方法及装置 | |
| Gherabati et al. | New engineering and geologic parameters to predict infill well performance in the Wolfcamp of the Delaware Basin | |
| Skjervheim et al. | Integrated workflow for consistent model building from depth conversion to flow simulation-North Sea field case | |
| CN104246538A (zh) | 用于地下储藏表征的系统和方法 | |
| CN107290803A (zh) | 建立基于改进Pride-Lee模型的岩石物理新模型的方法 | |
| CN113031056B (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 |