基于trap-3d软件寻找油气藏的方法
【专利摘要】本发明提供了一种基于TRAP-3D软件寻找油气藏的方法，该方法包括：根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体；将三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分；将三维岩性及断层数据立方体的深度切片顺序输入TRAP-3D软件进行油气藏评价。本发明提高了三维圈闭评价的精确性，有利于精确寻找油气藏，并且可绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并可得到勘探工区内油气藏的总圈闭量。
【专利说明】基于了以?-30软件寻找油气藏的方法
【技术领域】
[0001]本发明涉及油气勘探【技术领域】，尤其是涉及一种基于软件寻找油气藏的方法。
【背景技术】
[0002]适合于油气聚集、形成油气藏的场所称为圈闭(或油气圈闭目前圈闭概念是寻找石油和天然气的重要理论。没有圈闭，油气无法存在。但是长期以来，对于构造圈闭仍仅限于将构造隆起幅度和构造面积两项作为其圈闭评价的主要内容。但仅从上述两项指标难以定量衡量实际圈闭量的大小。而对于断块圈闭，通常仅根据断块的面积大小对圈闭面积进行估计，而对于该断块是否具备圈闭条件这一关键问题则缺乏精细研究。
[0003]目前断层封堵性、构造溢出点以及圈闭量估算研究的商业软件仍停留在定性分析或通过二维资料分析地下储层的阶段。其常规做法是在二维剖面上做断层两侧砂体构造叠合图，据此确定每个砂层的断层面溢出点，找出各封堵边界点；将上述封堵边界点投影至构造等高线平面并连线，以此计算封堵面积和圈闭量。而运用上述常规方法对由许多断层所形成的复杂断块的封堵性进行评价时，仅利用一系列断层两侧岩性二维对接剖面很难准确分析三维空间中砂层的连通及其总体的封堵情况，从而难以得到精确结果。
【发明内容】
[0004]本发明的目的在于提供一种基于I狀?-30软件寻找油气藏的方法，以提高三维圈闭评价的精确性。
[0005]为达到上述目的，本发明提供了一种基于I狀?-30软件寻找油气藏的方法，包括以下步骤:
[0006]根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体；
[0007]将所述三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分；
[0008]将所述三维岩性及断层数据立方体的深度切片顺序输入I狀?-30软件进行油气藏评价，其具体包括:
[0009]逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值；
[0010]根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值。
[0011]本发明的评价方法，所述根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体，具体包括:
[0012]将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。
[0013]本发明的评价方法，所述分别对每个深度切片进行独立砂体单元划分，具体为:
[0014](1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表不；
[0015](2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中所述每个砂点组的边界由泥岩和断层界定；
[0016](3)、将所述每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识；
[0017](4)、重复上述步骤(1)?(3)，直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
[0018]本发明的评价方法，所述预设的连通砂点搜索规则包括:
[0019]对于与断层不相邻的砂点，按照八方向搜索与该砂点相连通的相邻砂点；以及，
[0020]对于与断层相邻的砂点，按照四方向搜索与该砂点相连通的相邻砂点。
[0021]本发明的评价方法，在所述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括:
[0022]为每个深度切片内的各个独立砂体单元设置圈闭评价初值。
[0023]本发明的评价方法，在所述为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，包括:
[0024]如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则所述当前深度切片的某个独立砂体单元直接继承所述上一层深度切片的某个独立砂体单元的圈闭评价初值。
[0025]本发明的评价方法，所述逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括:
[0026](1)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点；
[0027](2)、根据公式3? =获取深度切片的圈闭量估值，其中，3?为第』个深度切片的圈闭量估值，为第』个深度切片的第1个独立砂体单元的空隙容积，
为第』个深度切片的第1个独立砂体单元的最终圈闭评价值；
[0028](3)、重复上述步骤(1)?(2)直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
[0029]本发明的评价方法，所述圈闭评价更新规则包括:
[0030]若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整所述各个独立砂体单元的圈闭评价值；
[0031]若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值:
[0032]= 00馭氺(^2?) /82^1/4
[0033]其中，I私1(为圈闭评价漏失值，00,为断层两侧独立砂体单元的圈闭评价值之差，1X2！！为断层两侧独立砂体单元的圈闭评价值岩性指数之和，8^1为漏失系数。
[0034]本发明首先把三维地震资料所获得的构造情况与测井资料获得的地下岩性的情况结合起来，加上断层的位置信息，综合成一个数据立方体；其次将数据立方体进行等深度划分，划分出若干等深度的深度切片，然后依次对每个深度切片进行圈闭评价，根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值，从而可以将工区储集油气的各种圈闭都找出来。由于本发明将数据立方体划分成若干等深度的深度切片，并依次对每个深度切片进行圈闭评价，从而通过这种更为精细的圈闭评价提高了三维圈闭评价的精确性，有利于精确寻找油气藏，并且本发明可绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并可得到勘探工区内油气藏的总圈闭量。
【专利附图】
【附图说明】
[0035]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0036]图1为本发明一个实施例的基于1狀？-30软件寻找油气藏的方法的流程图；
[0037]图2为本发明一个实施例中三维岩性及断层数据立方体的示意图；
[0038]图3为本发明一个实施例中一个深度切片的砂泥岩岩性划分及断层分布示意图；
[0039]图4为本发明一个实施例中一个深度切片的砂岩联通单元划分图；
[0040]图5为本发明一个实施例的平面甜点图；
[0041]图6为本发明一个实施例中自上而下的各深度切片内各砂体单元圈闭评价方法与评价传递方法示意图；
[0042]图7为本发明一个实施例中断层漏失量计算方法示意图；
[0043]图8为本发明一个实施例中经I狀?-30软件处理前后两个相邻切片的圈闭评价变化示意图；
[0044]图9为本发明一个实施例中经I狀?-30软件处理前后另外两个相邻垂直剖面的圈闭评价变化示意图；
[0045]图10为本发明一个实施例中深度域完成圈闭评价后输出的5条曲线示意图。
[0046]图1匕为本发明一个实施例中某工区均匀打采油井12 口的打井前储量分布示意图；
[0047]图116为本发明一个实施例中某工区均匀打采油井12 口的打井后，油气采完后的剩余储量分布示意图；
[0048]图123为本发明一个实施例中例油田的三维岩性数据立方体的一个切片的自然伽马值图；
[0049]图126和图133均为根据图12&中自然那伽马值推算的的砂泥岩岩性转换示意图；
[0050]图136为在图133中插入断层后的切片示意图；
[0051]图14为本发明一个实施例中例油田的一个经过I狀?-30软件分析后的垂直剖面的圈闭分析剖面图；
[0052]图153为本发明一个实施例中例油田按常规断块油田常规方法冒估的含油范围示意图；
[0053]图156为本发明一个实施例中例油田经I狀?-30软件处理后的含油范围示意图；
[0054]图16为本发明一个实施例中IX油田的深度域评价结果输出的5条统计曲线示意图。
【具体实施方式】
[0055]为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0056]下面结合附图，对本发明的【具体实施方式】作进一步的详细说明。
[0057]参考图1所示，本发明实施例的基于1狀？-30的油气藏评价方法包括以下步骤:
[0058]步骤3101、根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体。即将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。其中，三维岩性及断层数据立方体由按孔隙度大小区分的泥岩、粉砂岩、中等砂岩和好砂岩四种岩性以及断层面在三维等间距网格剖分空间中的展布信息组成，断层应为在三维空间等间距剖分网格中逐点连续展布的致密(无空缺点)断层面，例如相应的特征参数可以为:泥岩用0表示，粉砂岩用1表示，中等砂岩用2表示，好砂岩用3表示，断层用4表示，要求每个单元具有一个唯一的识别标识(10编号)，如图4所示。具体来说，根据项目需要，在已建立的构造模型的基础上，结合测井曲线((?)及声波曲线等对工区砕削岩进行岩性分类。解释的结果，以井间插值的方式进行了构造模型的岩性充填，并将沉积相解释结果作为边界条件对该模型进行了相控约束。用测井资料加上地震资料进行控制内插，将工区内按自然伽马值对砂泥岩岩性进行分类。根据60值的大小，把它转换为好砂岩3，中等砂岩2，粉砂岩1及泥岩0四种岩性。最终，建好的三维岩性及断层数据立方体，如图2所示。其中，I'狀?-30即三维圈闭。
[0059]步骤3102、将三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分。
[0060]如图3所示，为三维岩性及断层数据立方体的一个深度切片，其横坐标为水平方向X，纵坐标为水平方向V。图中白色为好砂岩，特征参数为3 ；浅灰色为中等砂岩，特征参数为2 ；深灰色为粉砂岩，特征参数为1 ；浅黑色为泥岩，特征参数为0 ；黑色为断层，特征参数为4。其中，分别对每个深度切片进行独立砂体单元划分，具体为:
[0061](1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表不。
[0062](2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中每个砂点组的边界由泥岩和断层界定；该预设的连通砂点搜索规则为:
[0063]对于与断层不相邻的砂点，按照八方向(例如0度、45度、90度、135度、180度、225度、270度和315度)搜索与该砂点相连通的相邻砂点；以及，
[0064]对于与断层相邻的砂点，按照四方向(例如0度、90度、180度和270度)搜索与该砂点相连通的相邻砂点。
[0065](3)、将每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识。
[0066](4)、重复上述步骤(1)?(3)，直至逐层顺序完成三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
[0067]举例来说，假设划分出的若干等深度的深度切片如图6所示。
[0068](^)切片[为特殊情况(因为它是第一张切片，上面没有评价控制):图中若有两个砂点，由于不知上方是否为泥岩，所以就给以最低评价IV = 16 ；
[0069](^)在切片吧里，泥岩里出现一个小砂点，它被上方附的泥岩所覆盖，因此，就给以最高评价靈=800 ；此外，令四周边水侵入。
[0070](0)在切片吧中，小砂点变成大砂饼，它们与上方是相连的，于是将上方吧中砂点的评价IV = 800传递给下它。
[0071]((1)在切片财中，大砂饼继续增大，但与边界相连了，于是边水侵入，圈闭破坏，评价降为最小，即靈=16。如果它与右方含水砂岩相连，IV也会降为最小。
[0072]步骤31033、逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值。其中，逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括:
[0073](1)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点；其中，圈闭评价更新规则为:
[0074]若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整各个独立砂体单元的圈闭评价值；
[0075]若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值:
[0076]= 00馭氺(^2?) /82^1/4
[0077]其中[拉汉为圈闭评价漏失值，0017为断层两侧独立砂体单元的圈闭评价值之差，1X2！！为断层两侧独立砂体单元的圈闭评价值岩性指数之和，3^1为漏失系数，8^1由用户指定，若32从=5就意味遇到对方5颗砂，才漏到和对方一样。
[0078]举例来说，例如图7所示，101好砂层在八80三处遇断层，0处断层对面是泥岩，8处断层对面是粉砂岩，八处对面是好砂岩，于是可以按照上述圈闭评价更新规则进行更新。
[0079](2)、根据公式3? =获取深度切片的圈闭量估值，其中，3?为第』个深度切片的圈闭量估值，为第』个深度切片的第1个独立砂体单元的空隙容积，
为第』个深度切片的第1个独立砂体单元的最终圈闭评价值。
[0080](3)、重复上述步骤(1)?(2)直至逐层顺序完成三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
[0081]步骤31031根据每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取勘探工区内储藏油气的总圈闭量估值。各个深度切片的圈闭量估值在深度2方向进行累积后绘出平面的总圈闭量，如图5所示，它反映着工区平面里每个点位处地下累计的圈闭量的大小，从该图可以看到，平面上哪里圈闭量最多，其含油气的可能性就大。就可以供布探井做抉择。
[0082]上述步骤3103^510313即为三维岩性及断层数据立方体的深度切片顺序输入丁狀?-30软件进行油气藏评价过程。此外，本发明实施例中，在上述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括:
[0083]为每个深度切片内的各个独立砂体单元设置圈闭评价初值。并且，在为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，还包括:
[0084]如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则当前深度切片的某个独立砂体单元直接继承上一层深度切片的某个独立砂体单元的圈闭评价初值。
[0085]本发明实施例中，如图8所示，假设两个相邻的切片41与42，图8左边分别是相邻的切片41与42的岩性指数的砂泥岩分布的显示，图8右边是经过I狀?-30软件进行圈闭分析后，切片41与42的圈闭评价值图。由右边两图可见由于边水渗入，也由于断层两侧砂岩见面，而使圈闭评价值逐步下降。轨迹如箭头所示，白色箭头表示边水通过断层联通而入侵，使圈闭评价降低的路径。
[0086]切片42的封闭含油范围比切片41更有所减少。
[0087]本发明实施例中，如图9所示，假设两个相邻的切片45与46为复杂岩性分布的复杂断块两条相邻的切片。图9左边分别是相邻的切片45与46的岩性指数仏的砂泥岩分布的显示，图9右边是经过I狀?-30软件进行圈闭分析后，切片45与45的圈闭评价值图，由此可以看到各砂层中的油水分布的情形(浅亮色是水其中，白色为圈闭评价极低的含水区，箭头所示处为油水界面。
[0088]此外，本发明实施例中，还可以在深度域绘出各深度切片的砂泥岩百分比及每个深度切片的圈闭量的对比图，如图10所示。这可更直观的看到在哪个深度上，哪个地层中圈闭体积最大。
[0089]本发明实施例中，如图1匕和图116所示，为打井12 口前后的储量分布及剩余未动用储量分布的示意图，在上述复杂断块模型的工区中均匀打井12 口，每口把所遇油层采干，结果还留下不能采到的“死油区”以及12 口井未钻遇的“剩余油”，从图116可以看到还有不少未动用储量可供开发。这些未动用储量也可以由I狀?-30软件计算出来。
[0090]下面以例油田为例，对本发明实施例进行举例说明:
[0091]例油田是一个复杂断块油田，本次研究面积约54平方千米，含油层系深度约1-2千米。工区内已打井8 口。见少量薄油层。本次切片分析的深度为1300111-1850111。
[0092]为了提供I狀?-30软件所需的输入数据格式，需要准备一个三维岩性及断层数据立方体。我们在8 口井测井数据的井控条件下作砂层内插，并加上地震解释的断层数据，形成三维岩性及断层数据立方体。假设三维岩性及断层数据立方体中每个单元点的大小为18111*18111*1111，深度方向按每1111为一片切片进行分析。
[0093](1)根据项目需要，在已建立的构造模型的基础上，结合测井曲线(⑶)及声波曲线等对工区砕削岩进行岩性分类。解释的结果，以井间插值的方式进行了构造模型的岩性充填，并将沉积相解释结果作为边界条件对该模型进行了相控约束。用测井资料加上地震资料进行控制内插，将工区内按自然伽马值对砂泥岩岩性进行分类。如图123所示，为例油田的某切片自然伽马⑶值彩色图。根据⑶值的大小，把它转换为好砂岩3，砂岩2，粉砂岩1及泥岩0，四种岩性，如图126和133图所示，为自然伽马值的转换图。
[0094]下一步再根据地震资料解释的断层数据融合进这些切片里去。
[0095]用通常的地震资料解释方法，对三维地震数据作过细的断层解释，从解释过的断层数据中导出断层文件，再将每个断层点见坐标插进解释好的三维岩性及断层数据立方体。给于每个断层点标上岩性指数为4。如图136所示，是砂泥岩划分的四种分布，右边是插入断层(黑色细线)后的情况。此图符合I'狀?-30程序的输入要求。
[0096]当每个切片都做完上述工作后，就完成了整个三维岩性及断层数据立方体的数据准备，就完成了适合于I'狀?-30软件分析的输入数据要求。
[0097](11)1^油田实际数据I'狀?-30试算结果
[0098]如图14所示，对于一条经过I'狀?-30软件处理的深度剖面，可看到砂层虽然很发育，但是只有少数地方具有圈闭。图14中右下方几条长条形黑线是断层，其它深黑色的是岩性圈闭，图14左方有一个背斜构造圈闭，如箭头所示。其中，灰色部分是泥岩，图中除下方两根垂向黑线是断层之外，其它黑色部分均为含油气的构造圈闭及岩性圈闭，白色为含水砂岩。
[0099]图153是例油田按断块油田常规方法冒估的含油范围，参考图156所示，例油田经I'狀?-30软件处理后的含油范围与上述冒估的含油范围基本吻合。由此可见，本申请的I'狀?-30软件更具有定量的概念，它包含了不少岩性圈闭，及北部的小背斜构造圈闭。
[0100]如图16所示，软件计算的按深度统计的各等深度切片中的砂泥岩百分比(曲线1及2)，根据311111(111(16X1计算的每个切片的空隙容积(曲线3)，表达储集性能的3皿|1”，(曲线4),以及最右边的曲线5表达着此深度切片中的圈闭量33 =811111 图16中，虽然深度16300处砂岩百分比高达60 %，但是圈闭量却不大，根据曲线4及5可看出深度1420-15200这一段砂岩虽然不那么多，但是圈闭量却较大，实际开采也证明这正是本油田的主要产层层段。可见I狀?-30软件给出了较好的定量的结果。
[0101](111)最后例油田的圈闭储量计算结果如下:
[0102]1)根据资料:平均孔隙度为28.5% ；平均含油饱和度为70% ；砂岩百分比为47.3% ；泥岩百分比为52.7% ；中部主体圈闭面积约10.5如12。
[0103]2)研究范围总圈闭储量为333，(即各切片全部加起来)。
[0104]888 = 2】X平均孔隙度^立方单兀体积/800
[0105]I'狀?-30软件计算结果:研究范围总圈闭储量约为0.53父108！:！3。
[0106]本发明实施例首先把三维地震资料所获得的构造情况与测井资料获得的地下岩性的情况结合起来，加上断层的位置信息，综合成一个数据立方体；其次将数据立方体进行等深度划分，划分出若干等深度的深度切片，然后依次对每个深度切片进行圈闭评价，从而可以将工区储集油气的各种圈闭都找出来。由于本发明实施例将数据立方体划分成若干等深度的深度切片，并依次对每个深度切片进行圈闭评价，从而通过这种更为精细的圈闭评价提高了三维圈闭评价的精确性。
[0107]以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种基于TRAP-3D软件寻找油气藏的方法，其特征在于，包括以下步骤: 根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体； 将所述三维岩性及断层数据立方体划分成若干等深度的深度切片，并分别对每个深度切片进行独立砂体单元划分； 将所述三维岩性及断层数据立方体的深度切片顺序输入TRAP-3D软件进行油气藏评价，其具体包括: 逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值； 根据所述每个深度切片内各个独立砂体单元的圈闭量估值绘出平面上的甜点图，以及展示纵向不同深度上的油气圈闭量，并获取所述勘探工区内储藏油气的总圈闭量估值。
2.根据权利要求1所述的方法，其特征在于，所述根据三维地震资料和测井资料建立勘探工区的三维岩性及断层数据立方体，具体包括: 将从三维地震资料中获得的构造信息和断层信息，与从测井资料中获得的岩性信息相结合建立勘探工区的三维岩性及断层数据立方体。
3.根据权利要求1所述的方法，其特征在于，所述分别对每个深度切片进行独立砂体单元划分，具体为: (1)、选择一个深度切片并将其内的岩性信息和断层信息分别用相应的特征参数表示； (2)、按照预设的连通砂点搜索规则搜索出该深度切片内所有的由若干相互连通的砂点构成的砂点组，其中所述每个砂点组的边界由泥岩和断层界定； (3)、将所述每个砂点组内所有砂点的特征参数合并，以作为一个独立砂体单元，并对每个独立砂体单元分别进行唯一性标识； (4)、重复上述步骤(I)?(3)，直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的独立砂体单元划分。
4.根据权利要求3所述的方法，其特征在于，所述预设的连通砂点搜索规则包括: 对于与断层不相邻的砂点，按照八方向搜索与该砂点相连通的相邻砂点；以及， 对于与断层相邻的砂点，按照四方向搜索与该砂点相连通的相邻砂点。
5.根据权利要求1所述的方法，其特征在于，在所述逐层对每个深度切片内的各个独立砂体单元分别进行圈闭评价之前还包括: 为每个深度切片内的各个独立砂体单元设置圈闭评价初值。
6.根据权利要求5所述的方法，其特征在于，在所述为每个深度切片内的各个独立砂体单元设置圈闭评价初值的过程中，包括: 如果当前深度切片的某个独立砂体单元内的砂点与其上一层深度切片的某个独立砂体单元内的砂点存在上下连通关系，则所述当前深度切片的某个独立砂体单元直接继承所述上一层深度切片的某个独立砂体单元的圈闭评价初值。
7.根据权利要求5所述的方法，其特征在于，所述逐层对每个深度切片内各个独立砂体单元分别进行圈闭评价，获取每个深度切片内的各个独立砂体单元的圈闭量估值，具体包括: (I)、选择一个深度切片，根据圈闭评价更新规则和该深度切片里的断层两侧的岩性的对置情况，分别对该深度切片内的各个独立砂体单元的圈闭评价初值进行更新，以获得该深度切片内的各个独立砂体单元的最终圈闭评价值，并将对应的最终圈闭评价值赋予对应独立砂体单元内的每个砂点； (2)、根据公式SS^= Sumiindexi^MVJ获取深度切片的圈闭量估值，其中，SS^为第j个深度切片的圈闭量估值，Indexi为第j个深度切片的第i个独立砂体单元的空隙容积，MVi为第j个深度切片的第i个独立砂体单元的最终圈闭评价值； (3)、重复上述步骤(I)?(2)直至逐层顺序完成所述三维岩性及断层数据立方体的所有深度切片的圈闭量估值的获取。
8.根据权利要求7所述的方法，其特征在于，所述圈闭评价更新规则包括: 若断层对侧为泥岩，则根据每个深度切片内的各个独立砂体单元的泥岩的涂抹封堵效应的不同对应调整所述各个独立砂体单元的圈闭评价值； 若断层对侧为砂岩，则当每个深度切片内的断层一侧的独立砂体单元的圈闭评价值高于该断层对侧的独立砂体单元的圈闭评价值时，按下式计算圈闭评价漏失值: LEAK = DDMV* (YXZH)/SEAL/4，其中，LEAK为圈闭评价漏失值，DDMV为断层两侧独立砂体单元的圈闭评价值之差，YXZH为断层两侧独立砂体单元的圈闭评价值岩性指数之和，SEAL为漏失系数。
【文档编号】G01V1/28GK104375179SQ201410635467
【公开日】2015年2月25日 申请日期:2014年11月5日 优先权日:2014年11月5日
【发明者】李庆忠, 张立彬, 张晓敏 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司