一种叠前方位P波各向异性裂缝检测方法及装置制造方法及图纸

公开了一种叠前方位P波各向异性裂缝检测方法及装置。所述裂缝检测方法包括以下步骤：1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；3)对步骤2)中投影得到的一组样本点进行椭圆拟合；4)求解出拟合的椭圆方程的最优椭圆参数解；5)基于所述最优椭圆参数解确定椭圆几何参数。本发明专利技术提出了一种叠前方位P波各向异性裂缝检测方法，提取P波叠前分方位振幅，求解椭圆方程并进行裂缝预测，可以为裂缝性油气藏的勘探开发提供重要的技术手段与评价工具。

A Pre-stack Azimuth P-wave Anisotropic Crack Detection Method and Device

A prestack azimuth P-wave anisotropic crack detection method and device are disclosed. The crack detection method includes the following steps: 1) extracting the azimuth-based seismic amplitude data of stacked gathers in the unit interval along the target layer; 2) projecting the azimuth amplitude information obtained in step 1 into polar coordinates; 3) fitting a group of sample points projected in step 2; 4) solving the optimal elliptic parameter solution of the fitting elliptic equation; 5) solving the optimal elliptic parameter solution based on the said optimum. The elliptic parameter solution determines the elliptic geometric parameters. The invention provides a prestack azimuth P-wave anisotropic fracture detection method, extracts the prestack azimuth amplitude of P-wave, solves elliptic equation and carries out fracture prediction, which can provide important technical means and evaluation tools for exploration and development of fractured reservoirs.

【技术实现步骤摘要】
一种叠前方位P波各向异性裂缝检测方法及装置
本专利技术涉及油气地球物理勘探领域，更具体地，涉及一种叠前方位P波各向异性裂缝检测方法及装置。
技术介绍
地震层位曲率可以用来描述地层的弯曲程度，进而可以表征地层所受的应力大小，因此可以用来预测受应力而形成的裂缝发育带。利用最大曲率、最小曲率、最大正曲率、最小负曲率等属性可以有效弥补相干属性在该工区预测裂缝时的不足。目前，利用地震资料来进行断裂检测的技术按照所用资料的不同分为两大类：一类是叠后预测技术，一类是叠前预测技术。常用叠后技术有相干分析技术、图像处理技术、曲率分析技术、倾角检测技术、应力分析技术等，它们通常以叠后地震资料为基础输入，研究的是整个网状断裂系统中那些较大尺度的断裂，从预测效果上来看预测的是网状断裂系统的骨架与轮廓，在平面上呈线性分布，数量较少；而对较小尺度、数量更庞大的小型或微断裂系统，叠后断裂预测技术根本无能为力。叠后断裂预测技术之所以失效，其根本原因在于叠后地震资料的信息量比较小，缺乏偏移距信息和方位角信息，失去了以方位各向异性来检测微断裂的理论依据和基础。因此，有必要开发一种以方位各向异性来检测微断裂的叠前预测方法。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
在碳酸盐岩、致密砂岩等致密地层中，油气储量及产量约占世界油气储量及产量的一半，碳酸盐岩、致密砂岩等致密地层中油气与裂缝关系密切，因此确定裂缝发育的区域、方位和密度在裂缝性油气勘探中具有重大意义。在垂直定向裂缝的物理模型实验及数值模拟中都证实了地震勘探过程中，P波在通过地质裂缝体时表现出很强的方位各向异性特征，P波反射振幅及旅行时与测线和裂缝方位有关，测线与裂缝平行时振幅最强、旅行时最短，随着测线与裂缝夹角的增大，振幅逐渐减弱、时间逐渐变长，至测线与裂缝方向垂直时，振幅最弱，时间最长。根据这些特征，在利用P波进行裂缝检测方面已经取得了比较完善的方法技术与成果，并且在实际资料的应用中取得了成功。另外，为了克服叠后裂缝预测技术的缺点，利用叠前地震资料具有数据量大、信息量大的特点，检测包含在叠前(观测)方位—偏移距(入射角)二维时空域中微断裂的地震响应信息来检测整个断裂体系，具有较强的理论依据和数据基础。因此，本专利技术提出了一种叠前方位P波各向异性裂缝检测方法，提取P波叠前分方位振幅，求解椭圆方程并进行裂缝预测，可以为裂缝性油气藏的勘探开发提供重要的技术手段与评价工具。根据本专利技术的一方面，提出一种叠前方位P波各向异性裂缝检测方法。所述方法可以包括以下步骤：1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；3)对步骤2)中投影得到的一组样本点进行椭圆拟合；4)求解出拟合的椭圆方程的最优椭圆参数解；5)基于所述最优椭圆参数解确定椭圆几何参数。优选地，基于椭圆一般方程Ax2+Bxy+Cy2+Dx+Ey+F＝0定义投影平面上一个样本点(x,y)到椭圆的几何距离为：d(a,x)＝Ax2+Bxy+Cy2+Dx+Ey+F(1)其中，a＝[ABCDEF]T，x＝[x2xyy2xy1]T。优选地，通过以下公式求得最小二乘意义下的最优椭圆参数解：优选地，通过将方程式(1)表示为广义特征值问题的形式：GTGa＝λHa(3)其中，λ为特征值，H为约束矩阵，a为特征向量a＝[ABCDEF]T，式中x1,x2,...,xn为第1、第2…第n点的横坐标，y1,y2,...,yn为第1、第2…第n点的纵坐标；基于椭圆参数的数值关系约束矩阵求解方程式(3)以得到特征向量a。优选地，通过以下公式计算椭圆几何参数：其中，(xe,ye)为椭圆中心相对坐标，a、b为椭圆长轴长度、短轴长度，θ为椭圆长轴与正北方向的夹角。根据本专利技术的另一方面，提出一种叠前方位P波各向异性裂缝检测装置。所述装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；3)对步骤2)中投影得到的一组样本点进行椭圆拟合；4)求解出拟合的椭圆方程的最优椭圆参数解；5)基于所述最优椭圆参数解确定椭圆几何参数。优选地，基于椭圆一般方程Ax2+Bxy+Cy2+Dx+Ey+F＝0定义投影平面上一个样本点(x,y)到椭圆的几何距离为：d(a,x)＝Ax2+Bxy+Cy2+Dx+Ey+F(1)其中，a＝[ABCDEF]T，x＝[x2xyy2xy1]T。优选地，通过以下公式求得最小二乘意义下的最优椭圆参数解：优选地，通过将方程式(1)表示为广义特征值问题的形式：GTGa＝λHa(3)其中，λ为特征值，H为约束矩阵，a为特征向量a＝[ABCDEF]T，式中x1,x2,...,xn为第1、第2…第n点的横坐标，y1,y2,...,yn为第1、第2…第n点的纵坐标；基于椭圆参数的数值关系约束矩阵求解方程式(3)以得到特征向量a。优选地，通过以下公式计算椭圆几何参数：其中，(xe,ye)为椭圆中心相对坐标，a、b为椭圆长轴长度、短轴长度，θ为椭圆长轴与正北方向的夹角。本专利技术涉及的叠前预测方法是一种基于椭圆方程的叠前方位P波各向异性裂缝检测方法。根据P波在通过裂缝体时表现出很强的方位各向异性特征：其叠前地震资料的方位振幅、NMO速度、频率、波阻抗以及振幅衰减等多个参数都随着方位角的不同呈椭圆形变化，通过对给定平面(超道集)范围上的一组样本点，以椭圆方程为模型进行拟合求解超定方程组，使单位平面内的椭圆方程尽量接近这些样本数据点，从而获取该椭圆方程的各个参数，达到了获取单元范围内裂缝体的发育方向与强度的目的。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1为根据本专利技术的示例性实施方案的叠前方位P波各向异性裂缝检测方法的流程图；图2为基于方位的叠加数据提取的结果；图3为方位振幅投影图；图4为椭圆拟合示例图；图5为方位振幅椭圆拟合示例图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。以下参考图1详细描述根据本专利技术示例性实施例的叠前方位P波各向异性裂缝检测方法。该方法主要包括：步骤1：沿目的层提取单元区间内基于方位的叠加道集地震振幅数据。如显示基于方位的叠加数据提取结果的图2所示，横轴代表方位角，单位为度(°)，图中所示为15°、45°、75°、105°、135°、165°等6个方位角的6道本文档来自技高网...

【技术保护点】
1.一种叠前方位P波各向异性裂缝检测方法，其特征在于，所述裂缝检测方法包括以下步骤：1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；3)对步骤2)中投影得到的一组样本点进行椭圆拟合；4)求解出拟合的椭圆方程的最优椭圆参数解；5)基于所述最优椭圆参数解确定椭圆几何参数。

【技术特征摘要】
1.一种叠前方位P波各向异性裂缝检测方法，其特征在于，所述裂缝检测方法包括以下步骤：1)沿目的层提取单元区间内基于方位的叠加道集地震振幅数据；2)将步骤1)获取的各方位振幅信息在极坐标中进行投影；3)对步骤2)中投影得到的一组样本点进行椭圆拟合；4)求解出拟合的椭圆方程的最优椭圆参数解；5)基于所述最优椭圆参数解确定椭圆几何参数。2.根据权利要求1所述的叠前方位P波各向异性裂缝检测方法，其特征在于，基于椭圆一般方程Ax2+Bxy+Cy2+Dx+Ey+F＝0定义投影平面上一个样本点(x,y)到椭圆的几何距离为：d(a,x)＝Ax2+Bxy+Cy2+Dx+Ey+F(1)其中，a＝[ABCDEF]T，x＝[x2xyy2xy1]T。3.根据权利要求2所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过以下公式求得最小二乘意义下的最优椭圆参数解：4.根据权利要求3所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过将方程式(1)表示为广义特征值问题的形式：GTGa＝λHa(3)其中，λ为特征值，H为约束矩阵，a为特征向量a＝[ABCDEF]T，式中x1,x2,...,xn为第1、第2…第n点的横坐标，y1,y2,...,yn为第1、第2…第n点的纵坐标；基于椭圆参数的数值关系约束矩阵求解方程式(3)以得到特征向量a。5.根据权利要求4所述的叠前方位P波各向异性裂缝检测方法，其特征在于，通过以下公式计算椭圆几何参数：其中，(xe,ye)为椭圆中心相对坐标，a、b为椭圆长轴长度、短轴长度，θ为椭圆长轴与正北方向的夹角。6.一种叠前方位P波各向异性裂缝检测装置，其特征...

【专利技术属性】
技术研发人员：陈楠，徐雷鸣，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11