韩元红，徐旭辉，陆建林，朱建辉，彭金宁，武英利，邱岐，王东燕

1.中国石油化工股份有限公司石油勘探开发研究院无锡石油地质研究所，江苏无锡 214151 2.国土资源部煤炭资源勘查与综合利用重点实验室，西安 710021

0 引言

苏北盆地是我国南方陆上重要的含油气盆地，古近系阜宁组是苏北盆地重要生油层系和勘探目的层[1]。对于苏北盆地油气勘探及研究工作早在上世纪五十年代已经展开，在近半个世纪的研究中取得显著进展和丰富成果[2-7]。但是，由于阜宁组湖盆边缘相缺失，加之阜宁期原型盆地是受控凹边界断裂所制约的箕状断陷，还是更广阔发育的统一陆相湖盆，一直存在争议，使得对于阜宁组物源的研究十分困难。目前，对于苏北盆地东台坳陷古近系碎屑岩物源一直存在争议。邱旭明等[8]认为来自北或北东区域。徐田武等[9-10]通过重矿法和碎屑岩类分析法指出提出主物源来自于鲁苏隆起和通杨隆起，次级物源来自于建湖隆起。周健等[11]通过对高邮凹陷古近系戴南组物源分析认为低凸起提供物源。王军等[12]通过碎屑岩岩屑类型分析，认为张八岭隆起为古近系碎屑岩主要物源区。

层序地层学相关理论的得出和计算机技术的不断发展为数值模拟研究奠定了基础。20世纪70年代以来，地层沉积演化数值模拟技术得到大力发展，2000年以后，出现了以法国石油研究院研发的DIONISOS和澳洲联邦科学与工业组织Griffiths团队研发的SEDISIM为代表的三维地层沉积演化模拟软件。三维地层沉积演化数值模拟技术的优势在于不仅可以重现一定时期地层沉积演化过程，精细刻画地层和岩相的时空展布，还可以验证地质模型及地质推论的合理性。

本次研究尝试运用三维地层沉积演化模拟技术，优化选取各类模拟参数，以目前各流派物源观点为依据，建立多组物源单因素对比模拟实验，验证各类物源观点的合理性，讨论苏北盆地东台坳陷阜宁组物源位置、方向及可能源区。

1 研究区概况

苏北盆地位于江苏省东北部，属苏北—南黄海盆地西部的陆上部分(图1a)，总体走向北东，延伸长度大于260 km，面积约38 000 km2。苏北盆地南北分别由东台和盐阜两大坳陷组成，西北临鲁苏隆起，南邻通扬隆起(图1b)，北部的滨海隆起与南黄海盆地的中部隆起相连。陆上的盐阜坳陷、东台坳陷向海域延伸与南黄海盆地南部坳陷相接。这些隆起、坳陷与海域的北部坳陷及陆上的建湖隆起共同组成了统一的苏北—南黄海盆地。控制盆地发育演化的断裂与隆起主要为嘉山响水断裂、杨村断裂、六合断裂、江都—吴堡断裂、海安—南港断裂和建湖隆起。

图1 苏北盆地构造位置与格架图Fig.1 Map showing regional location and tectonic framework of Subei Basin

阜宁组位于泰州组之上，戴南组之下，并进一步划分为四 个岩性段，自下而上分别为阜宁组一段(E1f1)、二段( E1f2) 、三段( E1f3) 和四段( E1f4) ，各岩性段之间大都为整合接触。阜宁组整体由下而上呈砂—泥—砂—泥的韵律变化。

2 模拟方法及参数

先进可靠模拟方法的使用和优化合理模拟参数的选取是搭建理想模型，获取可靠结果的基础。

2.1 模拟方法

本次研究借助DIONISOS展开对苏北盆地东台坳陷阜宁期的沉积演化模拟。 DIONOSOS是目前国际上最成熟的3D地层沉积演化正演模拟软件，可以按照勘探或地质评价的尺度，定量描述各类积环境下的碎屑岩系统盆地三维层序地层格架。它的优点在于在层序地层分析的过程中，可以采用实验的方法去考虑更多控制层序发育的因素，这有助于研究在地层概念模拟中体现不出来的或被假设的一些因素对层序的控制作用。DIONISOS中控制地层和岩相展布的因素主要包括三个方面，每一方面又由若干要素控制。第一方面为可容空间变化，主要受控于水平面变化和沉降。第二方面为沉积物供应，要素包括物源位置、物源供应速率、河流流量、碳酸盐产率和剥蚀率。第三方面沉积物搬运，包括河流搬运、波浪搬运和重力流搬运三方面。近年来，利用DIONISOS 定量分析控制地层和沉积相展布的主控因素的论著十分常见[13-18]。

2.2 参数获取

模拟成果的可靠性取决于地质模型的合理性以及模拟参数的可靠性，本次研究经过多次调试，多组对比模拟实验的调试，选取合适合理的数据参数搭建模型，获得最终模拟成果并进行讨论。数据参数主要包括可容空间、沉积物供应和沉积物搬运三方面。

本次研究选取了E1f2～E1f4为模拟层段。E1f1与E1f2之间的界面为构造界面，是从地震剖面明显识别所得的界面，而E1f2、E1f3、E1f4之间的界面为岩性界面，在地震剖面上没有明显可识别的界面，是由研究人员根据地质录井分层所得。因此，以E1f2～E1f4作为模拟研究层段，内部连续性更强，可使模拟结果更可靠。以E1f2、E1f3和E1f4为主要目的层系讨论物源，不影响结果分析。

2.2.1 可容空间

可容空间是指沉积物表面与沉积物基准面之间可供沉积物充填的所有空间。此次研究所用地层沉积演化模型中，可容空间主要参数包括起始水深、沉降和相对湖平面变化。

(1) 起始水深

起始水深决定着起始可容空间的大小，是模拟开始的关键。目前，对于古水深恢复大多是水深相对趋势或半定量水深的尝试。常用恢复方法：1)根据沉积相图[19]；2)根据地层厚度图[20]；3)根据古生物分布特征[21]；4)根据录井资料，如自然伽马[22]；5)元素地球化学方法[23]。本次研究模拟E1f2～E1f4的沉积演化过程，模拟起始水深即为E1f2沉积开始前湖泊水深。此次模拟所用的水深数据以地层厚度图、前人研究所得沉积相图(图2a)、地球化学指标(V/V+Ni)面上等值线图(图2b)为依据，设置了东深西浅的水深(图2c)。

(2) 阜宁组二段—阜宁组四段沉降

一定时期内，沉降是影响盆地可容空间变化最重要的因素，控制着盆地内沉积物供应和搬运，从而决定着地层构架和沉积演化。沉降包括构造沉降和压实沉降，本次模拟的层段为E1f2～E1f4，层内部为岩性界面，没有明显的不整合面(剥蚀界面)，因此用压实恢复过的原始地层厚度来表示盆地演化沉降。前期数据准备过程中，分别统计了E1f2(275口)、E1f3(514口)、E1f4(462口)钻井残余地层厚度(图3)，并根据压实曲线(图4)进行压实恢复作为沉降数据。

(3) 相对湖平面变化

相对湖平面变化是影响可容空间的另一重要因素。本次研究借鉴前人通过古生物等手段对湖平面升降和气候周期性变化的分析[8]：E1f2和E1f4是全盆相对湖平面最大上升期，E1f4发生的湖平面上升(亦称湖侵)的规模比E1f2大，E1f3中期具有局部性的或湖平面相对上升幅度较小的湖平面上升。模拟设置E1f2、E1f3和E1f4时间间隔分别是60.5～58 Ma，58～56 Ma和为56～55 Ma，设置相对湖平面变化如图5。

图2 模拟起始水深及其依据图件a.阜二段沉积相图[8]；b .V/V+Ni等值线图；c.模拟起始水深Fig.2 Sedimentary facies(a), isoline of V/V+Ni(b) and its indicating inital bathymetry(c) of E1f2 in Dongtai depression, Subei Basin

图3 苏北盆地东台坳陷E1f2、E1f3 和E1f4 残余地层厚度图Fig.3 Map showing stratum thickness of E1f2、E1f3 and E1f4 in Dongtai depression, Subei Basin

图4 苏北盆地东台坳陷阜宁组孔深曲线Fig.4 Relation between porosity and depth of Funing Formation in Dongtai depression, Subei Basin

图5 苏北盆地E1f2 ～ E1f4相对湖平面变化[8]Fig.5 Changes of lacustrine level during E1f2 to E1f4, Subei Basin[8]

2.2.2 沉积物供应

沉积物供应是沉积演化过程的物质基础，主要包括物源供应量、物源供应位置和物源砂体组成三部分内容。

(1) 物源供应位置

本次模拟物源位置最终设定是在文献资料的查阅、常规方法物源研究的基础上，还进行了对比模拟，优选选取，具体可见下文第三部分内容的分析。最终模型设置了北东区域和西南张八岭区域两支物源，且模拟过程中E1f2～E1f4整个过程中物源供应位置的设置保持不变。

(2) 物源供应量

本次模拟按照E1f2、E1f3和E1f4三个阶段进行物源输入设置。上文统计所得E1f2、E1f3、E1f4的地层厚度，并经过压实矫正，得出原始地层厚度，并经过“供应量=平均地层厚度×面积”计算，得出每一时期总供应量。同时，上文提到此次模拟设置物源两支(北偏东一支和南西一支)，物源总供应量按照常规方式研究所得扇体大小进行分配，西南支占35%，北东一支占65%。

(3) 物源岩性(砂体)组成

本次研究根据面上多口井阜宁组各层段岩性观察，得出阜二、阜四泥岩为主，含少量粉砂和细砂岩。阜三段以砂岩为主，包括粉砂、细砂和中砂，泥岩含量少于阜二、阜四段。因此，本次模拟研究设置了砂，粉砂和泥三种碎屑岩类型。同时，此次模拟运用的是碎屑岩和砂岩混合系统，除了物源设置的碎屑岩物源外，还有自生灰岩的生成。

碎屑岩组成随着时间的不同设置不同。整体而言，阜三段砂岩含量最高，阜二和阜四泥岩含量相对较高。整个阜二段起始砂岩含量最高，随后慢慢开始降低，阜四段仅在起始段含砂岩，后期基本全为泥岩。阜三段按照三分法的原则，分为三个旋回。

2.2.3 沉积物搬运

沉积物搬运由搬运系数、流量和坡度共同决定的。坡度是工区和地形的参数。流量是浓度的参数，在中国，浓度最小的河流在阿尔泰山脉，在旱季浓度为0.1 g/L。由于湖泊的扩散效应，水下分流河道浓度应低于0.1 g/L。因此，我们此次取值浓度为0.065 g/L。搬运系数由沉积物组成决定，沉积物组成参考各层段沉积相分析[8]。本次研究，物源供应速率为0.36 km3/Ma，因此，河流流量为3.6 m3/s。

3 结果及验证

在其他参数一定的条件下，通过不同物源方向的设定，模拟各流派物源方向条件的砂体展布特征，从而讨论物源方向。同时，建立了以北部滨海隆起为物源的全盆地层沉积演化模型，讨论中央建湖隆起是否为南部东台坳陷提供物源。最终，以物源位置单因素对比模拟和全盆滨海物源模拟结果为参考，建立了以北偏东和南西两支沉积物输入的三维地层沉积演化模型，并将模拟结果进行了验证。

3.1 物源位置单因素对比模拟

建立多组物源单因素对比模拟，模拟不同观点物源设置下的沉积相展布特征，并将其与常规方法做出的沉积相展布图对比分析，选择最正确的物源作为本次模拟的物源方向，验证物源。 以不同学者研究主张的不同物源方向为依据， 共设计了如表1所示的4组对比模拟，分别设定物源来自于N、SW、NE和NW方向的区域，且每组分别选用物源宽度为50 km和90 km进行模拟，各组对比模拟之间除物源方向有差异之外，其余参数均相同。

对比模拟得出砂体展布图(图6)，对应上述表1中各组模拟，自左至右分别为物源来自于N、SW、NE和NW模拟所得，第一列均为物源宽度为50 km模拟所得砂岩展布图，第二列则为物源宽度为90 km模拟所得砂岩展布图。

表1 对比模拟设计一览表

3.2 滨海隆起供源模拟

东台坳陷阜宁组源区的确定一直是本区研究的难点问题，有部分研究者认为阜宁期东台坳陷与盐阜凹陷为统一的盆地，物源来自于北、北东区域，滨海隆起可能提供物源[8]，也有研究者[9-10]则认为建湖隆起提供物源。因此，建湖隆起在古近系对苏北盆地盐阜坳陷和东台坳陷的分割作用大小，建湖隆起是否为其南部东台坳陷就近提供物源，尚未得到有利的解释。本研究设计了以滨海隆起作为物源区的全盆模拟，以期为这一争议的问题提供依据。

全盆模拟所用参数获取方式与东台坳陷所获参数一致， 模型设定建湖隆起和其他低凸起为水下隆起，设定物源来自于北部滨海隆起。

模拟得出若使得南部东台坳陷有砂体的展布，物源量设置须不断调试增大，使得足够物源供应。北部盐阜坳陷将演变为陆上部分，阜三段水深演变如图7所示。且砂体展布沿着湖盆北边缘展布，很难得出与沉积相图一致的砂体展布图。

图6 对比模拟所得苏北盆地东台坳陷阜三段砂岩展布Fig.6 Sandstone distribution of E1f3 in Dongtai depression, Subei Basin from the comparative simulation

图7 滨海隆起为物源模拟得出全盆水深演化Fig.7 Bathymetry evolution of Subei Basin from simulation when material sources from Binhai uplift

3.3 北偏东和南西双物源模拟

以物源位置单因素对比模拟和全盆滨海物源模拟结果为参考，模拟最终建立了以北偏东、南西两支物源供应的模型。最终模拟得到苏北东台坳陷E1f2～E1f4三维地层构架(图8)、剖面(图8)、单层岩相图(图9a)和单井岩性柱状图(图10)等。

3.4 模型验证

鉴于沉积演化模拟一些参数本身存在一定的不确定性，例如水深只能基于沉积环境和沉积相分析等给出的相对特征值。目前国际上对于模拟结果的验证多采用单井岩性对比和地层厚度而展开[19]。

选取自东至西6口井模拟所得单井岩性与实际录井岩性进行对比(图10)，得出单井岩性模拟所得与录井岩性较为一致，模拟设置水深特征与实际录井和测井特征匹配较好。拉平连井对比分析，得出地层厚度和旋回性也与实际地质模型十分一致。另外，通过对模拟与实测单井地层厚度对比，得出地层厚度误差值均在±10%以内(表2)，在误差范围以内。同时，模拟所得岩相图(图9a)与实测沉积相图(图9b)十分一致。综合验证了模拟的合理性。

表2 单井模拟与录井地层厚度对比

图8 模拟所得苏北东台坳陷E1f2～E1f4三维地层构架Fig.8 3D strata structure of E1f2～E1f4 in Dongtai depression, Subei Basin

图10 苏北东台坳陷E1f2～E1f4水深、岩性连井图Fig.10 Bathymetry and lithology features of wells of E1f2～E1f4 in Dongtai depression, Subei Basin

4 讨论

通过分析物源位置单因素对比模拟得出的砂体展布特征(图6)，得出同位置方向的两组模拟设置的物源宽度都足够大，使得砂体可以在垂直物源流向方向上扩散到最大，故宽度只在较小程度上影响砂体展布特征，砂体整体展布特征主要由物源位置和方向控制。同时，分别将四组模拟所得阜三段砂岩展布特征(图6)与前人常规分析研究所得阜三段沉积相图(图9b)对比分析，得出模拟所得第一(N)、三组(NE)砂体展布与沉积相图北部三角洲前缘相较为符合，第三组最为相似，西南区域砂体展布规模小于第二组(SW)模拟所得砂体展布。第四组(NW)砂体展布与沉积相图存在很大差异，故物源来自于北西区域的观点不合理。因此，推测这一时期沉积物源主要来自于北偏东或北区域，另有部分来自于西南区域。

全盆以北部滨海隆起为物源，将中央建湖隆起作为水下隆起模拟，得出北部盐阜坳陷逐渐被砂岩填充，演变为陆上部分。这与北部盐阜坳陷阜宁组地层厚度与南部东台坳陷较为一致，泥岩广泛展布的沉积事实不符合。因此反证推测沉积物可能就近来自于中央建湖隆起。

基于物源位置单因素模拟和滨海隆起作为物源的模拟结果分析，最终建立了北偏东和南西双物源模拟模型，得出模拟结果指示与地质模型地层沉积特征吻合度高，且模型验证显示良好。

综合上述三类模拟，推测苏北盆地东台坳陷阜宁组物源可能主要来自建湖隆起和西南区域张八岭隆起。东台坳陷阜三段砂岩沉积物粒度较细，整体以粉砂岩为主，可能与建湖隆起原始沉积物细粒有关。建湖隆起整体缺少阜宁组地层，以往对于其形态、规模及隆起时间等一无所知。由于苏北盆地阜宁组砂体展布显示出东部区域宽且偏南，西部窄且偏北特征，误导有人认为物源来自于湖盆北东部的南黄海区域，而结合本次模拟过程及结果认为砂体的这一展布特征可能与建湖隆起的位置、形态和规模有关，推测建湖隆起东部位置偏南，宽度较大，中西部位置偏北，宽度较窄，建湖隆起中西部高度较小，东部隆起高度远大于中西部。苏北盆地北部盐阜坳陷物源来自于北部滨海隆起，而此次研究认为南部东台坳陷物源来自于中央建湖隆起，推测当时苏北盆地的构造格局为建湖隆起分割南北，盐阜坳陷和东台坳陷在阜宁期相对独立，各自有着相对独立的物源供应系统。

5 结论

通过三维地层沉积演化模拟验证东台坳陷阜宁组物源，认为苏北盆地东台坳陷阜宁组物源来自于坳陷北部建湖隆起和西南区域张八岭隆起。苏北盆地在阜宁期被建湖隆起分割，盐阜坳陷和东台坳陷相对独立，各自有着相对独立的物源供应系统。