CN108318936B - 一种地层划分处理方法和装置 - Google Patents

一种地层划分处理方法和装置 Download PDF

Info

Publication number
CN108318936B
CN108318936B CN201711440328.5A CN201711440328A CN108318936B CN 108318936 B CN108318936 B CN 108318936B CN 201711440328 A CN201711440328 A CN 201711440328A CN 108318936 B CN108318936 B CN 108318936B
Authority
CN
China
Prior art keywords
data
density
target area
skeleton
cycle
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
Application number
CN201711440328.5A
Other languages
English (en)
Other versions
CN108318936A (zh
Inventor
李佳鸿
邓西里
蔚涛
胡亚斐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201711440328.5A priority Critical patent/CN108318936B/zh
Publication of CN108318936A publication Critical patent/CN108318936A/zh
Application granted granted Critical
Publication of CN108318936B publication Critical patent/CN108318936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例提供一种地层划分处理方法和装置。该方法包括：确定目标区的地层孔隙度数据；根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。利用本申请实施例提供的技术方案可以准确划分地层，为油田区域内油层的细分和对比、开发层系的划分等提供数据支持。

Description

一种地层划分处理方法和装置
技术领域
本发明涉及油气勘探技术领域，尤其涉及一种地层划分处理方法和装置。
背景技术
在油气勘探技术领域，油田区域内油层的细分和对比、开发层系的划分等是在地层划分对比的基础上进行的。因此，地层划分是石油地质研究的重要内容。正确的地层划分是试油勘探与开发工作的基础。
现有技术中，可以基于古生物化石对比法、岩石沉积特征对比法、岩层接触关系对比法和重矿物对比法进行地层划分处理。其中，应用比较广泛的是岩石沉积特征对比法中的旋回对比法。地壳周期性升降运动引起了地壳上海水进退和沉积环境相似重复的变化，导致岩石性质在纵向上有规律的变化，反映在地层剖面上，岩性呈有规律重复出现的现象就叫做沉积岩的旋回性，旋回由大到小可以划分为若干级。具体的，旋回对比法通过测井曲线来反映泥质含量，泥质含量的增加与减少反应沉积物能量的大小，进而划分正反旋回，然后，可以根据旋回的正反级次来划分地层。但现有技术中直接基于测井数据反映的泥质含量来进行正反旋回的划分，进而进行地层划分处理的只能适用于砂泥岩变化明显的地层的划分，对于冲积扇、扇三角洲等成因的快速堆积形成的块状砾岩、砂砾岩沉积体，泥质含量低于10％，测井响应区分度低，无法准确划分地层。因此，现有技术中亟需一种可以准确划分地层的地层划分处理方法。
发明内容
本申请的目的是提供一种地层划分处理方法和装置，可以准确划分地层，为油田区域内油层的细分和对比、开发层系的划分等提供数据支持。
本申请提供的地层划分处理方法和装置是这样实现的：
一种地层划分处理方法，所述方法包括：
确定目标区的地层孔隙度数据；
根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；
基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；
基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
在一个优选的实施例中，所述基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据包括：
基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段；
确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；
确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；
将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；
其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比。
在一个优选的实施例中，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度；
其中，所述密度测井数据包括体积密度数据和流体密度数据。
在一个优选的实施例中，所述确定目标区的地层孔隙度数据包括：
获取所述目标区的声波时差测井曲线；
利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
在一个优选的实施例中，所述确定目标区的地层孔隙度数据包括：
获取所述目标区的中值密度测井曲线；
利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
一种地层划分处理装置，所述装置包括：
地层孔隙度数据确定模块，用于确定目标区的地层孔隙度数据；
骨架密度数据计算模块，用于根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；
旋回划分处理模块，用于基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；
地层划分处理模块，用于基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
在一个优选的实施例中，所述旋回划分处理模块包括：
骨架密度段划分单元，用于基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段；
正旋回数据确定单元，用于确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；
反旋回数据确定单元，用于确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；
旋回数据确定单元，用于将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；
其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比。
在一个优选的实施例中，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度；
其中，所述密度测井数据包括体积密度数据和流体密度数据。
在一个优选的实施例中，所述地层孔隙度数据确定模块包括：
声波时差测井曲线获取单元，用于获取所述目标区的声波时差测井曲线；
第一计算单元，用于利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
在一个优选的实施例中，所述地层孔隙度数据确定模块包括：
中值密度测井曲线获取单元，用于获取所述目标区的中值密度测井曲线；
第二计算单元，用于利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
本申请提供的技术方案可以确定目标区的地层孔隙度数据；并根据目标区的孔隙度数据和密度测井数据计算出目标区的骨架密度数据；然后，基于骨架密度数据的数值变化趋势数据进行目标区的沉积岩的旋回划分处理，得到与目标区的沉积岩相对应的旋回数据；基于旋回数据对沉积岩所对应的地层进行划分处理，可以准确的对地层进行划分。与现有技术相比，利用本申请提供的技术方案可以准确划分地层，为油田区域内油层的细分和对比、开发层系的划分等提供数据支持。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的地层划分处理方法的一种实施例的流程图；
图2是本申请提供的测井曲线和骨架密度曲线对比示意图；
图3是本申请提供的基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据的一种实施例的流程示意图；
图4是本申请提供的测井的测井数据、骨架密度和旋回数据的一种实施例的部分示意图；
图5是本申请提供的地层划分处理装置的一种实施例中的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
以下首先介绍本申请一种地层划分处理方法的一种实施例。图1是本申请提供的地层划分处理方法的一种实施例的流程图，本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的如图1所示，所述方法可以包括：
S110：确定目标区的地层孔隙度数据。
具体的，本申请实施例中，确定目标区的地层孔隙度数据可以结合实际应用情况选取适合当前目标区的测井曲线来确定目标区的地层孔隙度数据。具体的，可以结合目标区的测井曲线与已知的地层孔隙度进行拟合来确定适合当前目标区的测井曲线。
在一个具体的实施例中，所述确定目标区的地层孔隙度数据可以包括：
获取所述目标区的声波时差测井曲线；
利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
具体的实施例中，所述利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据可以包括采用下述计算公式计算：
Φ＝(AC-a)/(b-a)
上式中，Φ表示地层孔隙度；AC表示目标区的声波时差数据，可以从钻井数据中获取，单位为μs/m；a表示岩石骨架声波时差，可以从钻井数据中获取，单位为μs/m，在一个具体的实施例中，a＝182，单位为μs/m；b表示岩石孔隙流体的声波时差，可以从钻井数据中获取，单位为μs/m；在一个具体的实施例中，b＝620。
在另一个具体的实施例，所述确定目标区的地层孔隙度数据可以包括：
获取所述目标区的中值密度测井曲线；
利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
此外，需要说明的是，本申请实施例中确定目标区的地层孔隙度数据并不仅限于上述的方法，在实际应用中还可以包括其他方法，本申请实施例并不以上述为限。
S120：根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据。
本申请实施例中，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据可以包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度。
具体的，在实际应用中，对地下储层(即本申请实施例所述地层)而言，流体密度可以取1g/m3。另外，在实际应用中，对于干燥的岩心样品，流体密度一般取值为0。所述密度测井数据可以在目标区的测井数据中获取，所述密度测井数据可以包括体积密度数据和流体密度数据。
在一个具体的实施例中，如图2所示，图2是本申请提供的测井曲线和骨架密度曲线对比示意图。图2中包括自然伽玛测井曲线GR，自然电位测井曲线SP，电阻率测井曲线(RXO、RI和RT)，密度测井曲线RHOB，声波时差测井曲线AC，中子测井曲线CNL，放射性钾测井曲线k，放射性钍测井曲线th，放射性铀测井曲线u以及骨架密度曲线。其中，伽马能谱测井显示该井在2770m处有一沉积间断面，间断面之下表现为高铀低钾，之上则表现为高钾低铀。沉积间断面由于暴露地表，通常风化、淋滤形成质地坚硬的风化壳，从而使骨架密度值呈现局部的极大值，因此，仅在骨架密度曲线上有响应，其他测井曲线对该沉积间断面无法识别。因此，本申请实施例利用骨架密度来进行旋回划分，可以更准确的进行地层划分。
S130：基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据。
本申请实施例中，在步骤S120得到骨架密度数据之后，可以基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据，具体的，图3是本申请提供的基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据的一种实施例的流程示意图，如图3所示，可以包括：
S131：基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段。
具体的，这里可以划分到多个而不同取值区间的骨架密度段。
S133：确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；
S135：确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；
S137：将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；
其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比。即，所述骨架密度段的取值区间越大，相应的旋回数据的越大；反之，所述骨架密度段的取值区间越小，相应的旋回数据的越小。
在一个具体的实施例中，如图4所示，图4是本申请提供的测井的测井数据、骨架密度和旋回数据的一种实施例的部分示意图。图4中干燥岩心骨架密度的曲线是流体密度取值为0时，根据上述计算骨架密度数据的计算公式计算确定的；而地层骨架密度的曲线则由密度测井曲线RHOB和声波时差测井曲线AC基于流体密度取值为1g/m3时，根据上述计算骨架密度数据的计算公式计算确定的。砾石成分、砂质成分和泥质含量是岩心薄片分析的结果。砾石含量与砂质和泥质含量呈现此消彼长的规律。干燥岩心骨架密度和地层骨架密度随砾石含量的降低、砂质和泥质含量的增加逐渐变大。依据骨架密度的数值变化趋势数据，可以将该地层划分为两个旋回，包括一个正旋回和一个反旋回。
S140：基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
本申请实施例中，在步骤S130得到多个旋回数据之后，可以基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。具体的可以基于经典层序地层学原理进行井间对比，最终建立地层格架。
由以上本申请一种地层划分处理方法的实施例可见，本申请通过确定目标区的地层孔隙度数据；并根据目标区的孔隙度数据和密度测井数据计算出目标区的骨架密度数据；然后，基于骨架密度数据的数值变化趋势数据进行目标区的沉积岩的旋回划分处理，得到与目标区的沉积岩相对应的旋回数据；基于旋回数据对沉积岩所对应的地层进行划分处理，可以准确的对地层进行划分。与现有技术相比，利用本申请提供的技术方案可以准确划分地层，为油田区域内油层的细分和对比、开发层系的划分等提供数据支持。
本申请另一方面还提供一种地层划分处理装置，图5是本申请提供的地层划分处理装置的一种实施例中的结构示意图；如图5所示，所述装置500可以包括：
地层孔隙度数据确定模块510，可以用于确定目标区的地层孔隙度数据；
骨架密度数据计算模块520，可以用于根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；
旋回划分处理模块530，可以用于基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；
地层划分处理模块540，可以用于基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
在一个优选的实施例中，所述旋回划分处理模块530可以包括：
骨架密度段划分单元，可以用于基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段；
正旋回数据确定单元，可以用于确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；
反旋回数据确定单元，可以用于确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；
旋回数据确定单元，可以用于将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；
其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比。
在一个优选的实施例中，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度；
其中，所述密度测井数据包括体积密度数据和流体密度数据。
在一个优选的实施例中，所述地层孔隙度数据确定模块510可以包括：
声波时差测井曲线获取单元，可以用于获取所述目标区的声波时差测井曲线；
第一计算单元，可以用于利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
在一个优选的实施例中，所述地层孔隙度数据确定模块510可以包括：
中值密度测井曲线获取单元，可以用于获取所述目标区的中值密度测井曲线；
第二计算单元，可以用于利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
由以上本申请一种地层划分处理方法或装置的实施例可见，本申请可以确定目标区的地层孔隙度数据；并根据目标区的孔隙度数据和密度测井数据计算出目标区的骨架密度数据；然后，基于骨架密度数据的数值变化趋势数据进行目标区的沉积岩的旋回划分处理，得到与目标区的沉积岩相对应的旋回数据；基于旋回数据对沉积岩所对应的地层进行划分处理，可以准确的对地层进行划分。与现有技术相比，利用本申请提供的技术方案可以准确划分地层，为油田区域内油层的细分和对比、开发层系的划分等提供数据支持。
本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (8)

1.一种地层划分处理方法，其特征在于，所述方法包括：
确定目标区的地层孔隙度数据；
根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；
基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；包括：基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段；确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比；
基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
2.根据权利要求1所述的方法，其特征在于，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度；
其中，所述密度测井数据包括体积密度数据和流体密度数据。
3.根据权利要求1所述的方法，其特征在于，所述确定目标区的地层孔隙度数据包括：
获取所述目标区的声波时差测井曲线；
利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
4.根据权利要求1所述的方法，其特征在于，所述确定目标区的地层孔隙度数据包括：
获取所述目标区的中值密度测井曲线；
利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
5.一种地层划分处理装置，其特征在于，所述装置包括：
地层孔隙度数据确定模块，用于确定目标区的地层孔隙度数据；
骨架密度数据计算模块，用于根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据；
旋回划分处理模块，用于基于所述骨架密度数据的数值变化趋势数据进行所述目标区的沉积岩的旋回划分处理，得到与所述目标区的沉积岩相对应的旋回数据；所述旋回划分处理模块包括：骨架密度段划分单元，用于基于所述骨架密度数据的数值变化趋势数据，将具有相同数值变化趋势数据的骨架密度数据划分一个骨架密度段；正旋回数据确定单元，用于确定数值变化趋势数据为依次变大的骨架密度段所对应的旋回数据为正旋回数据；反旋回数据确定单元，用于确定数值变化趋势数据为依次变小的骨架密度段所对应的旋回数据为反旋回数据；旋回数据确定单元，用于将所述正旋回数据和所述反旋回数据作为所述目标区的沉积岩相对应的旋回数据；其中，所述旋回数据包括多个不同级别的正旋回数据和反旋回数据，所述旋回数据的大小与相应的骨架密度段的取值区间的大小成正比；
地层划分处理模块，用于基于所述旋回数据对所述沉积岩所对应的地层进行划分处理，得到划分后的地层。
6.根据权利要求5所述的装置，其特征在于，所述根据所述目标区的孔隙度数据和密度测井数据计算所述目标区的骨架密度数据包括采用如下计算公式计算：
DENmatrix＝(DENvolume–DENfluid*Φ)/(1-Φ)
上式中，DENmatrix表示骨架密度，单位为g/m3；DENvolume表示体积密度，单位为g/m3；DENfluid表示流体密度，单位为g/m3；Φ表示地层孔隙度；
其中，所述密度测井数据包括体积密度数据和流体密度数据。
7.根据权利要求5所述的装置，其特征在于，所述地层孔隙度数据确定模块包括：
声波时差测井曲线获取单元，用于获取所述目标区的声波时差测井曲线；
第一计算单元，用于利用所述声波时差测井曲线计算所述目标区的地层孔隙度数据。
8.根据权利要求5所述的装置，其特征在于，所述地层孔隙度数据确定模块包括：
中值密度测井曲线获取单元，用于获取所述目标区的中值密度测井曲线；
第二计算单元，用于利用所述中值密度测井曲线计算所述目标区的地层孔隙度数据。
CN201711440328.5A 2017-12-27 2017-12-27 一种地层划分处理方法和装置 Active CN108318936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711440328.5A CN108318936B (zh) 2017-12-27 2017-12-27 一种地层划分处理方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711440328.5A CN108318936B (zh) 2017-12-27 2017-12-27 一种地层划分处理方法和装置

Publications (2)

Publication Number Publication Date
CN108318936A CN108318936A (zh) 2018-07-24
CN108318936B true CN108318936B (zh) 2020-02-14

Family

ID=62893604

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711440328.5A Active CN108318936B (zh) 2017-12-27 2017-12-27 一种地层划分处理方法和装置

Country Status (1)

Country Link
CN (1) CN108318936B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111399054B (zh) * 2019-01-03 2023-06-30 中国石油天然气股份有限公司 井震结合识别重复地层的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002046563A2 (fr) * 2000-12-04 2002-06-13 Innov-Pro Procede et dispositif pour determiner la qualite de reserve d'un puits de petrole
CN102313904A (zh) * 2011-05-05 2012-01-11 中国石油化工股份有限公司 利用地震信息自动划分地层旋回的方法
CN103645519A (zh) * 2013-12-17 2014-03-19 中国石油大学(华东) 一种基于火山岩风化壳识别及划分标准的油气勘探方法
CN105372716A (zh) * 2015-10-28 2016-03-02 中国石油大学(华东) 碳酸盐岩表生岩溶储层分布的评价方法
CN106707337A (zh) * 2016-11-15 2017-05-24 西安石油大学 基于分频地震反射能量地层切片进行沉积体系解释的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002046563A2 (fr) * 2000-12-04 2002-06-13 Innov-Pro Procede et dispositif pour determiner la qualite de reserve d'un puits de petrole
CN102313904A (zh) * 2011-05-05 2012-01-11 中国石油化工股份有限公司 利用地震信息自动划分地层旋回的方法
CN103645519A (zh) * 2013-12-17 2014-03-19 中国石油大学(华东) 一种基于火山岩风化壳识别及划分标准的油气勘探方法
CN105372716A (zh) * 2015-10-28 2016-03-02 中国石油大学(华东) 碳酸盐岩表生岩溶储层分布的评价方法
CN106707337A (zh) * 2016-11-15 2017-05-24 西安石油大学 基于分频地震反射能量地层切片进行沉积体系解释的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
测井曲线的沉积旋回划分方法及应用;李新虎 等;《西安科技大学学报》;20101130;第30卷(第6期);第702-705页 *

Also Published As

Publication number Publication date
CN108318936A (zh) 2018-07-24

Similar Documents

Publication Publication Date Title
CN109653725B (zh) 基于沉积微相和岩石相的混积储层水淹程度测井解释方法
US10048403B2 (en) Method and system for generation of upscaled mechanical stratigraphy from petrophysical measurements
CN110501744B (zh) 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质
Cheng et al. Field experiences with assisted and automatic history matching using streamline models
Lee et al. Subsidence analysis and visualization: for sedimentary basin analysis and modelling
Jolley et al. Structurally complex reservoirs: an introduction
CN110609319B (zh) 一种利用沉积正演模拟建立时间域层序地层剖面的方法
CN110632657B (zh) 泥岩涂抹型断层封闭性分析方法及装置
CN113050157A (zh) 一种基于露头资料的碳酸盐岩地震储层反演方法及系统
Wei et al. Porosity model and its application in tight gas sandstone reservoir in the southern part of West Depression, Liaohe Basin, China
CN108318936B (zh) 一种地层划分处理方法和装置
Leary et al. Power-law scaling of spatially correlated porosity and log (permeability) sequences from north-central North Sea Brae oilfield well core
CN108073744B (zh) 基于地质统计学的浊积扇推进距离预测方法
Li et al. Chemical weathering characteristics of the Late Cretaceous Nenjiang Formation from the Songliao Basin (Northeastern China) reveal prominent Milankovitch band variations
CN112649887B (zh) 基于钻井资料定量划分层序地层的方法及装置
CN110688728B (zh) 一种时间域内一维沉积物沉积特征及水体环境定量分析方法
Tebo et al. Use of volume-based 3-D seismic attribute analysis to characterize physical-property distribution: a case study to delineate sedimentologic heterogeneity at the Appleton field, southwestern Alabama, USA
Irofti et al. Characterization of a Tight Gas Reservoir Using the Integration of Electrofacies and Fracture Aperture, Ahnet, Algeria
Habrat et al. The concept of a computer system for interpretation of tight rocks using X-ray computed tomography results
Koch et al. Does fluvial channel-belt clustering predict net sand to gross rock volume? Architectural metrics and point-pattern analysis of a digital outcrop model
Nwafor et al. Geostatistical inversion of spectrally broadened seismic data for re-evaluation of oil reservoir continuity in Inas field, offshore Malay Basin
CN113138412A (zh) 深层页岩气孔隙度地震预测方法及装置
CN112305594A (zh) 非均质储层的油气分布确定方法及系统
Li et al. Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China
CN111708100B (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