伸展盆地构造演化的数值模拟的研究——以白云凹陷为例.pdfVIP

伸展盆地构造演化的数值模拟研究——以白云凹陷为例 廖杰L2，周蒂1，赵中贤L2，张云帆1,2，徐子英l，2 1．中国科学院边缘海地质重点实验室，中国科学院南海海洋研究所，广州510301 2．中国科学院研究生院，北京100039 摘要：在伸展盆地的数值模拟研究中，建立各种模型是基础。本文介绍了瞬时均匀拉伸模型、挠曲悬臂梁模型和多幕伸展模 型，特别强调各模型的摹本假设及适用条件。介绍了基于这些模型基础上发展出的模拟方法，如计算盆地沉降量的方法，计 算岩石圈伸展系数的传统方法以及为了弥补该方法的不足而发展出的新方法，即应用l叫剥反演和正演模拟通过模拟盆地的张 裂过程来获得岩石圈的伸展系数。在对盆地的数值模拟中，有时需要综合运用多种数值模型来突破单个模型中假设条件的约 束，如综合应用挠曲悬臂梁模型和多幕拉伸模型模拟盆地的多幕张裂。针对南海北部的白云凹陷，为了研究其裂后沉降特点， 应用文中介绍的回剥方法实测凹陷裂后沉降，应用正反演相结合的模拟方法模拟凹陷的多幕张裂过程，估测出沿整个剖面的 伸展系数并计算出理论热沉降，将实测与理论裂后沉降进行对比．结果表明，白云凹陷岩石圈的伸展系数大致呈钟形分布， 在凹陷中心处最大，大约为3．5；凹陷的实测裂后沉降远大于理论值，即存在裂后异常沉降，裂后期的异常沉降总量在凹陷 中心及北部裂后期达到2km以上。 关键词： 伸展盆地的数值模型构造演化数值模拟白云凹陷 伸展盆地是与在拉张应力作用下地壳和岩石圈伸展、减薄有关的一类张裂盆地，由陆内裂谷到被动大 陆边缘这一盆地演化序列所构成【11。伸展盆地作为沉积盆地的一种重要类型，其伸展模式及动力学机制一 直受到广泛重视。在伸展盆地的定量模拟研究中，建立各种数值模型是研究的基础。目前针对伸展盆地的 数值模拟，不同学者提出了不同的理论模型，这些模型在基本假设、应用条件及盆地演化的数学表达式上 不尽相同12】。在实际应用这些数值模型时，要分别考虑各个假设条件是否成立，有时根据实际情况需要综 合运用多种数值模型来突破单个模型中假设条件的限制。 本文介绍了瞬时均匀拉伸模型、挠曲悬臂梁模型和多幕伸展模型的假设条件和理论基础，以及基于这 些理论模型基础上发展出来的计算盆地沉降和岩石圈伸展系数的方法。盆地沉降包括应用回剥方法计算的 实测沉降和应用理论模型预测的理论沉降，回剥方法主要是指不考虑岩石圈强度的一维回剥13，4J和考虑挠 曲强度的二维同剥A61，通过除去沉积物荷载、海平面变化等非构造因素引起的沉降，得到构造沉降；理 时为了研究盆地的张裂特征及计算盆地的理论热沉降，需要知道岩石圈的伸展系数，而岩石圈伸展系数的 式计算出的理论热沉降曲线进行拟合对比，如果二者拟合较好，理论曲线对应的p值即为估测值。应用该 方法的前提条件是盆地的裂后沉降全部为热沉降，如果盆地裂后期尚存其他因素引起的沉降(即异常沉 降)，则不能应用该方法，针对这种情况本文介绍了综合应用回剥反演和正演模拟通过模拟盆地张裂过程 从而获得岩石圈伸展系数的方法，并以南海北部的白云凹陷为例，介绍了对这一方法的实际应用。 1．数值模型 1．1McKenzie瞬时均匀拉伸模型 模拟的理论基础。McKenzie模型最早将盆地演化划分成张裂期和裂后期两个阶段，并将张裂简化为在张 裂期的最后一刻瞬时发生，张裂期沉降S是由初始断层控制的沉降，取决于地壳厚度和伸展系数p；而裂 收稿FI期：修订日期： 基金项目：国家自然科学基金项目，863项目(2008AA092306)资助 作者简介：廖杰(1984一)，男，河北省衡水市人，硕士研究生，从事盆地分析研究．E-1ail：liaojie@scsio．aC．cn 488 后期的沉降S完全是由岩石圈冷却收缩引起，是岩石圈伸展系数p的函数。对于张裂及裂后沉降的计算， McKenzie【71给出了定量计算公式(式1，式2)： Si— (1) s，一E。等s；n芳【1一exp(一r／r)]cE。一二4jyll而pmaT,“，(2) 表1热沉降计算中的常量参数 上述初始沉降及热沉降的推导是基于空盆注水模犁，即盆地不含沉积物充填，在实际应用时，如果考 虑沉积物加载的盆地沉降，只需将公式中水的密度改换成沉积