CN111680386B - 断块圈闭的脆性盖层封闭性分析方法及装置 - Google Patents

断块圈闭的脆性盖层封闭性分析方法及装置 Download PDF

Info

Publication number
CN111680386B
CN111680386B CN202010304303.8A CN202010304303A CN111680386B CN 111680386 B CN111680386 B CN 111680386B CN 202010304303 A CN202010304303 A CN 202010304303A CN 111680386 B CN111680386 B CN 111680386B
Authority
CN
China
Prior art keywords
section
layer
interval
sandstone
factor
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
CN202010304303.8A
Other languages
English (en)
Other versions
CN111680386A (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 CN202010304303.8A priority Critical patent/CN111680386B/zh
Publication of CN111680386A publication Critical patent/CN111680386A/zh
Application granted granted Critical
Publication of CN111680386B publication Critical patent/CN111680386B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供了一种断块圈闭的脆性盖层封闭性分析方法及装置，该方法包括：在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。本发明可以在断面上分析脆性盖层封闭性，准确度高。

Description

断块圈闭的脆性盖层封闭性分析方法及装置
技术领域
本发明涉及地震勘探领域，尤其涉及一种断块圈闭的脆性盖层封闭性分析方法及装置。
背景技术
在断块油气藏勘探中，分析断层对盖层的破坏程度是重要的要素之一。断裂对不同脆韧性的盖层的破坏机理有着较大差异，导致断裂与盖层组合后垂向破坏盖层的能力不同，因此，在分析断层对盖层的破坏及其垂向调整油气能力时，要针对不同脆韧性盖层采取不同的评价方法。对于处于韧性阶段的盖层，其在断裂面上一般是具备封闭能力的；对于脆-韧性过渡阶段的盖层，采用泥岩涂抹因子来评价泥质岩盖层被断裂破坏的程度；对于脆性阶段的盖层，通常采用断接厚度法，但该方法的表示形式是基于大量钻井统计的图表或者是在某一地层平面上的展布，没有实现其在断层面上的展布显示，尤其是无法表现断面上不同位置盖层封闭能力的变化，使得封闭性分析准确度低。所以，上述方法无法直观的反应出脆性盖层在断面上的封闭差异，在应用时具有局限性。
发明内容
本发明实施例提出一种断块圈闭的脆性盖层封闭性分析方法，用以在断面上分析脆性盖层封闭性，准确度高，该方法包括：
在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。
本发明实施例提出一种断块圈闭的脆性盖层封闭性分析装置，用以在断面上分析脆性盖层封闭性，准确度高，该装置包括：
断层模型构建模块，用于在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
盖层封闭因子确定模块，用于根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
封闭性分析模块，用于根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。
本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述断块圈闭的脆性盖层封闭性分析方法。
本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述断块圈闭的脆性盖层封闭性分析方法的计算机程序。
在本发明实施例中，在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。在上述过程中，引入了预设的砂岩过滤窗尺寸来确定每一层段的盖层封闭因子的计算方式，体现了盖层段中存在的砂岩层对盖层封闭性的影响，使得最后得到的目标脆性盖层段的盖层封闭性更加准确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中断块圈闭的脆性盖层封闭性分析方法的流程图；
图2为本发明实施例中断块圈闭的脆性盖层封闭性分析方法的详细流程图；
图3为本发明实施例中是渤海湾盆地某一断面上地层和岩性的投影关系；
图4为本发明实施例中同时采用两种计算方式获得的盖层封闭因子的示意图；
图5为本发明实施例中只采用第二计算方式获得的盖层封闭因子的示意图；
图6为本发明实施例中断块圈闭的脆性盖层封闭性分析装置的示意图；
图7为本发明实施例中计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
在本说明书的描述中，所使用的“包含”、“包括”、“具有”、“含有”等，均为开放性的用语，即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施，其中的步骤顺序不作限定，可根据需要作适当调整。
图1为本发明实施例中断块圈闭的脆性盖层封闭性分析方法的流程图，如图1所示，该方法包括：
步骤101，在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
步骤102，根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
步骤103，根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。
在本发明实施例提出的方法中，引入了预设的砂岩过滤窗尺寸来确定每一层段的盖层封闭因子的计算方式，体现了盖层段中存在的砂岩层对盖层封闭性的影响，使得最后得到的目标脆性盖层段的盖层封闭性更加准确。
具体实施时，根据地震勘探数据和钻井数据，建立目标脆性盖层段的断层模型，上述断层模型包括若干个砂岩层段和/或泥岩层段这两种岩性属性，断层模型包括断面，有了断层模型可以获得断面上多个采样点的厚度和断距，从断层模型中，还可以确定目标脆性盖层的顶面、底面与断面的交线位置，用于后续的盖层封闭因子的计算。
在一实施例中，所述岩性属性包括砂岩层段和泥岩层段；
确定断层模型中每一层段的岩性属性，包括：
在断层模型中，通过标记砂岩层段在断面上的位置，确定岩性属性为砂岩层段的层段。
在上述实施例中，通过标记的方式可快速地确定砂岩层段，而非砂岩层段即为泥岩层段。
在一实施例中，根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子，包括：
对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子；否则确定采用第二计算方式计算该层段的盖层封闭因子。
在上述实施例中，预设的砂岩过滤窗尺寸SWL可以根据实际研究和应用情况确定，上述方式确定了两种计算方式，可以整理为：
该层段的岩性属性为砂岩层段时，且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子，以此方式是为了在断面上突出盖层中存在的较厚层砂岩段对盖层封闭性的影响，是对盖层中可能的渗漏区域在断面上的一种展示；否则确定采用第二计算方式计算该层段的盖层封闭因子。
该层段的岩性属性为泥岩层段时，确定采用第二计算方式计算该层段的盖层封闭因子。
上述分析考虑了砂岩过滤窗尺寸，并以此分开计算每一层段的盖层封闭因子，使得计算的盖层封闭因子更符合实际情况，准确度更高。一般而言，在断面上纵向上定义的第i个采样点的采样长度ΔHi要比砂岩过滤窗尺寸SWL小。
在一实施例中，所述第一计算方式为：在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，采用如下公式计算每一层段的盖层封闭因子：
Figure BDA0002455175680000051
其中，FJTi为断面上第i个采样点的盖层封闭因子；
N为断面在纵向上划分的采样点的数量；
SWFJTj为断面纵向上第j个连续砂岩层段的厚度；
Sk为断层模型的断面上第j个连续砂岩层段内第k个采样点的厚度；
m为第j个连续砂岩层段中纵向包含的采样点个数；
ΔHi为断面上每个层段中纵向第i个采样点的采样长度；
Di为断面纵向上每个层段中第i个采样点的断距。
在上述实施例中，
Figure BDA0002455175680000052
为单一连续砂岩砂岩层段的厚度，即在为第j个层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，采用公式(1)计算的数值将作为第j个层段的第k个采样点的盖层封闭因子。该算法考虑了盖层中砂岩层段对盖层封闭能力的影响，反应出脆性盖层在断面上的封闭差异性。
在一实施例中，根据该地层中每个采样点的厚度和断距，所述第二计算方式为：若该层段的岩性属性为砂岩层段且该层段的厚度不大于所述砂岩过滤窗尺寸，或该层段的岩性属性为泥岩层段，采用如下公式计算每个层段的盖层封闭因子：
Figure BDA0002455175680000053
其中，FJTi为断面上每个层段中纵向上第i个采样点的盖层封闭因子；
ΔHi为断面纵向上每个层段中第i个采样点的采样长度；
N为断面在纵向上划分的采样点的数量；
Di为断面纵向上每个层段中第i个采样点的断距。
在上述实施例中，第i个采样点的厚度和断距可以从断层模型中获得，也就是说在每个层段的岩性属性为砂岩层段且厚度不大于所述砂岩过滤窗尺寸，可以视为此砂岩层段被设置的过滤窗剔除，进而将其作为泥岩层段处理，计算盖层封闭因子，从而提高统计效率，满足工业界的应用需求。在层段为泥岩层时，也用上述公式计算。国内吕延防、付广等学者研究表明脆性盖层段的厚度与盖层的封闭性呈正相关关系，盖层厚度的减小，意味着盖层封闭性的降低，盖层一旦被断层错开以后，其横向的连续性和完整性遭到破坏。一旦有断层发育，对于盖层而言，起封闭作用的盖层厚度就不再是原盖层厚度，而是断层两盘盖层与盖层的对接厚度，用这一厚度作为盖层分析的方法被称为盖层的断接厚度分析方法，但该方法的表现形式是基于大量钻井统计的图表或者是在某一地层平面上的展布，没有实现其在断层面上的展布显示，尤其是无法表现断面上不同位置盖层封闭能力的变化，使得封闭性分析在断面上的准确度低，本实施例中公式(2)表达的恰是在断面上开展盖层分析的方法。
在一实施例中，根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析，包括：
对于每一层段的每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭。
在上述实施例中，可以在三维断面上将砂岩层段和泥岩层段的盖层封闭因子进行绘图显示，得到三维断面图，并根据三维断面图，确定目标脆性盖层段封闭因子阈值。
基于上述实施例，本发明提出如下一个实施例来说明断块圈闭的脆性盖层封闭性分析方法的详细流程，图2为本发明实施例中断块圈闭的脆性盖层封闭性分析方法的详细流程图，如图2所示，包括：
步骤201，在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型；
步骤202，在断层模型中，通过标记砂岩层段在断面上的位置，确定岩性属性为砂岩层段的层段；
步骤203，对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子；否则确定采用第二计算方式计算该层段的盖层封闭因子；基于所述计算方式，计算每一层段的盖层封闭因子；
步骤204，确定目标脆性盖层段封闭因子阈值；
步骤205，对于每一层段的每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭。
当然，可以理解的是，上述详细流程还可以有其他变化例，相关变化例均应落入本发明的保护范围。
下面给出一具体实施例，来说明断块圈闭的脆性盖层封闭性分析方法的具体应用。
本实施例以渤海湾盆地中的某拉张型正断层为例，首先根据地震勘探数据和钻井数据，建立渤海湾盆地某拉张型正断层的断层模型，并计算出断面上每个采样点的断距，还可以获得地层模型，与断层模型耦合，确定上述断层模型中盖层在断面上的空间投影窗口位置，包括盖层的顶面、底面与断面的交线位置。
将断层模型中砂岩层段标记出来，碎屑岩地层中砂泥岩互层地层是常见的，即使在厚层的泥岩层段中也会发育砂岩层段，图3为本发明实施例中是渤海湾盆地某1号断面上地层和岩性的投影关系，其中，图3中的(a)和(b)分别是不同剖面的示意图，B地层是以砂岩层段为主的储集层，A地层是以泥岩层段为主的盖层，在A地层内发育4层砂岩层段，分别为砂岩层段S1、砂岩层段S2、砂岩层段S3、砂岩层段S4，在平行于断面方向上A地层顶面距离B地层顶面的距离为360m。分析盖层封闭性时，需要考虑盖层段内砂岩层段的影响：一般而言，砂岩层段越厚，其对断面上盖层封闭性的影响越大，反之则越小。因此，在分析断层面上盖层时有必要设置一个砂岩过滤窗尺寸，对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所设定的砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子，突出了较厚的连续砂岩的影响，体现了断面上砂岩层处盖层的封闭能力与泥岩层的不同之处；否则，表示盖层中砂岩层段因厚度小，对盖层的影响很小，可以通过砂岩过滤窗滤除，可以视为泥岩统计，确定采用第二计算方式计算该层段的盖层封闭因子。基于所述计算方式，计算每一层段的盖层封闭因子。
表1为本发明实施例中图3中的断层模型的某一采样点的盖层封闭因子的示例。
表1断层模型的某一采样点的盖层封闭因子的示例
Figure BDA0002455175680000071
Figure BDA0002455175680000081
如表1所示，假设给定的砂岩过滤窗尺寸为12m，表1中砂岩层段S2的厚度为10m，这时候采用公式(2)计算，计算结果为180m，而其他砂岩层段的厚度均大于砂岩过滤窗尺寸12m，适用公式(1)来计算。
泥岩层段则直接采用公式(2)进行计算，此处的计算是在断面上盖层段内进行的。
在三维断面上将砂岩层段和泥岩层段的盖层封闭因子进行绘图显示，得到三维断面图，在本发明实施例中，以渤海湾某2号断层为例，该整体上该断面的上部以泥岩层段为主，下部以砂岩层段为主，图4为本发明实施例中同时采用两种计算方式获得的盖层封闭因子的示意图，此时砂岩过滤窗尺寸设置为10m时，图4中，整体表现出断面上左右两侧盖层封闭能力增强，中间减弱，上部强、下部弱的特征；而且在纵向细节上具备明显的分辨优势，突出了断面纵向上岩性对盖层封闭能力的影响，以三维图的方式展现出断面上不同位置盖层的封闭能力的差异性，符合地质认识。
为了宏观的显示整个断面上盖层的封闭能力，方便地质研究人员在顺着断层走向方向上快速判断其封闭性能，绘图时将纵向上岩性对盖层封闭能力的影响“隐藏”起来，只采用第二种计算方式计算盖层封闭因子，图5为本发明实施例中只采用第二方式获得的盖层封闭因子的示意图，其在横向上展示的规律也是左右两侧盖层封闭能力增强，中间减弱，只是纵向上无法呈现更多的细节。实际应用中，根据需求，上述两种绘图显示方式是互相参考的。
根据图4和图5，确定目标脆性盖层段封闭因子阈值，一般而言，目标脆性盖层段封闭因子阈值是根据区块不同而不同的，还可以根据已钻井数据确定出所统计的样本中盖层封闭因子的下限阈值，作为目标脆性盖层段封闭因子阈值，于每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭；且采样点的盖层封闭因子与目标脆性盖层段封闭因子阈值的差值越大，盖层在断面上的封闭性越强，反之越弱。
综上所述，在本发明实施例提出的方法中，在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。在上述过程中，引入了预设的砂岩过滤窗尺寸来确定每一层段的盖层封闭因子的计算方式，体现了盖层段中存在的砂岩层对盖层封闭性的影响，使得最后得到的目标脆性盖层段的盖层封闭性更加准确。
基于同样的原理，本发明实施例还给出一种断块圈闭的脆性盖层封闭性分析装置，其原理不再赘述，图6为本发明实施例中断块圈闭的脆性盖层封闭性分析装置的示意图，如图6所示，该装置包括：
断层模型构建模块601，用于在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
盖层封闭因子确定模块602，用于根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
封闭性分析模块603，用于根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。
在一实施例中，所述岩性属性包括砂岩层段和泥岩层段；
断层模型构建模块601具体用于：
在断层模型中，通过标记砂岩层段在断面上的位置，确定岩性属性为砂岩层段的层段。
在一实施例中，盖层封闭因子确定模块602具体用于：
对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子；否则确定采用第二计算方式计算该层段的盖层封闭因子。
在一实施例中，所述第一计算方式为：在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，采用如下公式计算每一层段的盖层封闭因子：
Figure BDA0002455175680000101
Figure BDA0002455175680000102
其中，FJTi为断面上第i个采样点的盖层封闭因子；
N为断面在纵向上划分的采样点的数量；
SWFJTj为断面纵向上第j个连续砂岩层段的厚度；
Sk为断层模型的断面上第j个连续砂岩层段内第k个采样点的厚度；
m为第j个连续砂岩层段中纵向包含的采样点个数；
ΔHi为断面上每个层段中纵向第i个采样点的采样长度；
Di为断面纵向上每个层段中第i个采样点的断距。
在一实施例中，所述第二计算方式为：若该层段的岩性属性为砂岩层段且该层段的厚度不大于所述砂岩过滤窗尺寸，或该层段的岩性属性为泥岩层段，采用如下公式计算每个层段的盖层封闭因子；
Figure BDA0002455175680000103
其中，FJTi为断面上每个层段中纵向上第i个采样点的盖层封闭因子；
ΔHi为断面纵向上每个层段中第i个采样点的采样长度；
N为断面在纵向上划分的采样点的数量；
Di为断面纵向上每个层段中第i个采样点的断距。
在一实施例中，封闭性分析模块603具体用于：
对于每一层段的每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭。
综上所述，在本发明实施例提出的装置中，在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析。在上述过程中，引入了预设的砂岩过滤窗尺寸来确定每一层段的盖层封闭因子的计算方式，体现了盖层段中存在的砂岩层对盖层封闭性的影响，使得最后得到的目标脆性盖层段的盖层封闭性更加准确。
本申请的实施例还提供一种计算机设备，图7为本发明实施例中计算机设备的示意图，该计算机设备能够实现上述实施例中的断块圈闭的脆性盖层封闭性分析方法中全部步骤，所述电子设备具体包括如下内容：
处理器(processor)701、存储器(memory)702、通信接口(CommunicationsInterface)703和总线704；
其中，所述处理器701、存储器702、通信接口703通过所述总线704完成相互间的通信；所述通信接口703用于实现服务器端设备、检测设备以及用户端设备等相关设备之间的信息传输；
所述处理器701用于调用所述存储器702中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的断块圈闭的脆性盖层封闭性分析方法中的全部步骤。
本申请的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的断块圈闭的脆性盖层封闭性分析方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的断块圈闭的脆性盖层封闭性分析方法的全部步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种断块圈闭的脆性盖层封闭性分析方法，其特征在于，包括：
在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析；
其中，根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，包括：
对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子；否则确定采用第二计算方式计算该层段的盖层封闭因子；
其中，所述第二计算方式为：若该层段的岩性属性为砂岩层段且该层段的厚度不大于所述砂岩过滤窗尺寸，或该层段的岩性属性为泥岩层段，采用如下公式计算每个层段的盖层封闭因子：
Figure FDA0003790450730000011
其中，FJTi为断面上每个层段中纵向上第i个采样点的盖层封闭因子；
ΔHi为断面纵向上每个层段中第i个采样点的采样长度；
N为断面在纵向上划分的采样点的数量；
Di为断面纵向上每个层段中第i个采样点的断距。
2.如权利要求1所述的断块圈闭的脆性盖层封闭性分析方法，其特征在于，所述岩性属性包括砂岩层段和泥岩层段；
确定断层模型中每一层段的岩性属性，包括：
在断层模型中，通过标记砂岩层段在断面上的位置，确定岩性属性为砂岩层段的层段。
3.如权利要求1所述的断块圈闭的脆性盖层封闭性分析方法，其特征在于，所述第一计算方式为：在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，采用如下公式计算每一层段的盖层封闭因子：
Figure FDA0003790450730000021
Figure FDA0003790450730000022
其中，FJTi为断面上第i个采样点的盖层封闭因子；
N为断面在纵向上划分的采样点的数量；
SWFJTj为断面纵向上第j个连续砂岩层段的厚度；
Sk为断层模型的断面上第j个连续砂岩层段内第k个采样点的厚度；
m为第j个连续砂岩层段中纵向包含的采样点个数；
ΔHi为断面上每个层段中纵向第i个采样点的采样长度；
Di为断面纵向上每个层段中第i个采样点的断距。
4.如权利要求1所述的断块圈闭的脆性盖层封闭性分析方法，其特征在于，根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析，包括：
对于每一层段的每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭。
5.一种断块圈闭的脆性盖层封闭性分析装置，其特征在于，包括：
断层模型构建模块，用于在砂泥岩互层地层的断块圈闭的脆性盖层段中，建立目标脆性盖层段的断层模型，确定断层模型中每一层段的岩性属性；
盖层封闭因子确定模块，用于根据预设的砂岩过滤窗尺寸、每一层段的岩性属性，确定每一层段的盖层封闭因子的计算方式，基于所述计算方式，计算每一层段的盖层封闭因子；
封闭性分析模块，用于根据每一层段的盖层封闭因子和目标脆性盖层段封闭因子阈值，在断面上实现目标脆性盖层段的盖层封闭性分析；
其中，所述盖层封闭因子确定模块具体用于：
对断层模型中的每一层段，在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，确定采用第一计算方式计算该层段的盖层封闭因子；否则确定采用第二计算方式计算该层段的盖层封闭因子；
其中，所述第二计算方式为：若该层段的岩性属性为砂岩层段且该层段的厚度不大于所述砂岩过滤窗尺寸，或该层段的岩性属性为泥岩层段，采用如下公式计算每个层段的盖层封闭因子；
Figure FDA0003790450730000031
其中，FJTi为断面上每个层段中纵向上第i个采样点的盖层封闭因子；
ΔHi为断面纵向上每个层段中第i个采样点的采样长度；
N为断面在纵向上划分的采样点的数量；
Di为断面纵向上每个层段中第i个采样点的断距。
6.如权利要求5所述的断块圈闭的脆性盖层封闭性分析装置，其特征在于，所述岩性属性包括砂岩层段和泥岩层段；
断层模型构建模块具体用于：
在断层模型中，通过标记砂岩层段在断面上的位置，确定岩性属性为砂岩层段的层段。
7.如权利要求5所述的断块圈闭的脆性盖层封闭性分析装置，其特征在于，所述第一计算方式为：在该层段的岩性属性为砂岩层段且该层段的厚度大于所述砂岩过滤窗尺寸时，采用如下公式计算每一层段的盖层封闭因子：
Figure FDA0003790450730000032
Figure FDA0003790450730000033
其中，FJTi为断面上第i个采样点的盖层封闭因子；
N为断面在纵向上划分的采样点的数量；
SWFJTj为断面纵向上第j个连续砂岩层段的厚度；
Sk为断层模型的断面上第j个连续砂岩层段内第k个采样点的厚度；
m为第j个连续砂岩层段中纵向包含的采样点个数；
ΔHi为断面上每个层段中纵向第i个采样点的采样长度；
Di为断面纵向上每个层段中第i个采样点的断距。
8.如权利要求5所述的断块圈闭的脆性盖层封闭性分析装置，其特征在于，封闭性分析模块具体用于：
对于每一层段的每一采样点，若该采样点的盖层封闭因子大于目标脆性盖层段封闭因子阈值，该采样点封闭，否则该采样点不封闭。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一项所述方法的计算机程序。
CN202010304303.8A 2020-04-17 2020-04-17 断块圈闭的脆性盖层封闭性分析方法及装置 Active CN111680386B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010304303.8A CN111680386B (zh) 2020-04-17 2020-04-17 断块圈闭的脆性盖层封闭性分析方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010304303.8A CN111680386B (zh) 2020-04-17 2020-04-17 断块圈闭的脆性盖层封闭性分析方法及装置

Publications (2)

Publication Number Publication Date
CN111680386A CN111680386A (zh) 2020-09-18
CN111680386B true CN111680386B (zh) 2022-11-04

Family

ID=72451488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010304303.8A Active CN111680386B (zh) 2020-04-17 2020-04-17 断块圈闭的脆性盖层封闭性分析方法及装置

Country Status (1)

Country Link
CN (1) CN111680386B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112965108B (zh) * 2021-02-05 2023-08-22 中国石油天然气股份有限公司 圈闭盖层垂向封闭性确定方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110632651A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断层封闭性定量分析方法及装置
CN110632655A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 非均质型断层的封闭性分析方法及装置
CN110632654A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断块圈闭含油边界确定方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110632651A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断层封闭性定量分析方法及装置
CN110632655A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 非均质型断层的封闭性分析方法及装置
CN110632654A (zh) * 2019-08-16 2019-12-31 中国石油天然气股份有限公司 断块圈闭含油边界确定方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
断块圈闭控藏模式――以辽河盆地西部凹陷为例;冷济高等;《油气地质与采收率》;20110325(第02期);全文 *
预测断层封堵性的FSPS软件系统;王晓红等;《石油勘探与开发》;20011030(第05期);全文 *

Also Published As

Publication number Publication date
CN111680386A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
AU2018340369B2 (en) Method and device for determining thin interlayer
CN106285660B (zh) 一种多层砂岩油藏低阻油层识别方法及装置
CN106324676B (zh) 一种确定断层封闭性的方法及装置
CN108122066B (zh) 储层岩性的确定方法和装置
CN111722277B (zh) 用于断溶体油气藏勘探的断溶体圈闭描述方法及装置
CN110632652B (zh) 断块油气藏的油气运移方式确定方法及装置
CN110632651B (zh) 断层封闭性定量分析方法及装置
CN106896406A (zh) 基于阻抗域属性预测优质储层的方法和装置
CN110632655B (zh) 非均质型断层的封闭性分析方法及装置
CN110632657B (zh) 泥岩涂抹型断层封闭性分析方法及装置
CN110632656B (zh) 断层的横向封闭性分析方法及装置
CN111680386B (zh) 断块圈闭的脆性盖层封闭性分析方法及装置
CN106338767A (zh) 多参数地震属性分析方法及系统
CN114063162B (zh) 一种基于小样本机器学习SVM Vp/Vs预测方法
CN106948812B (zh) 一种高渗带的渗透率下限值的确定方法和装置
CN115877464A (zh) 一种岩性识别方法、装置、计算机设备及存储介质
CN113534248A (zh) 圈闭断层的封闭性定量分析方法、装置及系统
Moradi et al. Learning from Behavioral Frac Maps: A Montney Case Study in Integration of Modern Microseismic and Production Data Analyses
CN105844082A (zh) 一种沉积相图中主导相的判断方法及系统
CN113031056B (zh) 构造约束下断层封闭性分析方法及装置
CN113466931B (zh) 一种地质时期断层封闭性评价方法及系统
Mendez et al. Rock Type Prediction Using Democratic Neural Networks Association DNNA: An Alternative Method to Describe Clastic Fillings in Ordovician Karsted Carbonates within a Region of the Tahe Oilfield, China
CN113495293B (zh) 油藏流体预测方法及装置
CN113552161B (zh) 页岩储层含油量确定方法、装置、设备及存储介质
Khodabakhshi et al. Field applications of a multi-scale multi-physics history matching approach

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