CN105182408B - 一种合成地震记录的制作方法和装置 - Google Patents
一种合成地震记录的制作方法和装置 Download PDFInfo
- Publication number
- CN105182408B CN105182408B CN201510542801.5A CN201510542801A CN105182408B CN 105182408 B CN105182408 B CN 105182408B CN 201510542801 A CN201510542801 A CN 201510542801A CN 105182408 B CN105182408 B CN 105182408B
- Authority
- CN
- China
- Prior art keywords
- data
- correction
- time
- curve
- time point
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例提供一种合成地震记录的制作方法和装置。该方法包括:获取包括时深关系数据和走廊叠加剖面数据的垂直地震数据;利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;利用深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;根据校正后的声波测井数据和校正后的密度测井数据计算得到反射系数序列;利用时深关系数据对反射系数序列进行时深转换得到时间域反射系数序列;根据走廊叠加剖面数据确定主频曲线;根据主频曲线生成雷克子波集;将雷克子波集和时间域反射系数序列进行褶积计算得到合成地震记录。本申请实施例提供的技术方案可获取符合真实地震反射规律的合成地震记录。
Description
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及一种合成地震记录的制作方法和装置。
背景技术
合成地震记录是地震数据和测井数据之间建立联系的一个重要桥梁。在地球物理勘探技术领域中合成地震记录在构造解释、储层分析、波阻抗反演和油藏描述等方面发挥着重要的作用。合成地震记录的制作精度将影响地震地质层位标定的准确性。
现有技术中一般利用地震子波和反射系数进行褶积制作合成地震记录。但采用的地震子波一般是给定子波或从井旁地震道提取的子波,这两种方法所利用的均为恒定子波,对于真实的地层反射情况描述存在着一定的偏差。同时,采用的反射系数是直接根据测井数据获得的,一般测井数据分辨率较低,导致得到的反射系数也无法准确反映地层反射情况。因此,这样制作的合成地震记录就不能准确的建立测井资料与地震资料之间的对应关系,无法有效的进行地震层位标定、储层解释等。
因此,亟需一种可以真实反映地震反射规律的合成地震记录的制作方法。
发明内容
本申请的目的是提供一种合成地震记录的制作方法和装置,以获取符合真实地震反射规律的合成地震记录,为地震层位标定、储层解释等提供支持。
为了实现上述目的,本申请提供了一种合成地震记录的制作方法,该方法包括:
利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;
利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;
根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列;
利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列;
根据走廊叠加剖面数据确定主频曲线;
根据所述主频曲线生成雷克子波集;
将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录。
在一个优选的实施例中,所述根据走廊叠加剖面数据确定主频曲线包括:
S1:设定扫描时窗;
S2:选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱;
S3:对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据;
S4:更换所述走廊叠加剖面中的时间点,循环步骤S2至S3,直到遍历完所述走廊叠加剖面中的时间点;
S5:将所述主频数据以时间关系排序确定主频曲线。
在一个优选的实施例中,所述根据所述主频曲线生成雷克子波集包括:
根据如下公式确定所述主频曲线中每个时间点的雷克子波:
上式中,Ai(t)代表i时刻的雷克子波振幅;fi代表i时刻的主频数据,单位为Hz;t代表i时刻所对应时窗中的时间点,单位为s;
根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集。
在一个优选的实施例中,所述根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列包括:
将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;
根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。
在一个优选的实施例中,所述将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录包括:
S11:将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算;
S12:获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据。
S13:更换所述雷克子波集中的时间点,循环步骤S11至S12,直到遍历完所述雷克子波中的时间点;
S14:将所述合成数据以时间关系排序确定合成地震记录。
本申请另一方面还提供一种合成地震记录的制作装置,该装置包括:
第一校正模块,用于利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;
第二校正模块,用于利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;
第一计算模块,用于根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列;
时深转换模块,用于利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列;
第一确定模块,用于根据走廊叠加剖面数据确定主频曲线;
雷克子波集生成模块,用于根据所述主频曲线生成雷克子波集;
第二计算模块,用于将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录。
在一个优选的实施例中,所述第一确定模块包括:
扫描时窗设定单元,用于设定扫描时窗;
频率谱获取单元,用于选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱;
主频数据获取单元,用于对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据;
第一数据处理单元,用于更换所述走廊叠加剖面中的时间点,并利用频率谱获取单元和主频数据获取单元获取所述时间对应时窗数据的主频数据,直到遍历完所述走廊叠加剖面中的时间点;
第二确定单元,用于将所述主频数据以时间关系排序确定主频曲线。
在一个优选的实施例中,所述第一计算模块包括:
第一确定单元,用于将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;
第一计算单元,用于根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。
在一个优选的实施例中,所述雷克子波集生成模块包括:
第二计算单元,用于根据如下公式确定所述主频曲线中每个时间点的雷克子波:
上式中,Ai(t)代表i时刻的雷克子波振幅;fi代表i时刻的主频数据,单位为Hz;t代表i时刻所对应时窗中的时间点,单位为s;
第三确定单元,用于根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集。
在一个优选的实施例中,所述第二计算模块包括:
第三计算单元,用于将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算;
第二数据处理单元,用于获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据;
第三数据处理单元,用于更换所述雷克子波集中的时间点,并利用所述第三计算单元和所述第二数据处理单元确定所述时间点对应合成记录数据,直到遍历完所述雷克子波中的时间点;
第四确定单元,用于将所述合成数据以时间关系排序确定合成地震记录。
由以上本申请实施例提供的技术方案可见,本申请实施例通过对测井数据井时深校正,提高了测井数据的分辨率,且根据校正后的测井数据计算得到的反射系数序列可以更准确反映地层反射情况,并将所述反射系数序列进行时深转换得到时间域反射系数序列。同时,根据垂直地震数据中的走廊叠加剖面数据确定的雷克子波集更具有描述地震地质反射特征的能力。最后,根据雷克子波集和时间域反射系数序列进行褶积可以计算得到符合真实地震反射规律的合成地震记录。与现有技术相比,本申请实施例制作的合成地震记录符合真实地震反射规律,为地震层位标定、储层解释等提供了有效支持。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种合成地震记录的制作方法的实施例的流程图;
图2是本申请实施例中所述校正后的声波测井数据对应的声波测井曲线部分示意图;
图3是本申请实施例中所述校正后的密度测井数据对应的度测井曲线部分示意图;
图4是本申请实施例中所述反射系数序列对应的深度域反射系数曲线部分示意图;
图5是本申请实施例中获取主频曲线的过程的示意图;
图6是本申请实施例中进行褶积计算得到合成地震记录过程的示意图;
图7是本申请实施例提供的合成地震记录镶嵌至地面地震中的部分示意图;
图8是本申请实施例提供的一种成地震记录的制作装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
在实际应用中,垂直地震方法是一种介于地震和测井方法之间的地球物理技术。该方法可以在地表设置震源,在井中设置检波器接收地震波,通过对所得资料经过静校正、能量补偿、反褶积、波场分离、同相叠加等处理,得到走廊叠加剖面。由于井中检波器排列接近于目的层位,所以走廊叠加剖面通常具有比地震勘探更高的纵向分辨率。因此,更具有描述地震地质反射特征的能力;同时垂直地震方法所获得的地震波时深关系也具有较高的可信度。
在一些实施例中,在制作合成地震记录之前需要先进行数据的准备。具体的,可以包括测井数据的准备和垂直地震数据的准备。所述测井数据的准备可以包括:将采集的测井数据进行去除重复深度、去除异常噪声、校正因不同时期采集造成的时差等处理工作,得到整理后的测井数据。所述测井数据可以至少包括下述之一:深度数据、声波时差、密度测井数据;且所述测井数据的深度起算面可以为地表,如果从井台起算的数据可以减去井台高度。所述垂直地震数据的准备可以包括:对采集的垂直地震数据进行静校正、能量补偿、反褶积、波场分离、同相叠加等处理工作,得到时深关系数据和走廊叠加剖面数据。且所述时深关系数据中的深度起算面可以为地表。此外,本申请实施例中所涉及深度、时间和密度单位为国际单位制。
以下介绍本申请一种合成地震记录的制作方法的实施例。结合附图1,该实施例包括:
S110:利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据。
在一些实施例中,所述时深关系数据可以包括时间域地震数据与深度域地震数据之间的对应关系。声波测井资料由于受地层、井径、井孔扩大、泥浆渗、入地层等诸多因素的影响,会有不同程度的失真。因此,可以利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据。所述深度校正数据可以包括所述时深关系数据在上述对声波测井数据进行校正处理过程中根据所述声波测井数据进行相应时差调整后的时深关系数据。
本申请利用时深关系数据对声波测井数据进行校正处理的过程可以包括将所述声波测井数据转换为时深数据,并计算所述时深数据的相对深度的时差;计算所述时深关系数据的相对深度的时差;然后,对上述两个相对深度的时差进行比较校正;将校正后的时深数据转换为声波测井数据。在一个具体的实施例中,可以详见申请号为201210418575.6,发明名称为一种利用VSP时深关系对声波测井资料进行校正的方法的专利文件。
S120:利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据。
在一些实施例中,在步骤S110获得所述时深校正数据之后,可以利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据。具体的,所述校正处理可以包括对所述密度测井数据中某一深度的密度值的校正处理,也可以包括对所述密度测井数据中某一段深度中的密度值校正处理。
S130:根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列。
在一些实施例中,可以根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列。具体的,可以将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;然后,根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。所述地层声波速度可以根据所述地层的层厚和相应的声波时差计算确定。反射系数可以通过如下公式计算获得:
上式中,R可以代表反射系数;ν2可以代表反射面下层的地震速度,单位为m/s;ν1可以代表反射面上层的地震速度,单位为m/s;ρ2可以代表反射面下层的岩石密度,单位为kg/m2;ρ1可以代表反射面上层的岩石密度,单位为kg/m3。
分别计算出每个方面的反射系数,并根据全部反射系数组成反射系数系列。
如图2所示是所述校正后的声波测井数据对应的声波测井曲线部分示意图,其中横坐标代表深度,单位为m,纵坐标代表声波时差,单位为m/s;如图3所示是所述校正后的密度测井数据对应的度测井曲线部分示意图,其中横坐标代表深度,单位为m,纵坐标代表密度,单位为kg/m3;如图4所示是所述反射系数序列对应的深度域反射系数曲线部分示意图,其中横坐标代表深度,单位为m,纵坐标代表反射系数。
S140:利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列。
在一些实施例中,在步骤S130获得反射系数序列之后,可以利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列。具体的,所述时深转换可以包括利用插值计算的方式,也可以包括利用商业测井处理软件的方式。此外,所述时深转换过程中的采样间隔可以为1ms。
此外,本申请实施例中所述时深转换并不仅限于上述的方式,在实际应用中,还可以根据具体的情况,选择其他时深转换的方式,本申请实施例并不仅限于此。
S150:根据走廊叠加剖面数据确定主频曲线。
在一些实施例中,可以根据走廊叠加剖面数据确定主频曲线。具体的,结合附图5,可以包括:
S151:设定扫描时窗。
在一些实施例中,为了保证频率计算结果稳定,扫描时窗一般可以至少应包含一个波长。由于垂直地震记录的主频通常为10-50Hz,因此扫描时窗可以确定为100ms。
S152:选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱。
在一些实施例,所述时频转换可以包括傅里叶变换,Gabor变换,连续小波变换等。
S153:对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据。
S154:更换所述走廊叠加剖面中的时间点,循环步骤S152至S153,直到遍历完所述走廊叠加剖面中的时间点。
S155:将所述主频数据以时间关系排序确定主频曲线。
此外,所述走廊叠加剖面数据的边界部分,如不能满足时窗长度,可以以0Hz填充。
S160:根据所述主频曲线生成雷克子波集。
在一些实施例中,在步骤S150之后,可以根据所述主频曲线生成雷克子波集。具体的,可以利用如下公式确定所述主频曲线中每个时间点对应的雷克子波:
上式中,Ai(t)可以代表i时刻的雷克子波振幅;fi可以代表i时刻的主频数据,单位为Hz;t可以代表i时刻所对应时窗中的时间点,单位为s。
根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集;所述雷克子波可以包括旁瓣振幅为主瓣振幅44.63%的正极性、零相位的子波。
S170:将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录。
在一些实施例中,可以将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录。具体的,结合附图6,可以包括:
S171:将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算。
S172:获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据。
S173:更换所述雷克子波集中的时间点,循环步骤S171至S172,直到遍历完所述雷克子波中的时间点。
S174:将所述合成数据以时间关系排序确定合成地震记录。
如图7所示的是利用本申请实施例技术方案制作的合成地震记录镶嵌至地面地震中的部分示意图,从图中可见,本申请实施例中制作的合成地震记录720在频率及波阻关系均与地面地震710对应良好。
由此可见,本申请一种合成地震记录的制作方法实施例提供的技术方案通过对测井数据井时深校正,提高了测井数据的分辨率,且根据校正后的测井数据计算得到的反射系数序列可以更准确反映地层反射情况,并将所述反射系数序列进行时深转换得到时间域反射系数序列。同时,根据垂直地震数据中的走廊叠加剖面数据确定的雷克子波集更具有描述地震地质反射特征的能力。最后,根据雷克子波集和时间域反射系数序列进行褶积可以计算得到符合真实地震反射规律的合成地震记录。与现有技术相比,本申请实施例制作的合成地震记录符合真实地震反射规律,为地震层位标定、储层解释等提供了有效支持。
本申请另一方面还提供一种合成地震记录的制作装置,结合附图8,该装置800包括:
第一校正模块810,可以用于利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;
第二校正模块820,可以用于利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;
第一计算模块830,可以用于根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列;
时深转换模块840,可以用于利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列;
第一确定模块850,可以用于根据走廊叠加剖面数据确定主频曲线;
雷克子波集生成模块860,可以用于根据所述主频曲线生成雷克子波集;
第二计算模块870,可以用于将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录。
在一个优选的实施例中,所述第一确定模块850可以包括:
扫描时窗设定单元,可以用于设定扫描时窗;
频率谱获取单元,可以用于选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱;
主频数据获取单元,可以用于对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据;
第一数据处理单元,可以用于更换所述走廊叠加剖面中的时间点,并利用频率谱获取单元和主频数据获取单元获取所述时间对应时窗数据的主频数据,直到遍历完所述走廊叠加剖面中的时间点;
第二确定单元,可以用于将所述主频数据以时间关系排序确定主频曲线。
在一个优选的实施例中,所述第一计算模块830可以包括:
第一确定单元,可以用于将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;
第一计算单元,可以用于根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。
在一个优选的实施例中,所述雷克子波集生成模块860可以包括:
第二计算单元,可以用于根据如下公式确定所述主频曲线中每个时间点的雷克子波:
上式中,Ai(t)代表i时刻的雷克子波振幅;fi代表i时刻的主频数据,单位为Hz;t代表i时刻所对应时窗中的时间点,单位为s;
第三确定单元,可以用于根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集。
在一个优选的实施例中,所述第二计算模块870可以包括:
第三计算单元,可以用于将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算;
第二数据处理单元,可以用于获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据;
第三数据处理单元,可以用于更换所述雷克子波集中的时间点,并利用所述第三计算单元和所述第二数据处理单元确定所述时间点对应合成记录数据,直到遍历完所述雷克子波中的时间点;
第四确定单元,可以用于将所述合成数据以时间关系排序确定合成地震记录。
由此可见,本申请一种合成地震记录的制作方法和装置实施例提供的技术方案通过对测井数据井时深校正,提高了测井数据的分辨率,且根据校正后的测井数据计算得到的反射系数序列可以更准确反映地层反射情况,并将所述反射系数序列进行时深转换得到时间域反射系数序列。同时,根据垂直地震数据中的走廊叠加剖面数据确定的雷克子波集更具有描述地震地质反射特征的能力。最后,根据雷克子波集和时间域反射系数序列进行褶积计算得到符合真实地震反射规律的合成地震记录。与现有技术相比,本申请实施例制作的合成地震记录符合真实地震反射规律,为地震层位标定、储层解释等提供了有效支持。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (8)
1.一种合成地震记录的制作方法,其特征在于,该方法包括:
利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;
利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;
根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列;
利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列;
根据走廊叠加剖面数据确定主频曲线;
根据所述主频曲线生成雷克子波集;
将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录;
其中,所述根据走廊叠加剖面数据确定主频曲线包括:
S1:设定扫描时窗;
S2:选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱;
S3:对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据;
S4:更换所述走廊叠加剖面中的时间点,循环步骤S2至S3,直到遍历完所述走廊叠加剖面中的时间点;
S5:将所述主频数据以时间关系排序确定主频曲线。
2.根据权利要求1所述的方法,其特征在于,所述根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列包括:
将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;
根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。
3.根据权利要求1所述的方法,其特征在于,所述根据所述主频曲线生成雷克子波集包括:
根据如下公式确定所述主频曲线中每个时间点的雷克子波:
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;f</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;f</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
</mrow>
上式中,Ai(t)代表i时刻的雷克子波振幅;fi代表i时刻的主频数据,单位为Hz;t代表i时刻所对应时窗中的时间点,单位为s;
根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集。
4.根据权利要求1所述的方法,其特征在于,所述将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录包括:
S11:将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算;
S12:获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据;
S13:更换所述雷克子波集中的时间点,循环步骤S11至S12,直到遍历完所述雷克子波中的时间点;
S14:将所述合成记录数据以时间关系排序确定合成地震记录。
5.一种合成地震记录的制作装置,其特征在于,该装置包括:
第一校正模块,用于利用时深关系数据对声波测井数据进行校正处理,获取校正后的声波测井数据及深度校正数据;
第二校正模块,用于利用所述深度校正数据对密度测井数据进行校正处理,获取校正后的密度测井数据;
第一计算模块,用于根据所述校正后的声波测井数据和所述校正后的密度测井数据计算得到反射系数序列;
时深转换模块,用于利用所述时深关系数据对所述反射系数序列进行时深转换得到时间域反射系数序列;
第一确定模块,用于根据走廊叠加剖面数据确定主频曲线;
雷克子波集生成模块,用于根据所述主频曲线生成雷克子波集;
第二计算模块,用于将所述雷克子波集和所述时间域反射系数序列进行褶积计算得到合成地震记录;
其中,所述第一确定模块包括:
扫描时窗设定单元,用于设定扫描时窗;
频率谱获取单元,用于选择所述走廊叠加剖面中时间点,以所述时间点为扫描时窗的中心选取对应的时窗数据,并对所述时窗数据进行时频转换得到频率谱;
主频数据获取单元,用于对所述频率谱进行包络拟合处理得到包络曲线,将所述包络曲线中的包络峰值作为所述时窗数据的主频数据;
第一数据处理单元,用于更换所述走廊叠加剖面中的时间点,并利用所述频率谱获取单元和所述主频数据获取单元获取所述时间对应时窗数据的主频数据,直到遍历完所述走廊叠加剖面中的时间点;
第二确定单元,用于将所述主频数据以时间关系排序确定主频曲线。
6.根据权利要求5所述的装置,其特征在于,所述第一计算模块包括:
第一确定单元,用于将根据所述校正后的声波测井数据确定地层声波速度作为地层地震速度;
第一计算单元,用于根据所述地层地震速度和所述校正后的密度测井数据计算得到反射系数序列。
7.根据权利要求5所述的装置,其特征在于,所述雷克子波集生成模块包括:
第二计算单元,用于根据如下公式确定所述主频曲线中每个时间点的雷克子波:
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;f</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&pi;f</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
</mrow>
上式中,Ai(t)代表i时刻的雷克子波振幅;fi代表i时刻的主频数据,单位为Hz;t代表i时刻所对应时窗中的时间点,单位为s;
第三确定单元,用于根据所述主频曲线中的每一个时间点对应的雷克子波确定雷克子波集。
8.根据权利要求5所述的装置,其特征在于,所述第二计算模块包括:
第三计算单元,用于将所述雷克子波集中时间点对应的雷克子波与所述时间域反射系数序列进行褶积计算;
第二数据处理单元,用于获取褶积结果中所述时间点对应的褶积值,将所述褶积值作为所述时间点的合成记录数据;
第三数据处理单元,用于更换所述雷克子波集中的时间点,并利用所述第三计算单元和所述第二数据处理单元确定所述时间点对应合成记录数据,直到遍历完所述雷克子波中的时间点;
第四确定单元,用于将所述合成记录数据以时间关系排序确定合成地震记录。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510542801.5A CN105182408B (zh) | 2015-08-28 | 2015-08-28 | 一种合成地震记录的制作方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510542801.5A CN105182408B (zh) | 2015-08-28 | 2015-08-28 | 一种合成地震记录的制作方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105182408A CN105182408A (zh) | 2015-12-23 |
CN105182408B true CN105182408B (zh) | 2017-09-01 |
Family
ID=54904604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510542801.5A Active CN105182408B (zh) | 2015-08-28 | 2015-08-28 | 一种合成地震记录的制作方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182408B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106094026B (zh) * | 2016-06-02 | 2018-06-01 | 中国石油天然气集团公司 | 一种获取垂直地震数据走廊叠加剖面的方法及装置 |
CN106772578B (zh) * | 2016-12-07 | 2018-11-09 | 中国矿业大学(北京) | 一种合成地震记录的方法和装置 |
CN107193040A (zh) * | 2017-06-27 | 2017-09-22 | 中国石油天然气股份有限公司 | 深度域合成地震记录的确定方法和装置 |
CN109212600A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 基于分辨能力的地震资料评价方法及系统 |
CN107561584B (zh) * | 2017-08-31 | 2020-11-24 | 中水北方勘测设计研究有限责任公司 | 一种岩体地震波、声波纵波速度对比换算的方法 |
CN108037532B (zh) * | 2017-12-06 | 2020-05-12 | 中国科学院地质与地球物理研究所 | 用于标定合成地震记录的方法、装置、系统和计算机可读介质 |
CN108344433A (zh) * | 2018-04-20 | 2018-07-31 | 吴伟 | 一种数据序列展示方法和装置 |
CN109061737B (zh) * | 2018-09-07 | 2021-04-09 | 中国石油化工股份有限公司 | 一种基于合成地震记录的储层预测方法及装置 |
CN109669223A (zh) * | 2019-01-07 | 2019-04-23 | 浙江大学 | 利用微测井获得浅地表声波测井与合成地震记录的方法 |
CN111060986B (zh) * | 2019-10-18 | 2022-01-28 | 中国石油化工股份有限公司 | 一种地层压力预测方法及岩性油藏评价方法 |
CN112946750B (zh) * | 2019-12-11 | 2024-01-30 | 中国石油天然气股份有限公司 | 井震标定方法及系统 |
CN111060970B (zh) * | 2019-12-30 | 2022-06-03 | 恒泰艾普(北京)能源科技研究院有限公司 | 一种储层约束反褶积的方法及应用 |
CN113376696A (zh) * | 2020-02-25 | 2021-09-10 | 中国石油天然气股份有限公司 | 基于多子波的井震标定方法及装置 |
CN112346129A (zh) * | 2020-10-29 | 2021-02-09 | 中国石油天然气集团有限公司 | 测井曲线提取及合成地震记录制作方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3961306A (en) * | 1971-10-28 | 1976-06-01 | Seiscom Delta Inc. | Method of forming color graphic displays from input data |
CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
CN103790569A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
CN104122588A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于谱分解提高叠后地震资料分辨率的方法 |
-
2015
- 2015-08-28 CN CN201510542801.5A patent/CN105182408B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3961306A (en) * | 1971-10-28 | 1976-06-01 | Seiscom Delta Inc. | Method of forming color graphic displays from input data |
CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
CN103790569A (zh) * | 2012-10-26 | 2014-05-14 | 中国石油化工股份有限公司 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
CN104122588A (zh) * | 2014-07-30 | 2014-10-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于谱分解提高叠后地震资料分辨率的方法 |
Non-Patent Citations (3)
Title |
---|
利用子波置换法提高地面地震资料的分辨率;何惺华;《油气地球物理》;20081031;第6卷(第4期);第4,15-18页 * |
基于数据驱动的合成地震记录技术及其应用研究;季敏 等;《石油物探》;20110731;第50卷(第4期);第373-377页 * |
高精度合成地震记录制作及层位精细标定;吴湘杰 等;《海洋地质前沿》;20111130;第27卷(第11期);第60-65页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105182408A (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105182408B (zh) | 一种合成地震记录的制作方法和装置 | |
Warner et al. | Anisotropic 3D full-waveform inversion | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
EP2920615B1 (en) | Process for creating image gathers | |
Ivandic et al. | Time-lapse analysis of sparse 3D seismic data from the CO2 storage pilot site at Ketzin, Germany | |
JP2020522699A (ja) | 地下構造の検出 | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
CN107817526B (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
WO2006036336A1 (en) | Method for depth migrating seismic data using pre-stack time migration, demigration, and post-stack depth migration | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
Gao et al. | Waveform tomography at a groundwater contamination site: VSP-surface data set | |
Holgate et al. | Constraining uncertainty in interpretation of seismically imaged clinoforms in deltaic reservoirs, Troll field, Norwegian North Sea: Insights from forward seismic models of outcrop analogs | |
Chapman et al. | On the geologic structure at the epicenter of the 1886 Charleston, South Carolina, earthquake | |
Colombo et al. | Near-surface full-waveform inversion in a transmission surface-consistent scheme | |
Kamei et al. | Application of waveform tomography to a crooked-line 2D land seismic data set | |
Pan et al. | Multichannel analysis of Love waves in a 3D seismic acquisition system | |
Smithyman et al. | Waveform tomography of field vibroseis data using an approximate 2D geometry leads to improved velocity models | |
US20100118655A1 (en) | Progressive 3D vertical seismic profiling method | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN109839659B (zh) | 一种对叠前深度偏移剖面进行迭代优化的方法 | |
CN104316959B (zh) | 一种基于等效流体声波阻抗的流体识别方法 | |
Zong et al. | Fractured basement imaging using random-space-shift reverse time migration: A vertical seismic profile survey in the Bohai Bay Basin, China | |
Caeiro et al. | Insights of processing early-outs from new high-quality 3D seismic data for exploration prospectivity: A land case example from Abu Dhabi | |
Gray et al. | Using advanced seismic imaging tools to see the invisible beneath foothills structures | |
O’Brien et al. | Multicomponent VSP imaging of tight-gas sands |
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 |