Reconstruction of Late Oligocene depositional systems in the Beikang Basin by seismic facies analysis
-
摘要: 依据位置（陆棚-陆架边缘-斜坡）、外部形态（丘状、楔状、板状等）、内幕结构（平行、斜交、发散等）和反射波属性 (振幅、频率、连续性等）4个方面的标识划分地震相的传统地震相分析法，在钻井较多的地区，可实现地震相-沉积相的准确转换；但在钻井稀少（无井）或地质条件复杂的地区，由于传统地震相分析的结果主要为环境相，在岩相判别上能力较弱，故在沉积相的转换上存在一定的缺陷和不足。通过新增地震反射波光滑性、地震反射整洁性和特殊波形3种标识，建立地震相-地震岩相-沉积相的分析方法，能有效提升岩相的辨别能力，结合沉积环境与水流体系，可实现准确的沉积体系重建，对油气勘探具重要价值。综合传统地震相分析与3种新增标识的各自优势，在北康盆地晚渐新世识别出6种典型地震岩相，即砂砾岩相、砂包泥相、砂泥岩互层相、泥包砂相、纯泥（页）岩相和火山岩相；划分出4类沉积体系，即河流三角洲、陆架浅海、陆坡浊积扇和深海盆地沉积体系；明确了两支来自婆罗洲地区的大型沉积物源；恢复了陆架三角洲-陆坡浊积体-深海盆地沉积体系以及火成岩体的空间展布，为进一步开展北康盆地油气调查提供基础支撑。Abstract: Traditional seismic facies analysis is usually carried out based on the features of topography (shelf, continental margin and slope), external morphology (dome-like, wedge-like, plate-like）、internal structure（parallel, oblique and radiate）and reflection waves (amplitude, frequency, and continuity). Lithofacies may be defined if there are enough drilling wells. However, in the areas only having few wells and complicated geological conditions, it is difficult to reach an appropriate conclusion of lithofacies without additional data support. In this paper, we will introduce three new marks to help making interpretation of seismic facies into lithofacies, namely, the smoothness of seismic reflection waves, the cleanliness of seismic reflection waves and special wave forms. Thus a new method for seismic interpretation from seismic facies to lithofacies to sedimentary facies is proposed, which is suitable for the reconstruction of sedimentary systems under any geological and data conditions and can effectively enhance the quality of lithofacies identification and sedimentary facies conversion, and it is of great significance to oil and gas exploration. Combining this method with the method for traditional seismic facies analysis, six typical seismic lithofacies are identified for the Late Oligocene of the Beikang Basin, namely, sandy conglomerate facies, sand-mud-encrusted facies, interbedded sand-mud facies, mud-encrusted facies, pure mud-shale facies and volcanic facies. Four kinds of depositional systems are recognized, namely, the fluvial delta, shelf shallow sea, slope turbidite fan and deep-sea basin depositional systems; Two large sediment sources from the Borneo region are recognized; The spatial distributions of shelf delta-continental slope turbidite, deep-sea basin depositional system and igneous body are clarified. The results have provided basic support for further exploration of oil and gas in the Beikang Basin.
-
Keywords:
- seismic facies /
- lithofacies /
- sedimentary facies /
- sedimentary system /
- southern South China Sea /
- Beikang Basin
-
天然气水合物是由水与天然气在高压低温条件下形成的结晶物质[1]，直观上与冰相似，广泛分布于海洋深水地层和高纬度极地冻土地层中[2]。目前普遍认为水合物中所储存的能量是所有传统化石能源的2倍[3]。同时由于其能量密度大、污染程度低等因素[4]，被认为是石油等传统化石能源的接替能源。目前水合物开发方向的研究主要包含3个方面：水合物成藏机理，水合物勘探开采以及相关地质问题研究。人们进行大量的水合物开采室内实验和数值模拟研究之后(例如Kurihara采用MH21-HYDRES[5]，White和McGrail采用STOMP-HYD[6]，Moridis采用TOUGH + HYDRATE[7, 8]等)，水合物分解所引起的地质问题开始逐渐得到各界的重视。研究认为，水合物开采分解所致的海底地质风险不仅会作用在开采范围，同时会诱发长期的巨大的地质灾害及环境风险，海底滑坡是其典型表现形式[9, 10]。目前关于水合物分解对海底边坡稳定性影响的科学研究尚处于理论探索阶段。
水合物分解会严重降低沉积层的强度[11, 12]，具体表现在以下2点：(1)水合物分解后其自身在水合物沉积层中的胶结、支撑作用的消失，其分布区域的力学强度大大降低；(2)注热开采或地层流动性较差的情况下，水合物分解产生的大量气体不能及时从地层中排出，会在孔隙内部形成异常高压，进而破坏储层结构，降低储层强度及承载能力。两方面作用使得储集层强度严重降低，进而引起地质失稳问题[13]。Sultana研究认为已知规模最大的海底滑坡Storgga滑坡的主要诱因是水合物分解，同时采用Comsol软件进行的水合物分解多场耦合的数值分析，进一步证实了水合物分解会引起大规模的海底滑坡[14]。Gupta等采用C++编程建立了系统的耦合热、化学、静水压力、岩石力学行为的水合物分解模型，论证了岩石力学行为在水合物分解过程中的重要性，并且结合数值模拟和实验得出结论：水合物分解会导致海底地层出现大于5%的体积应变[15]。Li等采用Flac3D建立了含水合物的海底边坡模型，论证水合物分解对海底边坡失稳的影响的同时分析了水合物分解时边坡各个位置向下运动的位移[16]。周丹分析认为海底滑坡规模远大于陆地滑坡且当边坡角度小于10°时，由于内部水合物层强度降低，仍有可能发生海底滑坡[17]。杨晓云等进一步探讨了水合物分解对海底沉积层岩石力学参数的影响机理及影响程度，并认为胶结强度的降低是海底边坡失稳的主导因素[18]。马云[19]、于桂林[20]、张振飞[21]分别采用ABAQUS软件建立了一系列含水合物边坡的数值模型，并评估了边坡稳定性和失稳范围及潜在的影响因素。
上述研究存在以下几方面的不足之处：(1)没有考虑水合物实际开采过程，采取单一水合物饱和度进行分析。水合物分解开采过程中，水合物自钻孔向两侧逐渐分解，即沿井眼径向上水合物饱和度不均匀；(2)本构模型及材料参数论证不足。不同水合物饱和度下沉积层强度参数直接影响整个边坡的力学行为，系统详实本构方程描述在相关分析中至关重要；(3)人工迭代试算的计算方法。上述诸多分析均采用人工计算的方式，使得所得结果不连续、非动态，且计算过程繁琐。本文在前人研究的基础上，采用基于ABAQUS的强度折减法，结合实际开采情况，选取连续的水合物沉积层本构模型，实现自动迭代折减计算，分析了开采过程中不同边坡角度、埋藏深度条件下水合物分解程度对边坡稳定性的影响情况。研究结果表明，水合物分解是边坡失稳的主导因素，在20°边坡条件下，井眼附近水合物地层强度低于沉积层强度时，会引起较大规模的海底滑坡；边坡角度和水合物分解程度共同作用下边坡会急剧失稳。同时，水合物带的位置及分解程度决定了边坡的稳定性及滑坡的规模。
1. 模型建立
1.1 水合物本构模型
水合物沉积层的本构模型建立是水合物分析中的重点难点部分。其具体表现为水合物对沉积层特有的3个作用:一是水合物对沉积物的胶结作用；二是水合物分解对沉积层结构的破坏；三是由于水合物强度高于单一沉积砂泥岩强度，使得含水合物的沉积层强度(在水合物未分解或轻度分解条件下)普遍高于不含水合物层[13, 22]。大量分析认为，影响水合物在地层中岩石力学行为特征的主要因素同样可以概括为3类：岩土粒径分布和岩性分布[23, 24]；水合物成藏原因[25-27]；水合物饱和度分布[28]。因此，本构模型的建立势必要引入水合物饱和度[29]。目前人们对水合物本构模型进行了大量的三轴实验并得到了一系列具有一定可信度的水合物沉积层本构模型[30]。其中应用较多的分别是基于摩尔—库伦模型和剑桥模型改进建立的水合物沉积层本构模型[31, 32]。综合分析各种已有模型，本文所选取的本构模型是基于非线性理论和摩尔—库伦模型所建立的，以水合物饱和度为基础变量，建模时离散沉积层和水合物层，并最终得到综合的抗剪强度参数及杨氏模量。该模型不仅可以较好地描述不同区块动态水合物饱和度下的沉积层应力应变关系[33]，同样与Miyazaki等进行的大量关于水合物的三轴实验结果相吻合，因此, 以该模型作为水合物沉积层力学特征的本构模型[34]。关系式如下所示[32]：
$$ c = {c_s} + {c_h} = {c_1}\left( {1 - {e^{ - \frac{{\sigma _3^\prime }}{{\sigma _t^\prime }}}}} \right) + {c_2}s_h^{{c_s}} $$ $$ \begin{array}{l}{\sin \varphi=\sin \varphi_{s}+\sin \varphi_{h}=} \\ {\sin \left[\delta_{1}-\delta_{2} \ln \left(\frac{\sigma^{\prime}_{3}}{1_{M P a}}\right)\right]+\delta_{3} s_{h}\left(1-e^{-\frac{\gamma}{\delta_{4}}}\right)}\end{array} $$ $$E=E_{s}+E_{h}=e_{1}\left(\frac{\sigma_{3}^{\prime}}{1_{M P a}}\right)^{\theta_{2}}+e_{3} s_{h} $$ 式中，c、cs、ch分别是水合物沉积层、砂和水合物的内聚力(kPa)；φ、φs、φh分别是水合物沉积层、砂和水合物的摩擦角(°)；E、Es、Eh分别是水合物沉积层、砂和水合物的弹性模量(kPa); sh是水合物饱和度。其他参数在具体某一区块为常数，σ′3为有效围压(kPa)；γ是内部剪切应变；δ1和δ2是与相对密度和固体材料相关的常数，δ3是与水合物的存在而引起的膨胀增加相关的常数，δ4是与脱粘/断裂相关的剪切应变相关的常数。c1、c2分别为所测沉积层与水合物的粘聚力(kPa)；c3为拟合所得的指数参数；e1、e2、e3是Pinkert等建立模型时拟合得出的常数。确定岩土强度的方法包括室内实验和现场原位测试两大类方法[2]，结合我国南海神狐海域水合物沉积层骨架强度，具体参数如表 1所示。
物理量 数值 单位 c1 293 kPa c2 1960 kPa c3 1.7 σ′3 1962 kPa σ′t 382 kPa δ1 0.4 δ2 0.2 δ3 0.22 δ4 0.04 e 2.72 e1 548000 kPa e2 265800 kPa e3 450000 kPa γ 1.2 基于ABAQUS的有限元强度折减法
传统的边坡稳定性分析方法主要包括极限分析、极限平衡和以有限元为主的数值计算方法[36]。随着计算机技术和有限元的发展，基于应力场的稳定分析方法和有限元强度折减法越来越得到普遍认可，并逐渐成为未来边坡稳定性分析的发展方向。陆上边坡稳定性分析方法同样适用于海底边坡稳定性分析[37, 38]。本文即采取有限元强度折减法作为水合物分解导致的边坡稳定性分析的指导方法。强度折减法由Zienkiewicz于1975年提出[39]，该方法的核心是抗剪强度折减系数，其定义为：在外载荷保持不变的情况下，边坡内土体所提供的最大抗剪强度与外载荷在边坡内产生的实际剪应力之比。临界破坏时的强度折减系数Fr即为边坡的安全系数Fs。
折减后的抗剪强度参数可分别表达为[40]：
$$ \begin{aligned} c_{m} &=c / F_{r} \\ \varphi_{m} &=\arctan \left(\tan \varphi / F_{r}\right) \end{aligned} $$ 式中：cm、φm是维持平衡所需要的或土体实际发挥的内聚力和摩擦角；Fr是强度折减系数。
计算中需要假定不同的强度折减系数Fr，使用折减后的强度参数进行有限元分析。在ABAQUS中，将折减系数定义为场变量，进而实现折减系数在计算过程中自动增加，并最终得到计算中止时的折减系数。在此基础上，根据下部分所选取的边坡失稳准则，确定边坡失稳的安全系数，并根据所得的计算结果进行分析。
1.3 边坡失稳准则
边坡失稳的判据同样直接影响对边坡稳定性的认识情况及强度折减法的准确性。目前判断土坡达到临界破坏的评价标准主要有如下3种[38, 40]：①形成连续的塑性贯通区；②边坡顶部拐点发生显著位移；③计算因不收敛而终止。
在实际计算过程中，塑性区一般在边坡底部开始形成，并逐渐向上发展。当塑性区最终发展形成连续的贯通区时，会导致边坡顶部拐点位置产生位移。随着岩土强度进一步折减，边坡顶部特征部位会产生一个严重的拐点，即说明此时边坡发生失稳，产生了较大变形。若再进一步折减岩土强度，就会引起计算的不收敛。即这3种判断依据是存在一个时间上的联系的。同时，根据数值模拟结果，发现这3种判据所确定的安全系数是逐渐增大的，但自塑性贯通区形成至计算中断，此过程折减系数变化范围不超过其数值的5%。
在实际分析过程中，只有以统一评价标准进行对比分析才可以得到可信的结果。而前2种评价标准只能根据计算过程中的云图分布或顶部节点位移变化曲线，人为地选取一折减系数作为安全系数，由于缺乏定量标准，会给不同情况下的对比带来主观误差。第3种判断方式具有统一的评价标准，如只需固定两迭代步之间的最小增量，即可形成统一评价标准。因此，本文以计算终止为边坡失稳的评价标准，以利于在同一标准下评价不同情况下的边坡稳定性问题。
2. 水合物边坡模型
2.1 数值模型建立
本模型模拟由于开发过程中水合物分解引起的水合物带强度特征变化对海底边坡稳定性的影响。随着时间推移，以井眼为中心，水合物带分解程度及分解范围逐渐扩大。通过模拟计算水合物特定分解范围、程度、时间点下边坡的安全系数，评估由水合物分解对边坡稳定性的影响。通过数值模拟塑性区分布云图及位移云图，确定水合物分解引起的边坡失稳规模及程度的变化。
对于理想的三维均质地层，通过基础弹塑性力学理论，可以简化为二维平面模型。对基本模型进行网格划分。为了保证计算精度，网格结构选择结构化的四边形单元，网格属性确定为平面应变网格。模型的基本尺寸如图 1所示。模型水平跨度1000m，总厚度400m，水合物埋深120m，边坡角度20°，坡顶水平跨度200m，1区域对称轴位置为井眼位置，将水合物带以井眼为中心，向两侧划分区域，以模拟不同的分解程度和分解范围[41]。本模型为基本模型，后续模型在本模型基础上进行适当修改。
2.2 边界条件和初始状态
结合典型边坡失稳分析，限制边坡左右边界在水平方向上的位移，限制边坡下部边界在各个方向的位移。同时，结合我国南海神狐海域水合物实际勘探情况，选取合适条件，确定数值模拟所需的基本参数。实际勘探数据与模型采用数据如表 2所示。
Table 2. Exploration data and numerical model selected data for hydrate formations in the Shenhu area参数 神狐海域勘探数据 模型选用数据 水深/m 500~2000 1300~1650 水合物埋深/m 7~350 50~150 海底坡脚/(°) 最高达30 10~25 颗粒密度/(g/cm3) 2.695~2.716 2.700 海水密度/(g/cm3) 1.040 1.040 水合物密度/(g/cm3) 0.980 0.980 水合物饱和度 0.45~0.92 0.8 沉积物孔隙度 0.27~0.63 0.4 整个模型只受来自海水的上覆压力及自身的重力影响[42-44]。采用地应力平衡分析对施加初始载荷的数值模型进行地应力平衡，在地应力平衡状态的基础上进行数值分析模拟。对于不同的水合物分解程度，可以计算出水合物残余饱和度，然后将其代入本构模型，可计算该条件下的强度参数。根据设定方案，计算所得各分解程度下材料强度如表 3所示。
表 3 模型材料强度参数Table 3. Strength parameters of numerical model材料 内聚力
/kPa摩擦角
/(°)杨氏模量
/MPa密度
/(g/cm3)沉积层 291.2 18.6 56.3 2.200 水合物层分解程度 未分解 1562.9 26.0 206.0 1.844 分解30% 504.2 22.7 89.7 1.794 分解60% 308.6 19.4 41.2 1.745 分解80% 89.2 17.3 23.7 1.711 分解100% 35.4 15.2 10.28 1.688 3. 模型计算分析
3.1 水合物分解程度对边坡稳定性影响
在如图 1所示模型上，设计7种不同的计算方案来分析水合物分解范围及程度逐渐扩大的过程。方案设计基本依据为：钻孔引起的水合物分解，随时间推移，水合物分解范围及分解程度不断变大。方案简介如下：作为对比方案，首先对不含天然气水合物层的海底边坡进行稳定性分析；之后对含水合物但未分解的海底边坡进行稳定性分析；随后将1区域水合物分解程度设为30%，其他区域未分解；将1区域分解程度设为60%，2区域分解程度设为30%，其他区域未分解；依次类推，具体计算方案如表 4所示。
表 4 模拟方案设计Table 4. Conceptual design of simulation方案 方案描述 方案1 不含水合物层的均质海底边坡模型 方案2 含水合物且水合物未分解的海底边坡模型 方案3 1区分解30%，其他区域未分解 方案4 1区分解60%，2区分解30%，其他区域未分解 方案5 1区分解80%，2区分解60%，3区分解30%，4区未分解 方案6 1区分解100%，2区分解80%，3区分解60%，4区分解30% 方案7 1、2区分解100%，3、4区分解80% 部分方案数值模拟如图 2所示。前3种方案下塑性区域发展情况基本相同：塑性区域首先出现在坡脚，逐渐发展至坡顶形成塑形贯通区，塑形贯通区及其所包含的区域是滑坡发生的范围。在ABAQUS中可以通过塑性区的位置以及边坡在水平方向上的位移来判断边坡失稳的范围与规模。对比方案1、2、3，说明水合物带的存在不会降低边坡的稳定程度，由于方案3下水合物带的强度是高于沉积层的强度的，因此，可以说明局部的强度增加对边坡整体的稳定性影响不显著。同样，对比方案2、3，可以认为局部强度的小幅度降低对边坡整体的稳定性影响不显著。随着水合物分解程度的增加，自方案4开始出现贯穿水合物带的塑形区，如图 3所示。
![]()
图 2 部分方案计算终止时塑性区分布云图Figure 2. Contours of plastic zone when ending calculation in several cases![]()
图 3 方案4下不同折减系数塑性区发展过程Figure 3. The development progress of plastic zone in reduction factors of case 4对比水合物带分解60%之后的水合物带强度与边坡均质土层强度，发现两者弹性模量、泊松比等强度参数特征相似，可以认为60%分解程度下的水合物带强度略低于边坡均质土层强度，两者强度关系对于分析塑性区发展过程十分重要。在此强度下的方案4计算过程中塑性区发展过程与之前所述方案明显不同。该方案下，如图 3a所示，塑性区首先在坡脚位置出现，紧接着在模型右侧和水合物带出现，而不是自坡脚逐渐向上发展。随着分析的进行，坡脚的塑性区逐渐向上发展至水合物带，水合物带的塑性区同时向上发展至坡顶，形成贯通的塑性区，如图 3b。此时塑形贯通区的形成并没有直接导致计算中止，折减系数进一步增大时，一个新的塑性带自坡脚发展至坡顶，并最终形成另一个较弱的塑性贯通区，即形成了2个塑性贯通区。如图 3c所示。随着折减系数继续增加，2个塑性区大部分合二为一，但仍存在小部分的间隔区域。
通过位移云图可以更加直观地认识边坡最终状态和具体失稳程度，图 3c与图 4c为同一折减系数，同为计算终止时的状态，可以发现2种云图的对应关系。同时，列举了几个方案下的位移云图，与图 2对比分析发现，水合物分解程度越高，滑坡区域的位移越大，即失稳程度越高。
![]()
图 4 部分方案下计算终止时位移云图Figure 4. Contours of displacement in horizon direction when ending calculation in several cases最后3种方案下，塑性区均是首先出现在水合物带，然后出现在坡脚位置。水合物带塑性区向上发展至坡顶，坡脚位置塑性区发展至水合物带，形成塑性贯通区。由于此3种方案下，水合物带强度已显著低于沉积层强度，因此，没有出现方案4相似的二次滑坡问题。前3种方案下，边坡失稳规模取决于其本身强度，后4种方案下，边坡失稳规模取决于水合物带所处的位置。因此，水合物分解直接影响边坡的稳定性及失稳规模。
各方案下安全系数如图 5所示。由于塑性区域总是发生在边坡强度最低或应力分布最集中的位置，前3个方案下均质沉积层的强度低于水合物带的强度，因此，安全系数保持不变。自方案4开始，水合物带分解后的强度低于沉积层的强度，安全系数取决于水合物带的强度。安全系数随水合物带强度的不断降低而降低，并最终小于1，这意味着随着水合物的分解，海底边坡会发生大规模滑坡。对比方案6及方案7的安全系数可以发现，边坡的安全系数不仅受边坡强度影响，还与边坡低强度带的规模相关，即局部小规模水合物带的高度分解对整个边坡的安全系数影响程度显著。与马云[19]、于桂林[20]、张振飞[21]等模拟所得结果进行对比，虽然各研究中地层强度参数、边坡模型尺寸、本构模型以及失稳准则等均有差异，但水合物分解前后所得安全系数不存在数量级之间的差异，且随水合物分解，安全系数的变化呈相同趋势。
3.2 不同上覆厚度下水合物带分解对边坡稳定性影响
在图 1所示的基本模型的基础上，只改变图 4模型中水合物带的竖直方向的位置，调整水合物带上覆沉积层的厚度，建立3个新模型，同时将基本模型纳入分析，每个模型根据水合物分解程度分别进行2次计算。具体计算方案如表 5所示。
表 5 水合物层不同上覆厚度下计算方案Table 5. Numerical procedure under different overlay thickness of hydrate计算模型 上覆沉积层厚度/m 水合物分解程度 模型1 50 未分解 1区域分解100%，2区域分解80%，3区域分解60%，4区域分解30% 模型2 80 模型3 120 模型4 150 当水合物带未分解时，水合物带相当于一个坚硬的夹层，边坡失稳的程度、规模以及安全系数取决于边坡自身的强度。但是其埋藏程度对边坡稳定性仍有一定影响。图 6反映了水合物未分解时的边坡失稳特征。当水合物带埋深低于边坡自身可能失稳时的滑坡深度时，具有较高强度的水合物带会限制边坡失稳的规模，但是对安全系数几乎没有影响。当水合物带埋深较深时，水合物带的存在对边坡自身稳定性几乎没有影响，但对边坡抵抗地质灾害的能力有所提升。
![]()
图 6 不同上覆厚度下水合物带未分解计算终止时位移云图Figure 6. Displacement contours of slope without hydrate dissociation under different overlay thickness由图 7可以看出，随着水合物带的分解，塑性区域总是会贯通整个水合物带，水合物带上部区域都是潜在的滑坡区域。边坡产生横向位移区域的下边界与塑性贯通区下边界是同一个边界。因此，横向位移云图不仅能反映边坡失稳的区域和规模，还能反映失稳的强度。对比各方案发现，塑性区总是可以贯通水合物带和坡脚坡顶，上覆沉积层厚度的不同直接导致滑坡规模的不同。随着上覆沉积层厚度的增加，边坡失稳区域随之增加。这种增加不仅仅表现为厚度的增加，坡顶失稳的范围也逐渐向内部发展。这是由于分解后的水合物带强度低于沉积土层所致。
![]()
图 7 不同上覆厚度下水合物带分解计算终止时横向位移云图Figure 7. Horizontal displacement contours of slope when ending calculation with different overlay thickness安全系数与水合物带埋深关系曲线如图 8所示，水合物未分解情况下，边坡的安全系数保持稳定；水合物带分解至低于沉积层强度时，随着上覆岩层厚度的增加，安全系数增加，且失稳的规模也增加。可能的原因是：水合物带分解后强度大幅降低后，塑形区总是会贯穿整个水合物层。上覆岩层厚度的增加，塑性贯通区以上所包含的均质沉积层的规模增大，所需能量越多，抵抗由于下部水合物带强度降低所引起的失稳的能力越大，导致安全系数逐渐上升。其具体原因有待进一步分析论证。根据图 8所示的两条曲线发展趋势，基本可以确定不同埋深条件下，水合物分解与未分解情况所得安全系数会存在交点。此时，水合物分解与否不会对边坡的稳定性产生影响，对边坡稳定性产生影响的只能是沉积层自身强度。同样，此规律尚无成型的数学预测模型，有待于进一步研究。
3.3 水合物带分解对不同角度边坡稳定性影响
在如图 1所示的基本模型的基础上，改变图 3-1模型中边坡角度，对应调整边坡左侧高度。本小节中共采用4个模型，其边坡角度分别是10°、15°、20°、25°，水合物带保持与边坡斜边平行，每个模型根据水合物分解程度分别进行2次计算。具体计算方案如表 6所示。
表 6 不同边坡角度的模型计算方案Table 6. Numerical procedure with different slope angel计算模型 边坡角度/(°) 水合物分解程度 模型5 10 未分解 1区域分解100%，2区域分解80%，3区域分解60%，4区域分解30% 模型6 15 模型7 20 模型8 25 水合物分解前后模拟结果如图 9、图 10所示，其中，边坡角度为20°时，水合物分解前、后云图分别见图 6c、图 7c。可以发现，边坡角度对边坡安全系数影响很大，但是对边坡可能的失稳规模影响不大。因为塑形区是连接整个水合物带和坡脚、坡顶的，在水合物带相对位置相同的情况下，塑形贯通区以上的失稳区域也大致相当。图 11所示为水合物分解前后不同角度下边坡的安全系数，随着边坡角度的增加，边坡的安全系数本身会有一个较大幅度的变化，但总是大于1(小于1意味着自然失稳)。正是由于水合物的分解作用，分解后水合物安全系数降至1以下，说明水合物分解对边坡稳定性影响显著, 同时说明在较大角度的海底斜坡上开采水合物风险很大，应充分评估其失稳可能性。
![]()
图 9 不同上覆厚度下水合物带未分解计算终止时位移云图Figure 9. Displacement contours of slope under different slope angle without hydrate dissociation![]()
图 10 不同边坡角度下水合物带未分解计算终止时横向位移分布云图Figure 10. Horizontal displacement contours when ending calculation under different slope angle4. 结论
(1) 塑性区域发展过程。正常情况下，塑性区自坡脚区域向上发展，形成塑形贯通区；水合物分解程度达到60%后，塑性区会先后出现在水合物带和坡脚区域，水合物带的塑性区向上发展至坡顶，坡脚区域的塑性区向上发展至水合物带，最终形成塑形贯通区。
(2) 水合物分解程度达到60%及以上，整个边坡的安全系数会降低，并最终失稳。分解程度越高，边坡失稳时的安全系数越小。
(3) 水合物未分解情况下，水合物带的存在有利于维持边坡稳定；水合物充分分解条件下，水合物埋深越深，整个边坡的安全系数越高。大角度边坡水合物开采风险较大。
-
![]()
图 2 北康盆地层序划分与构造演化综合柱状图
Figure 2. Integrated stratigraphic column and tectonic evolution of the Beikang Basin
![]()
图 6 钻井岩性与井旁地震反射特征
Figure 6. Well lithology and borehole seismic reflection characteristics
![]()
图 7 基于传统地震相分析和3种新增标识的综合地震相分析
Figure 7. Comprehensive seismic facies analysis based on traditional seismic facies analysis and the method using three new markers
![]()
图 8 北康盆地8种主要地震相与地震岩相特征
Figure 8. Characteristics of eight major seismic facies and seismic lithofacies in Beikang Basin
-
[1] 谢晓军, 赵志刚, 张功成, 等. 南海南部三大盆地油气地质条件差异性[J]. 地球科学, 2018, 43(3):802-811 XIE Xiaojun, ZHAO Zhigang, ZHANG Gongcheng, et al. Hydrocarbon geological differences of three basins in southern South China Sea [J]. Earth Science, 2018, 43(3): 802-811.
[2] 王龙, 谢晓军, 刘世翔, 等. 南海南部主要盆地油气分布规律及主控因素[J]. 天然气地球科学, 2017, 28(10):1546-1554 WANG Long, XIE Xiaojun, LIU Shixiang, et al. Analysis of hydrocarbon accumulation and diversity of the major basins in mid-southern part of the South China Sea [J]. Natural Gas Geoscience, 2017, 28(10): 1546-1554.
[3] 赵忠泉, 钟广见, 冯常茂, 等. 南海北部西沙海槽盆地新生代层序地层及地震相[J]. 海洋地质与第四纪地质, 2016, 36(1):15-26 ZHAO Zhongquan, ZHONG Guangjian, FENG Changmao, et al. Cenozoic sequence stratigraphy and seismic facies analysis of Xisha trough basin in northern South China Sea [J]. Marine Geology & Quaternary Geology, 2016, 36(1): 15-26.
[4] 吴时国, 袁圣强. 世界深水油气勘探进展与我国南海深水油气前景[J]. 天然气地球科学, 2005, 16(6):693-699,714 doi: 10.3969/j.issn.1672-1926.2005.06.002 WU Shiguo, YUAN Shengqiang. Advance of exploration and petroleum geological features of deep-water hydrocarbon in the world [J]. Natural Gas Geoscience, 2005, 16(6): 693-699,714. doi: 10.3969/j.issn.1672-1926.2005.06.002
[5] Vail P R, Mitchum R M Jr, Thompson S III. Seismic stratigraphy and global changes of sea level: Part 3. Relative changes of sea level from coastal Onlap: Section 2. Application of seismic reflection configuration to stratigraphic interpretation[M]//Seismic Stratigraphy-Applications to Hydrocarbon Exploration. AAPG Memoir, 1977, 26: 63-83.
[6] Mitchum R M Jr, Vail P R, Sangree J B. Seismic stratigraphy and global changes of sea level: Part 6. Stratigraphic interpretation of seismic reflection patterns in depositional sequences: Section 2. Application of seismic reflection configuration to stratigraphic interpretation[M]//Seismic Stratigraphy-Applications to Hydrocarbon Exploration. AAPG Memoir, 1977, 26: 117-135.
[7] Badley M E. Practical Seismic Interpretation[M]. Boston: IHRDC Press, 1985: 1-259.
[8] Sheriff R E. Aspects of seismic resolution: Chapter 1[M]//Seismic Stratigraphy II: An Integrated Approach to Hydrocarbon Exploration. AAPG Memoir, 1985, 39: 1-10.
[9] 鲜本忠, 牛花朋, 朱筱敏, 等. 准噶尔盆地西北缘下二叠统火山岩岩性、岩相及其与储层的关系[J]. 高校地质学报, 2013, 19(1):46-55 doi: 10.3969/j.issn.1006-7493.2013.01.008 XIAN Benzhong, NIU Huapeng, ZHU Xiaomin, et al. Early permian volcanic lithology, lithofacies and their relations to reservoir in northwestern margin of the Junggar basin [J]. Geological Journal of China Universities, 2013, 19(1): 46-55. doi: 10.3969/j.issn.1006-7493.2013.01.008
[10] 罗伟平, 李洪奇, 朱丽萍, 等. 地震与测井资料自动匹配的研究[J]. 石油地球物理勘探, 2014, 49(1):205-212 LUO Weiping, LI Hongqi, ZHU Liping, et al. Application of the cross-hole electromagnetic method (CHEM) in hydrocarbon reservoir monitoring [J]. Oil Geophysical Prospecting, 2014, 49(1): 205-212.
[11] 庞雄, 施和生, 朱明, 等. 再论白云深水区油气勘探前景[J]. 中国海上油气, 2014, 26(3):23-29 PANG Xiong, SHI Hesheng, ZHU Ming, et al. A further discussion on the hydrocarbon exploration potential in Baiyun deep water area [J]. China Offshore Oil and Gas, 2014, 26(3): 23-29.
[12] 张振波, 李东方, 轩义华, 等. 白云凹陷深水复杂构造区斜缆地震资料处理关键技术及应用[J]. 石油物探, 2014, 53(6):657-664 doi: 10.3969/j.issn.1000-1441.2014.06.005 ZHANG Zhenbo, LI Dongfang, XUAN Yihua, et al. Variable-depth streamer seismic data processing in deepwater complex structure area of Baiyun Sag [J]. Geophysical Prospecting for Petroleum, 2014, 53(6): 657-664. doi: 10.3969/j.issn.1000-1441.2014.06.005
[13] 郑胜, 吕成福, 陈国俊, 等. 南海西北部晚中新世红河海底扇储集层特征[J]. 中南大学学报:自然科学版, 2015, 46(5):1754-1762 ZHENG Sheng, LÜ Chengfu, CHEN Guojun, et al. Reservoir characteristics of Late Miocene Red River submarine fan, northwestern South China Sea [J]. Journal of Central South University:Science and Technology, 2015, 46(5): 1754-1762.
[14] 张文彪, 段太忠, 刘志强, 等. 深水浊积水道沉积构型模式及沉积演化: 以西非M油田为例[J]. 地球科学, 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 [J]. Earth Science, 2017, 42(2): 273-285.
[15] 徐怀大, 王世飞, 陈开远. 地震地层学解释基础[M]. 武汉: 中国地质大学出版社, 1990. XU Huaida, WANG Shifei, CHEN Kaiyuan. Fundamentals of Seismic Stratigraphic Interpretation[M]. Wuhan: China University of Geosciences Press, 1990.
[16] 刘彩燕, 潘树新, 梁苏娟. 松辽盆地西部地区地震沉积相研究[J]. 地球物理学进展, 2017, 32(5):2044-2050 doi: 10.6038/pg20170525 LIU Caiyan, PAN Shuxin, LIANG Sujuan. Seismic-sedimentologic facies of western area of Songliao basin [J]. Progress in Geophysics, 2017, 32(5): 2044-2050. doi: 10.6038/pg20170525
[17] 石巨业, 樊太亮, 周娇, 等. 南图尔盖盆地几种典型地震相分析及其地质意义[J]. 科学技术与工程, 2015, 15(34):133-138 doi: 10.3969/j.issn.1671-1815.2015.34.023 SHI Juye, FAN Tailiang, ZHOU Jiao, et al. Several typical seismic facies in South Turgay Basin and the geological meaning [J]. Science Technology and Engineering, 2015, 15(34): 133-138. doi: 10.3969/j.issn.1671-1815.2015.34.023
[18] 孔令辉, 凌涛, 叶青, 等. 地震相分析在沉积相研究中的应用[J]. 复杂油气藏, 2019, 12(2):36-40 KONG Linghui, LING Tao, YE Qing, et al. Application of seismic facies analysis in sedimentary facies research [J]. Complex Hydrocarbon Reservoirs, 2019, 12(2): 36-40.
[19] Berg O R. Seismic detection and evaluation of delta and turbidite sequences: their application to the exploration for the subtle trap[M]//The Deliberate Search for the Subtle Trap. AAPG Memoir, 1982, 32: 57-75.
[20] Rankey E C. Seismic architecture and seismic geomorphology of heterozoan carbonates: Eocene-Oligocene, Browse Basin, Northwest Shelf, Australia [J]. Marine and Petroleum Geology, 2017, 82: 424-443. doi: 10.1016/j.marpetgeo.2017.02.011
[21] Song C Y, Liu Z N, Wang Y J, et al. Multi-waveform classification for seismic facies analysis [J]. Computers & Geosciences, 2017, 101: 1-9.
[22] Posamentier H W. Seismic geomorphology and depositional systems of deep water environments; observations from offshore Nigeria, Gulf of Mexico, and Indonesia[C]//64th EAGE Conference & Exhibition. Florence, Italy: European Association of Geoscientists & Engineers, 2002.
[23] Zeng H L. Seismic geomorphology-based facies classification [J]. The Leading Edge, 2004, 23(7): 644-645, 688. doi: 10.1190/1.1776732
[24] 蔡全升, 胡明毅, 胡忠贵, 等. 松辽盆地徐家围子北部宋站地区沙河子组地震相与沉积相解释应用[J]. 西安石油大学学报:自然科学版, 2017, 32(4):1-10 CAI Quansheng, HU Mingyi, HU Zhonggui, et al. Seismic facies and sedimentary interpretation of Shahezi formation in Songzhan region, northern Xujiaweizi fault depresion, Songliao Basin [J]. Journal of Xi’an Shiyou University:Natural Science Edition, 2017, 32(4): 1-10.
[25] 尹青, 万朝大, 刘伟君, 等. 地震相分析及其在石油勘探中的应用[J]. 地质找矿论丛, 2011, 26(1):79-84 YIN Qing, WAN Chaoda, LIU Weijun, et al. Seismic phase analysis and the application to oil prospecting [J]. Contributions to Geology and Mineral Resources Research, 2011, 26(1): 79-84.
[26] 邱铁成, 王征, 纪中云, 等. 影响沉积盆地地震相分析结果的主要因素研究[J]. 地球物理学进展, 2014, 29(2):831-838 doi: 10.6038/pg20140248 QIU Tiecheng, WANG Zheng, JI Zhongyun, et al. The research on main influence factors of seismic facies analysis in sedimentary basin [J]. Progress in Geophysics, 2014, 29(2): 831-838. doi: 10.6038/pg20140248
[27] 高阳, 王春贤, 冯西会, 等. 地震相分析技术在煤田地震勘探中的应用[J]. 煤田地质与勘探, 2016, 44(1):107-111 doi: 10.3969/j.issn.1001-1986.2016.01.021 GAO Yang, WANG Chunxian, FENG Xihui, et al. Application of seismic facies analysis technology in coal seismic exploration [J]. Coal Geology & Exploration, 2016, 44(1): 107-111. doi: 10.3969/j.issn.1001-1986.2016.01.021
[28] 朱剑兵, 赵培坤. 国外地震相划分技术研究新进展[J]. 勘探地球物理进展, 2009, 32(3):167-171 ZHU Jianbing, ZHAO Peikun. Advances in seismic facies classification technology abroad [J]. Progress in Exploration Geophysics, 2009, 32(3): 167-171.
[29] 姚伯初, 刘振湖. 南沙海域沉积盆地及油气资源分布[J]. 中国海上油气, 2006, 18(3):150-160 doi: 10.3969/j.issn.1673-1506.2006.03.002 YAO Bochu, LIU Zhenhu. Sedimentary basins and petroleum resources in Nansha offshore area, South China Sea [J]. China Offshore Oil and Gas, 2006, 18(3): 150-160. doi: 10.3969/j.issn.1673-1506.2006.03.002
[30] 李三忠, 索艳慧, 刘鑫, 等. 南海的基本构造特征与成因模型: 问题与进展及论争[J]. 海洋地质与第四纪地质, 2012, 32(6):35-53 LI Sanzhong, SUO Yanhui, LIU Xin, et al. Basic structural pattern and tectonic models of the South China Sea: problems, advances and controversies [J]. Marine Geology & Quaternary Geology, 2012, 32(6): 35-53.
[31] 张功成, 屈红军, 刘世翔, 等. 边缘海构造旋回控制南海深水区油气成藏[J]. 石油学报, 2015, 36(5):533-545 doi: 10.7623/syxb201505002 ZHANG Gongcheng, QU Hongjun, LIU Shixiang, et al. Tectonic cycle of marginal sea controlled the hydrocarbon accumulation in deep-water areas of South China Sea [J]. Acta Petrolei Sinica, 2015, 36(5): 533-545. doi: 10.7623/syxb201505002
[32] 鲁宝亮, 王璞珺, 梁建设, 等. 古南海构造属性及其与特提斯和古太平洋构造域的关系[J]. 吉林大学学报:地球科学版, 2014, 44(5):1441-1450 LU Baoliang, WANG Pujun, LIANG Jianshe, et al. Structural properties of Paleo-South China Sea and their relationship with the Tethys and the Paleo-Pacific Tectonic Domain [J]. Journal of Jilin University:Earth Science Edition, 2014, 44(5): 1441-1450.
[33] 林长松, 初凤友, 高金耀, 等. 论南海新生代的构造运动[J]. 海洋学报, 2007, 29(4):87-97 LIN Changsong, CHU Fengyou, GAO Jinyao, et al. On tectonic movement in the South China Sea during the Cenozoic [J]. Acta Oceanologica Sinica, 2007, 29(4): 87-97.
[34] 钱坤, 闫义, 黄奇瑜, 等. 南海扩张过程及海陆变迁沉积记录[J]. 海洋地质前沿, 2016, 32(8):10-23 QIAN Kun, YAN Yi, HUANG Qiyu, et al. Sea floor spreading of South China Sea and its depositional records of sea and land changes [J]. Marine Geology Frontiers, 2016, 32(8): 10-23.
[35] 金庆焕, 李唐根. 南沙海域区域地质构造[J]. 海洋地质与第四纪地质, 2000, 20(1):1-8 JIN Qinghuan, LI Tanggen. Regional geologic tectonics of the Nansha sea area [J]. Marine Geology & Quaternary Geology, 2000, 20(1): 1-8.
[36] 姚永坚, 杨楚鹏, 李学杰, 等. 南海南部海域中中新世(T3界面)构造变革界面地震反射特征及构造含义[J]. 地球物理学报, 2013, 56(4):1274-1286 doi: 10.6038/cjg20130422 YAO Yongjian, YANG Chupeng, LI Xuejie, et al. The seismic reflection characteristics and tectonic significance of the tectonic revolutionary surface of mid-Miocene (T3 seismic interface) in the southern South China Sea [J]. Chinese Journal of Geophysics, 2013, 56(4): 1274-1286. doi: 10.6038/cjg20130422
[37] 张翀, 吴世敏, 丘学林. 南海南部海区前陆盆地形成与演化[J]. 海洋地质与第四纪地质, 2007, 27(1):61-70 ZHANG Chong, WU Shimin, QIU Xuelin. Formation of foreland basins in the south of the South China Sea [J]. Marine Geology & Quaternary Geology, 2007, 27(1): 61-70.
[38] 张厚和, 刘鹏, 廖宗宝, 等. 南沙海域主要盆地地质特征与油气分布[J]. 中国石油勘探, 2018, 23(1):62-70 doi: 10.3969/j.issn.1672-7703.2018.01.007 ZHANG Houhe, LIU Peng, LIAO Zongbao, et al. Geological characteristics and hydrocarbon distribution in major sedimentary basins in Nansha sea areas [J]. China Petroleum Exploration, 2018, 23(1): 62-70. doi: 10.3969/j.issn.1672-7703.2018.01.007
[39] 张强, 贺晓苏, 王彬, 等. 南海沉积盆地含油气系统分布特征及勘探潜力评价[J]. 中国海上油气, 2018, 30(1):40-49 ZHANG Qiang, HE Xiaosu, WANG Bin, et al. Petroleum system distributing characteristics and exploration potential assessment of sedimentary basins in South China Sea [J]. China Offshore Oil and Gas, 2018, 30(1): 40-49.
[40] 张功成, 李友川, 谢晓军, 等. 南海边缘海构造旋回控制深水区烃源岩有序分布[J]. 中国海上油气, 2016, 28(2):23-36 ZHANG Gongcheng, LI Youchuan, XIE Xiaojun, et al. Tectonic cycle of marginal sea controls the ordered distribution of source rocks of deep water areas in South China Sea [J]. China Offshore Oil and Gas, 2016, 28(2): 23-36.
[41] 赵长煜, 宋海斌, 杨振武, 等. 南海南部边缘沉积盆地构造-热演化历史[J]. 地球物理学报, 2014, 57(5):1543-1553 doi: 10.6038/cjg20140518 ZHAO Changyu, SONG Haibin, YANG Zhenwu, et al. Tectonic and thermal evolution modeling for the marginal basins of the southern South China Sea [J]. Chinese Journal of Geophysics, 2014, 57(5): 1543-1553. doi: 10.6038/cjg20140518
[42] 施秋华, 万志峰, 夏斌. 婆罗洲地质构造特征及其对南海南部盆地的影响[J]. 海洋地质前沿, 2013, 29(1):11-16 SHI Qiuhua, WAN Zhifeng, XIA Bin. Geology of Borneo Block and its influence on basins of southern South China Sea [J]. Marine Geology Frontiers, 2013, 29(1): 11-16.
[43] 周蒂, 孙珍, 杨少坤, 等. 南沙海区曾母盆地地层系统[J]. 地球科学-中国地质大学学报, 2011, 36(5):789-797 ZHOU Di, SUN Zhen, YANG Shaokun, et al. The stratigraphic system of the Zengmu Basin, Southern South China Sea [J]. Earth Science-Journal of China University of Geosciences, 2011, 36(5): 789-797.
[44] 孙珍, 赵中贤, 周蒂, 等. 南沙海域盆地的地层系统与沉积结构[J]. 地球科学-中国地质大学学报, 2011, 36(5):798-806 SUN Zhen, ZHAO Zhongxian, ZHOU Di, et al. The stratigraphy and the sequence architecture of the Basins in Nansha Region [J]. Earth Science-Journal of China University of Geosciences, 2011, 36(5): 798-806.
[45] N. A. 安蒂斯. 砂岩油气藏的地震勘探[M]. 石油工业出版社, 1987. N. A. Anstey. Seismic exploration of sand body reservoirs[M]. Petroleum Industry, 1987.
[46] Hutchison C S. Marginal basin evolution: the southern South China Sea [J]. Marine and Petroleum Geology, 2004, 21(9): 1129-1148. doi: 10.1016/j.marpetgeo.2004.07.002
[47] Madon M, Kim C L, Wong R. The structure and stratigraphy of deepwater Sarawak, Malaysia: implications for tectonic evolution [J]. Journal of Asian Earth Sciences, 2013, 76: 312-333. doi: 10.1016/j.jseaes.2013.04.040
-
期刊类型引用(1)
1. 骆帅兵，张莉，王笑雪，余秋华，雷振宇，帅庆伟. 南海南部南安盆地新生代以来构造-沉积演化特征及区域成藏模式. 海洋地质与第四纪地质. 2023(04): 116-128 . 
本站查看
其他类型引用(1)
下载:
计量
- 文章访问数: 1983
- HTML全文浏览量: 564
- PDF下载量: 56
- 被引次数: 2
