基于层序特征重构沉积条件的反演算法设计.pdfVIP

石油天然气学报 (江汉石油学院学报)2013年3月 第35卷 第3期 JournalofOilandGasTechnology(J．JPI)Mar．2013Vo1．35No．3 基于层序特征重构沉积条件的反演算法设计 张I潼’寿 ，长【江大学信息与数学学院，湖北荆州434023 、 水资源与水电科学国家重点实验室(武汉大学)，湖北武汉430o72』 陈忠 ，吕一兵 (长江大学信息与数学学院，湖北荆州434023) [摘要]设计了一种定量的迭代算法，对 已经确定的层序样式进行反演，该算法能够确定出在最大侵蚀速 率E一定的情况下与构成这种层序 的相应 的最大沉降速率B、最大沉积速率D组合，从而实现 了对地层层序 模 式的定量化研 究。 [关键词]基准面 ；地层旋 回；沉积旋 回 [中图分类号]TEl22 [文献标志码]A [文章编号]i000—9752(2013)03—0051一o3 层序地层学是一种划分、对 比和分析沉积岩的新方法。层序地层学已成为油气勘探中一种广泛应用 的技术 ，其主要作用就是在等时层序地层格架 中进行沉积体系的分析 ，预测生油层、储层和盖层的分布 及特征 。由于层序地层学理论的引入，对陆相含油气盆地地层及沉积研究起到极大 的推进作用。但是， 这一学科本身定性分析较多，在定量分析方面相对薄弱 ，从而造成技术操作的人为性和随意性。如何使 层序地层学的理论体系达到定量化是当前需要解决的一个重要问题。国内学者大都偏重于应用 ，对不同 地质背景下地层堆积样式及其控制因素的探索不够 ，对地层层序模式的定量化研究基本上是空 白[1]。 基于基准面旋 回为正弦波动的假定，在一个基准面旋 回内，沉积旋回显示为特定 的层序类型，而沉 积旋回是由不同的B(最大沉降速率)、D(最大沉积速率)、E(最大侵蚀速率)组合所形成 ，沉积旋回形成的 沉积地层都有一个相似的形态 ，这些形态都可以用标准形态进行描述。为此，笔者在提取了层序类型特征 的基础上 ]，通过一种反演算法确定了符合这一特征的相应的B、D、E组合 ，从而重构沉积条件 。 1 反演算法设计 l-1 算法思路 首先离散正弦基准面曲线的定义域区间，假设离散节点个数为 N(取 N 一 2000)，然后根据岩石面与 基准面的位置关系确定一个稳定 的沉积旋 回并计算 ： 1) o]，旋回的起始位置。 2) 1]，基准面位于沉积表面之上的分界点。 3) 2]，沉积面由低于基准面到上升到基准面位置 。 4) [3]，剥蚀作用结束 ，沉积作用开始的位置。 5)pr-4]，基准面与岩石面高度差最大的点 。 6) 5]，剥蚀开始的位置，即沉积间断开始的位置。 7)p[6-1，沉积面由原来高于基准面下降到与基准面相平的位置 。 8) [7]，岩石表面开始高于基准面的位置。 9) 8]，沉积面由低于基准面到上升到基准面位置。 [收稿日期]2012—12—04 [基金项目]国家自然科学基金项 目 61273179)；湖北省教育厅重点项目 (。 [作者简介]张涛 (1978一)，男，2002年大学毕业 ，讲师，博士生，现主要从事复杂系统建模与智能计算 。 石油天然气地质 2013年3月 1．2 具体算法 ． 算法 中参变量说明如下 ：Base为基准面 ；lace为岩石面；Sed为最大沉积速率 ；Erosion为最大侵蚀速 率；Sub为沉降速率 ；A为振幅。具体算法设计如下： 步骤 1 初始化 BaseEi3： Bas[]一Asin( +号) 一0，1，…，4000 如fPE4ooo~一BaseEo] Sub一 Sub Sed一 Sed hEli—BaseEi]一BaseEi一 1]+Sub hEi3一hEli—Sed h。[]一h。Ei一1]+h。Ei] 步骤 2 计算 Ce[i]。 ． 1)如果 知ceEi一1] BaseEi一1]，则 ceEi]一faceEi一1]+Seal—Sub；