塔中地区构造应力场数值模拟研究.pdfVIP

卷 (Volume)36，期 (Number)2，总 (SUM)133 页(Pages)168—175，2012，5(May，2012) 大 地 构 造 与成 矿 学 GeotectonicaetMetallogenia 塔 中地区构造应力场数值模拟研究 邱登峰，郑孟林，张瑜，张仲培 (中国石化 石油勘探开发研究院 构造与沉积储层实验 室，北京 100083) 摘 要 ：构造应力场研究对于油气运聚分析及储层特征描述具有重要的理论和实际意义。本文根据塔中隆起断裂 特征及其演化史，结合钻井资料 ，选择 中下奥陶统碳酸盐岩为主要研究对象，应用大型有限元分析软件 ANSYS，在 处理断层 问题时运用非连续的接触算法 ，模拟 了中奥陶世末、中泥盆世末 、三叠纪末 、新近纪末 的构造应力场 ，研究 了该区不 同时期的古应力分布特征。研究表 明，最大主压应力受东部车尔臣一星星峡走滑断裂影 响显著，在隆起 区沿断裂带呈条带状低值分布，在不 同时期塔 中 I号断裂带无一例外地位于最大主压应力的最低 区，反映 了塔 中 断裂带尤其是塔 中 I号断裂带是油气运移 的有利指 向区。最小主压应力沿断裂带呈条带状高值分布 ，在断裂 的上 盘、断裂走向发生变化和断裂的倾末端是张应力高值 区，为张裂缝发育的重要构造部位 。最小主压应力和最大剪 应力受车尔 臣 一星星峡走滑 断裂 的影 响明显减弱。 关键词 ：塔中；构造应力场 ；有限元 ；数值模拟 ；非连续 中图分类号 ：P554；TEl21．2 文献标志码 ：A 文章编号 ：1001—1552(2012)02-0168—008 构造应力场就是在一个空间范围内构造应力的 录等 ，2005；杨奎锋等，2006；刘翠荣 ，2002；刘聪等 ， 分布 (万天丰 ，1988)，各种构造迹象 ，包括矿物颗粒 2008；佟 彦 明，2007；张胜利 等，2007；桑广森等 ， 在三度空间排列方位的规律性 ，是岩石在构造应力 2010)。不过，由于计算条件和理论方法 的限制 ，以 场 中应力作用的反映(李 四光 ，1973)。它不仅控制 往的多数工作往往采用连续变形算法 ，在处理断层 着裂缝 的形成和演化 (Treagus，1981；单家增等 ， 问题时，常采用软弱带来模拟断层 的物理性质。这 2004)，也是油气运移 的重要动力之一 (Perrodon， 对断裂构造的模拟有很大的不足 (郑勇等 ，2007；雷 1992；李明诚 ，1994；王毅等 ，2005)。因此 ，开展古构 显权和陈运平 ，2011)，而断裂构造是造成地壳岩体 造应力场研究，尤其是古构造应力场的定量研究，对 中应力发生复杂变化 的主要 因素之一 (苏生瑞 ， 进行含油气盆地储集层特征描述及油气运聚分析具 2002)。不连续体概念和不连续变形理论 的逐渐成 有重要 的理论和实 际意义 (孙晓庆，2008)。构造应 熟使得利用数值模拟方法研究断层真实动力学状态 力场研究的基本方法为利用岩石中已存在的构造形 成为可能(石根华 ，1993；郑勇等 ，2007)。本文 以塔 变特征来反推形变作用发生前后 的构造应力状态 ， 中隆起为研究对象，根据塔 中隆起断裂特征及其演 除地质研究方法之外 ，还有模拟研究方法，即包括各 化史 ，结合钻井资料 ，应用非连续的接触模型模拟断 种物理模拟与数值模拟 (万天丰，1988)。有限单元 层 ，对中下奥陶统不 同时期 的构造应力场进行 了有 法是构造应力场数值模拟的重要方法之一。目前在 限元数值模拟研究，并对研究区内构造应力场与油 构造应力场数值模拟及其与油气运聚关系及裂缝预 气运聚关系进行了初步探讨。 测方面的研究取得了一批非常有价值 的成果 (蒋有 收稿 日期 ：2011—08—15；改回日期 ：2011—11—07 项 目资助 ：国家科技重大专项 “海相碳酸盐岩油气资源