一种确定沉积旋回类型的迭代算法.pdfVIP

学兔兔 石油天然气学报 (江汉石油学院学报) 2009年4月 第31卷 第2期 Journal of Oil and Gas Technology(J．JPI) Apr．2009 Vo1．31 No．2 一 种确定沉积旋回类型的迭代算法 张 涛，吕一兵，陈 忠 (长江大学信息与数学学院，湖北荆州434023) 张昌民，尹艳树 (长江大学地球科学学院，湖北荆州434023) [摘要]层序地层学中沉积旋回产生地层旋回，只有在充分了解沉积旋回的基础上才能准确解释地层旋回 样式，而沉积旋回是沉降速率B、沉积速率D和侵蚀速率E组合的结果．且旋回的形态受到8个基本点的控 制。利用一种迭代算法确定了这些基本点的位置，模拟出基准面、沉积面、沉积厚度、剩余可容纳空间、旋回 比例等参数在一个墓堆面旋回期间的变化规律，并将这些变化对应不同沉积阶段，对已经确定的层序样式 进行反演，反演的结果能够确定出在最大侵蚀速率E一定的情况下与构成这种层序的相应的B、D组合．从 而实现对地层层序模式的定量化研究。 [关键词]基准面；沉积旋回；地层旋回；迭代算法；层序样式 [中图分类号]TE121．34 [文献标识码]A [文章编号]1000—9752(2009)02—0025—05 Tipper基于基准面旋回正弦波动的假定，在考察基准面旋回中基准面与沉积面的相关关系的基础 上，提出了一种地层沉积模拟方法 (一维)DJ。其模拟的基本思路较简单，即在给定了沉降速率(B)、沉 积速率(D)和剥蚀速率(E)的情况下，在基准面波动过程中在特定点处可发生 5种地质作用：①基准面 约束下的沉积作用；②最大沉积速率约束下的沉积作用；③基准面约束下的剥蚀作用；④最大剥蚀速率 约束下的剥蚀作用；⑤沉积物的路过作用。通过改变基准面控制因素，确定地层叠加样式与各控制因素 之间的响应关系。在文献 E1]的模型中，假设地层基准面的波动完全符合正弦波特征；而构造沉降也 是固定的。对于沉积和侵蚀条件也进行了假设，即发生沉积作用时，具有最大沉积速率，当基准面位于 沉积面之上时，即使有足够的可容纳空间，沉积速率也不可能超过该速率；发生剥蚀作用时，具有最大 剥蚀速率，即使基准面远低于沉积面，也不能通过大于该速率的剥蚀速率来使沉积面下降至基准面的位 置，即剥蚀速率不能超过最大剥蚀速率。那么在一个基准面波动旋回内，基准面与沉积面之间的关系就 可以通过正弦曲线上的8个控制点 (P、c、P、q-，．、 、志、s)表征l2 (图1)。笔者通过一种直接迭代算法确定 了P、c、e、q、r、i、五、s这8个基本点的位置，进而通过这些点来确定层序旋回类型。 1 算法说明 下面根据这8个点计算出P(相滞后)、s(旋回比)、 (空效比)、H(间隙比)的值，进而通过P、s、 、H 来确定沉积旋回类型(层序类型的分类见图2)，再反演B、D、E组合。 刻画B、D、E的参数如下：A表示基准面波动幅度；w表示形成层序的时间；Loe 表示最大侵蚀速率； L 表示最大沉积物供给速率；B 表示最大沉降速率；M表示净流量。为方便模拟，将5个变量简化为3 个无量纲变量： B— B ×(W／2r：A)×100 D—Ls ×(Ⅳ／2uA)×1O0／M E—LDEF×(W／2uA)×lOO／M 式中，B为单位弧度内的沉降造成的相当于基准面波动振幅的百分比表示值，即沉降速率。 [收稿日期]2009—01—05 [基金项目]国家自然科学基金项目 (40572O78／Do206)；教育部重点实验室开放基金项目 (KLET0R06O8)；湖北省教育厅重点项 目 (D200512001)。 [作者简介]张涛 (1978一)，男，2002年大学毕业，硕士，、井师，现主要从事最优化理论与算法研究工作。 学兔兔 · 26· 石油天然气学报 (江汉石油学院学报