马上注册,下载丰富资料,享用更多功能,让你轻松玩转阳光石油论坛。
您需要 登录 才可以下载或查看,没有账号?欢迎注册
x
本帖最后由 刘萝卜锅 于 2020-12-19 17:17 编辑
品质源于技术 服务源于态度
这是阿什卡微信公众号第831篇原创文章 首发于2020年12月14日 ▽
1.概述
椭圆断层流构造模拟方法,是模拟和恢复具有变断距断层过程的独特方法,是MOVE软件开发的新技术,极大的丰富了MOVE软件可处理的断层类型。
此方法和其它构造模拟方法一样,都是根据野外地质露头实测数据统计分析得到的断层几何特征和断距之间的关系,预测地层的变形(Barnett等,1987;Gibson等,1989;Watterson1986)。
一般方法只能预测断层附近的断层相关褶皱,而不能预测变断距情况下的牵引褶皱,所以只能处理具有恒定断距或断距随深度加深而变化(生长断层)的断层(图1)。椭圆断层流方法的出现,解决了具有变断距的断层模拟和恢复的技术难题。
图1(a)使用斜剪切方法正断层发育模拟,断层上盘不同位置具有均一断距;(b)使用椭圆断层流方法,断层断距随远离断层中心逐渐减小。
2.椭圆断层流方法的原理
椭圆断层流方法使用如下理论计算断层周边的位移场:断层相关变形区域的空间范围;最大断层偏移量或滑距;两种公开关系来定义随着远离断层中心而减小的断距(Barnett等,1987;Watterson,1986)。该算法使用反映断层几何形状的椭圆形来定义断层相关变形的空间范围(图2),椭圆形的大小由最大偏移点和断层尖端(Q)之间的距离以及用户定义的半径(R)控制,该半径通常可以根据地震数据或现场观察得到。MOVE软件中,使用Movement表组和Displacement Radius窗口中的上、下、左和右选项来指定此椭圆的空间范围。
图2 使用椭圆断层流计算断距的核心参数:(a)直线正断层;(b)曲线正断层。黑色虚线内的区域表示断层相关形变的范围,彩色图显示断距,包括下盘的挠曲断距。
默认情况下,椭圆断层流方法中的断层偏移量或滑距(sn)根据Watterson (1986)发表的理想断层滑距场理论计算:
式中,Qn为断层上任一点与最大断距点的距离,根据最大断距点到断层端点(Q)的距离进行归一化。对于任何给定的点,Qn通过识别断层上最近的点,测量断层上最近的点到最大偏移点的距离(图2)来计算,然后将该值与Q进行归一化。
随着远离断层,断距逐渐减小,直至为0。所以断距为0的点连接在一起,即形成一个椭圆形,地层只在该椭圆内部有形变。在椭圆断层流方法内,断距递减梯度根据Barnett等人在1987年发表的方程计算,公式如下:
式中,dn为地层上的数据点距断层的距离(图2中的dj),相对于椭圆半径(R)进行归一化处理。Sn是归一化后断层断距,在断层表面的最近点根据方程1计算得到。
对于一个特定位置的点,方程1和方程2给出了归一化的位移值。在椭圆断层流方法中,可以进一步调整该值,用以考虑断层下盘由于地层向上挠曲而产生的位移(Resort 2008)。下盘位移的百分比(SFW)可以自定义(0%下盘无位移,100%上盘无位移),或者使用Gibson(1989)等人的方程根据断层倾角(θ)计算:
下盘中的任意一点,通过将自定义的最大断层断距乘以归一化的位移值和下盘变形百分比(SFW),来计算该点的真正位移大小;上盘中数据点,位移等于自定义的最大断层断距乘以归一化位移值和上盘位移百分比(1-SFW)。
椭圆断层流方法中,可使用位于2D Move-on-Fault和Horizons from Fault工具栏中Movement界面内的Maximum Deflection Angle参数来调整计算得到的位移方向。Maximum Deflection Angle参数的意义为模拟的运动方向与断层倾角之间最大的可能差异,如果数值为0°,地层将平行于断层滑动;如果数值>0,则将定义的角度值乘以归一化的偏差系数(图3)。
图3 不同正断层类型形的态归一化偏差系数和最大可能位移轨迹偏差
使用沿椭圆的弧长比(图3中的da和db)和据方程1、2、3计算得到的归一化位移场来计算归一化偏差系数。在椭圆的每四分之一处,当da和db相等时,最大偏差系数出现,且该值从该点开始线性递减。因此,偏差系数在断层平面和椭圆中轴处为零,并在最大断层偏移点处远离断层垂直延伸。然后,相对于归一化位移场按比例缩放计算出的偏转因子,并期望更大的位移对应更大的偏转。
以最大偏转角> 0°计算的位移矢量模拟了从许多地质力学方法(例如断层响应模型)获得的弯曲粒子路径。但是,实际上,最大偏转角为0°和> 0°时获得的结果相对相似,通常只能观察到很小的几何差异。
相反,调整椭圆的大小、最大断层断距点和断层形状,将会导致计算得到的磁层的几何形状存在明显差异(图4)。只要有可能,应通过分析基础数据来设置以上的参数;使用地层的几何拐点来定义椭圆的大小,并通过地震轴的断点来确定最大断距点和断层平面最可能的位置。
图4 具有相同的200米滑距和45°Maximum Deflection Angle参数、不同的R/Q椭圆比、最大断距点和断层形态的正断层发育模拟结果图:a)直线正断层,R/Q椭圆比0.75,最大断距点位于断层中心;b)直线正断层,R/Q椭圆比为1.5,最大断距点位于断层中心;c)直线正断层,R/Q椭圆比为0.75,最大断距点位于断层下部;d)直线正断层,R/Q椭圆比为0.75,最大断距点位于断层上部;e)铲式正断层,R/Q椭圆比为0.75,最大断距点位于断层中心;f)犁式正断层,R/Q椭圆比为0.75,最大断距点位于断层中心。
3.正演模拟
MOVE软件提供了多种不同构造类型发育模拟技术,来解决各种复杂程度不同的地质问题。
构造发育模拟根据解释好的断层,重现地层的变形过程,可用于验证和评估解释方案的合理性。如果模拟得到的结果与解释结果不匹配,则表明解释方案可能不准确,需要仔细检查。如果匹配,那么构造发育模拟可以提供区域的构造演化或特定油气藏运移和成藏的信息。
MOVE软件中用于构造发育模拟的Horizons from Fault和2D Move-on-Fault工具均提供了人机交互和预览功能,可以方便、快捷的进行构造发育模拟。对于椭圆断层流算法,用于定义变形空间范围的椭圆形态可以进行交互式调整,地层的构造形态可实时更新,从而可以轻松地确定最合适的参数,并建立与基础数据一致的模型。此外,在Horizons from Fault工具中,可以用鼠标拖动地层沿着断层运动,来改变地层的厚度和断层最大断距(图5)。
图5 使用断层椭圆流算法模拟墨西哥湾Baryant Canton地区一条孤立正断层,模拟结果与地震数据一致。
在MOVE软件中,可以模拟更复杂的构造,以便于更好地了解地区的构造演化史。使用2D Kinematic Modelling模块内的各种工具,在模拟过程中可考虑多条断层的影响、同生构造沉积、压实作用、均衡效应、裂谷期后的热沉降以及盐构造等。例如,使用椭圆断层流方法模拟正断层上盘褶皱核部由于逆牵引导致地层塌陷而形成的断裂系统(图6)。在建立此正演模拟模型时,用以表示前构造单元和同构造单元内地层变形程度的应变椭圆,已经加进了2D Move-on-Fault工具中Passive Objects参数栏,发育模拟运算完成后,通过测量整个剖面中变形后应变椭圆的长轴与短轴之比,来估算与背斜核部塌陷断裂系统发育相关的应变分布。当长轴和短轴的长度相似时,应变椭圆几乎没有变形,因此相对应的区域具有相对较低的应变强度。较高的应变导致圆的变形更大,长轴和短轴的长度之间存在较大差异。
图6 铲式正断层上盘核部塌陷断裂系统正演模拟结果,彩色显示的为根据应变椭圆变形得到的应变属性强度。
|