基于SSOM的流动单元划分方法及生产应用——以渤海湾盆地F油田古近系沙三中亚段湖底浊积水道为例.pdf

第 36 卷 第 2 期2024 年 3 月岩性油气藏LITHOLOGIC RESERVOIRSVol.36 No.2Mar.2024王亚等：基于 SSOM 的流动单元划分方法及生产应用收稿日期：2023-02-05；修回日期：2023-08-23；网络发表日期：2023-09-13基金项目：国家科技重大专项“低渗-致密油藏描述新方法与开发模式”（编号：2017ZX05009001）资助。第一作者：王亚（1994），男，博士，工程师，主要从事海上油气田开发地质研究工作。地址：（300452）天津市滨海新区海川路 2121 号海洋石油大厦 B 座。Email：。文章编号：1673-8926（2024）02-0160-10DOI：10.12108/yxyqc.20240215引用：王亚，刘宗宾，路研，等.基于 SSOM 的流动单元划分方法及生产应用：以渤海湾盆地 F 油田古近系沙三中亚段湖底浊积水道为例 J.岩性油气藏，2024，36（2）：160-169.Cite：WANG Ya，LIU Zongbin，LU Yan，et al.Flow unit division based on SSOM and its production application：A case study ofsublacustrine turbidity channels of middle Es3in F oilfield，Bohai Bay Basin J.Lithologic Reservoirs，2024，36（2）：160-169.基于 SSOM 的流动单元划分方法及生产应用以渤海湾盆地 F 油田古近系沙三中亚段湖底浊积水道为例王亚，刘宗宾，路研，王永平，刘超（中海石油（中国）有限公司天津分公司，天津 300452）摘 要：在储层构型级次划分的基础上，归纳总结了渤海湾盆地 F油田古近系沙三中亚段湖底浊积水道各级渗流屏障与构型界面的关系，并采用监督模式下的自组织映射神经网络算法开展流动单元定量评价，明确了构型模式控制下的流动单元分布规律。研究结果表明：渤海湾盆地F油田沙三中亚段储层可划分为、等4 类流动单元。基于 SSOM 算法的流动单元识别结果表现出较高的预测准确率，其中，256组训练样本的整体回判准确率为 82.81%，110 组测试样本的识别准确率为 80.91%，能够满足地质油藏研究的需求。垂向上，不同类型的单一水道内部发育的流动单元类型不同，造成流动单元垂向分布的差异性。类、类优质流动单元发育在浊积水道体系演化的中期，主要分布在二类单一水道；类、类流动单元发育在浊积水道体系演化的早期和晚期，其中，类流动单元分布广泛，在一类、二类、三类单一水道均有分布，类流动单元发育在一类、三类单一水道。侧向上，受水道体系不同沉积演化阶段的影响，流动单元的侧向分割方式不同。流动单元与非渗透层接触样式发育在浊积水道体系早期沉积旋回中，接触界面明显，属于一类水道沉积砂体；流动单元间的相互接触样式则发育在水道体系中晚期沉积旋回中，属于二类、三类水道沉积砂体。平面上，受单一水道侧向迁移和垂向加积作用的影响，复合水道内部流动单元平面展布差异性明显。类流动单元在单砂体内部广泛发育，井间连续性好，在浊积主水道、浊积水道及水道漫溢沉积砂体处均有发育；渗流能力较好的类、类流动单元分布局限、连续性差，仅在浊积水道主流线方向及主水道砂体处有分布，呈不连续的点状或带状分布；类流动单元呈环带状分布在类流动单元的外缘，在浊积水道漫溢沉积砂体处发育。关键词：SSOM 算法；浊积水道；储层构型；流动单元；自组织映射神经网络；监督模式；沙三中亚段；古近系；渤海湾盆地中图分类号：TE345；P618.3文献标志码：AFlow unit division based on SSOM and its production application：Acase study of sublacustrine turbidity channels of middle Es3inF oilfield，Bohai Bay BasinWANG 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王亚等：基于SSOM的流动单元划分方法及生产应用2024 年161Formation（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 controlof architecture models was clarified.The results show that：（1）The reservoirs of middle Es3in F oilfield of BohaiBay 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%for256 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 typeand flow units are developed in the middlestage 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 mainlydeveloped in categoryand 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 non-permeable 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 themiddle 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 flowunits in the composite channel was significantly different.The type flow units are developed widely in the singlesand body with good continuity between wells，and they are developed in the main turbidity channels，turbiditychannels and overflowing sand bodies.The distribution of typeand flow units with good seepage capacity islimited with poor continuity，and they are distributed only in the direction of the main stream line of turbiditychannels 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 areawhere the turbidity channels overflowed the sand bodies.Key words：SSOM algorithm；turbidity channel；reservoir architecture；flow unit；self-organizing map neuralnetwork；supervised mode；middle Es3；Paleogene；Bohai Bay Basin0引言浊积扇低渗透储层地下构型结构复杂，内部砂体形态、规模、延伸状况存在明显差异，砂体内广泛发育的隔层和夹层使得厚油层被分割成了多个连通或不连通的渗流单元。目前，已有不少学者开展了浊积扇储层构型模式的研究1-3，国外以 Mutti等4的五级划分方案应用最为广泛，而国内比较有代表性的是林煜等5和李志鹏等6的六级划分方案。这些浊积扇储层构型模式的研究为油田开发初期的井位部署提供了重要的指导作用，随着油田开发进入中后期，层间、层内的干扰导致开发效果变差，构型模式的认识难以对该阶段储层油水运动规律的表征形成指导作用，亟需在构型模式的基础上明确不同渗流单元的空间展布，以期为完善开发井网和注采结构调整提供重要依据。流动单元是表征储层渗流性能好坏的地质单元，自Hearn等71984 年首次提出流动单元的概念以来，被广泛应用于储层精细油藏描述8-9。流动单元的划分方法由最早的统计学方法，如线性判别分析、聚类分析、模糊数学逐渐向人工智能（支持向量机、人工神经网络等）领域迈进10-12，其评价参数也由单一的静态参数向多参数融合分析转变13。以往研究表明，2 种或多种方法的组合通常比单一方法具有更好的识别效果12-14。精细的储层构型模式研究是渗流屏障识别及连通体划分的基础15-17，对流动单元研究至关重要。一方面，储层构型模式研究可以找到识别和预测的成因单元，并在此结构内部依据岩性、物性及渗流特征等开展流动单元识别；另一方面，构型模式研究与渗流屏障识别关系密切，能够为渗流屏障识别和级别划分提供理论参考依据。然而，以往对基于储层构型模式的流动单元研究相对较少。以渤海湾盆地 F 油田沙三中亚段湖底浊积水道储层为例，综合运用岩心、测井及生产动态等资料，以取心层段定性划分的流动单元作为训练样本，采用监督模式下的自组织映射神经网络算法与因子分析法相结合的方式开展未取心层段的流动单元定量识别，并对构型模式控制下的流动单元分布规律进行分析，以期为后续剩余油挖潜、优化注162岩性油气藏第 36 卷第 2 期采井网及提高开发效果等提供依据。1地质概况F 油田位于渤海湾盆地东南部，构造断裂活动较弱，构造形态相对简单，是一被断层切割的穹隆背斜构造，东部由北东向和南东向的 2 条断层交叉遮挡，对砂体沉积及油气聚集起控制作用（图 1a，1b）。沙三中亚段 4 砂层组为其主要含油层系，垂向上根据隔夹层发育特征可划分出 4 个小层以满足油田开发。目的层油藏埋深为 2 8003 100 m，为典型的湖底浊积扇沉积，以中扇、外扇亚相为沉积主体（图 1c）。其中，浊积水道砂体发育，单砂层厚度为 6.613.7 m，平均为 10.7 m，垂向上水道砂体多期叠置，平面上呈大面积连片分布。受浊积水道垂向加积和侧向迁移变化快的影响，对于储层空间展布特征及连通性的研究难度较大。图 1渤海湾盆地 F 油田古近系沙三中亚段构造特征（a，b）及岩性地层综合柱状图（c）Fig.1Structural features（a，b）and stratigraphic column（c）of middle Es3in F oilfield，Bohai Bay Basin研究区沙三中亚段4 砂层组浊积扇沉积特征复杂、岩性多样，发育含砾砂岩、中粗砂岩、中细砂岩、粉砂岩、泥质粉砂岩及泥岩。储层岩石类型主要为岩屑长石砂岩，碎屑颗粒中石英体积分数占39.67%，长石体积分数占 33.60%，岩屑体积分数占26.73%，成分成熟度较低。砂岩分选中等，分选系数为 1.31.9，粒度中值为 0.120.16 mm。粒间填隙物以碳酸盐胶结物为主，呈孔隙式胶结，体积分数平均为 8.7%。受埋藏压实作用和碳酸盐胶结作用的影响，储层物性、渗流能力均较差，储层孔隙度平均为15.9%，渗透率平均为 11.6 mD，且不同小层间非均质性较强（图 2a）。开发进入中后期，受储层非均质性的影响，层间、层内干扰显著，注水开发效果变差，各层的采出程度差异大，造成水窜严重（图2b，2c）。图 2渤海湾盆地 F 油田古近系沙三中亚段各小层的渗透率（a）、相对吸水强度（b）和采出程度（c）对比Fig.2Comparison of permeability（a），relative water absorption strength（b）and recovery degree（c）of sublayers of middle Es3in F oilfield，Bohai Bay Basin渗透率/mD15105010.013.7（a）7.18.342434445小层相对吸水强度/%10503.98.1（b）4.17.242434445小层采出程度/%252015105013.1（c）42434445小层22.216.521.7地层系古近系组沙河街段沙三中砂层组4小层42434445GR/API10150深度/m录井岩性R25/（m）0.2200R4/（m）0.2200沉积相微相浊积水道水道间浊积水道水道间浊积水道水道间浊积水道亚相中扇相浊积扇2 8202 8402 860-2 800-2 840-2 880海拔/mF15-26F17-17F15-32F17-11F17-13F14-14F16-15F14-19F14-17F13-15F13-13F12-12F10F13-9F11-13F10-2F12-8F8-12F11-9F10-8F13-19F11-15F12-16F12-1F12-14F12-18F11-17F10-16F10-14F10-3F10-10F10-1F9-11F8-10F9-7F8-8F7-9F7-7F8-6F7-11F7-13F7-19F6-16F6-12F6-10F29F5-9F5-11F5-13F5-15F10-4F4-12F4-10F14-x102F8-14F291F8-22F8-16-2 860-2 850-2 840-2 830-2 820-2 810F8-18-2 800-2 7900250500 m0 10 km渤海湾盆地F8-18GRGRF8-14GRF8-511F8-509GRS42434445研究区井位断层构造等值线/m荧光油斑含砾细砂岩细砂岩粉砂岩钙质粉砂岩泥岩-2 800（c）（a）（b）王亚等：基于SSOM的流动单元划分方法及生产应用2024 年163因此亟需开展流动单元的相关研究，以明确不同渗流单元的空间展布，为完善开发井网、细分开发层系以及改善开发效果提供依据。2基于储层构型模式的流动单元分类2.1基于构型的渗流屏障及连通体分析根据Mutti等4和林煜等5针对浊积水道的构型划分方案，结合浊积水道储层的分布特征及其对油田开发的影响，建立了渤海湾盆地 F 油田沙三中亚段湖底浊积水道储层 8 个级次的构型界面划分方案：构型级次由高到低依次为湖底扇复合体单一湖底扇水道体系复合水道单一水道水道内部增生体或韵律层纹层组纹层。此外，依据浊积水道的成因类型和沉积模式可将研究区进一步划分为 3 类单一水道，分别为一类单一水道（以块状砂砾岩相及中粗砂岩相为主）、二类单一水道（以中细砂岩相为主）及三类单一水道（以泥质粉砂岩相为主），3 种类型的单一水道在垂向上具有不同的沉积模式（图 3）。一类单一水道形成于浊积水道体系演化的早期，是高密度浊流沉积的产物，自水道底部向顶部依次发育块状砂砾岩相、块状粗砂岩相、中粗砂图 3渤海湾盆地 F 油田古近系沙三中亚段单一水道类型及沉积模式Fig.3Single channel types and sedimentary model of middle Es3in F oilfield，Bohai Bay Basin岩夹薄层粉细砂岩相；鲍马序列的底部递变层理段最发育，顶层沉积不发育，序列模式为AB组合；自然伽马曲线（GR）和电阻率曲线（Rt）呈齿化箱形。二类单一水道形成于浊积水道体系演化的中期，属低密度浊流沉积，水道底部以块状砂砾岩相和块状中粗砂岩相为主，主体部位主要发育块状中细砂岩相及少量的薄层细砂岩，顶部主要发育泥质粉砂岩相，沉积物粒度自下而上逐渐变细，为典型的正韵律沉积；鲍马序列发育有ABCD 组合和 ABC组合；GR与Rt曲线呈钟形、箱形或箱形-钟形组合，内部存在回返。三类单一水道形成于浊积水道体系演化的晚期，属浊流消亡期的泥质浊流沉积，沉积物粒度较细，以泥质粉砂岩为主；鲍马序列发育最为完整，包括 ABCDE 组合和 ABCD 组合，以上平行纹层段的泥质粉砂岩和粉砂岩沉积为主，底部递变层理发育较少；GR与Rt曲线呈齿化严重的钟形，内部回返明显（图 3）。从水道体系、复合水道与单一水道 3 个级次开展渗流屏障分析。研究区为简单的背斜构造，钙质胶结物在储层中难以形成稳定分布的层间渗流屏障，以局部的层内渗流屏障存在，因此主要分析泥质渗流屏障。水道体系的划分界面为不同水道体系间区域稳定分布的泥质隔层，对应构型级次划分的六级（）界面，为一级渗流屏障；复合水道的泥岩屏障往往为发育广泛、分布相对稳定的垂向泥质隔层，对应构型级次划分的五级（）界面，为二级渗流屏障；单一水道间的渗流屏障为水道侧缘沉积或水道间的漫溢沉积，对应构型级次划分的四级（）界面，为三级渗流屏障；水道内部充填的泥质夹层和钙质夹层，对应构型级次划分的三级（）界A.递变层理段；B.下平行纹层段；C.流水沙纹层段；D.上平行纹层段；E.泥岩段。类型一类单一水道二类单一水道三类单一水道水道特征模式图鲍马序列模式典型岩心实例GRR25沉积成因类型高密度浊流沉积低密度浊流沉积浊流消亡期泥质浊流沉积BAAB 组合DCBACBAABCD 组合ABC 组合ABCD 组合ABCDE 组合DCBAEDCBADCECBABA5 cm5 cm5 cm5 cm5 cm5 cm164岩性油气藏第 36 卷第 2 期面，为四级渗流屏障。其中，一级、二级渗流屏障分布稳定，所限定的储层连通体之间没有流体流动；三级渗流屏障呈不连续分布，将连通体分割成了若干个连通或不连通的储集单元；四级渗流屏障所限定的储层连通体内部再无渗流屏障的存在，连通体内的渗流特征差异只与岩相充填特征有关（图 4）。图 4渤海湾盆地 F 油田古近系沙三中亚段渗流屏障级别与储层构型的关系Fig.4Relationship between seepage barriers and reservoirarchitecture of middle Es3in F oilfield，Bohai Bay Basin2.2流动单元划分基于浊积水道构型和渗流屏障的分析，明确浊积扇沉积体内部水道体系、复合水道与单一水道 3 个级次的构型要素及对应级次的渗流屏障分布，对单一水道开展流动单元划分及识别。采用流动分层指标法（FZI）法，结合储集物性、微观孔隙结构和生产动态等特征划分取心层段的流动单元类型，根据岩心物性分析数据计算的IFZ值（流动分层指标）的累计概率曲线确定 4 个主要的分布区间（图 5a）。以 FZI 法为主，结合储层的岩石物理特征和生产动态特征，可将取心层段的岩心样本划分为，等 4 类储层流动单元（图 5、表 1）。其中，类储层流动单元的储集物性和渗流能力最好，渗透率大于 12.1 mD，平均为 23.0 mD，初期产能最高，日产液量平均为 34.2 t。类储层流动单元储集物性和渗流能力相对较好，渗透率为 5.933.2 mD，平均为 12.9 mD，初期产能较高，日产液量平均为22.1 t。类储层流动单元的储集物性和渗流能力一般，渗透率为 3.628.6 mD，平均为 11.3 mD，初期产能较差，日产液量平均为 13.7 t。类储层流动单元的储集物性和渗流能力最差，渗透率平均为8.6 mD，初期产能最差，动用程度较低（表 1）。图 5渤海湾盆地 F 油田古近系沙三中亚段基于流动分层指标法（FZI）划分流动单元Fig.5Flow units division based on FZI method of middle Es3in F oilfield，Bohai Bay Basin表 1 渤海湾盆地 F 油田古近系沙三中亚段流动单元特征参数统计Table 1Characteristic parameters of different types of flow units of middle Es3in F oilfield，Bohai Bay Basin2.3单井流动单元单井流动单元能够反映储层内部的渗流特征分布，因此，在构型模式的指导下进行了单井流动单元研究。首先，在垂向上区分储层与非储层；其次，针对储层的物性特征、孔隙结构特征及流动特征差异开展单井流动单元识别。如图 6 所示，单井流动单元划分结果与初期产能分布具有较好的匹配关系。类流动单元与类流动单元初期产能注：最小值最大值（平均值）。（b）基于FZI法的孔渗交会图（a）FZI法累计概率曲线流动单元类型类类类类渗透率/mD12.153.4（23.0）5.933.2（12.9）3.628.6（11.3）3.116.7（8.6）孔隙度/%8.314.3（10.2）8.220.1（11.5）10.623.5（16.3）14.729.0（20.6）孔喉半径/m3.435.92（4.39）1.814.61（2.86）1.162.88（1.91）0.752.34（1.37）IFZ0.36（0.40）0.160.36（0.26）0.090.16（0.13）0.09（0.07）（Vsh）/%6.7553.55（26.65）4.3469.19（32.71）3.2977.35（34.15）7.5878.05（43.69）日产油量/t24.342.4（29.4）14.622.8（18.1）3.513.4（7.3）1.14.7（2.3）日产液量/t26.442.9（34.2）15.929.1（22.1）6.730.2（13.7）1.716.7（6.8）1008060402000.10.20.30.40.50.6IFZ/m孔隙度/%51015202530累计频率/%类类类类FZI 累积频率类流动单元类流动单元类流动单元类流动单元渗透率/mD100101R2=0.89R2=0.54R2=0.66R2=0.50类流动单元类流动单元类流动单元类流动单元王亚等：基于SSOM的流动单元划分方法及生产应用2024 年165较高，类流动单元初期产能一般。此外，由于同一单井的不同单一水道发育于不同的水道体系演化阶段，水道内部充填特征不同，导致流动单元类型及分布存在差异。对于不同单井的同一单层，由于水道发育位置不同，水道内部充填特征也存在差异，导致井间流动单元分布不同（图 6）。图 6渤海湾盆地 F 油田古近系沙三中亚段单井流动单元划分结果Fig.6Division results of flow units in single wells of middle Es3in F oilfield，Bohai Bay Basin3基于 SSOM 的流动单元测井定量识别3.1自组织映射神经网络算法原理自组织映射（SOM，Self-Organizing Maps）神经网络算法是 Kohonen18于1982年提出的一种无监督模式下的聚类分析算法，包括输入层和竞争层 2层神经网络结构（图 7a）。输入层的节点数即为流动单元输入样本的维数，当接收到流动单元样本的输入特征时，位于竞争层的神经元相互竞争对输入模式获取响应的机会。最终，获胜的神经元兴奋，失败的神经元受到抑制，获胜的神经元即为输入样本的流动单元分类。图 7自组织映射神经网络结构示意图19Fig.7Structure diagram showing self-organizingmap neural networkSOM 算法在岩性识别、流体预测及遥感图像分类等方面取得了较好的应用效果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 个储层评价参数。为避免维度过高导致运算困难，在进行神经网络运算之前，（b）有监督模式（a）无监督模式小层424344单砂体GR/API20110深度/mR25/（m）2200R4/（m）2200流动单元试油结论F14-17 井2 8402 8602 8802 900小层42434445单砂体GR/API20110深度/mR25/（m）2200R4/（m）2200CNL/%400DEN/（g cm-3）1.53.0流动单元试油结论F291 井2 8002 8202 8402 8602 880高产能中产能类流动单元类流动单元类流动单元类流动单元166岩性油气藏第 36 卷第 2 期采用因子分析算法，对输入参数进行了降维处理。当选取前 5 个主因子时，累计方差贡献率大于 80%（表 2）。因此，选取前 5 个主因子（F1，F2，F3，F4和F5）作为神经网络模型的输入特征；以 FZI 法获取的流动单元类型作为模型的输出类型，即可构建神经网络运算所必须的样本数据集。表 2 方差分析结果Table 2Rsults of variance analysis（2）网络结构初始化本次采用的神经网络为 1010 的对称网络，每个竞争层的神经网络节点代表一种流动单元储层参数响应特征，输出层的神经网络节点代表岩心样本的流动单元类别。对神经网络数据进行随机初始化，学习速率的初始值选 0.5，随着训练时间和迭代次数的增加逐渐减少至 0.000 1（学习率阈值）后结束训练；邻域距离的初始值设为 2，随着迭代的进行逐渐缩短邻域距离；设置输入层与竞争层之间的权值（Wij）以及竞争层与输出层之间的权值（Wjk）为 01 的相异随机数。（3）计算获胜神经元节点g本次研究改进了传统的最小欧式距离法，通过引入最大概率法以实现定量计算获胜神经元节点的目的。首先，计算输入样本 xi与所有竞争层神经元节点 j 的欧式距离DjDj=xi-wij,j=1,2,100（1）其次，根据欧式距离进行反距离加权，即可计算出神经元节点的获胜概率值P（j）P()g=maxP()j=max()Dj-j=1100()Dj-（2）式中：P（j）为单个神经元获胜的概率；P（g）为获胜神经元的获胜概率；为任意正实数（通常取 2）。（4）权系数调整训练获胜神经元邻域距离内的所有神经元，调整他们的权系数，包括获胜神经元节点邻域内所有节点与输入层的权值 Wij及其与输出层的权值 Wjk。当学习率小于阈值 0.000 1 时，训练结束。Wij=Wij+1()xi-Wij（3）Wjk=Wjk+2()Yk-Wjk（4）式中：Yk为取心井样本所属的流动单元类别；1，2均为学习率，随着训练时间及迭代次数的增加而线性降低。（5）预测模型输出与分析图 8 为反映流动单元分类结果的自组织映射拓扑图及 3D Sammon 投影图。拓扑图中保留了流动单元样本的输入特征及输出类别，每一个神经元节点都代表了一种储层参数响应特征，方格颜色代表该神经元节点所属的流动单元类别，且落在同一类流动单元内的数据点具有相同的岩石物理特征。3D Sammon 投影图是拓扑图在三维空间的投影，直观地反映了SSOM模型分类情况和各个竞争层神经元节点之间的距离。3.3流动单元识别结果与分析以取心井定性划分的 366 个流动单元样本数据作为模型训练和验证过程中的样本数据集，随机选取 70%的样本数据作为训练样本，其余 30%的样本作为测试样本对模型进行回判检验。基于SSOM 算法的流动单元识别结果表明，256 组训练样本仅出现 44 个错判，训练样本的整体回判准确率为 82.81%，且 110 组测试样本的识别准确率为80.91%，表现出较高的预测准确率（表 3）。测试样本中有 21 个样本发生错判，其中 19 个错判样本的流动单元类型被识别为相邻类别，主要是由于浊积水道非均质性强，仅靠 4 类流动单元类型难以全面表征其非均质性特征，而更加精细化的流动单元分类方案又需要大量岩心分析数据作为支撑，在有限的取心资料情况下，总体错判率能够满足地质油藏研究的需求。因此，监督模式下的自组织映射神经网络算法通过神经元竞争学习、相互监督的方式，能够输出最优的流动单元预测模型，为流动单元的研究提供了一种新的手段。因子F1F2F3F4F5初始特征值3.3842.8571.8151.5030.999方差贡献率/%26.02821.97413.95811.5597.683累积方差贡献率/%26.02848.00361.96173.52081.202提取平方和载入特征值3.3842.8571.8151.5030.999方差贡献率/%26.02821.97413.95811.5597.683累积方差贡献率/%26.02848.00361.96173.52081.202旋转平方和载入特征值3.3262.0621.8611.7831.525方差贡献率/%25.58215.86314.31213.71611.729累积方差贡献率/%25.58241.44555.75669.47381.202王亚等：基于SSOM的流动单元划分方法及生产应用2024 年167图 8自组织映射聚类网络拓扑图（a）及竞争神经元节点 3D Sammon 投影（b）Fig.8Self-organizing map topology（a）and 3D Sammon mapping projection of competing neuron nodes（b）表 3 渤海湾盆地 F 油田古近系沙三中亚段基于 SSOM 的流动单元预测精度统计Table 3Statistics of flow unit prediction accuracy based on SSOM of middle Es3in F oilfield，Bohai Bay Basin4储层构型控制下的流动单元分布规律渤海湾盆地 F 油田未取心井段流动单元的分类识别可以借助 SSOM 模型实现，在单井识别的基础上开展了井间剖面和平面流动单元分布预测24-25。4.1剖面分布规律受沉积物源、水动力条件、水道的迁移叠置等因素的影响，浊积水道体系不同演化时期发育的单一水道类型存在差异。自水道体系演化的早期晚期，单一水道类型表现出自一类向三类过渡的趋势。通过单一水道与流动单元剖面分布的对比分析发现，单一水道类型的差异性造成了流动单元剖面分布的差异性，即流动单元剖面展布特征受水道垂向演化控制作用显著，主要表现在以下 2 个方面。（1）垂向上，等 4 类流动单元在F油田4砂层组的分布差异明显。总体表现为类流动单元分布广泛，类次之，类、类流动单元发育较少。从水道体系级次来看，类、类优质流动单元发育在浊积水道体系演化的中期，类、类流动单元发育在浊积水道体系演化的早期和晚期。从单一水道来看，不同类型单一水道成因类型和沉积模式的差异性，导致不同类型单一水道内部发育的流动单元类型存在明显差异。类、类流动单元分布在二类单一水道；类流动单元分布广泛，在一类、二类、三类单一水道均有分布；类流动单元发育在一类、三类单一水道（图 9）。（2）侧向上，流动单元间呈 2 种接触样式，分别为流动单元与非渗透层（即隔层）直接接触和流动单元间的相互接触。前者发育在浊积水道体系早期沉积旋回中，接触界面明显，属于一类水道沉积砂体（如 44、45小层沉积砂体侧向上发育流动单元与非渗透层的直接接触）。后者接触界面不明显，通常显示出侧向上流动单元渗流特征的差异，发育在水道体系中晚期沉积旋回中，属于二类、三类水道沉积砂体（如 42、