CN105911586B - 一种垂直地震数据成像方法及装置 - Google Patents
一种垂直地震数据成像方法及装置 Download PDFInfo
- Publication number
- CN105911586B CN105911586B CN201610228403.0A CN201610228403A CN105911586B CN 105911586 B CN105911586 B CN 105911586B CN 201610228403 A CN201610228403 A CN 201610228403A CN 105911586 B CN105911586 B CN 105911586B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- vertical seismic
- data
- imaging
- parameter
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 claims abstract description 44
- 238000000605 extraction Methods 0.000 claims abstract description 31
- 238000007781 pre-processing Methods 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 22
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000000153 supplemental effect Effects 0.000 abstract 2
- 230000006870 function Effects 0.000 description 19
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001052 transient effect Effects 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
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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
技术领域
本申请涉及地球物理勘探技术领域,特别涉及一种垂直地震数据成像方法及装置。
背景技术
垂直地震剖面(VSP)方法是一种介于地震和测井方法之间的地球物理技术,其是一种地震观测方法。它与通常地面观测的地震剖面相对应。垂直地震剖面方法是在地表附近的一些点上激发地震波,在沿井孔不同深度布置的一些多级多分量的检波点上进行观测。在垂直地震剖面中,因为检波器通过井置于地层内部,所以不仅能接收到自下而上传播的上行纵波和上行转换波,也能接收到自上而下传播的下行纵波及下行转换波,甚至能接收到横波。因此该地震方法比常规地震方法具有更高的信噪比和分辨率。然而VSP成像技术却一直制约着该方法的发展,其中两个主要的原因分别是:该观测方式中不能利用均匀的多次覆盖改善成像效果,以及转换成像或偏移成像过程中产生了波场畸变。不能利用均匀的多次覆盖改善成像效果是由炮检位置关系决定,难以通过室内处理改善提高。转换成像或偏移成像过程中产生了波场畸变是由于VSP成像处理过程中的动校正拉伸或偏移距拉伸产生,可以通过一些保真处理手段进行改善。
为了改善波场畸变,现有的垂直地震数据成像方法通常包括:对垂直地震数据进行预处理操作,所述预处理操作可以包括:自动增益处理、道均衡处理或带通波处理,再将所述预处理后的垂直地震数据直接向地震共深度点道集(CDP)进行行数据转换,即VSP-CDP(垂直地震剖面-地震共深度点道集)成像。
发明人发现现有技术中至少存在如下问题:现有技术中的垂直地震数据成像方法,在利用VSP-CDP成像的过程中,先对垂直地震数据进行预处理操作,再将包括振幅、相位、频率的地震数据采用VSP-CDP技术进行成像操作,该成像过程可以保真振幅,但频率和相位仍然有较大的畸变,从而导致频率和相位的保真性较低。
发明内容
本申请实施例的目的是提供一种垂直地震数据成像方法及装置,以提高垂直地震数据的频率和相位的保真性。
为解决上述技术问题,本申请实施例提供一种垂直地震数据成像方法及装置是这样实现的:
一种垂直地震数据成像方法,包括:获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理;采用预设的参数提取方法提取所述预处理后的垂直地震数据的参数对应的数据;所述参数包括:振幅、频率和相位;采用第一成像方法将所述提取的每一种参数对应的数据分别转换为参数数据体;根据所述参数数据体,利用第二成像方法得到垂直地震数据成像。
优选方案中,所述预设的参数提取方法,具体包括:确定滑动时间窗长度和所述滑动时间窗的位置;在所述滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数;提取所述拟合函数中的所述参数对应的数据,将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据。
优选方案中,所述滑动时间窗长度大于或等于一个波长对应的时间长度。
优选方案中,所述滑动时间窗的长度为35毫秒~50毫秒。
优选方案中,所述预设的参数提取方法还包括:从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。
优选方案中,所述拟合函数的表达式为:
其中,x表示当前滑动时窗中拟合的函数;A表示振幅,f表示频率,t表示时间,表示相位;其中,A、f、为所述提取的参数对应的数据。
优选方案中,所述采用第一成像方法具体包括:采用VSP-CDP的成像方法,即垂直地震剖面-地震共深度点道集的成像方法。
优选方案中,根据所述参数数据体,利用第二成像方法得到垂直地震数据成像,包括:将相同位置处的三个参数数据体作为一组,带入正弦函数中计算得到新的地震数据,所有的地震数据的集合为垂直地震数据成像。
优选方案中,所述对采集到的原始垂直地震数据进行预处理,包括下述中的至少一种:带通滤波处理、相位转换处理、旋转定向处理、反褶积计算、波场分离处理。
一种垂直地震数据成像装置,包括:预处理模块、参数提取模块、参数分别转换模块和成像模块;其中,所述预处理模块,用于获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理;所述参数提取模块,用于采用预设的参数提取方法提取所述预处理模块预处理后的垂直地震数据的参数对应的数据;所述参数包括:振幅、频率和相位;所述参数分别转换模块,用于采用第一成像方法将所述参数提取模块提取的每一种参数对应的数据分别转换为参数数据体;所述成像模块,用于根据所述参数分别转换模块转换的参数数据体,利用第二成像方法得到垂直地震数据成像。
优选方案中,所述参数提取模块包括:确定模块、拟合模块和提取模块;其中,所述确定模块,用于确定滑动时间窗长度和所述滑动时间窗的位置;所述拟合模块,用于在所述确定模块,确定的滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数;所述提取模块,用于提取所述拟合模块,拟合的函数中的参数对应的数据,将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据;所述参数包括:振幅、频率和相位。
优选方案中,所述参数提取模块还包括:滑动模块,其用于从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。
由以上本申请实施例提供的技术方案可见,本申请实施例提供的垂直地震数据成像方法和装置,先提取垂直地震数据中的振幅、频率和相位三种参数,再基于VSP-CDP成像方法将三种参数的数据分别转换为数据体,再将同一位置的三个数据体进行合成,避免了VSP-CDP成像过程中数据拉伸和畸变影响,可以使成像数据的振幅、频率和相位得到最大限度的保护,从而提高垂直地震数据的频率和相位的保真性。同时,本申请对原始垂直地震数据进行预处理的过程中,没有采用破坏振幅、频率、相位相对关系的处理操作,可以进一步提高垂直地震数据的频率和相位的保真性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一个垂直地震数据成像方法实施例的流程图;
图2是本申请一个垂直地震数据成像装置实施例的模块图;
图3是本申请装置实施例中参数提取模块的一个组成图。
具体实施方式
本申请实施例提供一种垂直地震数据成像方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请一个垂直地震数据成像方法实施例的流程图。如图1所示,所述垂直地震数据成像方法可以包括:
S101:获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理。
可以获取采集到的原始垂直地震数据。所述原始垂直地震数据可以通过垂直地震剖面(VSP)方法获取到。
所述原始垂直地震数据中可以包括振幅、相位和频率。
可以对所述采集到的原始垂直地震数据进行预处理。
所述预处理可以包括下述操作中的至少一种:带通滤波处理、相位转换处理、旋转定向处理、反褶积计算、波场分离处理。
通过上述预处理操作,可以实现对原始垂直地震数据中的去噪处理,得到转换成像前的垂直地震数据。上述预处理的步骤一般不会破坏频率、振幅和相位的相对关系,也不会改变数据的连续性,因此可以较好地实现对原始垂直地震数据的保真。
S102:采用预设的参数提取方法提取所述预处理后的垂直地震数据的参数对应的数据。
可以采用预设的参数提取方法提取所述预处理后的垂直地震数据的参数对应的数据,其中,所述参数可以包括:振幅、频率和相位。那么,提取的数据可以分别是振幅数据、频率数据和相位数据。
所述预设的参数提取方法,具体可以包括:确定滑动时间窗长度和所述滑动时间窗的位置,在所述滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数,提取所述拟合函数中的所述参数对应的数据;将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据。
其中,所述滑动时间窗的长度可以根据所述垂直地震数据的频率来确定,具体地,所述滑动时间窗的长度可以大于或等于一个波长对应的时间长度。
例如,垂直地震数据的频率约为20-30HZ(赫兹),其对应的波长为33-50毫秒。那么所述滑动时间窗的长度可以为33-50毫秒,一般可以选取为50毫秒。
进一步地,所述拟合函数的表达式为:
其中,x表示当前滑动时窗中拟合的函数;A表示振幅,f表示频率,t表示时间,表示相位;其中,A、f、为所述提取的参数对应的数据。
需要说明的是,所述对滑动时间窗长度内所述垂直地震数据进行正弦拟合运算可以通过现有的数学计算方法实现,本申请对此并不作出限定
进一步地,所述预设的参数提取方法还可以包括:从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。所述单位时长可以是1毫秒。
通过逐步滑动所述滑动时间窗,可以确定每个时间窗中心位置处的参数数据。
所述确定的每个参数对应的数据的个数一般小于所述预处理后的垂直地震数据的参数对应的数据个数。由于滑动时间窗具备一定的时间窗长度,时间窗的中心位置不能到达所述预处理后的垂直地震数据的起始位置和终点位置,所以不能到达的位置处参数数据被舍弃。舍弃的参数数据相比于所有地震数据数量较小,因此,不会影响地震数据的真实性。
S103:采用第一成像方法将所述提取的每一种参数对应的数据分别转换为参数数据体。
可以采用第一成像方法,将所述提取的每一种参数对应的数据分别转换为参数数据体。
具体地,所述第一成像方法可以采用VSP-CDP的成像方法,可以实现振幅的保真性。
分别转换成像的过程中,所述第一成像方法还可以采用现有技术中的其他方法,本申请对此并不作出限定。
S104:根据所述参数数据体,利用第二成像方法得到垂直地震数据成像。
可以根据步骤S103中所述转换的参数数据体,利用第二成像方法得到垂直地震数据成像。所述第二成像方法具体可以包括:将相同位置处的三个参数数据体作为一组,带入正弦函数,即公式(1),中计算得到新的地震数据,所有的地震数据的集合为垂直地震数据成像。
上述实施例提供的垂直地震数据成像方法,先提取垂直地震数据中的振幅、频率和相位三种参数,再基于VSP-CDP成像方法将三种参数的数据分别转换为数据体,再将同一位置的三个数据体进行合成,避免了VSP-CDP成像过程中数据拉伸和畸变影响,可以使成像数据的振幅、频率和相位得到最大限度的保护,从而提高垂直地震数据的频率和相位的保真性。同时,本申请对原始垂直地震数据进行预处理的过程中,没有采用破坏振幅、频率、相位相对关系的处理操作,可以进一步提高垂直地震数据的频率和相位的保真性。
下面介绍本申请的垂直地震数据成像装置实施例。
图2是本申请一个垂直地震数据成像装置实施例的模块图。如图2所示,所述垂直地震数据成像装置可以包括:预处理模块201、参数提取模块202、参数分别转换模块203和成像模块204。其中,
所述预处理模块201,可以用于获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理。
所述参数提取模块202,可以用于采用预设的参数提取方法提取所述预处理模块201预处理后的垂直地震数据的参数对应的数据;所述参数包括:振幅、频率和相位。
所述参数分别转换模块203,可以用于采用第一成像方法将所述参数提取模块202提取的每一种参数对应的数据分别转换为参数数据体。
所述成像模块204,可以用于根据所述参数分别转换模块203转换的参数数据体,利用第二成像方法得到垂直地震数据成像。
图3是本申请装置实施例中参数提取模块的一个组成图。如图3所示,所述参数提取模块202可以包括:确定模块2021、拟合模块2022和提取模块2023。
所述确定模块2021,可以用于确定滑动时间窗长度和所述滑动时间窗的位置。所述滑动时间窗的长度可以根据所述垂直地震数据的频率来确定,具体地,所述滑动时间窗的长度可以大于或等于一个波长对应的时间长度。例如,垂直地震数据的频率约为20-30HZ(赫兹),其对应的波长为33-50毫秒。那么所述滑动时间窗的长度可以为33-50毫秒,一般可以选取为50毫秒。
所述拟合模块2022,可以用于在所述确定模块2021确定的滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数。
所述提取模块2023,可以用于提取所述拟合模块2022拟合的函数中的参数对应的数据,将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据。所述参数可以包括:振幅、频率和相位。
进一步地,所述参数提取模块202还可以包括:滑动模块2024,其可以用于从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。
上述实施例公开的垂直地震数据成像装置与本申请的垂直地震数据成像方法实施例相对应,可以实现本申请的方法实施例并取得方法实施例的技术效果。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种垂直地震数据成像方法,其特征在于,包括:
获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理;
采用预设的参数提取方法提取所述预处理后的垂直地震数据的参数对应的数据;所述参数包括:振幅、频率和相位;
采用第一成像方法将所述提取的每一种参数对应的数据分别转换为参数数据体;所述采用第一成像方法具体包括:采用VSP-CDP的成像方法,即垂直地震剖面-地震共深度点道集的成像方法;
根据所述参数数据体,利用第二成像方法得到垂直地震数据成像;具体包括:将相同位置处的三个参数数据体作为一组,带入正弦函数中计算得到新的地震数据,所有的地震数据的集合为垂直地震数据成像。
2.如权利要求1所述的一种垂直地震数据成像方法,其特征在于,所述预设的参数提取方法,具体包括:
确定滑动时间窗长度和所述滑动时间窗的位置;
在所述滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数;
提取所述拟合函数中的所述参数对应的数据,将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据。
3.如权利要求2所述的一种垂直地震数据成像方法,其特征在于,所述滑动时间窗长度大于或等于一个波长对应的时间长度。
4.如权利要求3所述的一种垂直地震数据成像方法,其特征在于,所述滑动时间窗的长度为35毫秒~50毫秒。
5.如权利要求2所述的一种垂直地震数据成像方法,其特征在于,所述预设的参数提取方法还包括:从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。
6.如权利要求2所述的一种垂直地震数据成像方法,其特征在于,所述拟合函数的表达式为:
其中,x表示当前滑动时窗中拟合的函数;A表示振幅,f表示频率,t表示时间,表示相位;其中,A、f、为所述提取的参数对应的数据。
7.如权利要求1所述的一种垂直地震数据成像方法,其特征在于,所述对采集到的原始垂直地震数据进行预处理,包括下述中的至少一种:带通滤波处理、相位转换处理、旋转定向处理、反褶积计算、波场分离处理。
8.一种垂直地震数据成像装置,其特征在于,包括:预处理模块、参数提取模块、参数分别转换模块和成像模块;其中,
所述预处理模块,用于获取采集到的原始垂直地震数据,对所述采集到的原始垂直地震数据进行预处理;
所述参数提取模块,用于采用预设的参数提取方法提取所述预处理模块预处理后的垂直地震数据的参数对应的数据;所述参数包括:振幅、频率和相位;
所述参数分别转换模块,用于采用第一成像方法将所述参数提取模块提取的每一种参数对应的数据分别转换为参数数据体;所述采用第一成像方法具体包括:采用VSP-CDP的成像方法,即垂直地震剖面-地震共深度点道集的成像方法;
所述成像模块,用于根据所述参数分别转换模块转换的参数数据体,利用第二成像方法得到垂直地震数据成像;具体包括:将相同位置处的三个参数数据体作为一组,带入正弦函数中计算得到新的地震数据,所有的地震数据的集合为垂直地震数据成像。
9.如权利要求8所述的一种垂直地震数据成像装置,其特征在于,所述参数提取模块包括:确定模块、拟合模块和提取模块;其中,
所述确定模块,用于确定滑动时间窗长度和所述滑动时间窗的位置;
所述拟合模块,用于在所述确定模块,确定的滑动时间窗长度内对所述垂直地震数据进行正弦拟合运算得到拟合函数;
所述提取模块,用于提取所述拟合模块,拟合的函数中的参数对应的数据,将所述提取的数据作为所述滑动时间窗中心位置处对应的参数数据;所述参数包括:振幅、频率和相位。
10.如权利要求9所述的一种垂直地震数据成像装置,其特征在于,所述参数提取模块还包括:滑动模块,其用于从所述预处理后的垂直地震数据的起始位置开始,每提取一次参数,滑动单位时长的所述滑动时间窗,直至所述预处理后的垂直地震数据的终点位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610228403.0A CN105911586B (zh) | 2016-04-13 | 2016-04-13 | 一种垂直地震数据成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610228403.0A CN105911586B (zh) | 2016-04-13 | 2016-04-13 | 一种垂直地震数据成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105911586A CN105911586A (zh) | 2016-08-31 |
CN105911586B true CN105911586B (zh) | 2018-10-16 |
Family
ID=56746788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610228403.0A Active CN105911586B (zh) | 2016-04-13 | 2016-04-13 | 一种垂直地震数据成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105911586B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
CN102213769A (zh) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 |
CN105093292A (zh) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种地震成像的数据处理方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2593687C2 (ru) * | 2012-04-02 | 2016-08-10 | Лэндмарк Графикс Корпорейшн | Система и способ вертикального сейсмического профилирования с представлением разведочных данных в виде комбинации параметризованных компрессионного, сдвигового и дисперсивного волновых полей |
-
2016
- 2016-04-13 CN CN201610228403.0A patent/CN105911586B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
CN102213769A (zh) * | 2010-04-07 | 2011-10-12 | 中国石油天然气集团公司 | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 |
CN105093292A (zh) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种地震成像的数据处理方法和装置 |
Non-Patent Citations (2)
Title |
---|
"复杂构造地区零井源距VSP成像方法研究";蔡志东 等;《石油物探》;20150531;第54卷(第3期);第309-316页 * |
"大阵列3D-VSP技术在大庆油田的应用";李彦鹏 等;《石油地球物理勘探》;20110430;第46卷(第2期);第311-316页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105911586A (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106054251B (zh) | 一种初至波拾取方法及装置 | |
US11086034B2 (en) | Method and apparatus for determining permeability of reservoir | |
CN107193045B (zh) | 一种地震数据处理方法及装置 | |
CN107563297B (zh) | 一种波形匹配方法及装置 | |
CN106338763B (zh) | 一种确定地震剖面显示数据的方法及装置 | |
CN106646611A (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN106842304B (zh) | 一种叠前深度偏移方法及装置 | |
CN108897036B (zh) | 一种地震数据处理方法及装置 | |
CN106353816B (zh) | 一种地震采集脚印噪音压制方法及系统 | |
CN107144878B (zh) | 一种断层识别方法及装置 | |
CN105044773A (zh) | 转换波井震标定方法及装置 | |
CN108181644A (zh) | 一种采样方法及装置 | |
CN108035708B (zh) | 一种去除地层界面反射波的方法及装置 | |
CN111060961A (zh) | 基于多信息约束反演的品质因子确定方法、装置及系统 | |
CN105911586B (zh) | 一种垂直地震数据成像方法及装置 | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
CN106199693B (zh) | 地震数据速度谱自动拾取方法和装置 | |
CN111596358B (zh) | 一种多次波的压制方法、设备以及系统 | |
CN106772596B (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN111239832B (zh) | 3d3c-vsp成像处理方法、装置及设备 | |
CN107561588A (zh) | 一种地震数据噪声压制方法及装置 | |
CN106094026B (zh) | 一种获取垂直地震数据走廊叠加剖面的方法及装置 | |
CN108693558A (zh) | 地震数据处理方法和装置 | |
CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
CN106054253A (zh) | 一种道集同相轴校正的方法及装置 |
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 |