CN104267434B - 一种三维多分量地震观测系统炮检距分布获取方法及装置 - Google Patents
一种三维多分量地震观测系统炮检距分布获取方法及装置 Download PDFInfo
- Publication number
- CN104267434B CN104267434B CN201410466484.9A CN201410466484A CN104267434B CN 104267434 B CN104267434 B CN 104267434B CN 201410466484 A CN201410466484 A CN 201410466484A CN 104267434 B CN104267434 B CN 104267434B
- Authority
- CN
- China
- Prior art keywords
- observation system
- geophone
- geophone offset
- dimensional multi
- bin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供一种三维多分量地震观测系统炮检距分布获取方法及装置,所述方法包括:获取三维多分量地震观测系统数据;针对某一深度的目标地层,按照面元位置,对三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;根据重新排序的三维多分量地震观测系统数据中的震源点与检波点位置信息,计算目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:根据目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。本发明可在地震观测系统设计阶段实现对三维多分量地震观测系统的炮检距均匀性效果的评价。
Description
技术领域
本发明涉及三维多分量地震观测系统分析评价技术,尤其涉及一种三维多分量地震观测系统炮检距分布获取方法及装置。
背景技术
三维地震勘探是石油与天然气勘探的主要工具。人们很早就认识到地震波是一种弹性波,这种弹性波包含纵波和横波信息,即多分量地震波信息。多分量地震勘探具有许多优势,比如可以得到更高分辨率的地震资料和更多的地下岩石物性参数等。但长期以来,由于技术和成本的限制,地震勘探实际上大多只利用纵波信息,多分量地震勘探仍然面临许多的困难。由于横波震源的施工难度极大,如今的多分量地震一般采用纵波激发,纵波和横波同时接收的勘探方式,其主要过程包括:(1)数据资料采集。在陆上地震勘探数据采集现场一般进行下述三个工作:地震观测系统设计、根据设计方案在野外布设震源和检波器,地震波的激发和接收。首先在室内进行地震观测系统设计,以确定震源点和检波点的最佳摆放位置。然后根据设计方案,在野外布设纵波震源、纵波检波器和横波检波器,陆上地震勘探的震源点一般采用炸药震源,并沿地震测线等间距布置多个检波器来接收地震波信号,现代地震勘探中检波器的数量多在1000或10000个。震源点在爆炸后产生出地震波,地震波遇岩层界面反射回来被检波器接收并传到仪器车,仪器车将检波器传来的信号记录下来,这就获得了用以研究地下油气埋藏情况的地震记录。(2)地震数据处理。地震数据处理是把第一步骤采集到的地震数据资料(包括纵波地震数据和横波地震数据)输入专用电子计算机,按不同要求用一系列功能不同的程序进行处理运算,把数据进行归类编排,突出有效的,除去无效和干扰的,最后把经过各种处理的数据进行叠加和偏移,最终分别得到二维或三维的纵波和横波地震数据体文件。(3)资料解释。资料解释是把经过处理的地震信息变成地质成果的过程,包括运用波动理论和地质知识,综合地质、钻井、测井等各项资料,作出构造解释、地层解释、岩性和烃类检测解释及综合解释,绘出有关成果图件,对勘探区域作出含油气评价,提出钻探井位置等。
为了摆脱人工判断三维地震观测系统采样均衡性所带来的不确定性,中国专利申请号200610114254.1提供了一种地震观测系统定量分析方法,通过定量的分析三维地震观测系统采样均衡性,使得判断结果更准确。中国专利申请号200610114254.1提供的定量分析方法中考虑了四个象限内炮检距分布的均衡性,没有从更加精细的网格尺度上来分析面元的属性,同时也没有综合考虑各个面元所有属性因素(包括覆盖次数、炮检距、方位角等)给出该面元的整体均衡度值(需分别依据五个参数进行判断),因此在应用过程中会产生判断结果不确定的问题。中国专利申请号201010569364.3提供了一种三维地震观测系统方案的整体均衡度定量分析方法,通过引入信息论中熵的概念分析三维地震观测系统整体采样的均衡性。但这些专利技术均基于传统的纵波勘探,未考虑三维多分量地震观测系统的特殊性,及其对炮检距分布均匀性的影响。
对于三维多分量地震观测系统,由于转换波的下行波是纵波,上行波是横波,故其射线路径是不对称的,进而造成转换波炮检距分布的不均匀。这种炮检距的不规则分布并非地震观测系统本身的分布特征造成,而是由于转换波的不对称特性所带来的假象,并可直接影响后续地震成像的效果。但一直以来缺乏适当的方法对三维多分量地震观测系统转换波炮检距分布的均匀性进行有效的定量评估。
发明内容
本发明实施例提供一种三维多分量地震观测系统炮检距分布获取方法及装置,以对三维多分量地震观测系统转换波炮检距分布的均匀性进行有效的定量评估提供了一种技术方案。
一方面,本发明实施例提供了一种三维多分量地震观测系统炮检距分布获取方法,所述三维多分量地震观测系统炮检距分布获取方法包括:
获取三维多分量地震观测系统数据;
针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;
根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。
另一方面,本发明实施例提供了一种三维多分量地震观测系统炮检距分布获取装置,所述三维多分量地震观测系统炮检距分布获取装置包括:
获取单元,用于获取三维多分量地震观测系统数据;
排序单元,用于针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;
非均匀性系数计算单元,用于根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
非均匀性系数分布图绘制单元,用于根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。
上述技术方案具有如下有益效果:因为采用所述三维多分量地震观测系统炮检距分布获取方法包括:获取三维多分量地震观测系统数据;针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图的技术手段,所以达到了如下的技术效果:克服了现有炮检距均匀性定量分析方法无法被应用于多分量地震观测系统转换波偏移噪声分析的缺陷,提出了一种三维多分量地震观测系统转换波炮检距均匀性的分析方法。利用本发明所述技术方案,可将观测系统炮检距均匀性定量分析算法推广至多分量地震勘探的转换波分析,在地震观测系统设计阶段实现对多分量三维地震观测系统的炮检距均匀性效果的评价,有力地促进了多分量地震勘探的实际应用。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种三维多分量地震观测系统炮检距分布获取方法流程图;
图2为本发明实施例一种三维多分量地震观测系统炮检距分布获取装置结构示意图;
图3为本发明实施例非均匀性系数计算单元结构示意图;
图4为本发明应用实例三维多分量地震观测系统转换波传播过程的路径示意图;
图5为本发明应用实例某三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,为本发明实施例一种三维多分量地震观测系统炮检距分布获取方法流程图,所述三维多分量地震观测系统炮检距分布获取方法包括:
101、获取三维多分量地震观测系统数据;
102、针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;
103、根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
104、根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。
优选的,所述三维多分量地震观测系统数据包括:震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系;其中,所述震源点的三维空间坐标包括:高程;所述检波点三维空间坐标包括:高程;所述面元为棋盘状水平网格。
优选的,所述根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,包括:根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距;根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数。
优选的,所述根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距,包括:
利用如下公式计算所述目标地层每一个面元位置处所对应的转换波炮检距:
其中,震源点为(xs,ys,zs),目标点为(xo,yo,zo),检波点为(xr,yr,zr)。
优选的,所述根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数,包括:
假设一个面元所对应的炮检距为R1,R2,……,Rn,则利用如下公式计算该面元炮检距分布的非均匀系数:
C代表了该面元处炮检距分布的非均匀性,C数值越小则炮检距分布越均匀。
如图2所示,为本发明实施例一种三维多分量地震观测系统炮检距分布获取装置结构示意图,所述三维多分量地震观测系统炮检距分布获取装置包括:
获取单元21,用于获取三维多分量地震观测系统数据;
排序单元22,用于针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;
非均匀性系数计算单元23,用于根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
非均匀性系数分布图绘制单元24,用于根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图。
优选的,所述三维多分量地震观测系统数据包括:震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系;其中,所述震源点的三维空间坐标包括:高程;所述检波点三维空间坐标包括:高程;所述面元为棋盘状水平网格。
优选的,如图3所示,为本发明实施例非均匀性系数计算单元结构示意图,非均匀性系数计算单元23包括:
转换波炮检距计算模块231,用于根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距;
非均匀性系数计算模块232,用于根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数。
优选的,所述转换波炮检距计算模块231,进一步具体用于利用如下公式计算所述目标地层每一个面元位置处所对应的转换波炮检距:
其中,震源点为(xs,ys,zs),目标点为(xo,yo,zo),检波点为(xr,yr,zr)。
优选的,所述非均匀性系数计算模块232,进一步具体用于假设一个面元所对应的炮检距为R1,R2,……,Rn,则利用如下公式计算该面元炮检距分布的非均匀系数:
C代表了该面元处炮检距分布的非均匀性,C数值越小则炮检距分布越均匀。
本发明实施例上述技术方案具有如下有益效果:因为采用所述三维多分量地震观测系统炮检距分布获取方法包括:获取三维多分量地震观测系统数据;针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序;根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图的技术手段,所以达到了如下的技术效果:克服了现有炮检距均匀性定量分析方法无法被应用于多分量地震观测系统转换波偏移噪声分析的缺陷,提出了一种三维多分量地震观测系统转换波炮检距均匀性的分析方法。利用本发明所述技术方案,可将观测系统炮检距均匀性定量分析算法推广至多分量地震勘探的转换波分析,在地震观测系统设计阶段实现对多分量三维地震观测系统的炮检距均匀性效果的评价,有力地促进了多分量地震勘探的实际应用。
以下结合应用实例对本发明实施例进行详细说明:
如图4所示,为本发明应用实例三维多分量地震观测系统转换波传播过程的路径示意图。由于转换波的下行波是纵波,上行波是横波,故其射线路径是不对称的。其反射点位于x/(1+vs/vp)处(x为震源点到检波点的距离,vs和vp分别为横波和纵波的速度),与纵波的共反射点位于x/2处有较大的差别,且对于不同深度、不同速度比的地层,其反射点位置各不相同。总体而言,反射点比较靠近接收点的一边,且随着地层由深至浅更加靠近接收点。因此,对于转换波地震观测系统设计而言,应该首先确定目标地层的深度,再以此深度为基础进行后续分析工作。
依据地下地震速度模型与三维地震观测系统设计方案,直接计算三维地震观测系统的预期偏移噪声。由于采用离散退化方法求解波动方程,该方法对于地下强对比度介质条件与广角成像范围也具有很高的计算精度。由于输入地震速度模型为三维网格模型,通过合理设计网格间距,该方法可适用于任意复杂的介质模型。本发明应用实例的技术方案如下:
1)计算多分量地震转换波的射线路径
地震波总是从震源点(xs,ys,zs)开始下行,在目标点(xo,yo,zo)发生等角度反射,然后上行到达检波点(xr,yr,zr)。因此,可以针对某一深度的目标地层,根据入射波与反射波角度相同的原理得到目标点(xo,yo,zo)的位置,进而计算多分量地震转换波的射线路径。
2)计算多分量地震转换波的炮检距
将震源点(xs,ys,zs)与目标点(xo,yo,zo)的距离与目标点(xo,yo,zo)和检波点(xr,yr,zr)的距离的和作为多分量地震转换波的炮检距
其中,震源点为(xs,ys,zs),目标点为(xo,yo,zo),检波点为(xr,yr,zr)。
3)计算一个面元所对应的所有炮检对(一个震源点和其对应的一个检波点)。面元为棋盘状水平网格,比如200*200个,每个网格大小可为12.5米*12.5米或25米*25米等。当震源点和检波点的水平坐标中点位于面元内时,我们认为该震源点和检波点对信息属于该面元。计算同一个面元内所有炮检对的炮检距差的倒数和。假设一个面元所对应的炮检距为R1,R2,……,Rn,该面元炮检距非均匀系数为:
C代表了该面元处炮检距分布的非均匀性,C数值越小则炮检距分布越均匀。
1)获取三维多分量地震观测系统数据
三维多分量地震观测系统数据包括检波点的三维空间坐标(含高程)、震源点三维空间坐标(含高程),和震源点与检波点之间的对应关系。一个地震观测系统设计方案包含大量的震源点和检波点的三维坐标信息。这里的震源点指地震勘探普遍采用的人工炸药震源。检波点指地震勘探中在地表布设的检波器,一般至少有上千个,用于接收震源发出经地下反射后回到地表的地震波。
2)针对某一深度的目标地层,按照面元位置,对多分量观测系统数据的震源点与检波点位置信息进行重新排序。
面元为棋盘状水平网格,比如200*200个,每个网格大小可为12.5米*12.5米或25米*25米等。炮检对为地震观测系统的最小数据单位,包含震源点三维空间坐标信息和其所对应的一个检波点三维空间坐标信息。每个炮检对需要用6个参数进行描述:震源点横坐标xs、震源点纵坐标ys、震源点高程zs、检波点横坐标xr、检波点纵坐标yr和检波点高程zr。如果一个炮检对所对应的目标点(xo,yo,zo)位于某个面元位置,则认为该炮检对的信息属于该面元位置。值得注意的是,由于转换波的下行波是纵波,上行波是横波,其射线路径不对称,故目标点的水平位置一般不位于炮检对的中点。重新排序前,观测系统数据一般按照震源点位置排序,即按照震源点的信息和其所对应的所有检波点的信息排序。基于面元位置的重新排列过程,将每个面元位置所对应的观测系统信息整理在一起,为下一步的基于单个面元的分析计算进行数据的准备。
3)计算一个面元所对应的转换波炮检距分布非均匀性系数:
为了剥离地下介质因素的影响,独立的分析复杂地表对面元振幅能量均衡性的影响,我们对地下介质情况进行简化,假设其为均匀介质,即介质速度为一恒定值,而且仅考虑水平目标层位情况。
地震转换波传播示意图如图4所示。利用公式(1)计算一个面元所对应所有转换波的炮检距,然后利用公式(2)计算该面元处转换波炮检距分布的非均匀性。
4)分别对某一目的层位的每个面元,分别计算各个面元位置处转换波炮检距分布的非均匀性系数,并绘制转换波炮检距分布的非均匀性系数分布图:
如图5所示,为本发明应用实例某三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图,横纵坐标分别表示X和Y坐标,单位米。图5中的颜色深度值代表每个面元的波炮检距分布的非均匀性系数。每个面元的值越小,该面元的炮检距均匀性越好。该技术可以实现多分量地震勘探的炮检距分布均匀性的定量评估,克服了以往定性方法的缺点。应用本发明可对转换波炮检距分布进行定量评价,得到的转换波炮检距分布图。利用该技术可实现对三维多分量地震观测系统转换波炮检距分布的定量分析。
本发明实施例及应用实例从地震勘探采集的源头实现了对转换波所引起的地震波炮检距均匀性的定量评估,对提高多分量地震采集资料保真度提供了保障,进而为多分量地震偏移成像、油气储层预测、油藏描述的可靠性奠定了基础,具有重要应用价值。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrative logical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrative components),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种三维多分量地震观测系统炮检距分布获取方法,其特征在于,所述三维多分量地震观测系统炮检距分布获取方法包括:
获取三维多分量地震观测系统数据;
针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序,当所述震源点和检波点的水平坐标中点位于所述面元内时,则所述震源点和检波点对信息属于所述面元;
根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图;
所述根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距,包括:
利用如下公式计算所述目标地层每一个面元位置处所对应的转换波炮检距:
其中,震源点为(xs,ys,zs),目标点为(xo,yo,zo),检波点为(xr,yr,zr);
所述根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数,包括:
假设一个面元所对应的炮检距为R1,R2,……,Rn,则利用如下公式计算该面元炮检距分布的非均匀系数:
C代表了该面元处炮检距分布的非均匀性,C数值越小则炮检距分布越均匀。
2.如权利要求1所述三维多分量地震观测系统炮检距分布获取方法,其特征在于,所述三维多分量地震观测系统数据包括:震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系;其中,所述震源点的三维空间坐标包括:高程;所述检波点三维空间坐标包括:高程;所述面元为棋盘状水平网格。
3.如权利要求1所述三维多分量地震观测系统炮检距分布获取方法,其特征在于,所述根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,包括:
根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距;
根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数。
4.一种三维多分量地震观测系统炮检距分布获取装置,其特征在于,所述三维多分量地震观测系统炮检距分布获取装置包括:
获取单元,用于获取三维多分量地震观测系统数据;
排序单元,用于针对某一深度的目标地层,按照面元位置,对所述三维多分量地震观测系统数据中的震源点与检波点位置信息进行重新排序,当所述震源点和检波点的水平坐标中点位于所述面元内时,则所述震源点和检波点对信息属于所述面元;
非均匀性系数计算单元,用于根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数:
非均匀性系数分布图绘制单元,用于根据所述目标地层每一个面元位置处所对应的转换波炮检距分布的非均匀性系数,获取所述三维多分量地震观测系统的转换波炮检距分布的非均匀性系数分布图;
所述转换波炮检距计算模块,进一步具体用于利用如下公式计算所述目标地层每一个面元位置处所对应的转换波炮检距:
其中,震源点为(xs,ys,zs),目标点为(xo,yo,zo),检波点为(xr,yr,zr);
所述非均匀性系数计算模块,进一步具体用于假设一个面元所对应的炮检距为R1,R2,……,Rn,则利用如下公式计算该面元炮检距分布的非均匀系数:
C代表了该面元处炮检距分布的非均匀性,C数值越小则炮检距分布越均匀。
5.如权利要求4所述三维多分量地震观测系统炮检距分布获取装置,其特征在于,所述三维多分量地震观测系统数据包括:震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系;其中,所述震源点的三维空间坐标包括:高程;所述检波点三维空间坐标包括:高程;所述面元为棋盘状水平网格。
6.如权利要求4所述三维多分量地震观测系统炮检距分布获取装置,其特征在于,所述非均匀性系数计算单元包括:
转换波炮检距计算模块,用于根据重新排序的所述三维多分量地震观测系统数据中的震源点与检波点位置信息,计算所述目标地层每一个面元位置处所对应的转换波炮检距;
非均匀性系数计算模块,用于根据所述目标地层每一个面元位置处所对应的转换波炮检距,计算所述转换波炮检距分布的非均匀性系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466484.9A CN104267434B (zh) | 2014-09-12 | 2014-09-12 | 一种三维多分量地震观测系统炮检距分布获取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410466484.9A CN104267434B (zh) | 2014-09-12 | 2014-09-12 | 一种三维多分量地震观测系统炮检距分布获取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104267434A CN104267434A (zh) | 2015-01-07 |
CN104267434B true CN104267434B (zh) | 2017-01-11 |
Family
ID=52158973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410466484.9A Expired - Fee Related CN104267434B (zh) | 2014-09-12 | 2014-09-12 | 一种三维多分量地震观测系统炮检距分布获取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104267434B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109738944B (zh) * | 2019-03-05 | 2020-05-08 | 中国石油大学(北京) | 基于广角反射的地震采集参数确定方法及装置 |
CN112444843B (zh) * | 2019-08-30 | 2024-06-18 | 中国石油化工股份有限公司 | 一种观测系统的调整方法 |
CN110579799B (zh) * | 2019-09-04 | 2020-10-27 | 中国科学院地质与地球物理研究所 | 一种等旅行时间间隔的地震采集观测方法及系统 |
CN112748458B (zh) * | 2019-10-30 | 2024-05-28 | 中国石油天然气集团有限公司 | 用于现场精细变观的在线实时面元获取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2378507B (en) * | 2001-06-08 | 2005-08-17 | Pgs Americas Inc | Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN102928878A (zh) * | 2012-11-05 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 三维地震观测系统振幅均衡性定量评估方法 |
-
2014
- 2014-09-12 CN CN201410466484.9A patent/CN104267434B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2378507B (en) * | 2001-06-08 | 2005-08-17 | Pgs Americas Inc | Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion |
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN102928878A (zh) * | 2012-11-05 | 2013-02-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 三维地震观测系统振幅均衡性定量评估方法 |
Non-Patent Citations (2)
Title |
---|
三维观测系统属性分析研究;黄迎雪,等;《工程地球物理学报》;20120131;第9卷(第1期);25-28 * |
炮检距属性的非均匀性系数分析;赵虎,等;《石油地球物理勘探》;20110228;第46卷(第1期);22-27 * |
Also Published As
Publication number | Publication date |
---|---|
CN104267434A (zh) | 2015-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bortoli et al. | Constraining stochastic images to seismic data: Stochastic simulation of synthetic seismograms | |
CN106094032B (zh) | 一种构建地层速度模型的方法 | |
Boullenger et al. | Studying CO 2 storage with ambient-noise seismic interferometry: A combined numerical feasibility study and field-data example for Ketzin, Germany | |
CN104793247B (zh) | 近地表速度建模方法及建模装置 | |
CN103645503B (zh) | 一种三维时间域照明分析及振幅补偿方法 | |
CN104155701B (zh) | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 | |
CN102262240B (zh) | 三维观测系统满覆盖区域自动布设方法 | |
CN104267434B (zh) | 一种三维多分量地震观测系统炮检距分布获取方法及装置 | |
CN106405644B (zh) | 裂缝确定方法和装置 | |
Janson et al. | 3-D forward seismic model of an outcrop-based geocellular model | |
CN108957548A (zh) | 一种多波多分量联合观测地震页岩气富集区预测技术 | |
CN105137479B (zh) | 一种面元覆盖次数的计算方法及装置 | |
CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
CN103645460B (zh) | 基于目标层可视性分析的激发位置确定方法 | |
US6147929A (en) | Method for predicting continuous and discontinuous waveguide targets using interwell seismic signature characteristics | |
US20220236435A1 (en) | Low-Frequency Seismic Survey Design | |
US20100118655A1 (en) | Progressive 3D vertical seismic profiling method | |
Karpekin et al. | Borehole acoustic reflection survey in horizontal wells: High resolution reservoir structure to guide properties distribution | |
Zhang | Assessing Attenuation, Fractures, and Anisotropy Using Logs, Vertical Seismic Profile, and Three-component Seismic Data: Heavy Oilfield and Potach Mining Examples | |
CN109991663A (zh) | 工区地震速度体校正方法和装置 | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
CN113614578A (zh) | 使用双向虚拟源重建的地震勘测 | |
Dutta et al. | A novel approach to fracture characterization utilizing borehole seismic data | |
CN112180443A (zh) | 页岩气二维地震甜点区优选方法及装置 | |
Botter et al. | Seismic attribute analysis of a fault zone in the Thebe field, Northwest shelf, Australia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170111 Termination date: 20190912 |
|
CF01 | Termination of patent right due to non-payment of annual fee |