[发明专利]确定层状各向同性介质转换波叠加速度的方法和装置

说明书

技术领域

本发明涉及地震勘探技术，尤其涉及确定层状各向同性介质转换波叠加速度的方法和装置。

背景技术

转换波地震勘探开始于二十世纪八十年代。从那时起，国内外的地球物理学家投入了很多精力来试图建立一个标准的转换波处理流程。与纵波地震数据处理不同的是，直到目前为止，也没有一个公认的转换波地震数据的标准处理流程，这主要是因为转换波传播路径的非对称性所致。

转换波地震数据处理中的一个难题是转换波地震数据的正常时差校正。只有先求取层状各向同性介质转换波的叠加速度，才能实现转换波地震数据正常时差校正、叠加与偏移成像等数据处理工作。

转换波是下行的地震纵波在地下地层界面发生反射、由入射的下行地震纵波转换成的出射的上行地震横波。在转换波地震勘探中，下行波为地震纵波，上行波为地震横波，因为在同一介质中，纵波和横波的传播速度不同，因此即使在水平层状介质中，下行地震纵波和上行地震横波的传播路径也不对称。由于转换波传播路径的非对称性，导致描述转换波的反射旅行时间的时距曲线方程也表现为非双曲形式，因此适用于地震纵波数据的双曲形式的正常时差校正方法不再适用于转换波数据处理。

发明内容

本发明就是鉴于现有技术中的上述问题而提出的。本发明提供了一种新的确定层状各向同性介质转换波叠加速度的方法。

本发明的发明人提出通过利用两个控制参数来描述层状各向同性介质中转换波反射旅行时间，据此对转换波地震数据进行双参数扫描，可以获得转换波双参数叠加速度模型，用于转换波地震数据的正常时差校正处理。

在本发明的一个方面中，提出了一种确定层状各向同性介质转换波叠加速度的方法，该方法包括以下步骤：

步骤A：获取所述研究区域的纵波叠加速度；

步骤B：根据所述纵波叠加速度确定所述转换波叠加速度的取值范围，根据所述目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围，并确定纵波叠加速度和转换波叠加速度的平方比的取值范围；

步骤C：将所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值，其中，M、N和K为自然数；

步骤D：从所述研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据，并抽取所提取的转换波地震数据的渐近共转换点道集；

步骤E：对于所述渐近共转换点道集中的每一个渐近共转换点道集执行以下步骤：

步骤E1：利用所述M个转换波叠加速度值、所述N个纵波叠加速度和转换波叠加速度的平方比值和所述K个转换波垂直双程反射传播时间值，计算该渐近共转换点道集中的每一个地震道的转换波反射旅行时间，并将该渐近共转换点道集中的每一个地震道的转换波反射旅行时间前后2L+1个采样点的数据采样值叠加后求平方并累加求和，以得到该渐近共转换点道集的扫描计算能量数据；

步骤E2：对应转换波垂直双程反射传播时间的每个取值，确定所得到的扫描计算能量数据{E_m，n，k}_{m＝1…M，n＝1…N}中的极大值的位置，提取该极大值的位置处所对应的取值作为所述转换波叠加速度以及所述纵波叠加速度和转换波叠加速度的平方比。

在本发明的另一个方面中，提出了一种确定层状各向同性介质转换波叠加速度的装置，该装置包括：

存储部，该存储部存储研究区域的纵波叠加速度以及转换波地震数据；

取值范围确定部，该取值范围确定部根据所述存储部中存储的所述纵波叠加速度确定转换波叠加速度的取值范围，根据目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围，并且确定纵波叠加速度和转换波叠加速度的平方比的取值范围；

离散部，该离散部将由所述取值范围确定部确定的所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值，其中，M、N和K为自然数；

渐近共转换点道集抽取部，该渐近共转换点道集抽取部从所述存储部中存储的研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据，并抽取所提取的转换波地震数据的渐近共转换点道集；