CN109100792A - 基于台站与三维地震联合采集资料的速度反演方法 - Google Patents
基于台站与三维地震联合采集资料的速度反演方法 Download PDFInfo
- Publication number
- CN109100792A CN109100792A CN201811292083.0A CN201811292083A CN109100792A CN 109100792 A CN109100792 A CN 109100792A CN 201811292083 A CN201811292083 A CN 201811292083A CN 109100792 A CN109100792 A CN 109100792A
- Authority
- CN
- China
- Prior art keywords
- station
- data
- seismics
- wave
- velocity
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000004587 chromatography analysis Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 6
- 238000009933 burial Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 abstract 1
- 239000009671 shengli Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于台站与三维地震联合采集资料的速度反演方法,该基于台站与三维地震联合采集资料的速度反演方法包括台站布设位置设计、台站与三维地震数据联合采集、台站与三维地震观测系统坐标匹配、台站数据截取、台站与三维地震数据匹配、台站初至波走时拾取、三维地震反射波走时拾取、初始速度模型建立、初至波与反射波走时联合反演、速度模型输出。该基于台站与三维地震联合采集资料的速度反演方法可以减少深层速度反演的多解性,提高深层速度反演的精度,改善深层地震资料的成像质量。
Description
技术领域
本发明涉及油气地震勘探领域,特别是涉及到一种基于台站与三维地震联合采集资料的速度反演方法。
背景技术
自从在亚特兰大召开的第54届地球物理勘探学家协会年会上设置了地震速度成像研究内容的专题之后,经过以Daily,Somerstein,Bishop,Dyer和Worthington等人为代表的大量研究,地震波速度反演的理论、方法和技术以数值模拟得以迅速的发展。就多种波旅行时联合速度反演这一重要分支而言,其发展现状如下:1992年,Zelt提出了一种基于任意类型的地震体波的同时确定二维地壳的速度及界面形态的地震旅行时反演方法。1995年,华标龙推导出了二维层状介质中适用于反射波和透射波的地震波旅行时对界面偏导数的公式,并将其用于二维层状介质中界面和速度的联合反演。1997年,McCaughey利用透射波、广角反射以及折射旅行时来同时反演速度和界面深度,该方法通过对系数矩阵内不同量纲的值加以权衡值,并在反演中加入正则化来稳定反演的解。1998年,Zhang Jie提出了一种非线性折射波和反射波旅行时反演方法,该方法利用折射波和反射波旅行时同时进行速度反演。2003年,James利用透射、广角反射、折射以及多次波旅行时对三维速度及界面进行同时速度,在反演过程中使数据误差和模型粗糙度趋向于最小来寻找一个合适的层速度—界面模型。2003年,Zhou提出了基于反演界面形态的可变形层析成像方法。由于绝大多数沉积盆地都具有层状速度结构,解决当前基于网格和节点参数化的层析方法难以高精度重构层状速度结构的问题。2006年,Zhou在模型参数化时,又引入多尺度参数化,发展出初至多尺度可变形层析成像方法和反射尺度可变形层析成像方法,与单尺度层析相比,多尺度层析反演结果更加稳定可靠。2006年,周龙泉依据波逆行原理推导了三维层状介质中地震波旅行时对界面偏导数的公式,得到了复杂界面情况下,反演时所需的旅行时对界面偏导数系数矩阵,通过模型试算,证实了该方法对三维速度结构和界面重构是有效的。2011年,白超英等人将不同震相种类数据的权系数以及不同种类参数的归一化因子引入到反演算法之中,实现了三维复杂层状介质中多震相走时联合反演成像。
关于多震相联合反演问题的研究由来已久,但是面对的都是盆地规模的宏观速度反演问题,而且现有的技术都是基于同一种观测系统和同一种检波器,解决的是速度反演精度的问题。如何利用好初至波和反射波各自优点提高油气田深层速度模型反演的精度,未见公开发表的技术成果,这是因为受限于经济成本和施工可行性,现有的油气地震勘探观测系统的最大偏移距受到了限制,无法接收到来自深层的地震波信息。
胜利油田近几年深层油气储量显著增长,石油与天然气地质储量丰富,但整体勘探程度较低,存在构造落实难、有效储层描述困难等问题。二次采集的高精度三维资料虽然大幅度提高了深层地震资料的成像质量,但对于地质和油藏描述要求,还存在一定的多解性,主要矛盾就是深层速度模型精度的问题。为此我们发明了一种新的基于台站与三维地震联合采集资料的速度反演方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种能够有效提高深层速度建模精度的基于台站与三维地震联合采集资料的速度反演方法。
本发明的目的可通过如下技术措施来实现:基于台站与三维地震联合采集资料的速度反演方法,该基于台站与三维地震联合采集资料的速度反演方法包括:步骤1:根据三维地震观测系统,结合实际地表情况,在三维地震施工前布设好台站,进行连续记录;步骤2:对台站和三维地震坐标系统进行匹配,使两种观测系统的坐标统一;步骤3:从连续记录的台站数据中,截取放炮时间段内的数据;步骤4:将截取的台站记录数据与三维地震记录的数据进行匹配,使两种数据做到时间上的同步;步骤5:在匹配后的台站数据上进行初至波到时拾取;步骤6:在三维地震数据上进行反射波到时拾取;步骤7:建立地下初始速度模型;步骤8:利用拾取的台站初至波走时和三维地震反射波走时进行多尺度联合层析反演;步骤9:当反演结果收敛时,将结果输出作为最终反演的速度模型。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,台站的布设采用的是平行交错式图形。
在步骤1中,采用的台站可连续记录至少4天,采集的数据自动存储在自身携带的存储卡上,更换备用电池后可继续工作。
在步骤1中,台站采用挖坑掩埋的方法,保证台站距地表最少60cm。
在步骤2中,两种观测系统的坐标匹配方法是利用台站与检波器重合点的检波器坐标求取未知的台站坐标,通过求解如下方程组求取与检波器坐标一致的台站坐标:
其中,(x,y)为待求台站的坐标,(x1,y1)、(x2,y2)、(x3,y3)为待求台站附近三个台站的坐标,d1、d2、d3分别为待求台站到附近三个台站的距离,该距离通过手持GPS设备测量的坐标计算得到。
在步骤3中,从连续记录的台站数据中截取放炮时间段内的数据,首先抽样统计台站计时系统与三维地震计时系统的大概时差范围,然后以三维地震记录的每一炮的野外放炮时间为中心,前后截取时差范围内的数据输出。
在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配前,需要进行两种记录波形的匹配,使两种记录的波形尽量保持一致。
在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配前,需要对波形匹配后的两种记录进行滤波和归一化处理。
在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配,是通过计算对应的两种记录互相关叠加值,通过最大互相关叠加值确定台站数据上准确的放炮时刻,并按照设定的数据长度截取输出台站数据。
在步骤5中,台站数据上进行初至波到时拾取是在台站数据的共检波点数据集内完成的,即将共炮点台站数据抽取成共检波点台站数据,并按照偏移距远近重新排列,然后利用长短时窗能量比方法自动拾取初至。
在步骤6中,在三维地震数据上进行反射波到时拾取是在共炮点数据集完成的,拾取时选择信噪比较高的反射同相轴作为标准反射层,通过自动和人工的方式进行反射波到时拾取。
在步骤7中,初始速度模型是利用台站数据的初至波走时进行多尺度网格层析反演得到的层状速度模型,其过程是以不同大小的网格离散化实际地质模型,同时对各个尺度模型进行反演。
在步骤8中,台站初至波走时和三维地震反射波走时联合层析反演,是将地下速度结构离散为层状模型,反演速度界面的起伏和层内的速度值。
在步骤8中,台站初至波走时和三维地震反射波走时联合反演,是通过计算界面多尺度分解后各个尺度上的节点位置变化对地震波到时的扰动量,同时使用初至波和反射波走时反演速度模型。
在步骤8中,台站初至波走时和三维地震反射波走时联合反演,是以不同大小的网格离散化实际地质模型,在反演过程中同时对各个尺度模型进行反演。
在步骤8中,台站初至波走时和三维地震反射波走时联合层析反演方程中,同时使用了初至波和反射波到时进行速度模型的反演。
本发明中的基于台站与三维地震联合采集资料的速度反演方法,改变了常规三维地震采用单一观测系统的方法,在基本不增加采集成本的条件下,通过布设稀疏台站获得超大偏移距初至波信息,结合三维地震反射信息进行初至波与反射波联合层析反演,能够有效提高深层速度建模精度。利用该方法可以减少深层速度反演的多解性,提高深层速度反演的精度,改善深层地震资料的成像质量。
附图说明
图1是胜利油田工区A理论层状速度模型及射线追踪结果的示意图;
图2是胜利油田工区A初至波反演结果的示意图;
图3是胜利油田工区A反射波反演结果的示意图;
图4是胜利油田工区A初至波与反射波联合反演结果的示意图;
图5是胜利油田工区B稀疏台站位置(三角形)、三维地震观测系统图;
图6是胜利油田工区B截取的一段台站数据(共炮点数据集)的示意图;
图7是胜利油田工区B匹配完成后的共台站数据集按偏移距显示的示意图;
图8是胜利油田工区B某一条测线初至波与反射波多尺度联合层析反演结果的示意图;
图9为本发明的基于台站与三维地震联合采集资料的速度反演方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图9所示,图9为本发明的基于台站与三维地震联合采集资料的速度反演方法的流程图。
步骤101:根据三维地震观测系统,结合实际地表情况,在三维地震施工前布设好台站,进行连续记录;台站的布设采用的是平行交错式图形。采用的台站可连续记录至少4天,采集的数据可自动存储在自身携带的存储卡上,更换备用电池后可继续工作。台站采用挖坑掩埋的方法,保证台站距地表最少60cm。
步骤102:对台站和三维地震坐标系统进行匹配,使两种观测系统的坐标统一;两种观测系统的坐标匹配方法是利用台站与检波器重合点的检波器坐标求取未知的台站坐标,可通过求解如下方程组求取与检波器坐标一致的台站坐标:
其中,(x,y)为待求台站的坐标,(x1,y1)、(x2,y2)、(x3,y3)为待求台站附近三个台站的坐标,d1、d2、d3分别为待求台站到附近三个台站的距离,该距离通过手持GPS设备测量的坐标计算得到。
步骤103:从连续记录的台站数据中,截取放炮时间段内的数据;从连续记录的台站数据中截取放炮时间段内的数据,首先抽样统计台站计时系统与三维地震计时系统的大概时差范围,然后以三维地震记录的每一炮的野外放炮时间为中心,前后截取时差范围内的数据输出。
步骤104:将截取的台站记录数据与三维地震记录的数据进行匹配,使两种数据做到时间上的同步;将截取的台站记录数据与三维地震记录的数据进行匹配前,需要进行两种记录波形的匹配,使两种记录的波形尽量保持一致。将截取的台站记录数据与三维地震记录的数据进行匹配前,需要对波形匹配后的两种记录进行滤波和归一化处理。将截取的台站记录数据与三维地震记录的数据进行匹配,是通过计算对应的两种记录互相关叠加值,通过最大互相关叠加值确定台站数据上准确的放炮时刻,并按照设定的数据长度截取输出台站数据。
步骤105:在匹配后的台站数据上进行初至波到时拾取;台站数据上进行初至波到时拾取是在台站数据的共检波点数据集内完成的,即将共炮点台站数据抽取成共检波点台站数据,并按照偏移距远近重新排列,然后利用长短时窗能量比方法自动拾取初至。
步骤106:在三维地震数据上进行反射波到时拾取;在三维地震数据上进行反射波到时拾取是在共炮点数据集完成的,拾取时选择信噪比较高的反射同相轴作为标准反射层,通过自动和人工的方式进行反射波到时拾取。
步骤107:利用初至波多尺度层析反演结果建立地下层状初始速度模型;使用台站数据的初至波走时,进行多尺度网格层析反演,利用反演结果建立层状初始速度模型。
步骤108:利用拾取的台站初至波走时和三维地震反射波走时进行联合层析反演;台站初至波走时和三维地震反射波走时联合层析反演,是将地下速度结构离散为层状模型,反演速度界面的起伏和层内的速度值。台站初至波走时和三维地震反射波走时联合层析反演中,通过计算界面多尺度分解后各个尺度上的节点位置变化对地震波到时的扰动量,同时使用初至波和反射波走时反演速度模型。台站初至波走时和三维地震反射波走时联合层析反演方程中,以不同大小(尺度)的网格离散化实际地质模型,在反演过程中同时对各个尺度模型进行反演。
步骤109:当反演结果收敛时,将结果输出作为最终反演的速度模型。
以下为应用本发明的具体实施例,实施例1来源于中石化胜利油田一个理论速度模型的反演结果,通过该实施例可验证方法的有效性,具体实施方式为:
(1)根据工区地震资料解释成果,建立二维层状速度模型,如图1所示;
(2)利用射线追踪计算得到初至波与反射波走时,图1中的灰色折线为初至波射线追踪结果,图1中黑色折线为反射波射线追踪结果;
(3)只利用初至波走时反演速度模型,结果如图2,结果与真实速度模型相差较大;
(4)只利用反射波走时反演速度模型,结果如图3,结果优于单独利用初至波反演;
(5)利用初至波与反射波走时进行联合层析反演,结果如图4,联合反演结果明显优于单独的初至波或反射波反演结果。
实施例2来源于中石化胜利油田一个台站与三维地震联合采集项目,具体实施方式为:
(1)根据三维地震观测系统,结合实际地表情况,在三维地震施工前布设好台站,进行连续记录,台站位置与三维地震观测系统如图5所示。
(2)对台站和三维地震坐标系统进行匹配,使两种观测系统的坐标统一。
(3)从连续记录的台站数据中,截取放炮时间段内的数据,图6为截取的某一炮的台站数据。
(4)将截取的台站记录数据与三维地震记录的数据进行匹配,使两种数据做到时间上的同步,图7为匹配完成后的某一台站的共台站数据集按偏移距显示结果。
(5)在匹配后的台站数据上进行初至波到时拾取。
(6)在三维地震数据上进行反射波到时拾取。
(7)利用初至波多尺度网格层析反演结果建立地下层状初始速度模型。
(8)利用拾取的台站初至波走时和三维地震反射波走时进行联合层析反演,图8为台站初至走时与三维地震反射走时联合层析反演结果。
(9)当反演结果收敛时,将结果输出作为最终反演的速度模型。
Claims (16)
1.基于台站与三维地震联合采集资料的速度反演方法,其特征在于,该基于台站与三维地震联合采集资料的速度反演方法包括:
步骤1:根据三维地震观测系统,结合实际地表情况,在三维地震施工前布设好台站,进行连续记录;
步骤2:对台站和三维地震坐标系统进行匹配,使两种观测系统的坐标统一;
步骤3:从连续记录的台站数据中,截取放炮时间段内的数据;
步骤4:将截取的台站记录数据与三维地震记录的数据进行匹配,使两种数据做到时间上的同步;
步骤5:在匹配后的台站数据上进行初至波到时拾取;
步骤6:在三维地震数据上进行反射波到时拾取;
步骤7:建立地下初始速度模型;
步骤8:利用拾取的台站初至波走时和三维地震反射波走时进行多尺度联合层析反演;
步骤9:当反演结果收敛时,将结果输出作为最终反演的速度模型。
2.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤1中,台站的布设采用的是平行交错式图形。
3.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤1中,采用的台站可连续记录至少4天,采集的数据自动存储在自身携带的存储卡上,更换备用电池后可继续工作。
4.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤1中,台站采用挖坑掩埋的方法,保证台站距地表最少60cm。
5.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤2中,两种观测系统的坐标匹配方法是利用台站与检波器重合点的检波器坐标求取未知的台站坐标,通过求解如下方程组求取与检波器坐标一致的台站坐标:
其中,(x,y)为待求台站的坐标,(x1,y1)、(x2,y2)、(x3,y3)为待求台站附近三个台站的坐标,d1、d2、d3分别为待求台站到附近三个台站的距离,该距离通过手持GPS设备测量的坐标计算得到。
6.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤3中,从连续记录的台站数据中截取放炮时间段内的数据,首先抽样统计台站计时系统与三维地震计时系统的大概时差范围,然后以三维地震记录的每一炮的野外放炮时间为中心,前后截取时差范围内的数据输出。
7.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配前,需要进行两种记录波形的匹配,使两种记录的波形尽量保持一致。
8.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配前,需要对波形匹配后的两种记录进行滤波和归一化处理。
9.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤4中,将截取的台站记录数据与三维地震记录的数据进行匹配,是通过计算对应的两种记录互相关叠加值,通过最大互相关叠加值确定台站数据上准确的放炮时刻,并按照设定的数据长度截取输出台站数据。
10.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤5中,台站数据上进行初至波到时拾取是在台站数据的共检波点数据集内完成的,即将共炮点台站数据抽取成共检波点台站数据,并按照偏移距远近重新排列,然后利用长短时窗能量比方法自动拾取初至。
11.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤6中,在三维地震数据上进行反射波到时拾取是在共炮点数据集完成的,拾取时选择信噪比较高的反射同相轴作为标准反射层,通过自动和人工的方式进行反射波到时拾取。
12.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤7中,初始速度模型是利用台站数据的初至波走时进行多尺度网格层析反演得到的层状速度模型,其过程是以不同大小的网格离散化实际地质模型,同时对各个尺度模型进行反演。
13.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤8中,台站初至波走时和三维地震反射波走时联合层析反演,是将地下速度结构离散为层状模型,反演速度界面的起伏和层内的速度值。
14.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤8中,台站初至波走时和三维地震反射波走时联合反演,是通过计算界面多尺度分解后各个尺度上的节点位置变化对地震波到时的扰动量,同时使用初至波和反射波走时反演速度模型。
15.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤8中,台站初至波走时和三维地震反射波走时联合反演,是以不同大小的网格离散化实际地质模型,在反演过程中同时对各个尺度模型进行反演。
16.根据权利要求1所述的基于台站与三维地震联合采集资料的速度反演方法,其特征在于,在步骤8中,台站初至波走时和三维地震反射波走时联合层析反演方程中,同时使用了初至波和反射波到时进行速度模型的反演。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811292083.0A CN109100792B (zh) | 2018-10-31 | 2018-10-31 | 基于台站与三维地震联合采集资料的速度反演方法 |
PCT/CN2019/000145 WO2020087767A1 (zh) | 2018-10-31 | 2019-07-18 | 基于台站与三维地震联合采集资料的速度反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811292083.0A CN109100792B (zh) | 2018-10-31 | 2018-10-31 | 基于台站与三维地震联合采集资料的速度反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100792A true CN109100792A (zh) | 2018-12-28 |
CN109100792B CN109100792B (zh) | 2021-06-22 |
Family
ID=64870029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811292083.0A Active CN109100792B (zh) | 2018-10-31 | 2018-10-31 | 基于台站与三维地震联合采集资料的速度反演方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109100792B (zh) |
WO (1) | WO2020087767A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020087767A1 (zh) * | 2018-10-31 | 2020-05-07 | 中国石油化工股份有限公司 | 基于台站与三维地震联合采集资料的速度反演方法 |
CN112036008A (zh) * | 2020-07-28 | 2020-12-04 | 中国石油天然气集团有限公司 | 一种利用混合模型记录反演速度的方法及系统 |
CN112824939A (zh) * | 2019-11-20 | 2021-05-21 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的剩余静校正方法、装置、设备及可读存储介质 |
CN113917524A (zh) * | 2020-07-10 | 2022-01-11 | 中国石油化工股份有限公司 | 基于无缆节点台站的近地表精细探测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
CN102053263A (zh) * | 2009-10-30 | 2011-05-11 | 中国石油化工集团公司 | 调查表层结构的方法 |
CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
CN102841373A (zh) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于方位角约束的微地震事件定位方法 |
CN102937721A (zh) * | 2012-11-07 | 2013-02-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用初至波走时的有限频层析成像方法 |
US8406081B2 (en) * | 2009-09-25 | 2013-03-26 | Landmark Graphics Corporation | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers |
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
CN107843922A (zh) * | 2017-12-25 | 2018-03-27 | 中国海洋大学 | 一种基于地震初至波和反射波走时联合的层析成像方法 |
CN108375789A (zh) * | 2018-02-11 | 2018-08-07 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN110261899A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 地震数据z字形干扰波去除方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109100792B (zh) * | 2018-10-31 | 2021-06-22 | 中国石油化工股份有限公司 | 基于台站与三维地震联合采集资料的速度反演方法 |
-
2018
- 2018-10-31 CN CN201811292083.0A patent/CN109100792B/zh active Active
-
2019
- 2019-07-18 WO PCT/CN2019/000145 patent/WO2020087767A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
US8406081B2 (en) * | 2009-09-25 | 2013-03-26 | Landmark Graphics Corporation | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers |
CN102053263A (zh) * | 2009-10-30 | 2011-05-11 | 中国石油化工集团公司 | 调查表层结构的方法 |
CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
CN102841373A (zh) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于方位角约束的微地震事件定位方法 |
CN102937721A (zh) * | 2012-11-07 | 2013-02-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用初至波走时的有限频层析成像方法 |
CN104391323A (zh) * | 2014-11-21 | 2015-03-04 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
CN107843922A (zh) * | 2017-12-25 | 2018-03-27 | 中国海洋大学 | 一种基于地震初至波和反射波走时联合的层析成像方法 |
CN108375789A (zh) * | 2018-02-11 | 2018-08-07 | 中国石油化工股份有限公司 | 联合采集地震数据的同步匹配方法 |
CN110261899A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 地震数据z字形干扰波去除方法 |
Non-Patent Citations (3)
Title |
---|
叶佩: "地震初至波和反射波旅行时联合层析成像", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
彭愿: "地震数据的可视化技术", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
魏剑平: "高精度初至波和反射波联合速度反演", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020087767A1 (zh) * | 2018-10-31 | 2020-05-07 | 中国石油化工股份有限公司 | 基于台站与三维地震联合采集资料的速度反演方法 |
CN112824939A (zh) * | 2019-11-20 | 2021-05-21 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的剩余静校正方法、装置、设备及可读存储介质 |
CN113917524A (zh) * | 2020-07-10 | 2022-01-11 | 中国石油化工股份有限公司 | 基于无缆节点台站的近地表精细探测方法 |
CN113917524B (zh) * | 2020-07-10 | 2024-03-29 | 中国石油化工股份有限公司 | 基于无缆节点台站的近地表精细探测方法 |
CN112036008A (zh) * | 2020-07-28 | 2020-12-04 | 中国石油天然气集团有限公司 | 一种利用混合模型记录反演速度的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
WO2020087767A1 (zh) | 2020-05-07 |
CN109100792B (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100487488C (zh) | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 | |
CN109100792A (zh) | 基于台站与三维地震联合采集资料的速度反演方法 | |
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN106094032B (zh) | 一种构建地层速度模型的方法 | |
CN104597494B (zh) | 地震地层体分析方法及装置 | |
CN103645503B (zh) | 一种三维时间域照明分析及振幅补偿方法 | |
Nakanishi et al. | Three-dimensional plate geometry and P-wave velocity models of the subduction zone in SW Japan: Implications for seismogenesis | |
CN102841375A (zh) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 | |
CN104133245A (zh) | 一种地震资料的静校正方法及系统 | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
CN104793247A (zh) | 近地表速度建模方法及建模装置 | |
Zelt et al. | 3D simultaneous seismic refraction and reflection tomography of wide‐angle data from the central Chilean margin | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN109188506A (zh) | 一种适用于高铁隧底地震ct的纯地表三维观测系统 | |
Wolf et al. | Through thick and thin: A new model for the Colorado Plateau from seismic refraction data from Pacific to Arizona crustal experiment | |
CN106415321A (zh) | 用于储层建模的基于瞬时等时属性的地质体识别 | |
CN104166159B (zh) | 四维微地震监测的裂缝形态处理方法和系统 | |
CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
CN107817516A (zh) | 基于初至波信息的近地表建模方法及系统 | |
CN102866422B (zh) | 一种深度域地质实体模型生成方法 | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
CN109100798B (zh) | 实现折射多次波层析反演的方法、装置及处理终端 | |
CN102096099A (zh) | 折射波射线层析成像方法 | |
CN111650641B (zh) | 一种地壳三维结构模型融合方法及装置 | |
CN109143398A (zh) | 一种自动网格层析深度域速度的建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |