确定层状各向同性介质转换波叠加速度的方法和装置制造方法及图纸

技术编号:8734813 阅读:269 留言:0更新日期:2013-05-26 11:38
本发明专利技术提供了确定层状各向同性介质转换波叠加速度的方法和装置,该方法包括:确定目标地层的深度;获取研究区域的纵波叠加速度;确定转换波叠加速度、转换波垂直双程反射传播时间、纵波叠加速度和转换波叠加速度的平方比的取值范围;分别对这些取值范围进行离散采样;从转换波地震数据中提取近、中炮检距的转换波地震数据,并抽取研究区域转换波地震数据的所有渐近共转换点道集;对于每一个渐近共转换点道集:利用上述经离散采样的值,计算该渐近共转换点道集中每一个地震道的转换波反射旅行时间,并得到该渐近共转换点道集的扫描计算能量数据;对应于每个转换波垂直双程反射传播时间的离散采样值,确定扫描计算能量数据中的极大值的位置,提取该极大值的位置处所对应的转换波叠加速度、纵波叠加速度与转换波叠加速度的平方比。

【技术实现步骤摘要】

本专利技术涉及地震勘探技术,尤其涉及确定层状各向同性介质转换波叠加速度的方法和装置
技术介绍
转换波地震勘探开始于二十世纪八十年代。从那时起,国内外的地球物理学家投入了很多精力来试图建立一个标准的转换波处理流程。与纵波地震数据处理不同的是,直到目前为止,也没有一个公认的转换波地震数据的标准处理流程,这主要是因为转换波传播路径的非对称性所致。转换波地震数据处理中的一个难题是转换波地震数据的正常时差校正。只有先求取层状各向同性介质转换波的叠加速度,才能实现转换波地震数据正常时差校正、叠加与偏移成像等数据处理工作。转换波是下行的地震纵波在地下地层界面发生反射、由入射的下行地震纵波转换成的出射的上行地震横波。在转换波地震勘探中,下行波为地震纵波,上行波为地震横波,因为在同一介质中,纵波和横波的传播速度不同,因此即使在水平层状介质中,下行地震纵波和上行地震横波的传播路径也不对称。由于转换波传播路径的非对称性,导致描述转换波的反射旅行时间的时距曲线方程也表现为非双曲形式,因此适用于地震纵波数据的双曲形式的正常时差校正方法不再适用于转换波数据处理。
技术实现思路
本专利技术就是鉴于现有技术中的上述问题而提出的。本专利技术提供了一种新的确定层状各向同性介质转换波叠加速度的方法。本专利技术的专利技术人提出通过利用两个控制参数来描述层状各向同性介质中转换波反射旅行时间,据此对转换波地震数据进行双参数扫描,可以获得转换波双参数叠加速度模型,用于转换波地震数据的正常时差校正处理。在本专利技术的一个方面中,提出了一种确定层状各向同性介质转换波叠加速度的方法,该方法包括以下步骤:步骤A:获取所述研究区域的纵波叠加速度;步骤B:根据所述纵波叠加速度确定所述转换波叠加速度的取值范围,根据所述目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围,并确定纵波叠加速度和转换波叠加速度的平方比的取值范围;步骤C:将所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值,其中,M、N和K为自然数;步骤D:从所述研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据,并抽取所提取的转换波地震数据的渐近共转换点道集;步骤E:对于所述渐近共转换点道集中的每一个渐近共转换点道集执行以下步骤:步骤El:利用所述M个转换波叠加速度值、所述N个纵波叠加速度和转换波叠加速度的平方比值和所述K个转换波垂直双程反射传播时间值,计算该渐近共转换点道集中的每一个地震道的转换波反射旅行时间,并将该渐近共转换点道集中的每一个地震道的转换波反射旅行时间前后2L+1个采样点的数据采样值叠加后求平方并累加求和,以得到该渐近共转换点道集的扫描计算能量数据;步骤E2:对应转换波垂直双程反射传播时间的每个取值,确定所得到的扫描计算能量数据{Em,n,Jm=L..^n=L..,中的极大值的位置,提取该极大值的位置处所对应的取值作为所述转换波叠加速度以及所述纵波叠加速度和转换波叠加速度的平方比。在本专利技术的另一个方面中,提出了一种确定层状各向同性介质转换波叠加速度的装置,该装置包括:存储部,该存储部存储研究区域的纵波叠加速度以及转换波地震数据;取值范围确定部,该取值范围确定部根据所述存储部中存储的所述纵波叠加速度确定转换波叠加速度的取值范围,根据目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围,并且确定纵波叠加速度和转换波叠加速度的平方比的取值范围;离散部,该离散部将由所述取值范围确定部确定的所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取 值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值,其中,Μ、N和K为自然数;渐近共转换点道集抽取部,该渐近共转换点道集抽取部从所述存储部中存储的研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据,并抽取所提取的转换波地震数据的渐近共转换点道集;扫描部,该扫描部利用所述M个转换波叠加速度值、所述N个纵波叠加速度和转换波叠加速度的平方比值和所述K个转换波垂直双程反射传播时间值,计算由所述渐近共转换点道集抽取部抽取的所述渐近共转换点道集中的每一个渐近共转换点道集的各个地震道的转换波反射旅行时间,并将所述每一个渐近共转换点道集中的各个地震道的转换波反射旅行时间前后2L+1个采样点的数据采样值叠加后求平方并累加求和,以得到所述每一个渐近共转换点道集的扫描计算能量数据,其中,所述扫描部确定所得到的扫描计算能量数据中的极大值的位置,提取该极大值的位置处所对应的取值作为所述转换波叠加速度以及所述纵波叠加速度和转换波叠加速度的平方比。根据本专利技术,采用两个控制参数来描述层状各向同性介质中转换波反射旅行时间以进行转换波叠加速度分析,同时考虑了非对称路径对转换波反射旅行时间的影响,对近偏移距、中等偏移距转换波地震数据渐近共转换点(ACCP)道集进行双参数扫描,对扫描结果进行分析能够得到准确的转换波叠加速度模型。本专利技术对于油气勘探中转换波地震数据处理具有十分重要的应用价值。附图说明图1示出了本专利技术第一实施方式的确定层状各向同性介质转换波叠加速度的方法的流程图;图2示出了理论合成五层水平层状介质的转换波道集;图3示出了近炮检距和中炮检距范围的理论合成五层水平层状介质的转换波道集;图4示出了对第一 层Tc0=800ms的转换波双参数扫描结果;图5示出了对第二层Tc0=1426ms的转换波双参数扫描结果;图6示出了对第三层Tc0=1954ms的转换波双参数扫描结果;图7示出了对第四层Tc0=2431ms的转换波双参数扫描结果;图8示出了对第五层Tc0=2875ms的转换波双参数扫描结果;图9示出了本专利技术第二实施方式的确定层状各向同性介质转换波叠加速度的装置的框图。具体实施例方式为使本专利技术的目的、技术方案和优点更加清楚,下面结合附图对本专利技术的具体实施方式进行详细说明。专利技术人通过包含(ν ,Y is。)的双参数模型来对层状各向同性介质中转换波的反射波旅行时间进行量化描述如下:权利要求1.一种确定层状各向同性介质转换波叠加速度的方法,该方法包括以下步骤: 步骤A:获取所述研究区域的纵波叠加速度; 步骤B:根据所述纵波叠加速度确定所述转换波叠加速度的取值范围,根据目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围,并确定纵波叠加速度和转换波叠加速度的平方比的取值范围; 步骤C:将所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值,其中,M、N和K为自然数; 步骤D:从所述研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据,并抽取所提取的转换波地震数据的渐近共转换点道集; 步骤E:对本文档来自技高网...

【技术保护点】
一种确定层状各向同性介质转换波叠加速度的方法,该方法包括以下步骤:步骤A:获取所述研究区域的纵波叠加速度;步骤B:根据所述纵波叠加速度确定所述转换波叠加速度的取值范围,根据目标地层的深度和所述纵波叠加速度确定转换波垂直双程反射传播时间的取值范围,并确定纵波叠加速度和转换波叠加速度的平方比的取值范围;步骤C:将所述转换波叠加速度的取值范围、所述纵波叠加速度和转换波叠加速度的平方比的取值范围、以及所述转换波垂直双程反射传播时间的取值范围分别离散采样为M个转换波叠加速度值、N个纵波叠加速度和转换波叠加速度的平方比值和K个转换波垂直双程反射传播时间值,其中,M、N和K为自然数;步骤D:从所述研究区域的转换波地震数据中提取近炮检距和中炮检距的转换波地震数据,并抽取所提取的转换波地震数据的渐近共转换点道集;步骤E:对于所述研究区域的所有渐近共转换点道集中的每一个渐近共转换点道集执行以下步骤:步骤E1:利用所述M个转换波叠加速度值、所述N个纵波叠加速度和转换波叠加速度的平方比值和所述K个转换波垂直双程反射传播时间值,计算该渐近共转换点道集中的每一个地震道的转换波反射旅行时间,并将该渐近共转换点道集中的每一个地震道的转换波反射旅行时间前后2L+1个采样点的数据采样值叠加后求平方并累加求和,以得到该渐近共转换点道集的扫描计算能量数据;步骤E2:对应转换波垂直双程反射传播时间的每个取值,确定所得到的扫描计算能量数据中的极大值的位置,提取该极大值的位置处所对应的取值作为所述转换波叠加速度以及所述纵波叠加速度和转换波叠加速度的平方比。...

【技术特征摘要】

【专利技术属性】
技术研发人员:李向阳陈双全
申请(专利权)人:中国石油天然气集团公司中国石油大学北京
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1