一种地震数据处理VDA双参数分析方法与流程

一种地震数据处理vda双参数分析方法
技术领域
[0001]
本发明涉及地震数据处理技术领域，特别涉及一种地震数据处理vda双参数分析方法。


背景技术：

[0002]
在地震数据处理过程中，有多种零偏移距成像和速度分析方法。
[0003]
在几乎所有地震数据处理系统赖以为基础的共中心点(cmp)方法中，cmp速度v
cmp
其实只是一个信号叠加参数，没有地质意义，这是因为所有实际数据远不是理想的水平层状均匀介质，而只有在水平层状均匀介质的情况下才能得到合理的均方根速度v
rms
。反射界面的稍微倾斜或者出现速度非均质性，都会导致v
cmp
的急剧变化。
[0004]
还有其它一些得到零偏移距成像剖面的方法，如倾角时差校正(dmo)叠加、共反射面元(crs)叠加和多聚焦(mf)成像等，这些方法都是经过各自不同的变换进行时距曲线叠加。但是，时距曲线的反射信号属于反射界面上的不同点。
[0005]
利用不同方法求得的零偏移距成像速度，它们往往都与地下反射界面的倾角和偏移距有关，这就不可避免地导致在叠加时得到的是等效速度，与真正的均方根速度v
rms
存在差别，有些差别很大。众所周知，根据等效速度只有在一种情况下能够得到正确的层速度，那就是理想水平层状均匀介质。对于具有倾向倾角各异的若干反射界面等复杂地质情况来说，不管用哪种方法，得到的速度的准确性更低，这个速度往往只是一个叠加参数，没有实际地质意义，这势必影响成像的准确性，无法实现真正意义上的共反射点零偏移距成像，获得真正的零偏移距时间剖面。


技术实现要素：

[0006]
本发明的目的旨在至少解决所述技术缺陷之一。
[0007]
为此，本发明的目的在于提出一种地震数据处理vda双参数分析方法。
[0008]
为了实现上述目的，本发明的实施例提供一种地震数据处理vda双参数分析方法，包括如下步骤：
[0009]
步骤s1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：l为炮点s与检波点r之间的距离，v为地震波传播速度，l0表示反射点o处反射界面的法线在地面出射点c与炮点的距离，β为入射波射线和反射波射线之间夹角的一半，θ为反射点o处反射界面的地层倾角；
[0010]
步骤s2，通过炮点s、检波点r和反射点o构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角2β固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点；
[0011]
相等的角β将圆弧sr分成相等的两段弧，该圆的圆心位于炮点s-检波点r线段的中垂线上；令t表示沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时；其中反射
点o处法线与炮检点的连线sr相交于点c，c即为反射点o对应的成像道的地表位置，t0表示中间法向射线oc段的双程旅行时，令表示中间法向射线延长到p点即op段对应的双程旅行时，v表示地震波速度，计算的vda双参数分析方法的变换方程为：
[0012][0013]
步骤s3，通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
[0014]
进一步，所述θ进一步为法向射线与垂向线的夹角。
[0015]
进一步，所述计算vda双参数分析方法的变换方程，包括如下步骤：
[0016]
已知如下数学关系：
[0017][0018]
由圆的性质，有如下关系式：
[0019][0020]
把(2)中的关系代入(1),得到:
[0021][0022]
通过计算可以求取圆的直径d：
[0023]
d＝l/(2sin2β)
ꢀꢀꢀꢀ
(4)
[0024]
进而可以推出如下关系式，
[0025][0026]
利用公式(1)-(4)，可以推出
[0027][0028]
结合椭圆展开变换方程，有
[0029][0030]
公式(6)和(7)联立,可以推出
[0031][0032]
地层倾角θ与射线参数的关系式为：
[0033][0034]
进而变换得出关于速度v和地层倾角θ的新方程
[0035][0036]
进一步，所述步骤s3，还包括：根据所述vda双参数分析方法的变换方程，生成综合地层倾角和速度的双参数谱，并估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，以获得更理想的成像剖面。
[0037]
根据本发明实施例的地震数据处理vda双参数分析方法，基于vda双参数成像原理，估算制作用于双参数分析的叠加能量剖面以及速度和地层倾角双参数谱，通过分析获得参数速度和地层倾角，最终实现双参数共反射点成像。该方法没有地下反射层为水平的假设前提条件，地下反射层可以倾斜或弯曲。该方法可以估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，考虑此双参数的叠加成像方法就消除了传统单一速度分析方法估算的叠加速度与实际地震波射线速度不同导致的成像计算误差，修正仅使用速度一个参数进行叠加时的不聚焦现象，提高了成像质量。
[0038]
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0039]
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
[0040]
图1为根据本发明实施例的地震数据处理vda双参数分析方法的流程图；
[0041]
图2为根据本发明实施例的地震波传播路径示意图；
[0042]
图3为根据本发明实施例的时vda双参数成像及双参数分析原理图；
[0043]
图4为根据本发明实施例的vda双参数分析界面图；
[0044]
图5为根据本发明实施例的vda双参数分析获得的地层倾角剖面图。
具体实施方式
[0045]
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0046]
本发明提供一种地震数据处理vda双参数分析方法，考虑地层倾角和速度两个参数进行零偏移距成像，修正仅使用速度一个参数进行零偏移距成像时的不聚焦现象，提高成像质量。
[0047]
如图1所示，本发明实施例的地震数据处理vda(velocity andformationdipangle)双参数分析方法，包括如下步骤：
[0048]
步骤s1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：l为炮点s与检波点r之间的距离，v为地震波传播速度，l0表示反射点o处反射界面的法线在地面出射点c与炮点的距离，β为入射波射线和反射波射线之间夹角的一半，θ为反射点o处反射界面的地层倾角。θ进一步为法向射线与垂向线的夹角。
[0049]
反射波速度与入射波速度相等时,地震波传播路径如图2所示。
[0050]
步骤s2，通过炮点s、检波点r和反射点o构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角2β固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点。
[0051]
相等的角β将圆弧sr分成相等的两段弧，该圆的圆心位于炮点s-检波点r线段的中垂线上；令t表示沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时；其中反射点o处法线与炮检点的连线sr相交于点c，c即为反射点o对应的成像道的地表位置，t0表示中间法向射线oc段的双程旅行时，令表示中间法向射线延长到p点即op段对应的双程旅行时，v表示地震波速度，计算的vda双参数分析方法的变换方程为：
[0052][0053]
图3展示了vda双参数分析方法原理。如图3所示，设置炮点s、检波点r，其位置坐标已知。它们对应地下某任意反射地层上的一个反射点o，三者满足地震波反射定律。o处相应地层倾角为θ。下面推导用于vda双参数分析的变换方程。炮点s、检波点r和反射点o三点确定一个圆。根据反射定律，在o点处作反射界面的法线oc，其在炮检点连线sr上的出露点为c，与圆的交点为p。oc即入射波射线so和反射波射线or的夹角的角平分线，将角∠sor分为入射角和反射角，二者相等，记为β。c点即为反射点o对应的零偏移距成像道的地表位置。
[0054]
该圆具有很好的几何性质：该圆可以看作是一组可能反射点的轨迹，对于该圆上的每一个反射点，入射线和反射线之间的夹角固定不变；在每个反射点处作角平分线，该角平分线即是各个反射点处反射界面的法线，所有的法线都恰巧相交于点p，p点也是圆和纵向垂直的直径的交点，称之为极点；该圆的圆心位于sr线段的中垂线上。
[0055]
令t表示地震波沿入射波射线的旅行时t1和反射波射线的旅行时t2的总旅行时。t0表示法向射线oc段的双程旅行时。令表示法向射线延长到p点即op段对应的双程旅行时，
v表示地震波速度。l表示炮点s与检波点r之间的距离，l0表示反射点o处反射界面的法线在sr上的出射点c与炮点s的距离。
[0056]
计算vda双参数分析方法的变换方程，包括如下步骤：
[0057]
已知如下数学关系：
[0058][0059]
由圆的性质，有如下关系式：
[0060][0061]
把(2)中的关系代入(1),得到:
[0062][0063]
通过计算可以求取圆的直径d：
[0064]
d＝l/(2sin2β)
ꢀꢀꢀꢀ
(4)
[0065]
进而可以推出如下关系式，
[0066][0067]
利用公式(1)-(4)，可以推出
[0068][0069]
结合椭圆展开变换方程，有
[0070][0071]
公式(6)和(7)联立,可以推出
[0072][0073]
地层倾角θ与射线参数的关系式为：
[0074][0075]
进而变换得出关于速度v和地层倾角θ的新方程
[0076]
[0077]
公式(10)是均匀介质时的vda双参数分析的变换方程，该方程中含有速度和地层倾角两个参数，对于均匀介质来说，它绝对准确，通过该方程，可以实现速度和地层倾角双参数分析，估算拾取速度和地层倾角，用于零偏移距成像。步骤s3，通过上述变换方程，可以把t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
[0078]
具体的，公式(10)是vda双参数分析方法的变换方程，该方程中含有速度和地层倾角两个参数，绝对准确，通过该方程，可以实现速度和地层倾角双参数分析。
[0079]
步骤s3，通过上述变换方程，计算生成速度和地层倾角双参数分析界面，包括(1)生成双参数谱，用于双参数分析估算速度和地层倾角双参数值；(2)将t转换为t0，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，辅助双参数分析。
[0080]
图4和图5为某数据的vda双参数分析界面。
[0081]
图4中左图为叠加能量剖面与零偏移距成像剖面的叠合显示，从浅至深，无论是断面、低倾角地层、陡倾角地层、地质异常体和背斜构造，二者吻合一致，各反射层能量集中，干涉噪音小，地层接触关系清晰。右图为cp1380处的双参数谱的时间-速度剖面，用于显示速度变化趋势和选择要人工分析的双参数谱的时间。中间小图为cp1380处t0为1328ms时刻的双参数谱的速度-地层倾角剖面，横向尺度为速度，纵向尺度为地层倾角。通过该谱可以实现双参数分析，估算获得双参数的值。
[0082]
此外，步骤s3还包括：根据所述vda双参数分析方法的变换方程，生成综合地层倾角和速度的双参数谱，并估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，以获得更理想的成像剖面。具体的，图5为本发明自动估算的地层倾角剖面与零偏移距成像剖面的叠合显示。从浅至深，无论是低倾角地层还是陡倾角地层，估算的地层倾角的变化趋势与实际零偏移距成像剖面上反射层的形态基本一致，较好的说明了该方法的有效性。右下角小图为cp1380处t0为1328ms时刻的双参数谱速度-地层倾角剖面，横向尺度为速度，纵向尺度为地层倾角。通过该谱可以实现双参数分析，估算获得不同位置、不同零偏移距成像时间处相应的速度和地层倾角双参数的值，以获得更理想的零偏移距成像剖面。
[0083]
根据本发明实施例的地震数据处理vda双参数分析方法，基于vda双参数成像原理，估算制作用于双参数分析的叠加能量剖面以及速度和地层倾角双参数谱，通过分析获得参数速度和地层倾角，最终实现双参数共反射点成像。该方法没有地下反射层为水平的假设前提条件，地下反射层可以倾斜或弯曲。该方法可以估算不同位置、不同成像时间处相应的速度和地层倾角双参数的值，考虑此双参数的叠加成像方法就消除了传统单一速度分析方法估算的叠加速度与实际地震波射线速度不同导致的成像计算误差，修正仅使用速度一个参数进行叠加时的不聚焦现象，提高了成像质量。
[0084]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0085]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例
性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。