用方位属性识别小断层.docVIP

倾角和方位角属性在小断层识别中的应用实例 摘要:当隐蔽断层和其它不连续界面在3D地震资料上显示出微小变化时，常规地震解释不容易找出测线之间的相关性，因此，通常很难绘出隐蔽断层和其它不连续界面的图像。为了绘制这些变化的图像，Mike Bahorich和Steve Farmer （1995年）介绍了一种利用3D叠加数据进行不连续界面成像的处理方法，该方法将反射振幅3D数据体转换成定义地震波形区域相似性的相干数据体。除了从3D叠加地震资料直接提相干性之外，可能抽提出另外两种代表时间-倾向值和倾斜方位角的属性。为了导出这些属性，使用了印度石油天然气委员会（ONGC）开发的一项尖端技术。本文介绍了该项技术在印度Assam盆地Sonari地区3D地震体的应用及其结果。 　　一、引 言 　　在地震资料解释中，最重要的问题之一是断层及断裂带的检测和识别。传统方法的基础是在叠加地震剖面上用肉眼识别反射层的不连续性。当解释者的地质知识能够控制时，这样的方法是最有效的，然而，这样的方法主要取决于解释者的经验和前人解释的地质信息。因为在处理过程的反射同向轴的定线中可能平滑或者完全消除小的不连续层，所以小幅度断层的多解性特别严重。此外，不仅断层在地震剖面上产生这样的不连续层，静态问题、速度误差等也可能产生不连续层。 　　为了从地震资料中获得可能的地质构造不连续层独立的、无系统误差的信息，期望用一些只依靠地震资料的自动处理方法。这样的方法必须利用在地震剖面上与存在断层有关的特定构造知识。然后，可以用全部有效的地质资料分析和验证这一信息（Evgeny Landa等，1998年）。 　　为了适应这种需要，ONGC的J．V．S．N．Murty博士和G．V．Rao发明了这项技术；即众所周知的DATS（倾角和方位角转换软件）。该项技术接收西部编码或SEGY格式的3D叠加地震数据，并将其转换为称作时间-倾向和倾斜方位角数据体的不同数据体。这项技术介绍了如何产生属性时间-倾向和倾斜方位角的3D数据体，目前还没有这方面出版物或者专利。该项专利被认为是计算时间倾向和倾斜方位角数据体的技术发明。这些转换的3D数据体重点突出了不同种类的原始地震数据，并且有助于了解构造形态和其它的不连续层，而只使用叠加的3D地震数据很难解释出这些特征。用该项技术得出的倾斜方位角并需要预先解释数据。为了研究这些结果，应用该项技术大约转换了75.6平方公里的3D数据体。主测线262条，16～278号，测线间距为100m；连络测线444条，每条测线上的地震道距为25m。 　　研究区位于Assam - Arkan盆地之中，它是印度的主要石油产区之一。该盆地包括上Assam、 Dhansiri谷、Schuppen 断裂带、Manipur、Mizoram、Cachar 和 Tripura地区，总面积为10万平方公里。因为白垩纪早期的大陆裂谷作用，该盆地得以继承发展。在始新世至渐新世期间，印度板块与缅甸板块碰撞导致了部分印度洋的完全闭合。这一碰撞形成了一系列迭瓦状逆掩断层和隆起，以及古老沉积物的叠置，它们最终叠置到Assam陆架上，这一逆掩断裂带被称为Schuppe断裂带。这个开始于始新世的逆掩断裂作用到中新世中期更加强烈。 　　印度国家石油公司和ONGC已经发现了许多石油和天然气。其中一些油田是Nahorkatia、Digboi、Moran、Ruthrasagar、Lakwa、 Gelki、 Demulgaon、Borholla、Khoraghat和Badarpur，估计在该盆地最深处的沉积层厚度超过11km。 　　大多数油田是渐新统和中新统的Barail组和Tipam组内的构造圈闭。然而，在Sylhet 组、Tura组和裂缝基底也发现了油气。Kopili页岩是该盆地的主要生油岩。所有的构造形态分布图都是借助于地震反射资料完成的。 　　二、理论和方法 　　考虑一个体积元： 　　dv ＝ dX ． dY． dT （1） 　　将样品置于3D叠加数据中心，该数据的dX 和 dY 具有100m的数量级，dT 具有子波长度的数量级。如果计算了体积元内包含的使用地震波形的归一化能量的倾角和方位角，那么，沿着该层的倾角和方位角（如果存在）定位，这个体积元归一化能量达到最大。因此，能够将样品的倾向数量级计算成微秒/米，样品的倾斜方位角与网格的北向成比例。为了进行计算，用两个互相垂直方向上的不同倾角计算了归一化能量， 即X （主测线）方向、Y （连络测线）方向，得到的时间-倾角和倾斜方位角为： 　　时间-倾角 ＝ sqrt（ dipX * dipX ＋ dipY * dipY） （2） 　　倾斜方位角 ＝ atan（ dipX / dipY） （3） 　　式中： dipX和dipY 分别是X和Y方向上的归一化能量最