岩性油气藏  2024, Vol. 36 Issue (2): 160-169       PDF    
×
基于SSOM的流动单元划分方法及生产应用——以渤海湾盆地F油田古近系沙三中亚段湖底浊积水道为例
王亚, 刘宗宾, 路研, 王永平, 刘超     
中海石油(中国)有限公司天津分公司, 天津 300452
摘要: 在储层构型级次划分的基础上,归纳总结了渤海湾盆地F油田古近系沙三中亚段湖底浊积水道各级渗流屏障与构型界面的关系,并采用监督模式下的自组织映射神经网络算法开展流动单元定量评价,明确了构型模式控制下的流动单元分布规律。研究结果表明:①渤海湾盆地F油田沙三中亚段储层可划分为Ⅰ、Ⅱ、Ⅲ、Ⅳ等4类流动单元。②基于SSOM算法的流动单元识别结果表现出较高的预测准确率,其中,256组训练样本的整体回判准确率为82.81%,110组测试样本的识别准确率为80.91%,能够满足地质油藏研究的需求。③垂向上,不同类型的单一水道内部发育的流动单元类型不同,造成流动单元垂向分布的差异性。Ⅰ类、Ⅱ类优质流动单元发育在浊积水道体系演化的中期,主要分布在二类单一水道;Ⅲ类、Ⅳ类流动单元发育在浊积水道体系演化的早期和晚期,其中,Ⅲ类流动单元分布广泛,在一类、二类、三类单一水道均有分布,Ⅳ类流动单元发育在一类、三类单一水道。侧向上,受水道体系不同沉积演化阶段的影响,流动单元的侧向分割方式不同。流动单元与非渗透层接触样式发育在浊积水道体系早期沉积旋回中,接触界面明显,属于一类水道沉积砂体;流动单元间的相互接触样式则发育在水道体系中晚期沉积旋回中,属于二类、三类水道沉积砂体。④平面上,受单一水道侧向迁移和垂向加积作用的影响,复合水道内部流动单元平面展布差异性明显。Ⅲ类流动单元在单砂体内部广泛发育,井间连续性好,在浊积主水道、浊积水道及水道漫溢沉积砂体处均有发育;渗流能力较好的Ⅰ类、Ⅱ类流动单元分布局限、连续性差,仅在浊积水道主流线方向及主水道砂体处有分布,呈不连续的点状或带状分布;Ⅳ类流动单元呈环带状分布在Ⅲ类流动单元的外缘,在浊积水道漫溢沉积砂体处发育。
关键词: SSOM算法    浊积水道    储层构型    流动单元    自组织映射神经网络    监督模式    沙三中亚段    古近系    渤海湾盆地    
Flow unit division based on SSOM and its production application: A case study of sublacustrine turbidity channels of middle Es3in F oilfield, Bohai Bay Basin
WANG Ya, LIU Zongbin, LU Yan, WANG Yongping, LIU Chao     
Tianjin Branch of CNOOC China Limited, Tianjin 300452, China
Abstract: Based on the classification of reservoir architecture, the relationship between seepage barriers and architecture interfaces of the sublacustrine turbidity channels of the middle third member of Paleogene Shahejie Formation(Es3)in F oilfield of Bohai Bay Basin was summarized, a supervised self-organizing map neural network algorithm was used to quantitatively evaluate flow units, and the distribution of flow units under the control of architecture models was clarified. The results show that: (1)The reservoirs of middle Es3 in F oilfield of Bohai Bay Basin can be divided into four types of flow units, namely type Ⅰ, Ⅱ, Ⅲ and Ⅳ.(2)The flow unit recognition results based on SSOM algorithm show high prediction accuracy, with an overall accuracy of 82.81% for 256 training samples and 80.91% for 110 testing samples, which can meet the needs of geological reservoir research. (3)Vertically, the flow units developed in different types of single channels may vary greatly, resulting in difference of vertical distribution of flow units. The favorable typeⅠand Ⅱ flow units are developed in the middle stage of turbidity channel system evolution, mainly distributed in the category Ⅱ single channels. Type Ⅲ and Ⅳ flow units are developed in the early and late stages of turbidity channel system evolution. Among them, type Ⅲ flow units are widely distributed in category Ⅰ, Ⅱ and Ⅲ single channels, while type Ⅳ flow units are mainly developed in categoryⅠand Ⅲ single channels. Laterally, influenced by the sedimentary evolution stages of channels systems, the lateral division of flow units is different. The contact pattern between the flow units and the nonpermeable layers develops in the early sedimentary cycle of the turbidity channel system, with obvious contact interfaces, which belongs to category Ⅰ channel sand bodies. The contact pattern between the flow units develops in the middle and late sedimentary cycles of the channel system, which belongs to category Ⅱ and Ⅲ channel sand bodies. (4)Horizontally, due to the lateral migration and vertical accretion of single channels, the distribution of flow units in the composite channel was significantly different. The type Ⅲ flow units are developed widely in the single sand body with good continuity between wells, and they are developed in the main turbidity channels, turbidity channels and overflowing sand bodies. The distribution of typeⅠand Ⅱ flow units with good seepage capacity is limited with poor continuity, and they are distributed only in the direction of the main stream line of turbidity channels and the sand bodies of the main channels, showing discontinuous point or ribbon distribution. The Type Ⅳ flow units are distributed in a ring band along the outer edge of type Ⅲ flow units and developed in the area where the turbidity channels overflowed the sand bodies.
Key words: SSOM algorithm    turbidity channel    reservoir architecture    flow unit    self-organizing map neural network    supervised mode    middle Es3    Paleogene    Bohai Bay Basin    
0 引言

浊积扇低渗透储层地下构型结构复杂,内部砂体形态、规模、延伸状况存在明显差异,砂体内广泛发育的隔层和夹层使得厚油层被分割成了多个连通或不连通的渗流单元。目前,已有不少学者开展了浊积扇储层构型模式的研究[1-3],国外以Mutti等[4]的五级划分方案应用最为广泛,而国内比较有代表性的是林煜等[5]和李志鹏等[6]的六级划分方案。这些浊积扇储层构型模式的研究为油田开发初期的井位部署提供了重要的指导作用,随着油田开发进入中后期,层间、层内的干扰导致开发效果变差,构型模式的认识难以对该阶段储层油水运动规律的表征形成指导作用,亟需在构型模式的基础上明确不同渗流单元的空间展布,以期为完善开发井网和注采结构调整提供重要依据。

流动单元是表征储层渗流性能好坏的地质单元,自Hearn等[7]1984年首次提出流动单元的概念以来,被广泛应用于储层精细油藏描述[8-9]。流动单元的划分方法由最早的统计学方法,如线性判别分析、聚类分析、模糊数学逐渐向人工智能(支持向量机、人工神经网络等)领域迈进[10-12],其评价参数也由单一的静态参数向多参数融合分析转变[13]。以往研究表明,2种或多种方法的组合通常比单一方法具有更好的识别效果[12-14]。精细的储层构型模式研究是渗流屏障识别及连通体划分的基础[15-17],对流动单元研究至关重要。一方面,储层构型模式研究可以找到识别和预测的成因单元,并在此结构内部依据岩性、物性及渗流特征等开展流动单元识别;另一方面,构型模式研究与渗流屏障识别关系密切,能够为渗流屏障识别和级别划分提供理论参考依据。然而,以往对基于储层构型模式的流动单元研究相对较少。

以渤海湾盆地F油田沙三中亚段湖底浊积水道储层为例,综合运用岩心、测井及生产动态等资料,以取心层段定性划分的流动单元作为训练样本,采用监督模式下的自组织映射神经网络算法与因子分析法相结合的方式开展未取心层段的流动单元定量识别,并对构型模式控制下的流动单元分布规律进行分析,以期为后续剩余油挖潜、优化注采井网及提高开发效果等提供依据。

1 地质概况

F油田位于渤海湾盆地东南部,构造断裂活动较弱,构造形态相对简单,是一被断层切割的穹隆背斜构造,东部由北东向和南东向的2条断层交叉遮挡,对砂体沉积及油气聚集起控制作用(图 1a1b)。沙三中亚段4砂层组为其主要含油层系,垂向上根据隔夹层发育特征可划分出4个小层以满足油田开发。目的层油藏埋深为2 800~3 100 m,为典型的湖底浊积扇沉积,以中扇、外扇亚相为沉积主体(图 1c)。其中,浊积水道砂体发育,单砂层厚度为6.6~13.7 m,平均为10.7 m,垂向上水道砂体多期叠置,平面上呈大面积连片分布。受浊积水道垂向加积和侧向迁移变化快的影响,对于储层空间展布特征及连通性的研究难度较大。

下载原图 图 1 渤海湾盆地F油田古近系沙三中亚段构造特征(a,b)及岩性地层综合柱状图(c) Fig. 1 Structural features(a, b)and stratigraphic column(c)of middle Es3 in F oilfield, Bohai Bay Basin

研究区沙三中亚段4砂层组浊积扇沉积特征复杂、岩性多样,发育含砾砂岩、中—粗砂岩、中—细砂岩、粉砂岩、泥质粉砂岩及泥岩。储层岩石类型主要为岩屑长石砂岩,碎屑颗粒中石英体积分数占39.67%,长石体积分数占33.60%,岩屑体积分数占26.73%,成分成熟度较低。砂岩分选中等,分选系数为1.3~1.9,粒度中值为0.12~0.16 mm。粒间填隙物以碳酸盐胶结物为主,呈孔隙式胶结,体积分数平均为8.7%。受埋藏压实作用和碳酸盐胶结作用的影响,储层物性、渗流能力均较差,储层孔隙度平均为15.9%,渗透率平均为11.6 mD,且不同小层间非均质性较强(图 2a)。开发进入中后期,受储层非均质性的影响,层间、层内干扰显著,注水开发效果变差,各层的采出程度差异大,造成水窜严重(图 2b2c)。因此亟需开展流动单元的相关研究,以明确不同渗流单元的空间展布,为完善开发井网、细分开发层系以及改善开发效果提供依据。

下载原图 图 2 渤海湾盆地F油田古近系沙三中亚段各小层的渗透率(a)、相对吸水强度(b)和采出程度(c)对比 Fig. 2 Comparison of permeability(a), relative water absorption strength(b)and recovery degree(c)of sublayers of middle Es3 in F oilfield, Bohai Bay Basin
2 基于储层构型模式的流动单元分类 2.1 基于构型的渗流屏障及连通体分析

根据Mutti等[4]和林煜等[5]针对浊积水道的构型划分方案,结合浊积水道储层的分布特征及其对油田开发的影响,建立了渤海湾盆地F油田沙三中亚段湖底浊积水道储层8个级次的构型界面划分方案:构型级次由高到低依次为湖底扇复合体—单一湖底扇—水道体系—复合水道—单一水道—水道内部增生体或韵律层—纹层组—纹层。此外,依据浊积水道的成因类型和沉积模式可将研究区进一步划分为3类单一水道,分别为一类单一水道(以块状砂砾岩相及中—粗砂岩相为主)、二类单一水道(以中—细砂岩相为主)及三类单一水道(以泥质粉砂岩相为主),3种类型的单一水道在垂向上具有不同的沉积模式(图 3)。

下载原图 图 3 渤海湾盆地F油田古近系沙三中亚段单一水道类型及沉积模式 A. 递变层理段;B. 下平行纹层段;C. 流水沙纹层段;D. 上平行纹层段;E. 泥岩段。 Fig. 3 Single channel types and sedimentary model of middle Es3 in F oilfield, Bohai Bay Basin

一类单一水道形成于浊积水道体系演化的早期,是高密度浊流沉积的产物,自水道底部向顶部依次发育块状砂砾岩相、块状粗砂岩相、中—粗砂岩夹薄层粉—细砂岩相;鲍马序列的底部递变层理段最发育,顶层沉积不发育,序列模式为AB组合;自然伽马曲线(GR)和电阻率曲线(Rt)呈齿化箱形。二类单一水道形成于浊积水道体系演化的中期,属低密度浊流沉积,水道底部以块状砂砾岩相和块状中—粗砂岩相为主,主体部位主要发育块状中—细砂岩相及少量的薄层细砂岩,顶部主要发育泥质粉砂岩相,沉积物粒度自下而上逐渐变细,为典型的正韵律沉积;鲍马序列发育有ABCD组合和ABC组合;GRRt曲线呈钟形、箱形或箱形-钟形组合,内部存在回返。三类单一水道形成于浊积水道体系演化的晚期,属浊流消亡期的泥质浊流沉积,沉积物粒度较细,以泥质粉砂岩为主;鲍马序列发育最为完整,包括ABCDE组合和ABCD组合,以上平行纹层段的泥质粉砂岩和粉砂岩沉积为主,底部递变层理发育较少;GRRt曲线呈齿化严重的钟形,内部回返明显(图 3)。

从水道体系、复合水道与单一水道3个级次开展渗流屏障分析。研究区为简单的背斜构造,钙质胶结物在储层中难以形成稳定分布的层间渗流屏障,以局部的层内渗流屏障存在,因此主要分析泥质渗流屏障。水道体系的划分界面为不同水道体系间区域稳定分布的泥质隔层,对应构型级次划分的六级(⑥)界面,为一级渗流屏障;复合水道的泥岩屏障往往为发育广泛、分布相对稳定的垂向泥质隔层,对应构型级次划分的五级(⑤)界面,为二级渗流屏障;单一水道间的渗流屏障为水道侧缘沉积或水道间的漫溢沉积,对应构型级次划分的四级(④)界面,为三级渗流屏障;水道内部充填的泥质夹层和钙质夹层,对应构型级次划分的三级(③)界面,为四级渗流屏障。其中,一级、二级渗流屏障分布稳定,所限定的储层连通体之间没有流体流动;三级渗流屏障呈不连续分布,将连通体分割成了若干个连通或不连通的储集单元;四级渗流屏障所限定的储层连通体内部再无渗流屏障的存在,连通体内的渗流特征差异只与岩相充填特征有关(图 4)。

下载原图 图 4 渤海湾盆地F油田古近系沙三中亚段渗流屏障级别与储层构型的关系 Fig. 4 Relationship between seepage barriers and reservoir architecture of middle Es 3 in F oilfield, Bohai Bay Basin
2.2 流动单元划分

基于浊积水道构型和渗流屏障的分析,明确浊积扇沉积体内部水道体系、复合水道与单一水道3个级次的构型要素及对应级次的渗流屏障分布,对单一水道开展流动单元划分及识别。采用流动分层指标法(FZI)法,结合储集物性、微观孔隙结构和生产动态等特征划分取心层段的流动单元类型,根据岩心物性分析数据计算的IFZ值(流动分层指标)的累计概率曲线确定4个主要的分布区间(图 5a)。以FZI法为主,结合储层的岩石物理特征和生产动态特征,可将取心层段的岩心样本划分为Ⅰ,Ⅱ,Ⅲ,Ⅳ等4类储层流动单元(图 5表 1)。其中,Ⅰ类储层流动单元的储集物性和渗流能力最好,渗透率大于12.1 mD,平均为23.0 mD,初期产能最高,日产液量平均为34.2 t。Ⅱ类储层流动单元储集物性和渗流能力相对较好,渗透率为5.9~33.2 mD,平均为12.9 mD,初期产能较高,日产液量平均为22.1 t。Ⅲ类储层流动单元的储集物性和渗流能力一般,渗透率为3.6~28.6 mD,平均为11.3 mD,初期产能较差,日产液量平均为13.7 t。Ⅳ类储层流动单元的储集物性和渗流能力最差,渗透率平均为8.6 mD,初期产能最差,动用程度较低(表 1)。

下载原图 图 5 渤海湾盆地F油田古近系沙三中亚段基于流动分层指标法(FZI)划分流动单元 Fig. 5 Flow units division based on FZI method of middle Es3 in F oilfield, Bohai Bay Basin
下载CSV 表 1 渤海湾盆地F油田古近系沙三中亚段流动单元特征参数统计据 Table 1 Characteristic parameters of different types of flow units of middle Es3 in F oilfield, Bohai Bay Basin
2.3 单井流动单元

单井流动单元能够反映储层内部的渗流特征分布,因此,在构型模式的指导下进行了单井流动单元研究。首先,在垂向上区分储层与非储层;其次,针对储层的物性特征、孔隙结构特征及流动特征差异开展单井流动单元识别。如图 6所示,单井流动单元划分结果与初期产能分布具有较好的匹配关系。Ⅰ类流动单元与Ⅱ类流动单元初期产能较高,Ⅲ类流动单元初期产能一般。此外,由于同一单井的不同单一水道发育于不同的水道体系演化阶段,水道内部充填特征不同,导致流动单元类型及分布存在差异。对于不同单井的同一单层,由于水道发育位置不同,水道内部充填特征也存在差异,导致井间流动单元分布不同(图 6)。

下载原图 图 6 渤海湾盆地F油田古近系沙三中亚段单井流动单元划分结果 Fig. 6 Division results of flow units in single wells of middle Es3 in F oilfield, Bohai Bay Basin
3 基于SSOM的流动单元测井定量识别 3.1 自组织映射神经网络算法原理

自组织映射(SOM,Self-Organizing Maps)神经网络算法是Kohonen[18]于1982年提出的一种无监督模式下的聚类分析算法,包括输入层和竞争层2层神经网络结构(图 7a)。输入层的节点数即为流动单元输入样本的维数,当接收到流动单元样本的输入特征时,位于竞争层的神经元相互竞争对输入模式获取响应的机会。最终,获胜的神经元兴奋,失败的神经元受到抑制,获胜的神经元即为输入样本的流动单元分类。

下载原图 图 7 自组织映射神经网络结构示意图[19] Fig. 7 Structure diagram showing self-organizing map neural network

SOM算法在岩性识别、流体预测及遥感图像分类等方面取得了较好的应用效果[19-20],但在解决非线性关系强烈的模式分类问题时,存在分类准确度不高、样本统计困难、难以挖掘出非显著信息等缺点[21]。因此,国内外学者对其进行了改进,提出了一种监督模式下的自组织神经网络(SSOM)算法[21-23]。相比SOM算法,SSOM算法增加了输出层,其神经网络结构如图 7b所示。监督模式下的学习训练过程,通过调整获胜神经元节点邻域内所有节点与输入层和输出层节点之间的权值完成竞争层与输入层、输出层之间的权值校正。其算法的实质是通过分析流动单元输入样本的预测类别和岩心定性分析流动单元类别之间的吻合程度,进行权值调整,从而建立流动单元预测模型,对未取心井的流动单元进行预测。

3.2 基于SSOM神经网络的流动单元预测

(1)特征曲线的选取及数据预处理

输入特征的选取直接影响着训练结果的准确性,本研究选取了影响岩石渗流能力的4类参数,包括储集物性参数:孔隙度(φ)、渗透率(k)、泥质含量(Vsh)、含油饱和度(So)、储层品质指数(IRQ)、渗透率突进系数(Tk)、渗透率级差(Jk);沉积特征参数:沉积微相、砂岩厚度(H)、夹层厚度(hd)及夹层密度(ρd);微观孔隙结构参数:孔喉半径(R35)及流动分层指标(IFZ)等在内的13个储层评价参数。为避免维度过高导致运算困难,在进行神经网络运算之前,采用因子分析算法,对输入参数进行了降维处理。当选取前5个主因子时,累计方差贡献率大于80%(表 2)。因此,选取前5个主因子(F1F2F3F4F5)作为神经网络模型的输入特征;以FZI法获取的流动单元类型作为模型的输出类型,即可构建神经网络运算所必须的样本数据集。

下载CSV 表 2 方差分析结果 Table 2 Rsults of variance analysis

(2)网络结构初始化

本次采用的神经网络为10×10的对称网络,每个竞争层的神经网络节点代表一种流动单元储层参数响应特征,输出层的神经网络节点代表岩心样本的流动单元类别。对神经网络数据进行随机初始化,学习速率η的初始值选0.5,随着训练时间和迭代次数的增加逐渐减少至0.000 1(学习率阈值)后结束训练;邻域距离的初始值设为2,随着迭代的进行逐渐缩短邻域距离;设置输入层与竞争层之间的权值(Wij)以及竞争层与输出层之间的权值(Wjk)为0~1的相异随机数。

(3)计算获胜神经元节点g

本次研究改进了传统的最小欧式距离法,通过引入最大概率法以实现定量计算获胜神经元节点的目的。首先,计算输入样本xi与所有竞争层神经元节点j的欧式距离Dj

$ {D_j} = \left\| {{x_i} - {w_{ij}}} \right\|, j = {\rm{ }}1, {\rm{ }}2, {\rm{ }} \cdots , {\rm{ }}100 $ (1)

其次,根据欧式距离进行反距离加权,即可计算出神经元节点的获胜概率值Pj

$ {P(g) = {\rm{ max}}P(j) = {\rm{ max}}\frac{{({D_j}) - \mu }}{{\sum\limits_{j = 1}^{100} {({D_j}) - \mu } }}} $ (2)

式中:Pj)为单个神经元获胜的概率;Pg)为获胜神经元的获胜概率;μ为任意正实数(通常取2)。

(4)权系数调整

训练获胜神经元邻域距离内的所有神经元,调整他们的权系数,包括获胜神经元节点邻域内所有节点与输入层的权值Wij及其与输出层的权值Wjk。当学习率η小于阈值0.000 1时,训练结束。

$ {{W_{ij}} = {W_{ij}} + {\eta _1}({x_i} - {W_{ij}})} $ (3)

$ {{W_{jk}} = {W_{jk}} + {\eta _2}({Y_k} - {W_{jk}})} $ (4)

式中:Yk为取心井样本所属的流动单元类别;η1η2均为学习率,随着训练时间及迭代次数的增加而线性降低。

(5)预测模型输出与分析

图 8为反映流动单元分类结果的自组织映射拓扑图及3D Sammon投影图。拓扑图中保留了流动单元样本的输入特征及输出类别,每一个神经元节点都代表了一种储层参数响应特征,方格颜色代表该神经元节点所属的流动单元类别,且落在同一类流动单元内的数据点具有相同的岩石物理特征。3D Sammon投影图是拓扑图在三维空间的投影,直观地反映了SSOM模型分类情况和各个竞争层神经元节点之间的距离。

下载原图 图 8 自组织映射聚类网络拓扑图(a)及竞争神经元节点3D Sammon投影(b) Fig. 8 Self-organizing map topology(a)and 3D Sammon mapping projection of competing neuron nodes(b)
3.3 流动单元识别结果与分析

以取心井定性划分的366个流动单元样本数据作为模型训练和验证过程中的样本数据集,随机选取70% 的样本数据作为训练样本,其余30% 的样本作为测试样本对模型进行回判检验。基于SSOM算法的流动单元识别结果表明,256组训练样本仅出现44个错判,训练样本的整体回判准确率为82.81%,且110组测试样本的识别准确率为80.91%,表现出较高的预测准确率(表 3)。测试样本中有21个样本发生错判,其中19个错判样本的流动单元类型被识别为相邻类别,主要是由于浊积水道非均质性强,仅靠4类流动单元类型难以全面表征其非均质性特征,而更加精细化的流动单元分类方案又需要大量岩心分析数据作为支撑,在有限的取心资料情况下,总体错判率能够满足地质油藏研究的需求。因此,监督模式下的自组织映射神经网络算法通过神经元竞争学习、相互监督的方式,能够输出最优的流动单元预测模型,为流动单元的研究提供了一种新的手段。

下载CSV 表 3 渤海湾盆地F油田古近系沙三中亚段基于SSOM的流动单元预测精度统计 Table 3 Statistics of flow unit prediction accuracy based on SSOM of middle Es3 in F oilfield, Bohai Bay Basin
4 储层构型控制下的流动单元分布规律

渤海湾盆地F油田未取心井段流动单元的分类识别可以借助SSOM模型实现,在单井识别的基础上开展了井间剖面和平面流动单元分布预测[24-25]

4.1 剖面分布规律

受沉积物源、水动力条件、水道的迁移叠置等因素的影响,浊积水道体系不同演化时期发育的单一水道类型存在差异。自水道体系演化的早期—晚期,单一水道类型表现出自一类向三类过渡的趋势。通过单一水道与流动单元剖面分布的对比分析发现,单一水道类型的差异性造成了流动单元剖面分布的差异性,即流动单元剖面展布特征受水道垂向演化控制作用显著,主要表现在以下2个方面。

(1)垂向上,Ⅰ,Ⅱ,Ⅲ,Ⅳ等4类流动单元在F油田4砂层组的分布差异明显。总体表现为Ⅲ类流动单元分布广泛,Ⅱ类次之,Ⅰ类、Ⅳ类流动单元发育较少。从水道体系级次来看,Ⅰ类、Ⅱ类优质流动单元发育在浊积水道体系演化的中期,Ⅲ类、Ⅳ类流动单元发育在浊积水道体系演化的早期和晚期。从单一水道来看,不同类型单一水道成因类型和沉积模式的差异性,导致不同类型单一水道内部发育的流动单元类型存在明显差异。Ⅰ类、Ⅱ类流动单元分布在二类单一水道;Ⅲ类流动单元分布广泛,在一类、二类、三类单一水道均有分布;Ⅳ类流动单元发育在一类、三类单一水道(图 9)。

下载原图 图 9 渤海湾盆地F油田古近系沙三中亚段流动单元与单一水道剖面分布特征 Fig. 9 Distribution of flow units and single channels of middle Es3 in F oilfield, Bohai Bay Basin

(2)侧向上,流动单元间呈2种接触样式,分别为流动单元与非渗透层(即隔层)直接接触和流动单元间的相互接触。前者发育在浊积水道体系早期沉积旋回中,接触界面明显,属于一类水道沉积砂体(如44、45小层沉积砂体侧向上发育流动单元与非渗透层的直接接触)。后者接触界面不明显,通常显示出侧向上流动单元渗流特征的差异,发育在水道体系中晚期沉积旋回中,属于二类、三类水道沉积砂体(如42、43小层沉积砂体侧向上多发育流动单元间的相互接触)(图 9)。

4.2 平面分布规律

平面上以43-2单砂体为例,分析复合水道级次流动单元的平面分布规律。受单一水道侧向迁移和垂向加积作用的影响,复合水道内部泥质渗流屏障较为发育,水道内部渗流差异明显,造成了复合水道内部流动单元平面展布的差异性。研究结果表明,43-2单砂体内部主要发育Ⅲ类流动单元,井间连续性好,在浊积主水道、浊积水道及水道漫溢沉积砂体处均有发育;渗流能力较好的Ⅰ类、Ⅱ类流动单元分布局限、连续性差,仅在浊积水道主流线方向及主水道砂体处有分布,呈不连续的点状或带状分布;Ⅳ类流动单元呈环带状分布在Ⅲ类流动单元的外缘,在浊积水道漫溢沉积砂体处发育。总体来说,43-2单砂体平面非均质性较强,流动单元平面分布差异性明显(图 10)。

下载原图 图 10 渤海湾盆地F油田古近系沙三中亚段43-2单砂体流动单元分布与构型单元平面分布对比 Fig. 10 Planar distribution of flow units and architecture elements of 43-2 single sand layer of middle Es3 in F oilfield, Bohai Bay Basin
5 结论

(1)研究区湖底浊积水道砂体构型级次由高到低依次为湖底扇复合体—单一湖底扇—水道体系—复合水道—单一水道—水道内部增生体或韵律层—纹层组—纹层。浊积水道砂体可进一步划分为一类单一水道(以块状砂砾岩相及中—粗砂岩相为主)、二类单一水道(以中—细砂岩相为主)及三类单一水道(以泥质粉砂岩相为主)。

(2)在取心层段流动单元划分的基础上,采用SSOM算法开展未取心层段的流动单元定量识别,测试样本正判率达80.91%,具有较高的预测精度和显著的分类效果。

(3)流动单元垂向展布受水道垂向演化控制作用明显。从水道体系来看,优质流动单元主要发育在水道体系演化的中期,而渗流能力较差的流动单元多发育在水道体系演化的早期和晚期,且不同演化阶段的流动单元侧向分割方式不同;从单一水道来看,不同类型的单一水道内部流动单元类型不同。

(4)复合水道级次的流动单元的平面展布差异性明显,优质流动单元仅发育在浊积水道主流线方向及主水道砂体处,渗流能力较差的流动单元多发育在水道漫溢沉积砂体处。

参考文献
[1]
满晓, 胡德胜, 吴洁, 等. 北部湾盆地涠西南凹陷始新统流一段湖底扇发育特征及成藏模式. 岩性油气藏, 2023, 35(4): 137-144.
MAN Xiao, HU Desheng, WU Jie, et al. Development characteristics and accumulation model of sublacustrine fans of the first member of Eocene Liushagang Formation in Weixinan Sag, Beibuwan Basin. Lithologic Reservoirs, 2023, 35(4): 137-144.
[2]
WYNN R B, CRONIN B T, PEAKALL J. Sinuous deep-water channels: Genesis, geometry and architecture. Marine and Petroleum Geology, 2007, 24(6/7/8/9): 341-387.
[3]
张文彪, 段太忠, 刘志强, 等. 深水浊积水道沉积构型模式及沉积演化: 以西非M油田为例. 地球科学, 2017, 42(2): 273-285.
ZHANG Wenbiao, DUAN Taizhong, LIU Zhiqiang, et al. Architecture model and sedimentary evolution of deepwater turbidity channel: A case study of M oilfield in West Africa. Earth Science, 2017, 42(2): 273-285.
[4]
MUTTI E, NORMARK W R. Comparing examples of modern and ancient turbidite systems: Problems and concepts. Marine Clastic Sedimentology, 1987, 1: 1-38.
[5]
林煜, 吴胜和, 王星, 等. 深水浊积水道体系构型模式研究: 以西非尼日尔三角洲盆地某深水研究区为例. 地质论评, 2013, 59(3): 510-520.
LIN Yu, WU Shenghe, WANG Xing, et al. Research on architecture model of deep water turbidity channel system: A case study of a deepwater research area in Niger Delta Basin, West Africa. Geological Review, 2013, 59(3): 510-520.
[6]
李志鹏, 杨勇, 侯加根, 等. 渤南五区浊积水道砂体储层构型研究. 中国石油大学学报(自然科学版), 2015, 39(5): 36-42.
LI Zhipeng, YANG Yong, HOU Jiagen, et al. Reservoir architecture of turbidity channels in the 5 th area of Bonan Oilfield. Journal of China University of Petroleum(Edition of Natural Science), 2015, 39(5): 36-42. DOI:10.3969/j.issn.1673-5005.2015.05.005
[7]
HEARN C L, EBANKS W J, TYE R S. Geological factors influencing reservoir performance of the Hartzog Draw Field, Wyoming. Journal of Petroleum Technology, 1984, 36(8): 1335-1344. DOI:10.2118/12016-PA
[8]
袁丙龙, 张辉, 叶青, 等. 基于三角洲复合砂体构型的流动单元划分及剩余油分布模式. 沉积学报, 2021, 39(5): 1253-1263.
YUAN Binglong, ZHANG Hui, YE Qing, et al. Flow-unit classification based on compound sand-body architecture of delta and distribution pattern of remaining oil. Acta Sedimentologica Sinica, 2021, 39(5): 1253-1263.
[9]
NOORUDDIN H A, HOSSAIN M E. Modified Kozeny-Carmen correlation for enhanced hydraulic flow unit characterization. Journal of Petroleum Science and Engineering, 2012, 80(1): 107-115.
[10]
周游, 李治平, 景成, 等. 基于"岩石物理相-流动单元"测井响应定量评价特低渗透油藏优质储层: 以延长油田东部油区长6油层组为例. 岩性油气藏, 2017, 29(1): 116-123.
ZHOU You, LI Zhiping, JING Cheng, et al. Quantitative evaluation of favorable reservoir in ultra-low permeable reservoir based on "petrophysical facies-flow unit" log response: A case study of Chang 6 oil reservoir set in Yanchang Oilfield. Lithologic Reservoirs, 2017, 29(1): 116-123.
[11]
HATAMPOUR A, SCHFFI M, JAFARI S. Hydraulic flow units, depositional facies and pore type of Kangan and Dalan formations, south pars gas field, Iran. Journal of Nature Gas Science Engineering, 2015, 23: 171-183. DOI:10.1016/j.jngse.2015.01.036
[12]
SHAN L, CAO L, GUO B. Identification of flow units using the joint of WT and LSSVM based on FZI in a heterogeneous carbonate reservoir. Journal of Petroleum Science and Engineering, 2018, 161: 219-230. DOI:10.1016/j.petrol.2017.11.015
[13]
陈哲, 陆军, 赵耀辉, 等. 绥靖油田杨19区延91油藏流动单元划分合理性研究. 石油与天然气地质, 2015, 36(3): 497-503.
CHEN Zhe, LU Jun, ZHAO Yaohui, et al. Research on flow unit division rationality of Yan 91 oil reservoir in Yang 19 block of Suijing oilfield, Ordos Basin. Oil & Gas Geology, 2015, 36(3): 497-503.
[14]
GHIASI-FREEZ J, KADKHODAIE-ILKHCHI A, ZIAII M. Improving the accuracy of flow units prediction through two committee machine models: An example from the south pars gas field, Persian gulf basin, Iran. Computers & Geosciences, 2012, 46: 10-23.
[15]
任梦怡, 胡光义, 范廷恩, 等. 秦皇岛32-6油田北区新近系明化镇组下段复合砂体构型及控制因素. 岩性油气藏, 2022, 34(6): 141-151.
REN Mengyi, HU Guangyi, FAN Tingen, et al. Composite sand body architecture and controlling factors of the lower Minghuazhen Formation of Neogene in northern Qinhuangdao 32-6 Oilfield. Lithologic Reservoirs, 2022, 34(6): 141-151.
[16]
张皓宇, 李茂, 康永梅, 等. 鄂尔多斯盆地镇北油田长3油层组储层构型及剩余油精细表征. 岩性油气藏, 2021, 33(6): 177-188.
ZHANG Haoyu, LI Mao, KANG Yongmei, et al. Reservoir architecture and fine characterization of remaining oil of Chang 3 reservoir in Zhenbei Oilfield, Ordos Basin. Lithologic Reservoirs, 2021, 33(6): 177-188.
[17]
王石, 万琼华, 陈玉琨, 等. 基于辫状河储层构型的流动单元划分及其分布规律. 油气地质与采收率, 2015, 22(5): 47-51.
WANG Shi, WAN Qionghua, CHEN Yukun, et al. Flow units division and their distribution law based on braided river reservoir architecture. Petroleum Geology and Recovery Efficiency, 2015, 22(5): 47-51.
[18]
KOHONEN T. Self-organized formation of topologically correct feature maps. Biological Cybernetics, 1982, 43(1): 59-69. DOI:10.1007/BF00337288
[19]
任培罡, 夏存银, 李媛, 等. 自组织神经网络在测井储层评价中的应用. 地质科技情报, 2010, 29(3): 114-118.
REN Peigang, XIA Cunyin, LI Yuan, et al. Application of self-organizing neural network to logging reservoir evaluation. Geological Science and Technology Information, 2010, 29(3): 114-118. DOI:10.3969/j.issn.1000-7849.2010.03.018
[20]
仲鸿儒, 成育红, 林孟雄, 等. 基于SOM和模糊识别的复杂碳酸盐岩岩性识别. 岩性油气藏, 2019, 31(5): 84-91.
ZHONG Hongru, CHENG Yuhong, LIN Mengxiong, et al. Lithology identification of complex carbonate based on SOM and fuzzy recognition. Lithologic Reservoirs, 2019, 31(5): 84-91.
[21]
王亚, 杨少春, 路研, 等. 基于测井岩石物理相识别的低渗透储层评价方法: 以东营凹陷高青地区蒙阴组上段为例. 中国矿业大学学报, 2018, 47(6): 1264-1275.
WANG Ya, YANG Shaochun, LU Yan, et al. Evaluation method of low permeability reservoirs based on logging petrophysical facies identification: A case study of the upper member of Mengyin formation in Gaoqing area, Dongying depression. Journal of China University of Mining & Technology, 2018, 47(6): 1264-1275.
[22]
许建华, 蔡瑞. 有监督SOM神经网络在油气预测中的应用. 石油物探, 1998, 37(1): 71-76.
XU Jianhua, CAI Rui. Application of the supervised SOM neural network to oil and gas prediction. Geophysical Prospecting for Petroleum, 1998, 37(1): 71-76.
[23]
WANG Ya, LU Yan. Diagenetic facies prediction using a LDAassisted SSOM method for the Eocene beach-bar sandstones of Dongying Depression, East China. Journal of Petroleum Science and Engineering, 2021, 196: 108040. DOI:10.1016/j.petrol.2020.108040
[24]
宁正福, 赵洋, 程林松. 基于因子分析的流动单元研究. 中国石油大学学报(自然科学版), 2012, 36(4): 107-111.
NING Zhengfu, ZHAO Yang, CHENG Linsong. Study on flow units based on factor analysis. Journal of China University of Petroleum(Edition of Natural Science), 2012, 36(4): 107-111.
[25]
毛雪莲, 朱继田, 姚哲, 等. 琼东南盆地深水区中央峡谷砂体成因与展布规律. 岩性油气藏, 2017, 29(6): 60-68.
MAO Xuelian, ZHU Jitian, YAO Zhe, et al. Sandbody genesis and distribution regularity of central canyon in deepwater area of Qiongdongnan Basin. Lithologic Reservoirs, 2017, 29(6): 60-68.