一种基于层位模型的倾角计算方法技术

本发明专利技术提供了一种基于层位模型的倾角计算方法。所述方法包括：在地震深度剖面上进行层位拾取，并构建层位模型；在层位模型上选定一深度位置信息，找到与该深度位置信息的道号对应的共中心点道集，对所述共中心点道集所对应的多个振幅进行相关运算，得到最大能量；以所述最大能量对应的深度信息替换所述深度位置信息的深度信息；以替换后的深度位置信息计算所述选定深度位置的反射面的倾角属性。本发明专利技术能够获得更加准确的地震剖面倾角属性；不用求解复杂的计算方程组，更加简便高效。

【技术实现步骤摘要】

本专利技术涉及地震几何属性计算
，具体来讲，涉及一种地震倾角属性的计算方法，更具体来讲，涉及一种用于石油地震勘探中地质模型速度建模中网格层析的关键方法——地震倾角属性计算方法。
技术介绍
通常，在石油地震勘探领域中，断层和裂缝作为十分重要的地质结构往往是找气找油的关键。地震的几何属性揭示了地震数据在空间上的几何结构特征，通过描述相邻道的关系，计算出地层反射界面的几何特征，能够真实的再现地下构造形态，成为了描述刻画地下断层裂缝的重要工具。地震倾角属性作为地震几何属性的一员，在地震勘探资料处理与解释中也被广泛的使用。倾角属性是地震解释的一个重要工具，它能够确定局部反射面，处理人员能够通过倾角属性来判断反射面是否连续并提取起连续有用的部分用于网格层析速度分析等后续处理。然而，由于常规的基于层位模型的倾角计算方法通常通过人工拾取层位的深度信息来计算反射面倾角数据，所以通常难以得到比较准确的倾角属性。
技术实现思路
本专利技术的目的在于解决现有技术存在的上述不足中的至少一项。专利技术人经过研究发现，常规的基于层位模型的倾角计算方法之所以存在难以得到比较准确的倾角属性的问题，其主要原因是速度畸变以及人工拾取层位与真实的地震数据往往存在误差(甚至于较大的偏差)。为了解决计算倾角属性时由于速度畸变以及人工拾取层位带来的计算误差问题，本专利技术提供了一种基于层位模型的倾角计算方法。所述方法包括：在地震深度剖面上进行层位拾取，并构建层位模型；在层位模型上选定一深
度位置信息，找到与该深度位置信息的道号对应的共中心点道集，对所述共中心点道集所对应的多个振幅进行相关运算，得到最大能量；以所述最大能量对应的深度信息替换所述深度位置信息的深度信息；以替换后的深度位置信息计算所述选定深度位置的反射面的倾角属性。在本专利技术的一个示例性实施例中，所述相关运算通过相关函数来进行。在本专利技术的一个示例性实施例中，所述计算选定深度位置的反射面的倾角属性的步骤可通过下式计算，ax＝cosθcosφay＝cosθsinφaz＝sinφ，其中，ax，ay，az分别表示反射面倾角单位向量沿反射面的平面元素点的x，y，z方向的分向量，倾角θ为反射面与由x、y构成的水平面之间的夹角，方位角φ为法平面与y轴的夹角。与现有技术相比，本专利技术的有益效果包括：能够获得更加准确的地震剖面倾角属性；不用求解复杂的计算方程组，更加简便高效。附图说明图1示出了根据常规方法所构建的层位模型的示意图。图2示出了根据本专利技术的方法更新层位信息的示意图。图3示出了根据本专利技术的方法计算倾角的示意图。具体实施方式在下文中，将结合附图和示例性实施例来详细说明本专利技术的基于层位模型的倾角计算方法。在本专利技术的一个示例性实施例中，基于层位模型的倾角计算方法可按照如下步骤进行：首先，可根据常规方法在地震深度剖面上进行人工的层位拾取，并构建层位模型(例如，可为封闭的层位模型)。在地震剖面上所拾取的数据包含多个深度位置的信息(例如，多个点的信息)，每个深度位置的信息由深度和道号构成。如图1所示，构建的层位模型中横坐标表示道号，纵坐标表示深度。接下来，如图1所示，在层位模型上选定某一深度位置，例如，可选定待确定其倾角属性的深度位置。该深度位置的信息包括道号信息和深度信息。然后，根据选定深度位置信息，找到与该深度位置信息的道号对应的共中心点(CMP)道集，对CMP道集所对应的多个振幅进行常规的相关运算。例如，可通过下面的式(1)进行相关运算。式(1)的Rf(τ)为相关函数，其反映的是f(τ)的两个函数在不同的相对位置上互相匹配的程度。另外，式(1)也能够描述随机信号f(t)在不同的深度t1，t2的平均值作为延迟深度t的函数，是随机信号与延迟信号f(τ)的相似性度量，当延迟信号为零时，则成为随机信号的均方值，此时能量最大。 R f ( τ ) = f ( τ ) * f * ( - τ ) = ∫ - ∞ ∞ f ( t + τ ) f * ( t ) d t = ∫ - ∞ ∞ f ( t ) f * ( t - τ ) d t - - - ( 1 ) ]]>其中，“*”是卷积算符；(·)*为取共轭；f(t)为地震信号，也就是振幅信号；τ为f(t)的某一深度，属于两个深度(如，t1、t2)区间之内。另外，式(1)也能够描述随机信号f(t)在不同的深度t1，t2的平均值(τ即为不同的深度t1，t2的平均值)作为延迟深度t的函数，是随机信号与延迟信号f(τ)的相似性度量，当延迟信号为零时，则成为随机信号的均方值，此时能量最大。接下来，通过对CMP道集所对应的多个振幅进行相关运算，得到最大能量。以该最大能量换算得到深度值，并用该深度值更新上述选定深度位置信息的深度信息。图2示出了根据本专利技术的方法更新层位信息的示意图。随后，运用更新后的层位深度信息计算该选定深度位置的反射面的倾角属性。倾角属性的计算方式可如图3所示，具体描述如下：a、由层位上更新后的深度信息找到该层位深度处的反射界面；b、设反射面的平面元素为点X＝(x，y，z)，令反射面的法向量为n＝(nx，ny，nz)，其中nx，ny和nz分别是沿x，y和z轴上的分量，反射面倾角单位向量为a＝(ax，ay，az)；c、其中，水平面与反射面的夹角为倾角θ，法平面与y轴的夹角则为方位角Φ，于是可得出下面的式(2)：ax＝cosθcosφay＝cosθsinφaz＝sinφ (2)即可算出倾角与方位角。综上所述，本专利技术的方法能够获得更加准确的地震剖面倾角属性；不用求解复杂的计算方程组，更加简便高效。尽管上面已经结合附图和示例性实施例描述了本专利技术，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。本文档来自技高网...

【技术保护点】
一种基于层位模型的倾角计算方法，所述方法包括：在地震深度剖面上进行层位拾取，并构建层位模型，其特征在于，所述方法还包括以下步骤：在层位模型上选定一深度位置，该深度位置的信息由道号和深度构成，找到与该深度位置信息的道号对应的共中心点道集，对所述共中心点道集所对应的多个振幅进行相关运算，得到最大能量；以所述最大能量计算得到的深度替换所述深度位置信息的深度；以替换后的深度位置信息计算所述选定深度位置的反射面的倾角属性。

【技术特征摘要】
1.一种基于层位模型的倾角计算方法，所述方法包括：在地震深度剖面上进行层位拾取，并构建层位模型，其特征在于，所述方法还包括以下步骤：在层位模型上选定一深度位置，该深度位置的信息由道号和深度构成，找到与该深度位置信息的道号对应的共中心点道集，对所述共中心点道集所对应的多个振幅进行相关运算，得到最大能量；以所述最大能量计算得到的深度替换所述深度位置信息的深度；以替换后的深度位置信息计算所述选定深度位置的反射面的倾角属性。2...

【专利技术属性】
技术研发人员：陈三平，唐虎，徐春春，张恩嘉，段鹏飞，熊晶璇，陈立，王珑，孙伟鹏，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51