浅水三角洲前缘砂体广泛发育于湖盆环境,是陆相盆地中极具勘探潜力的重要油气储集层类型(李元昊等,
2009;张昌民等,
2010;邹才能等,
2015;Zhu
et al.,
2017),但其内部构型复杂,给储集层精细预测与井位部署带来巨大挑战(王于恒等,
2024;Li
et al.,
2025)。由于浅水湖盆湖平面频繁升降,相比于正常三角洲,浅水三角洲砂体往往经历多期沉积叠置,分流河道频繁迁移改造,砂体纵向多期叠置、平面砂体相互交织,导致内部非均质性强、展布形态极为复杂(吴胜和等,
2021;Zhang
et al.,
2022;徐振华等,
2024)。前人对于浅水三角洲构型的表征主要依靠野外露头、现代沉积以及密井网等手段(纪友亮等,
2016;林承焰等,
2017;张莉等,
2017;张宪国等,
2019;Cai
et al.,
2022;Wang
et al.,
2023)。然而,在缺乏足够井控约束的稀疏井网研究区,不同成因单砂体边界难以准确厘定(邓文龙等,
2024; 孙以德等,
2024),砂体间连通性难以判别(于海波等,
2024),导致常规地质建模与地震属性预测精度显著降低,严重制约了浅水三角洲储集层的高效勘探开发。
面对沉积构型复杂、储集层非均质性强等问题,地球物理方法在薄互层,尤其是致密储集层的识别中同样存在显著技术瓶颈(Anees
et al.,
2022;曹思远等,
2023)。一方面,浅水三角洲河道砂体普遍具有层厚薄、泥砂互层频繁的特征,目标储集层厚度远低于地震资料1/4波长的垂向分辨率,传统反演方法受限于分辨率难以清晰揭示薄储集层(翟亮等,
2023;焦社宝等,
2025);另一方面,致密砂岩普遍存在孔隙度和渗透率普遍较低、波阻抗与围岩对比度低等问题,进一步削弱了地震响应特征,使储集层的准确识别和预测变得更加困难(史浩等,
2024;王武荣等,
2024)。在此背景下,近年来发展起来的波形指示反演(Waveform Indication Inversion,WII)技术凭借其对地震波形横向变化特征的敏感识别能力,已逐步成为揭示砂体展布与沉积构型的重要工具(顾雯等,
2016;王于恒等,
2024;童姜楠等,
2025; 朱绍鹏等,
2025)。WII方法通过对三维地震波形的动态分析,建立波形模式与沉积相类型之间的映射关系,有效提高了薄层储集层的分辨率与相控识别能力,尤其在稀疏井网区表现出良好的适应性(陈彦虎等,
2020;Wang
et al.,
2025a)。与此同时,自动化机器学习(AutoML)技术近年来发展迅速,作为一种新兴的数据驱动机器学习方法,有效降低了技术门槛,助力开发人员高效部署高性能模型(Shen
et al.,
2018; He
et al.,
2021;Zöller and Huber,
2021)。以AutoGluon为代表的先进框架,具备自动参数优化、多模型融合与特征选择等功能,能够从大规模、多维度数据中高效提取关键特征组合,减少人为干预,提高模型稳定性和预测精度(Erickson
et al.,
2020)。然而,尽管AutoML在其他领域应用广泛,其在地震地质领域的研究仍处于起步阶段,尤其在浅水三角洲薄层砂体的识别与定量预测方面尚属空白。
针对稀疏井网条件下浅水三角洲致密薄层砂体地震预测精度低的难题,作者提出一种面向薄层砂体识别与定量预测的新方法,融合基于自动化机器学习(AutoML)的分频多属性融合技术与波形指示反演(WII)技术,集成了自动特征选择、多模型优化与地震波形模式识别等优势,显著提升了地震对薄层砂体的响应能力与分辨率。以川中地区天府气田侏罗系沙溪庙组一段(沙一段,J2s1)浅水三角洲为研究对象,采用新方法实现了对浅水三角洲致密砂岩储集层的精细刻画,明确了单砂体的形态、叠置样式与空间规模,构建了浅水三角洲沉积构型模式,为复杂构型浅水三角洲储集层的智能预测与高效勘探开发提供了新的技术路径与理论支撑。
1 区域地质概况
四川盆地是扬子板块的次级构造单元,属于中新生代形成的大型稳定克拉通盆地,具有稳定背景下的多旋回沉积与构造演化特征,主要受基底结构控制(Yue
et al.,
2018;王翘楚等,
2025)。二叠纪前整体为碳酸盐岩地台,受元古代基底运动影响,形成川西低陡带、川北低缓带、川东高陡断褶带、川西南低褶带、川南低陡穹形带及川中平缓带6大构造单元(Wang
et al.,
2025b)(
图1-a)。天府气田简阳区块位于四川省简阳市境内,构造上位于川中平缓构造带西侧,构造形态相对简单,整体为南东向北西下倾的单斜,地形平缓(
图1-b)。整体构造作用较弱,研究区断层、裂缝发育较少(杨雨等,
2023)。
四川盆地相继沉积了中、下三叠统及其以下的以碳酸盐岩为主的海相地层和上三叠统、侏罗系以砂泥岩为主的陆相地层(杨雨等,
2023),整体地层无沉积间断(
图2-a)。其中中侏罗统沙溪庙组作为简阳区块主力产气层段,发育1套巨厚的陆相碎屑岩,岩性主要以紫红色泥岩夹细—中粒砂岩为主,厚度介于800~2200 m之间(
图2-a)。沙溪庙组内部以“叶肢介页岩”为界,分为沙一段和沙二段。根据岩性组合、沉积旋回特征,沙一段可以划分为3个亚段(J
2${s}_{1}^{1}$、J
2${s}_{2}^{1}$、J
2${s}_{3}^{1}$),根据河道砂体河道演化特征和沉积旋回特征可将沙一段进一步细分为6个小层(J
2${s}_{1-1}^{1}$、J
2${s}_{1-2}^{1}$、J
2${s}_{1-3}^{1}$、J
2${s}_{2-1}^{1}$、J
2${s}_{2-2}^{1}$、J
2${s}_{3}^{1}$)(肖富森等,
2020)(
图2-b)。
川中地区沙溪庙组整体处于浅水三角洲—湖泊—河流沉积体系(肖富森等,
2020;杨跃明等,
2022)。在沉积期间,气候呈现出干湿交替的特点,水体以淡水为主,但在局部地区也存在微咸—咸水的特征(余文强等,
2024)。
其中,沙一段沉积期的气候环境为干湿交替,整体处于温暖湿润与半干旱—半湿润之间,物源充足,发育典型的浅水三角洲内前缘沉积(冯潇飞等,
2024;赵晓明等,
2024)。相比于正常三角洲,浅水三角洲内前缘的分流河道水体能量强,河道频繁改道使前期发育的河口坝多被后期的分流河道冲刷破坏,较少被保存下来,仅在少量钻井岩心上见到残余的河口坝沉积物(杨跃明等,
2022)。
2 构型单元类型及特征
天府气田简阳地区因靠近川南物源,湖盆水体较浅且动荡,地形相对宽缓,主要发育浅水三角洲内前缘沉积,基于对取心井岩心及测井资料的分析并借鉴前人研究成果,进一步识别出水下分流河道、河口坝、溢岸和分流间湾等构型单元(
表1),本研究河口坝构型单元在内前缘指受冲刷改造后的残存坝体。
天府气田简阳地区发育的水下分流河道构成了湖盆砂体的骨架,控制着沉积体系的空间展布与储集层发育格局。该类河道为分流河道向湖盆内部水体延伸的产物。其岩性以中砂岩和细砂岩为主,常见楔状交错层理、槽状交错层理和平行层理,底部多为不规则冲刷面,局部可见定向排列的泥砾,反映出较强的水动力环境。垂向上呈粒度向上变细的正韵律,自然伽马曲线呈钟形或箱形,底界多为突变接触,顶部则为突变或渐变接触。河口坝为分流河道携带沉积物在流入湖盆水体过程中因流速骤减而沉降形成的沉积体,岩性以细砂岩和中—细砂岩为主,发育槽状交错层理、砂纹层理等,具有向上粒度变粗、泥质含量减少的反韵律特征,自然伽马曲线表现为齿化的漏斗形或箱形,顶部为突变接触,底部为渐变接触。需要指出的是,在本区仅发育局部残存的河口坝沉积,分布零星,形态保存不完整。溢岸沉积包括天然堤与泛滥期形成于水道间低洼带的细粒沉积,岩性多为细砂岩与粉细砂岩,发育砂纹层理、小型交错层理,韵律性不明显,伽马曲线表现为泥岩背景上的中—低幅齿状或指状曲线。分流间湾为分布广泛的相对静水条件下的细粒沉积,岩性以粉砂质泥岩或泥岩为主,常含碳屑及植物根茎化石,发育块状层理、水平层理,自然伽马曲线平直微齿化,缺乏明显韵律性。
3 复合砂体级次构型表征
3.1 基于自动化机器学习的分频属性融合砂体预测
在稀疏井网条件下,地震属性分析仍是预测储集层砂体分布的关键手段之一,但受调谐效应影响,振幅与频率信息存在多解性,显著降低了地震属性对薄层砂体的识别能力,尤其在构型复杂的浅水三角洲体系中更为突出(Li
et al.,
2019;岳大力等,
2025)。研究区面积约1200 km
2,共部署37口井,其中直井/斜井27口、水平井10口。井网稀疏且分布不均,平均井距约5.1 km,最小井距1.45 km,按面积折算的井网密度0.031口/km
2,显著低于成熟开发区块常见的合理井网密度(约3口/km
2)及其对应的开发井距(约0.4~0.8 km)(何东博等,
2013)。
尽管研究区三维地震资料全区覆盖,但地震资料主频约为39 Hz,频宽14~64 Hz(
图3-a)。研究区小层层厚平均约30 m,而其中单砂体厚度整体较薄(10~20 m),按主频估算的调谐厚度约30.37 m,高于研究区目标层单砂体厚度,与地震垂向分辨率处于同一量级,易受调谐效应与多解性影响,增加了薄层砂体识别与定量预测的难度。导致传统地震属性与砂体厚度的相关性较低,其中最大相关系数仅为0.64(
表2),仅能反映小层级次的宏观展布趋势,难以实现对储集层的精细刻画。鉴于地震属性与储集层参数之间普遍存在复杂的非线性关系,传统线性方法已难以满足薄层储集层定量预测的需求。
3.1.1 分频预处理与属性优选
单一频段地震属性包含地震信息有限,且受调谐现象影响,对砂体预测存在多解性,因此对原始地震数据进行分频,通过融合多个频段地震属性以提高砂体预测精度(Li
et al.,
2019)。分频技术是利用短时傅里叶变换、小波变换和S变换等算法,将地震数据由时间域转至频率域进行分解。小波变换对于高频分量采用窄时间窗提升时间分辨率,对于低频分量使用宽时间窗增强频率分辨率,通过自适应的时频窗口可以有效捕捉地震波传播中的非平稳特征(高静怀等,
1996)。因此,采用基于Morlet小波的时频变换技术,对原始地震资料在有效带宽内以10 Hz频段间距划分9个频段,结合调谐厚度与目的层5~20 m砂体厚度,优选主频为23.4 Hz(频宽14.8~32.0 Hz)、41.7 Hz(频宽26.2~57.2 Hz)、52.2 Hz(频宽34.8~69.6 Hz)的低、中、高频数据体(
图3-b),优选的分频数据体总频宽可覆盖原始地震数据,既保留了原始数据的有效信息,又减少了重复信息。
不同地震属性能够从不同角度表征储集层特征,其中振幅类属性可以有效指示岩性、刻画砂体厚度; 相位类属性有助于识别储集层边界、断层等; 频率类属性更适合反映储集层在垂向上的岩性变化(岳大力等,
2022)。本研究选取测井资料齐全的27口井,通过前期地震层位解释结果约束目的层层面属性提取时窗,分别提取各类属性,包括振幅类、相位类、能量类及频率类等21种常见的、无残缺信息的地震属性。
分析各属性与砂体厚度相关性,共得到10种与砂体厚度相关系数大于0.4的地震属性,由于部分地震属性间共线性较强(较强的相关性),因而在地震多属性融合前利用无监督聚类技术将10种地震属性进一步划分为3类(
图4),结合各类地震属性与砂体厚度的相关系数,最终优选出与砂体厚度相关性最好的3个属性进行属性融合,即标准差振幅(Standard deviation of amplitude)、平均正振幅(Average positive amplitude)和平均波形长度(Interval average: arithmetic)。其中标准差振幅主要反映局部振幅变化程度,高值常对应砂泥界面密集分布区,能够有效识别互层发育的复合砂体边界,而平均正振幅增强了与砂体充填有关的能量响应,有助于识别砂体主轴部位的强反射特征,对连续河道识别有明显优势,最后平均波形长度本质上反映局部波形稳定性,长波形常与连续泥质段对应,而波形突变区往往为砂体边界,有助于提高砂体边缘刻画精度。本研究优选3类属性具备清晰的地质意义,适用于浅水三角洲薄互层与复杂叠置构型的识别。
3.1.2 砂体定量预测及可靠性验证
本研究针对定量储集层预测与地震属性融合的需求,基于AutoGluon自动化机器学习框架构建了融合分频属性的非线性回归模型,并应用于浅水三角洲致密砂体的智能预测。整体流程包括数据预处理、模型训练与评估3个阶段(
图5)。在优选出频段和属性的基础上,分别提取标准差振幅、平均正振幅和平均波形长度3类关键属性,作为输入特征的核心变量。考虑到主频差异对应不同厚度与边界特征的响应能力,该多频多属性组合有效提升了对薄层、互层及复合砂体的适应性。模型训练阶段选取研究区27口测井资料齐全的井,其中按约3︰1比例划分为训练样本20口井与独立测试样本7口井,计算井点周围(半径25 m)属性的平均值作为输入特征,并以小层的测井解释砂厚为监督标签。
在完成井震标定、归一化与特征编码等预处理步骤后,利用AutoGluon框架进行多模型自动调参与集成学习,充分挖掘属性间的非线性耦合关系,相较于传统线性拟合方法更具泛化能力与预测精度。AutoGluon的集成机制可自动组合多种回归模型(如RF、XGBoost、NN等),并动态优化其权重与结构,显著减少人为干预,提高了建模效率和适应性。在模型评估阶段,通过设置精度阈值与交叉验证策略对模型性能进行量化,若结果不达预期,可自动回溯优化输入属性组合与超参数配置,确保预测结果稳定可靠。最终,将训练完成的模型推广至全区,得到小层砂体厚度的融合属性预测图,实现了在井网稀疏背景下对致密砂岩储集层的定量刻画。
传统地震属性预测方法多基于主频地震体,提取单一属性与实测砂厚建立线性关系。然而,在浅水三角洲沉积背景下,砂体厚度普遍小于调谐厚度,易受频散与调谐效应干扰,导致预测精度较低,难以实现薄互层砂体的精细刻画。以J
2${s}_{1-2}^{1}$小层为例,采用传统主频体提取的原始属性中,相关性最高的为标准差振幅属性(
图6-b),其与实测砂厚的相关系数仅为0.638(
图7-a),预测能力有限,尤其对薄层或窄河道砂体的响应效果较差。相较之下,基于自动化机器学习的分频属性融合方法在空间识别精度和砂体边界清晰度方面表现更优,特别是在J
2${s}_{1-2}^{1}$小层中的窄条带状河道砂体区域(
图6-a、
6-b红色箭头处),融合后河道砂体边界更加清晰,砂体更加连续。同时对于砂体厚度的预测与单井解释更为符合,融合结果更能体现砂体厚度的差异,砂体厚度预测精度显著提高,预测结果与测井解释砂厚的相关系数由原始的0.638提高到0.820(
图7-a,7-b),7口测试井的相关性也高达0.816(
图7-c),进一步说明基于自动化机器学习的分频属性融合预测砂体的可靠性。
3.2 复合砂体级次构型分布特征
天府气田简阳区块井资料较少、井网稀疏,无法采用传统的密井网解剖的思路。因此,本研究综合考虑研究区沉积特征,以单井构型解释结果及分频属性融合砂体预测结果作为约束,采用“平剖互动、井震结合”的思路,表征了复合砂体级次构型的分布特征。研究区沙一段小层内部发育单一河道沉积,在相对湿润气候条件下,受浅水湖泊的改造,河道频繁分叉改道,河道对河口坝的切叠改造较为明显,分流河道与下伏早期河口坝沉积呈冲刷接触,仅局部发育残余河口坝,导致形成以河道为骨干砂体的网状浅水三角洲,发育较多分支窄河道(
图8)。
简阳地区受龙门山断裂带南部的物源控制,整体上研究区复合砂体呈西南—东北流向,西南部河道较东北部宽,西南部河道宽度普遍在600 m以上,最大可达2000 m,而东北部河道宽度较小(约300~600 m)。垂向上,沉积旋回对砂体叠置样式具有明显控制作用。长期旋回阶段湖平面上升、可容空间增大,河道砂体呈孤立条带状,弯曲度逐渐增加。而高频的小幅周期波动(短期旋回)则控制了砂体的空间分布和叠置样式。以J
2${s}_{1}^{1}$为例,J
2${s}_{1-1}^{1}$与J
2${s}_{1-2}^{1}$小层处于短期上升旋回,分流河道规模较小,发育较多分支窄河道(
图8-a,8-b),而J
2${s}_{1-3}^{1}$小层处于短期下降旋回,河道侧向切叠明显,发育宽度较大的复合砂体叠置体(
图8-c)。
4 单砂体级次构型表征
4.1 基于地震波形指示反演的单砂体叠置样式
地震波形指示反演基于沉积学理论,通过关联地震资料横向波形变化与钻井测井曲线特征,定量反映储集层空间沉积相变规律,并通过高频成分估计预测储集层(陈彦虎等,
2020)。与传统反演方法不同的是,其在样本优选统计分析中,按空间距离排序样本点,选关键井作初始模型,进行高频成分无偏最优估计(顾雯等,
2016)。该方法利用储集层相变在地震波形上的响应特征差异,以地震波形横向变化为相控约束,协同模拟反演钻井数据与地震波形,有效降低了反演结果随机性,提升了储集层反演精度。尤其有利于复杂相带中砂体叠置关系及规模较小河道的精细刻画(Wang
et al.,
2025a)。
4.1.1 反演参数优化及可靠性分析
地震反演的核心在于构建储集层特征曲线与地震属性之间的高精度映射关系,其关键在于测井参数的优选与优化。针对沙一段致密化河道砂体纵横向叠置交织的特征,传统波阻抗参数对于砂泥岩识别存在重叠区,而无铀伽马(KTH)能够更好地区分砂岩、泥岩(
图9),故在测井曲线标准化的基础上,选取标准化后的无铀伽马曲线作为反演学习的目标曲线。
在研究区分别采用常规波阻抗地质统计学反演和地震波形指示反演方法,得到反演数据体,随机选取未参与反演学习的yq213井与yq214井提取井旁道对反演结果进行对比验证(
图10-a,10-b)。就整体而言,地震波形指示反演成果不仅与地震波形响应吻合度高,反映了反演过程充分考虑了地震信息,而且在细节上,相比于常规波阻抗地质统计学反演,地震波形指示反演分辨率更高,且与单井砂体解释吻合度高,对10 m以上砂体有稳定响应,垂向上可清晰区分砂泥薄互层,同时部分4 m左右的孤立砂体也出现了地震反演响应,可以有效区分单砂体,识别砂体叠置样式,满足稀疏井网下构型单元精细刻画的要求。
4.1.2 构型单元组合样式的地震反演响应图版
基于地震波形指示反演结果,结合砂体垂向演化规律、成因类型及单砂体侧向接触关系,建立了沙一段浅水三角洲前缘构型单元组合样式的反演图版。分析识别出6种构型组合样式,其中垂向叠置样式4种,侧向切叠样式2种(
图11)。
垂向叠置样式包括: 多期水下分流河道连续叠置型、多期水下分流河道间隔叠置型、水下分流河道—河口坝叠置型和水下分流河道—溢岸叠置型。其中,多期水下分流河道连续叠置型、水下分流河道—河口坝叠置型和水下分流河道—溢岸叠置型多形成于沉积物供给速率(A)大于可容空间增长速率的环境(即A/S>1),后期河道对前期砂体改造强,砂体叠置紧密,泥质夹层少,垂向连通性好; 而多期水下分流河道间隔叠置型常发育于A/S<1的环境,砂体间保留细粒沉积,测井曲线可见明显回返,泥质夹层发育且较厚,连通性较差。
侧向切叠样式表现为不同期次砂体在平面上迁移、切入和拼接,形成横向连接的复合砂体,主要包括水下分流河道切叠型和水下分流河道—溢岸对接型。前者发育于A/S>1的条件下,河道迁移频繁,侧向侵蚀作用强,后期砂体可大面积切叠早期砂体,形成大接触面,砂体间泥质夹层少,连通性好;后者多见于A/S≈1的环境,河道侧向侵蚀能力弱,叠合方式以对接为主,砂体间接触面积小,缺乏明显冲刷侵蚀界面,连通性较差。
4.2 单砂体级次构型定量表征
单砂体是指垂向上和平面上都连续,且与上、下砂体间有泥岩或不渗透夹层分隔的砂体(何文祥等,
2005;张庆国等,
2008)。单砂体是本研究最低的构型解剖级次,准确识别单砂体边界是构型表征的基础。在复合砂体识别与预测的基础上,进一步识别单一河道的分布边界是实现精细地质建模与储集层连通性评价的关键环节。本研究结合单井沉积相分析成果与地震波形指示反演数据,对研究区沙一段浅水三角洲前缘单一河道的垂向期次划分及侧向展布边界进行精细刻画。
4.2.1 单砂体级次构型河道期次划分
在垂向上,单一河道期次的识别主要依赖单井岩心及测井资料所反映的沉积间断面,判断砂体上下叠置关系。考虑地震资料的分辨率限制,仅识别可由分流间湾泥岩清晰区分的河道单元,即以泥岩隔层为界,划分叠置的单期河道砂体。该方法可有效剔除高频次小尺度的沉积变化干扰,保障了地震反演与单井资料之间的尺度匹配性。
在侧向上,单一河道的展布边界主要通过3种地震响应特征进行判识: (1)厚度差异特征: 河道中心部位表现为砂体厚度最大,向边部过渡过程中砂体厚度逐渐减薄,波形指示反演成果颜色从中心的亮红色逐渐过渡到浅绿色,反映出典型的河道充填形态(
图12-a,12-b中①);(2)高程差异特征: 在同一沉积期内,多个河道单元发育于不同的地貌高程上,揭示了河道相对位置的差异(
图12-a,12-b中②);(3)侧向相变特征: 在相邻河道之间可见残存的河口坝或薄层溢岸砂体充填,波形指示反演结果颜色由亮红逐渐过渡到黄色、绿色,相变过渡特征明显,为判定河道期次及单砂体边界提供了可靠依据(
图12-a,12-b中③)。
在融合多属性识别复合砂体的基础上,为进一步厘清单一河道的沉积期次,在等时地层格架约束下,采用层拉平技术,以切物源方向的波形指示反演剖面为基础,开展单期河道的期次划分分析。以J
2${s}_{1-2}^{1}$小层为例,识别出单砂体级次构型下的4期分流河道,依次为1号、2号、3号和4号河道(
图13-a),表明该小层内河道体系经历了多阶段的演化过程。在精细划分过程中,考虑到河道频繁改道和分支,同一沉积期内主干河道与其分支仍存在内部期次之分。通过分析各河道顶部至J
2${s}_{1-2}^{1}$小层顶面的垂向距离,结合地震剖面叠置关系进行期次判别。以3号河道及其分支为例,3
2分支河道顶部高于3
3分支河道,3
3分支河道又明显切叠于4号河道之上,反映3号河道早于4号河道(
图13-b)。进一步结合BB’(
图13-c)与CC’剖面(
图13-d),对各期河道的上下叠置关系进行验证,可明确1号河道早于2号,2号早于3号,3号早于4号。该方法有效厘清了复杂河道体系的沉积期次关系,提供了关键的期次划分依据。
4.2.2 单砂体级次构型定量规模
河道宽度是衡量河道规模及河道识别的重要参数,开展单砂体定量规模研究有助于地下储集层连通性预测,更好地服务于油气勘探开发(孙以德等,
2024)。在沙一段单一河道砂体级次构型表征的基础上,对各个小层的河道厚度与宽度进行了统计分析,结果表明沙一段各小层河道厚度与宽度具有良好的正相关关系(
图14),由此建立了单一河道定量构型模式,应用该经验公式可以根据沉积条件下活动河道厚度对河道宽度进行预测,这对地下单一河道识别起到了规模控制作用。
不同小层河道宽厚比差异明显,垂向上旋回对单一河道构型定量规模的控制作用明显,各个亚段表现出相似的沉积特征,底部河道沉积处于可容空间扩张阶段,河道砂体宽度厚度比值较低,地层记录上表现为向上变细的退积沉积旋回序列。随着湖盆水位下降,可容空间逐渐转换到压缩阶段,顶部河道下切能力增强,厚度逐渐增大,平面上分布广泛,整体宽厚比较大,地层沉积多表现为向上变粗的进积或加积序列。这一演化过程说明,在川中缓坡稳定的背景下,微小的构造沉降变化或湖平面振荡即可引起可容空间的阶段性增减,从而引发河道沉积形态在孤立型窄河道和叠加型宽河道之间反复循环演化(
表3)。
5 浅水三角洲三维构型模式
在结合前人研究成果的基础上建立了浅水三角洲不同亚相单砂体构型沉积演化模式。在四川盆地沙溪庙组浅水三角洲体系中,三维构型呈现出明显的顺物源分异特征。由近源向远源依次可划分为三角洲平原、内前缘、外前缘、前三角洲和滨浅湖5个亚相带(杨跃明等,
2022;邓文龙等,
2024)(
图15-a)。天府气田简阳地区因靠近川南物源,湖盆水体较浅且动荡,地形相对宽缓,主要发育浅水三角洲内前缘沉积。其中三角洲平原发育厚层分流河道和洪泛泥质沉积,内前缘以水下分流河道为主,河道频繁侧向迁移,呈“网状”展布,砂体以多期叠置、相互切叠为主; 外前缘水动力减弱,发育窄幅分流河道和舌状至朵状河口坝,砂体之间以垂向叠加和侧向对接为主; 前三角洲主要沉积粉砂质泥岩,滨浅湖发育低能泥质沉积。整体构型由近源强动力的侵蚀—搬运—堆积,过渡至远源低能的沉降—沉积控制模式(
图15-b)。
垂向上,在切物源剖面上,分流河道砂体多呈“顶平底凸”特征,表明其基底侵蚀强烈,而河口坝砂体则呈“底平顶凸”形态,反映其在静水环境中堆积发育,形成典型河道直接叠置于河口坝之上的“坝上河”构型模式。沙一段沉积时期,研究区主要发育浅水三角洲内前缘沉积,沉积物供给速率(A)和可容空间的变化(S)直接影响内前缘的各成因单砂体发育规模和叠置关系。内前缘靠近物源,沉积物供给充足,A/S>1,河道频繁侧向迁移,在平面上呈现“网状”的形态,此时河道作用强,河口坝被河道侵蚀切割,难以形成大规模河口坝,仅以残余河口坝形式存在,水下分流河道为骨架砂体,各单砂体间以垂向侧向相互切叠接触为主。随水下分流河道不断向湖盆推进,到达外前缘沉积中心,远源区沉积物供给速率下降,搬运能力减弱,A/S<1,水下分流河道规模逐渐减小,在平面呈现窄条带状,靠近湖盆中心的砂体受湖浪改造较强,河道侵蚀作用减弱,河口坝与席状砂体开始大规模发育,河口坝与席状砂体多呈舌状分布,部分地区河口坝交织成朵状坝,显示出明显的空间演化规律(
图15-c,15-d)。
6 结论
1)在稀疏井网条件下,采用基于自动化机器学习(AutoML)的分频多属性融合方法,结合地震波形指示反演(WII)技术,实现了致密砂岩薄层的精细表征,显著提高了浅水三角洲薄层砂体的预测精度,但受限于地震分辨率及调谐效应影响,陆相致密薄层砂体仍存在显著的不确定性与多解性挑战,未来将进一步探索新的正则化和约束机制来提高预测可靠度,将储集层非均质性特征和沉积地质模式融入深度学习框架,以期进一步提高稀疏井网下的复杂浅水三角洲致密气藏的储集层构型预测精度。
2)在四川盆地天府气田简阳区块沙一段主要发育浅水三角洲内前缘沉积,共识别出水下分流河道、河口坝、溢岸和分流间湾4种构型单元。基于地震波形指示反演结果,结合砂体垂向演化规律、成因类型及单砂体侧向接触关系,分析识别出多期水下分流河道连续叠置型、多期水下分流河道间隔叠置型、水下分流河道—河口坝叠置型、水下分流河道—溢岸叠置型、水下分流河道切叠型和水下分流河道—溢岸对接型共6种构型组合样式。在沙一段单一河道砂体级次构型表征的基础上,对各个小层的河道厚度与宽度进行了统计分析,结果表明沙一段各小层河道厚度与宽度具有良好的正相关关系,并通过数据统计说明了在川中缓坡稳定的背景下,微小的构造沉降变化或湖平面振荡即可引起可容空间的阶段性增减,从而引发河道沉积形态在孤立型窄河道和叠加型宽河道之间反复循环演化。
3)建立了浅水三角洲沉积构型模式,明确构型具有明显的空间分异特征,由近源强动力的侵蚀—搬运—堆积,逐渐过渡至远源低能的沉降—沉积控制。其中近源三角洲平原为厚层分流河道与洪泛泥; 内前缘A/S>1,水下分流河道呈网状、多期切叠,仅发育残余河口坝; 外前缘A/S<1,水下河道逐渐变窄,河口坝及席状砂呈舌—朵状; 前三角洲与滨浅湖广泛沉积泥质。
*国家自然科学基金项目(42272186)
国家自然科学基金项目(42302128)
国家自然科学基金项目(42412179)
中国石油大学(北京)科研基金(2462023YJRC039)