CN104375179B - 基于trap‑3d软件寻找油气藏的方法 - Google Patents
基于trap‑3d软件寻找油气藏的方法 Download PDFInfo
- Publication number
- CN104375179B CN104375179B CN201410635467.3A CN201410635467A CN104375179B CN 104375179 B CN104375179 B CN 104375179B CN 201410635467 A CN201410635467 A CN 201410635467A CN 104375179 B CN104375179 B CN 104375179B
- Authority
- CN
- China
- Prior art keywords
- trap
- depth
- section
- evaluation
- body unit
- 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 27
- 239000004576 sand Substances 0.000 claims abstract description 158
- 238000011156 evaluation Methods 0.000 claims abstract description 87
- 235000021185 dessert Nutrition 0.000 claims abstract description 7
- 238000003325 tomography Methods 0.000 claims description 50
- 239000004575 stone Substances 0.000 claims description 19
- 239000004215 Carbon black (E152) Substances 0.000 claims description 8
- 229930195733 hydrocarbon Natural products 0.000 claims description 8
- 150000002430 hydrocarbons Chemical class 0.000 claims description 8
- 239000011800 void material Substances 0.000 claims description 4
- 239000007789 gas Substances 0.000 description 18
- 238000010586 diagram Methods 0.000 description 14
- 239000011435 rock Substances 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000009826 distribution Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000006467 substitution reaction Methods 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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/644—Connectivity, e.g. for fluid movement
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/645—Fluid contacts
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
本发明提供了一种基于TRAP‑3D软件寻找油气藏的方法，该方法包括：根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体；将三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分；将三维岩性及断层数据立方体的深度切片顺序输入TRAP‑3D软件进行油气藏评价。本发明提高了三维圈闭评价的精确性，有利于精确寻找油气藏，并且可绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并可得到勘探工区内油气藏的总圈闭量。
Description
技术领域
本发明涉及油气勘探技术领域，尤其是涉及一种基于TRAP-3D软件寻找油气藏的方法。
背景技术
适合于油气聚集、形成油气藏的场所称为圈闭(或油气圈闭)。目前圈闭概念是寻找石油和天然气的重要理论。没有圈闭，油气无法存在。但是长期以来，对于构造圈闭仍仅限于将构造***幅度和构造面积两项作为其圈闭评价的主要内容。但仅从上述两项指标难以定量衡量实际圈闭量的大小。而对于断块圈闭，通常仅根据断块的面积大小对圈闭面积进行估计，而对于该断块是否具备圈闭条件这一关键问题则缺乏精细研究。
目前断层封堵性、构造溢出点以及圈闭量估算研究的商业软件仍停留在定性分析或通过二维资料分析地下储层的阶段。其常规做法是在二维剖面上做断层两侧砂体构造叠合图，据此确定每个砂层的断层面溢出点，找出各封堵边界点；将上述封堵边界点投影至构造等高线平面并连线，以此计算封堵面积和圈闭量。而运用上述常规方法对由许多断层所形成的复杂断块的封堵性进行评价时，仅利用一系列断层两侧岩性二维对接剖面很难准确分析三维空间中砂层的连通及其总体的封堵情况，从而难以得到精确结果。
发明内容
本发明的目的在于提供一种基于TRAP-3D软件寻找油气藏的方法，以提高三维圈闭评价的精确性。
为达到上述目的，本发明提供了一种基于TRAP-3D软件寻找油气藏的方法，包括以下步骤：
根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体；
将所述三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分；
将所述三维岩性及断层数据立方体的深度切片顺序输入TRAP-3D软件进行油气藏评价，其具体包括：
逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值；
根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值。
本发明的评价方法，所述根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体，具体包括：
将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。
本发明的评价方法，所述分别对每个深度切片进行独立砂体单元划分，具体为：
(1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表示；
(2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中所述每个砂点组的边界由泥岩和断层界定；
(3)、将所述每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识；
(4)、重复上述步骤(1)～(3)，直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
本发明的评价方法，所述预设的连通砂点搜索规则包括：
对于与断层不相邻的砂点，按照八方向搜索与该砂点相连通的相邻砂点；以及，
对于与断层相邻的砂点，按照四方向搜索与该砂点相连通的相邻砂点。
本发明的评价方法，在所述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括：
为每个深度切片内的各个独立砂体单元设置圈闭评价初值。
本发明的评价方法，在所述为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，包括：
如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则所述当前深度切片的某个独立砂体单元直接继承所述上一层深度切片的某个独立砂体单元的圈闭评价初值。
本发明的评价方法，所述逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括：
(1)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点；
(2)、根据公式SSj＝Sum{indexi*MVi}获取深度切片的圈闭量估值，其中，SSj为第j个深度切片的圈闭量估值，indexi为第j个深度切片的第i个独立砂体单元的空隙容积，MVi为第j个深度切片的第i个独立砂体单元的最终圈闭评价值；
(3)、重复上述步骤(1)～(2)直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
本发明的评价方法，所述圈闭评价更新规则包括：
若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整所述各个独立砂体单元的圈闭评价值；
若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值：
LEAK＝DDMV*(YXZH)/SEAL/4
其中，LEAK为圈闭评价漏失值，DDMV为断层两侧独立砂体单元的圈闭评价值之差，YXZH为断层两侧独立砂体单元的圈闭评价值岩性指数之和，SEAL为漏失系数。
本发明首先把三维地震资料所获得的构造情况与测井资料获得的地下岩性的情况结合起来，加上断层的位置信息，综合成一个数据立方体；其次将数据立方体进行等深度划分，划分出若干等深度的深度切片，然后依次对每个深度切片进行圈闭评价，根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值，从而可以将工区储集油气的各种圈闭都找出来。由于本发明将数据立方体划分成若干等深度的深度切片，并依次对每个深度切片进行圈闭评价，从而通过这种更为精细的圈闭评价提高了三维圈闭评价的精确性，有利于精确寻找油气藏，并且本发明可绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并可得到勘探工区内油气藏的总圈闭量。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明一个实施例的基于TRAP-3D软件寻找油气藏的方法的流程图；
图2为本发明一个实施例中三维岩性及断层数据立方体的示意图；
图3为本发明一个实施例中一个深度切片的砂泥岩岩性划分及断层分布示意图；
图4为本发明一个实施例中一个深度切片的砂岩连通单元划分图；
图5为本发明一个实施例的平面甜点图；
图6为本发明一个实施例中自上而下的各深度切片内各砂体单元圈闭评价方法与评价传递方法示意图；
图7为本发明一个实施例中断层漏失量计算方法示意图；
图8为本发明一个实施例中经TRAP-3D软件处理前后两个相邻切片的圈闭评价变化示意图；
图9为本发明一个实施例中经TRAP-3D软件处理前后另外两个相邻垂直剖面的圈闭评价变化示意图；
图10为本发明一个实施例中深度域完成圈闭评价后输出的5条曲线示意图。
图11a为本发明一个实施例中某工区均匀打采油井12口的打井前储量分布示意图；
图11b为本发明一个实施例中某工区均匀打采油井12口的打井后，油气采完后的剩余储量分布示意图；
图12a为本发明一个实施例中TN油田的三维岩性数据立方体的一个切片的自然伽马值图；
图12b和图13a均为根据图12a中自然那伽马值推算的砂泥岩岩性转换示意图；
图13b为在图13a中***断层后的切片示意图；
图14为本发明一个实施例中TN油田的一个经过TRAP-3D软件分析后的垂直剖面的圈闭分析剖面图；
图15a为本发明一个实施例中TN油田按常规断块油田常规方法冒估的含油范围示意图；
图15b为本发明一个实施例中TN油田经TRAP-3D软件处理后的含油范围示意图；
图16为本发明一个实施例中TN油田的深度域评价结果输出的5条统计曲线示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
下面结合附图，对本发明的具体实施方式作进一步的详细说明。
参考图1所示，本发明实施例的基于TRAP-3D的油气藏评价方法包括以下步骤：
步骤S101、根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体。即将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。其中，三维岩性及断层数据立方体由按孔隙度大小区分的泥岩、粉砂岩、中等砂岩和好砂岩四种岩性以及断层面在三维等间距网格剖分空间中的展布信息组成，断层应为在三维空间等间距剖分网格中逐点连续展布的致密(无空缺点)断层面，例如相应的特征参数可以为：泥岩用0表示，粉砂岩用1表示，中等砂岩用2表示，好砂岩用3表示，断层用4表示，要求每个单元具有一个唯一的识别标识(ID编号)，如图4所示。具体来说，根据项目需要，在已建立的构造模型的基础上，结合测井曲线(GR)及声波曲线等对工区碎屑岩进行岩性分类。解释的结果，以井间插值的方式进行了构造模型的岩性充填，并将沉积相解释结果作为边界条件对该模型进行了相控约束。用测井资料加上地震资料进行控制内插，将工区内按自然伽马值对砂泥岩岩性进行分类。根据GR值的大小，把它转换为好砂岩3，中等砂岩2，粉砂岩1及泥岩0四种岩性。最终，建好的三维岩性及断层数据立方体，如图2所示。其中，TRAP-3D即三维圈闭。
步骤S102、将三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分。
如图3所示，为三维岩性及断层数据立方体的一个深度切片，其横坐标为水平方向X，纵坐标为水平方向Y。图中白色为好砂岩，特征参数为3；浅灰色为中等砂岩，特征参数为2；深灰色为粉砂岩，特征参数为1；浅黑色为泥岩，特征参数为0；黑色为断层，特征参数为4。其中，分别对每个深度切片进行独立砂体单元划分，具体为：
(1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表示。
(2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中每个砂点组的边界由泥岩和断层界定；该预设的连通砂点搜索规则为：
对于与断层不相邻的砂点，按照八方向(例如0度、45度、90度、135度、180度、225度、270度和315度)搜索与该砂点相连通的相邻砂点；以及，
对于与断层相邻的砂点，按照四方向(例如0度、90度、180度和270度)搜索与该砂点相连通的相邻砂点。
(3)、将每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识。
(4)、重复上述步骤(1)～(3)，直至逐层顺序完成三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
举例来说，假设划分出的若干等深度的深度切片如图6所示。
(a)切片N1为特殊情况(因为它是第一张切片，上面没有评价控制):图中若有两个砂点，由于不知上方是否为泥岩,所以就给以最低评价MV＝16；
(b)在切片N2里，泥岩里出现一个小砂点，它被上方N1的泥岩所覆盖，因此，就给以最高评价MV＝800；此外，令四周边水侵入。
(c)在切片N3中，小砂点变成大砂饼，它们与上方是相连的，于是将上方N2中砂点的评价MV＝800传递给下它。
(d)在切片N4中,大砂饼继续增大，但与边界相连了，于是边水侵入，圈闭破坏，评价降为最小，即MV＝16。如果它与右方含水砂岩相连，MV也会降为最小。
步骤S103a、逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值。其中，逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括：
(1)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点；其中，圈闭评价更新规则为：
若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整各个独立砂体单元的圈闭评价值；
若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值：
LEAK＝DDMV*(YXZH)/SEAL/4
其中LEAK为圈闭评价漏失值，DDMV为断层两侧独立砂体单元的圈闭评价值之差，YXZH为断层两侧独立砂体单元的圈闭评价值岩性指数之和，SEAL为漏失系数，SEAL由用户指定，若SEAL＝5就意味遇到对方5颗砂，才漏到和对方一样。
举例来说，例如图7所示，101好砂层在ABC三处遇断层，C处断层对面是泥岩，B处断层对面是粉砂岩，A处对面是好砂岩，于是可以按照上述圈闭评价更新规则进行更新。
(2)、根据公式SSj＝Sum{indexi*MVi}获取深度切片的圈闭量估值，其中，SSj为第j个深度切片的圈闭量估值，indexi为第j个深度切片的第i个独立砂体单元的空隙容积，MVi为第j个深度切片的第i个独立砂体单元的最终圈闭评价值。
(3)、重复上述步骤(1)～(2)直至逐层顺序完成三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
步骤S103b、根据每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取勘探工区内储藏油气的总圈闭量估值。各个深度切片的圈闭量估值在深度Z方向进行累积后绘出平面的总圈闭量，如图5所示，它反映着工区平面里每个点位处地下累计的圈闭量的大小，从该图可以看到，平面上哪里圈闭量最多，其含油气的可能性就大。就可以供布探井做抉择。
上述步骤S103a-S103b即为三维岩性及断层数据立方体的深度切片顺序输入TRAP-3D软件进行油气藏评价过程。此外，本发明实施例中，在上述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括：
为每个深度切片内的各个独立砂体单元设置圈闭评价初值。并且，在为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，还包括：
如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则当前深度切片的某个独立砂体单元直接继承上一层深度切片的某个独立砂体单元的圈闭评价初值。
本发明实施例中，如图8所示，假设两个相邻的切片41与42，图8左边分别是相邻的切片41与42的岩性指数(Litho-Index)的砂泥岩分布的显示，图8右边是经过TRAP-3D软件进行圈闭分析后，切片41与42的圈闭评价值图。由右边两图可见由于边水渗入，也由于断层两侧砂岩见面，而使圈闭评价值逐步下降。轨迹如箭头所示，白色箭头表示边水通过断层连通而入侵，使圈闭评价降低的路径。
切片42的封闭含油范围比切片41更有所减少。
本发明实施例中，如图9所示，假设两个相邻的切片45与46为复杂岩性分布的复杂断块两条相邻的切片。图9左边分别是相邻的切片45与46的岩性指数(Litho-Index)的砂泥岩分布的显示，图9右边是经过TRAP-3D软件进行圈闭分析后，切片45与45的圈闭评价值图，由此可以看到各砂层中的油水分布的情形(浅亮色是水)。其中，白色为圈闭评价极低的含水区，箭头所示处为油水界面。
此外，本发明实施例中，还可以在深度域绘出各深度切片的砂泥岩百分比及每个深度切片的圈闭量的对比图，如图10所示。这可更直观的看到在哪个深度上，哪个地层中圈闭体积最大。
本发明实施例中，如图11a和图11b所示，为打井12口前后的储量分布及剩余未动用储量分布的示意图，在上述复杂断块模型的工区中均匀打井12口，每口把所遇油层采干，结果还留下不能采到的“死油区”以及12口井未钻遇的“剩余油”，从图11b可以看到还有不少未动用储量可供开发。这些未动用储量也可以由TRAP-3D软件计算出来。
下面以TN油田为例，对本发明实施例进行举例说明：
TN油田是一个复杂断块油田，本次研究面积约54平方千米，含油层系深度约1-2千米。工区内已打井8口。见少量薄油层。本次切片分析的深度为1300m-1850m。
为了提供TRAP-3D软件所需的输入数据格式，需要准备一个三维岩性及断层数据立方体。我们在8口井测井数据的井控条件下作砂层内插，并加上地震解释的断层数据，形成三维岩性及断层数据立方体。假设三维岩性及断层数据立方体中每个单元点的大小为18m*18m*1m，深度方向按每1m为一片切片进行分析。
(i)根据项目需要，在已建立的构造模型的基础上，结合测井曲线(GR)及声波曲线等对工区碎屑岩进行岩性分类。解释的结果，以井间插值的方式进行了构造模型的岩性充填，并将沉积相解释结果作为边界条件对该模型进行了相控约束。用测井资料加上地震资料进行控制内插，将工区内按自然伽马值对砂泥岩岩性进行分类。如图12a所示，为TN油田的某切片自然伽马GR值彩色图。根据GR值的大小，把它转换为好砂岩3，砂岩2，粉砂岩1及泥岩0，四种岩性，如图12b和13a图所示，为自然伽马值的转换图。
下一步再根据地震资料解释的断层数据融合进这些切片里去。
用通常的地震资料解释方法，对三维地震数据作过细的断层解释，从解释过的断层数据中导出断层文件，再将每个断层点XYZ坐标插进解释好的三维岩性及断层数据立方体。给于每个断层点标上岩性指数为4。如图13b所示，是砂泥岩划分的四种分布，右边是***断层(黑色细线)后的情况。此图符合TRAP-3D程序的输入要求。
当每个切片都做完上述工作后，就完成了整个三维岩性及断层数据立方体的数据准备，就完成了适合于TRAP-3D软件分析的输入数据要求。
(ii)TN油田实际数据TRAP-3D试算结果
如图14所示，对于一条经过TRAP-3D软件处理的深度剖面，可看到砂层虽然很发育，但是只有少数地方具有圈闭。图14中右下方几条长条形黑线是断层，其它深黑色的是岩性圈闭，图14左方有一个背斜构造圈闭，如箭头所示。其中，灰色部分是泥岩，图中除下方两根垂向黑线是断层之外，其它黑色部分均为含油气的构造圈闭及岩性圈闭，白色为含水砂岩。
图15a是TN油田按断块油田常规方法冒估的含油范围，参考图15b所示，TN油田经TRAP-3D软件处理后的含油范围与上述冒估的含油范围基本吻合。由此可见，本申请的TRAP-3D软件更具有定量的概念，它包含了不少岩性圈闭，及北部的小背斜构造圈闭。
如图16所示，TRAP-3D软件计算的按深度统计的各等深度切片中的砂泥岩百分比(曲线1及2)，根据Sum{index}计算的每个切片的空隙容积(曲线3)，表达储集性能的Sum{MV}，(曲线4)，以及最右边的曲线5表达着此深度切片中的圈闭量SS＝Sum{index*MV}图16中，虽然深度1630m处砂岩百分比高达60％，但是圈闭量却不大，根据曲线4及5可看出深度1420-1520m这一段砂岩虽然不那么多，但是圈闭量却较大，实际开采也证明这正是本油田的主要产层层段。可见TRAP-3D软件给出了较好的定量的结果。
(iii)最后TN油田的圈闭储量计算结果如下：
1)根据资料:平均孔隙度为28.5％；平均含油饱和度为70％；砂岩百分比为47.3％；泥岩百分比为52.7％；中部主体圈闭面积约10.5km2。
2)研究范围总圈闭储量为SSS，(即各切片全部加起来)。
SSS＝Σ【Sum{index*MV}】×平均孔隙度*立方单元体积/800
TRAP-3D软件计算结果:研究范围总圈闭储量约为0.53×108m3。
本发明实施例首先把三维地震资料所获得的构造情况与测井资料获得的地下岩性的情况结合起来，加上断层的位置信息，综合成一个数据立方体；其次将数据立方体进行等深度划分，划分出若干等深度的深度切片，然后依次对每个深度切片进行圈闭评价，从而可以将工区储集油气的各种圈闭都找出来。由于本发明实施例将数据立方体划分成若干等深度的深度切片，并依次对每个深度切片进行圈闭评价，从而通过这种更为精细的圈闭评价提高了三维圈闭评价的精确性。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于TRAP-3D软件寻找油气藏的方法，其特征在于，包括以下步骤：
根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体；
将所述三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分；
将所述三维岩性及断层数据立方体的深度切片顺序输入TRAP-3D软件进行油气藏评价，其具体包括：
逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值；
根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值。
2.根据权利要求1所述的方法，其特征在于，所述根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体，具体包括：
将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。
3.根据权利要求1所述的方法，其特征在于，所述分别对每个深度切片进行独立砂体单元划分，具体为：
(1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表示；
(2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中所述每个砂点组的边界由泥岩和断层界定；
(3)、将所述每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识；
(4)、重复上述步骤(1)～(3)，直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
4.根据权利要求3所述的方法，其特征在于，所述预设的连通砂点搜索规则包括：
对于与断层不相邻的砂点，按照八方向搜索与该砂点相连通的相邻砂点；以及，
对于与断层相邻的砂点，按照四方向搜索与该砂点相连通的相邻砂点。
5.根据权利要求1所述的方法，其特征在于，在所述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括：
为每个深度切片内的各个独立砂体单元设置圈闭评价初值。
6.根据权利要求5所述的方法，其特征在于，在所述为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，包括：
如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则所述当前深度切片的某个独立砂体单元直接继承所述上一层深度切片的某个独立砂体单元的圈闭评价初值。
7.根据权利要求5所述的方法，其特征在于，所述逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括：
(1)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点；
(2)、根据公式SSj＝Sum{indexi*MVi}获取深度切片的圈闭量估值，其中，SSj为第j个深度切片的圈闭量估值，indexi为第j个深度切片的第i个独立砂体单元的空隙容积，MVi为第j个深度切片的第i个独立砂体单元的最终圈闭评价值；
(3)、重复上述步骤(1)～(2)直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
8.根据权利要求7所述的方法，其特征在于，所述圈闭评价更新规则包括：
若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整所述各个独立砂体单元的圈闭评价值；
若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值：
LEAK＝DDMV*(YXZH)/SEAL/4，其中，LEAK为圈闭评价漏失值，DDMV为断层两侧独立砂体单元的圈闭评价值之差，YXZH为断层两侧独立砂体单元的圈闭评价值岩性指数之和，SEAL为漏失系数。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410635467.3A CN104375179B (zh) | 2014-11-05 | 2014-11-05 | 基于trap‑3d软件寻找油气藏的方法 |
| CA2912626A CA2912626C (en) | 2014-11-05 | 2015-04-23 | 3d trap evaluation method of searching for oil-gas reservoir |
| PCT/CN2015/077283 WO2016070587A1 (zh) | 2014-11-05 | 2015-04-23 | 基于trap-3d软件寻找油气藏的方法 |
| EP15794444.8A EP3217192B1 (en) | 2014-11-05 | 2015-04-23 | Method for looking for oil and gas reservoir on the basis of trap-3d software |
| US14/946,466 US10228478B2 (en) | 2014-11-05 | 2015-11-19 | Method of searching for oil-gas reservoir based on trap-3D software |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410635467.3A CN104375179B (zh) | 2014-11-05 | 2014-11-05 | 基于trap‑3d软件寻找油气藏的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104375179A CN104375179A (zh) | 2015-02-25 |
| CN104375179B true CN104375179B (zh) | 2017-02-15 |
Family
ID=52554207
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410635467.3A Active CN104375179B (zh) | 2014-11-05 | 2014-11-05 | 基于trap‑3d软件寻找油气藏的方法 |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP3217192B1 (zh) |
| CN (1) | CN104375179B (zh) |
| WO (1) | WO2016070587A1 (zh) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2912626C (en) | 2014-11-05 | 2018-02-13 | China National Petroleum Corporation | 3d trap evaluation method of searching for oil-gas reservoir |
| CN104375179B (zh) * | 2014-11-05 | 2017-02-15 | 中国石油天然气集团公司 | 基于trap‑3d软件寻找油气藏的方法 |
| CN106443779B (zh) * | 2016-08-17 | 2018-08-07 | 中国石油化工股份有限公司 | 一种中小型砂砾岩体的储层预测方法 |
| CN109343115B (zh) | 2018-11-21 | 2019-12-03 | 成都理工大学 | 一种基于测井约束的含气储层刻画方法 |
| CN109738947B (zh) * | 2018-12-12 | 2020-11-20 | 核工业北京地质研究院 | 一种圈定砂岩型铀矿找矿远景区的物化探组合方法 |
| CN109597129B (zh) * | 2019-01-07 | 2019-12-17 | 中国地质大学（武汉） | 基于目标检测的缝洞型油藏串珠状反射特征识别方法 |
| CN109917455B (zh) * | 2019-02-19 | 2021-03-30 | 中国石油天然气集团有限公司 | 一种地层圈闭的刻画方法及*** |
| CN113589371B (zh) * | 2020-04-30 | 2023-03-21 | 中国石油化工股份有限公司 | 碳酸盐岩内断层的封闭性评价方法 |
| CN111505720B (zh) * | 2020-05-21 | 2023-03-28 | 中国海洋石油集团有限公司 | 一种岩性圈闭的刻画方法 |
| EP3926368B1 (en) | 2020-06-15 | 2024-01-10 | TotalEnergies OneTech | Analyzing a hydrocarbon trap |
| CN112034513B (zh) * | 2020-09-03 | 2023-03-14 | 中海石油(中国)有限公司 | 一种已开发油田周边潜力资源量的定量评价方法 |
| CN112255676B (zh) * | 2020-10-15 | 2023-12-19 | 中国海洋石油集团有限公司 | 旋向相反的共轭走滑断层共轭区有利圈闭定量寻找方法 |
| CN112360444B (zh) * | 2020-10-16 | 2023-08-22 | 中国石油天然气股份有限公司 | 圈闭目标的量化排序处理方法及装置 |
| CN114594518B (zh) * | 2020-12-04 | 2023-04-07 | 中国石油化工股份有限公司 | 基于井震交替的开发后期复杂断块精细地层对比方法 |
| CN112883554B (zh) * | 2021-01-22 | 2024-03-29 | 中国科学技术大学 | 一种构造圈闭识别方法及*** |
| CN112983401B (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种边水气藏水侵边界确定方法 |
| CN113514875B (zh) * | 2021-06-03 | 2024-04-02 | 德仕能源科技集团股份有限公司 | 一种基于大数据的岩性油气藏的勘测和描述方法及*** |
| CN113514883B (zh) * | 2021-06-18 | 2023-03-17 | 中国石油化工股份有限公司 | 一种断层-岩性油藏刻画方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4008459A (en) * | 1975-07-18 | 1977-02-15 | Texaco Inc. | Method of seismic surveying for stratigraphic traps |
| US5563513A (en) * | 1993-12-09 | 1996-10-08 | Stratasearch Corp. | Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps |
| CN101520517A (zh) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | 一种能准确评价碎屑岩盆地含油气目标的方法 |
| CN103604916A (zh) * | 2013-11-21 | 2014-02-26 | 中国石油大学(北京) | 一种连续型致密砂岩气藏含气范围的测量方法及*** |
| CN103792592A (zh) * | 2013-03-14 | 2014-05-14 | 中国石油大学(北京) | 一种确定陆相断陷盆地油气藏分布最有利范围的方法 |
| CN104662445A (zh) * | 2012-10-26 | 2015-05-27 | 雪佛龙美国公司 | 圈闭完整性分析的***和方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5563949A (en) * | 1994-12-12 | 1996-10-08 | Amoco Corporation | Method of seismic signal processing and exploration |
| FR2909200B1 (fr) * | 2006-11-23 | 2009-07-03 | Total Sa | Procede, programme et systeme informatique de simulation de chenaux |
| GB0900687D0 (en) * | 2009-01-16 | 2009-02-25 | Rock Deformation Res Ltd | Fault analysis system |
| US9194968B2 (en) * | 2010-05-28 | 2015-11-24 | Exxonmobil Upstream Research Company | Method for seismic hydrocarbon system analysis |
| CN104375179B (zh) * | 2014-11-05 | 2017-02-15 | 中国石油天然气集团公司 | 基于trap‑3d软件寻找油气藏的方法 |
-
2014
- 2014-11-05 CN CN201410635467.3A patent/CN104375179B/zh active Active
-
2015
- 2015-04-23 WO PCT/CN2015/077283 patent/WO2016070587A1/zh active Application Filing
- 2015-04-23 EP EP15794444.8A patent/EP3217192B1/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4008459A (en) * | 1975-07-18 | 1977-02-15 | Texaco Inc. | Method of seismic surveying for stratigraphic traps |
| US5563513A (en) * | 1993-12-09 | 1996-10-08 | Stratasearch Corp. | Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps |
| CN101520517A (zh) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | 一种能准确评价碎屑岩盆地含油气目标的方法 |
| CN104662445A (zh) * | 2012-10-26 | 2015-05-27 | 雪佛龙美国公司 | 圈闭完整性分析的***和方法 |
| CN103792592A (zh) * | 2013-03-14 | 2014-05-14 | 中国石油大学(北京) | 一种确定陆相断陷盆地油气藏分布最有利范围的方法 |
| CN103604916A (zh) * | 2013-11-21 | 2014-02-26 | 中国石油大学(北京) | 一种连续型致密砂岩气藏含气范围的测量方法及*** |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104375179A (zh) | 2015-02-25 |
| EP3217192A1 (en) | 2017-09-13 |
| WO2016070587A1 (zh) | 2016-05-12 |
| EP3217192A4 (en) | 2018-06-27 |
| EP3217192B1 (en) | 2021-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104375179B (zh) | 基于trap‑3d软件寻找油气藏的方法 | |
| CN106934858B (zh) | 一种矿集区尺度区域三维地质建模方法及*** | |
| Ringrose et al. | Reservoir model design | |
| Cannon | Reservoir modelling: A practical guide | |
| US10228478B2 (en) | Method of searching for oil-gas reservoir based on trap-3D software | |
| CN105469443B (zh) | 基于地质路线（prb）过程双重建模生成三维地质图的方法 | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| CN107133670A (zh) | 一种基于决策树数据挖掘算法的复杂岩性识别方法及*** | |
| Mahmud et al. | Petrophysical interpretations of subsurface stratigraphic correlations, Baram Delta, Sarawak, Malaysia | |
| CN108663714A (zh) | 一种沉积微相刻画方法 | |
| Lyu et al. | Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China | |
| Majeed et al. | Constructing 3D Geological Model for Tertiary Reservoir in Khabaz Oil Field by using Petrel software. | |
| Naidu et al. | Calibrated basin modelling to understand hydrocarbon distribution in Barmer Basin, India | |
| RU2681250C1 (ru) | Способ квазитрехмерного моделирования эффективных газо- и нефтенасыщенных толщин залежей углеводородов | |
| Masoud et al. | Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya | |
| Freeman et al. | A New Facies Classification Scheme Using Gamma Ray and Bulk Density Logs, With Multiple Practical Applications in North Kuwait Heavy Oil Fields | |
| Yang | Realistic and theoretical 3D modelling of the sedimentation, burial, thermal and tectonic history of the gippsland rift basin | |
| Ba et al. | Applying multi-point statistical methods to build the facies model for Oligocene formation, X oil field, Cuu Long basin | |
| Chapin et al. | Integrated static and dynamic modeling of the Pinedale tight gas field, Wyoming | |
| Abdul Majeed et al. | Constructing 3D Geological Model for Tertiary Reservoir in Khabaz Oil Field by using Petrel software. | |
| Airen et al. | Reservoir Hydrocarbon Volumetric Analysis of Sapele Deep Field, Niger Delta, Southern Nigeria | |
| Taura et al. | Structural and petrophysical controls on the remaining fluid distribution in the reservoirs of Gharif Formation before abandonment | |
| Tellez | Integrated characterization of tight siliciclastic reservoirs: examples from the Cretaceous Burro Canyon Formation, Colorado, and Mississippian Meramec Strata, Oklahoma | |
| CN108510586B (zh) | 一种构建地质导向地层模型的方法 | |
| Kurianova et al. | Comprehensive Analysis of the Geological and Geophysical Data in the Study of the Upper Miocene Turbidite Systems of the Nam Con Son Basin, Vietnam |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |