CN108828664B - 一种多次波识别方法及装置 - Google Patents
一种多次波识别方法及装置 Download PDFInfo
- Publication number
- CN108828664B CN108828664B CN201810579839.3A CN201810579839A CN108828664B CN 108828664 B CN108828664 B CN 108828664B CN 201810579839 A CN201810579839 A CN 201810579839A CN 108828664 B CN108828664 B CN 108828664B
- Authority
- CN
- China
- Prior art keywords
- reflection coefficient
- determining
- record
- primary
- seismic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 13
- 238000002310 reflectometry Methods 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 2
- 238000000605 extraction Methods 0.000 claims 1
- 238000005755 formation reaction Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 239000011435 rock Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000032798 delamination Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 206010013883 Dwarfism Diseases 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- JPMIIZHYYWMHDT-UHFFFAOYSA-N octhilinone Chemical compound CCCCCCCCN1SC=CC1=O JPMIIZHYYWMHDT-UHFFFAOYSA-N 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000011426 transformation method 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
- G01V1/30—Analysis
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
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
技术领域
本发明涉及石油勘探技术,特别是关于石油勘探的数据处理技术,具体的讲是一种多次波识别方法及装置。
背景技术
在陆上地震勘探中,地震波通过地下岩层传播,当地下有如基岩面、不整合面、火成岩面等强反射界面时,会产生层间多次波。这类多次波与深层一次波振幅相当,动校正时差较小,与深层一次波相干涉,难以识别,影响地震数据偏移质量,给后续的解释工作带来困扰。
目前层间多次波压制方法以叠前为主,多次覆盖技术以来,人们认识到地震数据中的多次波可以根据其时距曲线特征加以识别。一般情况下,地震波传播速度随着深度的增加而增加,与多次波相比,一次波传播速度要高。共中心点(CMP)道集中,多次波与一次波同相轴之间出现倾角或时间差异,NMO之后,一次波是平的,而多次波由于NMO校正时间不足是下拉的,通常可以采用一个数学变换的方法将地震数据变换到新的数据域,从而将一次波与多次波投影到可分离的不同区域。但是,如果地层速度发生反转,与浅层高速体有关的多次波速度会接近该高速体以下的一次波速度,使得两者干涉,难以分离。在很多地区,在叠前阶段经过多次波压制后,仍然会有较强能量、频带宽的多次波残留。
在叠后阶段,可以通过叠加地震剖面上强反射界面和强反射界面之下地震反射的形态和时间识别多次波。当地下反射界面倾角较小时,如果在强反射界面下方存在另一组与其形态一致,且出现时间小于其2倍的反射同相轴时,就可以判断该组同相轴为多次反射波。对于倾斜地层,叠加剖面上二次反射波的倾角约为其一次反射波的2倍。
现有技术中,叠后阶段识别多次波的方法相对单一,主要依靠地震剖面上地震反射的外观、形态判断,主观性较强,识别的结果说服力不足。当强反射界面较多、多次波来源可能性较多,多次波波场特征复杂时,无法难以判断其来源层位。
发明内容
针对当前多次波识别方法的不足,本发明实施例提供一种多次波识别方法,包括:
获取目标区域的叠后地震数据和测井数据;
根据所述测井数据确定各反射界面的一次反射系数;
将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面;
根据确定的强反射界面和一次反射系数确定多次反射系数;
根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录;所述合成地震记录包括:一次波的正演记录和包含多次波的正演记录;
在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。
本发明实施例中,所述的测井数据包括:声波时差曲线、密度测井曲线以及深度数据。
本发明实施例中,所述的根据所述测井数据确定目标区域的一次反射系数包括:
根据声波时差曲线确定速度数据;
根据速度数据和深度数据确定双程旅行时间;
根据速度和密度曲线确定波阻抗;
对所述波阻抗按等时时间间隔dt方波化处理确定方波化的波阻抗;
根据确定的方波化阻抗和下式确定一次反射系数;
其中,Rj为一次反射系数;Z1为方波化的波阻抗;j为采样点序号,Z1j+1为第j+1个采样点的方波化阻抗,Z1j第j个采样点的方波化阻抗。
本发明实施例中,所述的根据确定的强反射界面和一次反射系数确定多次反射系数包括:
以等时时间间隔dt划分地层;
利用反射率法确定划分的各地层的上行波能量和下行波能量,将各地层的上行波能量和下行波能量的比值作为全井段包含多次波的反射系数;
利用所述强反射界面筛选所述全井段包含多次波的反射系数确定多次反射系数。
本发明实施例中,所述的根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录包括:
从所述叠后地震数据中提取统计子波;
统计子波与反射系数序列褶积确定正演记录;其中,所述反射系数序列为:一次反射系数、全井段包含多次波的反射系数以及多次反射系数组成的反射系数序列;
将所述正演记录进行归一化处理生成合成地震记录,所述合成地震记录包括:一次波正演记录Sy,包含多次波的正演记录Syn,其中n=0,1...n。
本发明实施例中,所述的在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波包括:
在合成地震记录中创建时间窗口;
确定各时间窗口中,实际地震道的正演记录分别与一次波正演记录Sy、包含多次波的正演记录Sy0的相关系数;
实际地震道的正演记录与正演记录Sy0的相关系数大于实际地震道的正演记录与正演记录Sy时,根据多次波的正演记录Syn中波组特征的变化确定产生多次波的地层。
本发明提供一种剥层法分析和识别多次波的方法和装置,利用反射率法正演模拟,了解地震波传播特点、识别多次波波场信息。通过此方法可以直接针对某些复杂层位,识别出产生的多次波的地层,在常规的叠加剖面上进一步识别出残留的多次波。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明公开的多次波识别方法的流程图;
图2为本发明实施例中的流程图;
图3a为过A井0~4s的叠后地震剖面;
图3b为本发明实施例中A井0~5800m的声波时差曲线AC;
图3c为本发明实施例中A井0~5800m的密度曲线Density;
图4为本发明实施例中地震道及地震记录曲线;
图5为本发明实施例中Rj层的地震波反射路径;
图6为本发明实施例中Rj层的地震波反射路径;
图7本发明实施例中Syn中1.2s到2.4s的地震记录。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供一种多次波识别方法,如图1所示,该方法包括:
步骤S101,获取目标区域的叠后地震数据和测井数据;
步骤S102,根据所述测井数据确定各反射界面的一次反射系数;
步骤S103,将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面;
步骤S104,根据确定的强反射界面和一次反射系数确定多次反射系数;
步骤S105,根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录;所述合成地震记录包括:一次波的正演记录和包含多次波的正演记录;
步骤S106,在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波。
本发明在于通过剥层法识别和分析多次波,通过正演剥层前后含有多次波的正演记录,与一次波的正演记录、实际地震道对比分析,识别该层可能产生的多次波,最后识别出实际地震数据中的残留多次波和多次波来源。
本发明实施例识别多次波的步骤,主要包括:
获取地震数据、声波时差、密度测井曲线。
根据声波时差求取倒数得到速度,并对测井曲线进行深时转换;速度和密度相乘得到波阻抗;对波阻抗激进行等时间隔方波化处理;由反射系数公式计算得到一次反射系数序列。
根据一次反射系数序列的大小排序筛选出强反射地层。
利用获得的一次反射系数序列通过反射率法计算得到全井段包含多次波的反射系数序列;结合筛选的强反射界面,依次去除强反射层及其以上的反射系数,重新通过反射率法计算得到强反射层以下井段包含多次波的反射系数序列。
一次反射系数序列、包含多次波的反射系数序列与子波褶积得到正演记录;从叠后剖面提取测井位置处的实际地震道;对实际地震道、正演记录振幅归一化处理。
创建时窗,分析时窗内实际地震道和正演记录之间的波组特征和互相关系数识别出实际叠后剖面中的残留多次波和来源。
图2为本发明实施例的流程图,如图2所示,该方法包括:
S201,获取地震数据和测井数据,具体为:
1)获取叠后地震数据
2)获取测井资料,包括声波时差和密度曲线
如图3a为过A井0~4s的叠后地震剖面,图3b、图3c分别为A井0~5800m的声波时差曲线AC和密度曲线Density。
S202,计算一次反射系数:
1)测井曲线深时转换:
声波时差曲线AC求取倒数得到速度Vp,并利用深度Depth和速度Vp关系求得双程旅行时间。
过程如下:
load Depth.txt;%加载深度数据
load AC.txt;%加载声波时差数据
Vp=1000000./AC;
n=length(Depth);
Time(1)=2*Depth(1)/Vp(1);
for i=2:n
Time(i)=Time(i-1)+4*(Depth(i)-Depth(i-1))/(Vp(i)+Vp(i-1));
end
其中,n记录深度Depth数据个数,Time记录双程旅行时间数据。
2)计算方波化波阻抗曲线:
由于声波和密度测井测量方法和仪器所限,测井曲线会包含较多异常值,存在较多高频噪声。因此需要对测井曲线进行方波化处理,消除薄层和异常值的干扰,保持厚层界面信息,使测井曲线频带接近地震频带。
速度Vp和密度Density相乘计算得到波阻抗。
Z=Vp×Density
对波阻抗Z按等时时间间隔dt进行方波化处理得到的波阻抗Z1,并得到方波化的时间Time1。
其中N为方波化时间间隔dt对应的方波化长度,Z1为方波化波阻抗;DT、ρ分别为地层声波时差和密度;N为方波化间隔dt对应的方波化长度。
如图4所示,a为波阻抗Z,b为方波化的波阻抗Z1。
本实例中的方波化间隔dt为0.004s,过程如下:
其中,m记录方波化个数,mean为方波化函数。
3)计算一次反射系数:
图4中的c为由公式计算得到的一次反射系数R,Z1为方波化的波阻抗;j为采样点序号。
程序实现过程如下:
for i=1:m-2
R(i,1)=(Z1(i+1)-Z1(i))/(Z1(i+1)+Z1(i));
end
S203,筛选强反射界面:
地层一次反射系数R从大到小排序并记录对应方波化的时间值Time1。本实例中的排序利用Matlab中的sort函数。
[RR index]=sort(R(:,1));
Time1(index);
其中RR为排序后的反射系数,Time1(index)为对应的时间。
如表1所示,本实施例中A井地层反射系数从大到小依次为0.2432、0.2211、0.2024、0.1918、0.1337,……,具体的,对应的强反射源界面Time1时间值分别为1.9160s对应奥陶系底O1、1.7800s对应上二叠底P2l、1.1240s对应上三叠系须家河底T3x1、0.8000s对应侏罗系底J1、1.5040s对应下三叠系飞四段底T1f4、2.1400s对应寒武系筇竹寺底E1q。
表1
序号 | R | Time1(s) |
1 | 0.2432 | 1.9160 |
2 | 0.2211 | 1.7800 |
3 | 0.2024 | 1.1240 |
4 | 0.1918 | 0.8000 |
5 | 0.1662 | 1.5040 |
6 | 0.1337 | 2.1400 |
…… | …… | …… |
本发明实施例中,筛选出大于2/3最大反射系数的界面作为强反射界面,从而得到了5个强反射界面,按Time1从小到大排序分别是:T1=0.8000s,T2=1.1240s,T3=1.5040s,T4=1.7800s,T5=1.9160s。
S204,计算多次反射系数:
以方波化时间间隔dt为间隔,将地层划分为m层(方波化个数),每层用Time1表示。
图5、图6为Rj层(对应时间Time1(j)层)的地震波反射路径,Rj表示反射界面的一次波下行反射系数,利用反射率法求多次波的反射系数即求取图中Rj层的地表的上行波能量u。
一般地,对于任意反射系数为Rj的反射层,下行波能量的公式为:
dj+1=-Rj*u+(1-Rj)*dj
上行波能量的公式为:
u=Rj*dj+(1+Rj)*u
由此得到每个Time1层的上行波能量u,将上行波能量与下行波能量的比值作为多次反射系数。
如图4的d到i所示,根据一次反射系数R由反射率法计算得到全井段包含多次波的反射系数R0;结合筛选强反射界面T1,去除T1层及其以上的一次反射系数R,重新由反射率法计算得到T1层以下井段包含多次波的反射系数R1,同理依次去除T2、T3、T4和T5层及其以上的反射系数分别得到R2,R3,R4,R5。
具体实现过程如下:
测井数据:深度声波时差密度都是数组;
声波时差DT求倒数得到速度速度Vp乘以密度ρ得到:
波阻抗
速度Vp和深度Depth得到时间
每个时间点i都对应一个波阻抗i=1~n;
方波化:因为不是等时间隔采样,Time1表示的是每层的包括所有层,以dt为时间间隔重采样得到:
每个时间点j都对应一个波阻抗j=1~m;
由公式,通过Z1计算一次反射系数得到
进一步,关于多次反射系数计算:
更新上行波公式:u=Rj*dj+(1+Rj)*u;
更新下行波公式:dj+1=-Rj*u+(1-Rj)*dj;
初始条件:下行波能量为1,上行波能量为0,即d1=1,u=0;
更新上行波公式:u=R1*d1+(1+R1)*u=R1
依次计算下去得到整个
多次反射系数R1的计算:比如本实施例中,一次反射系数中R010为最大值,是最大的反射界面,去掉R010及以上的数据得到:这时只有m-10个数据了。
多次反射系数
初始条件:下行波能量为1,上行波能量为0,即d1=1,u=0;
更新上行波公式:u=R11*d1+(1+R11)*u=R11
依次计算下去得到整个
代码过程如下:
其中,startT为起始计算时间。
S205,计算合成地震记录:
1)提取地震子波,合成地震记录:
以方波化时间间隔dt为采样间隔,从叠后地震数据中提取A井处实际地震道S,并从地震道S中提取出统计子波W;反射系数序列R、R0、R1、R2、R3、R4、R5与地震子波W褶积得到的正演记录Sy、Sy0、Sy1、Sy2、Sy3、Sy4、Sy5。
Syn=Rn*W(n=,0,1,2,3,4,5)
Syn为测井合成地震记录;W为地震子波;Rn为反射系数。
2)地震道能量归一化处理:
将实际地震道S及正演记录Sy、Sy0、Sy1、Sy2、Sy3、Sy4、Sy5的振幅分别除以该道振幅绝对值的最大值,更新原记录。
Syn=[S,Syn]
Syn=Syn./max(abs(Syn))
本实例中,提取实际地震道、统计子波借用了地震反演软件HRS的Strata模块,褶积运算借用了Matlab中的conv函数,过程如下:
load W.txt;%加载统计子波
for i=1:length(startT)+1
Syn(:,i)=conv(R(:,i),W);
end
load S.txt;%加载实际地震道
Syn=[S,Syn];
Syn=Syn./max(abs(Syn));
如图7所示,为Syn中1.2s到2.4s的地震记录;
S106,识别多次波:
1)确定多次波时窗:
在地震记录Syn中,创建时间窗口,对比时间窗口中一次波的正演记录Sy和实际地震道S的波组特征,并求互相关系数。理论上,一次波的正演记录Sy和实际地震道S二者波组特征越相近,互相关系数越大,地震道多次波干扰越少。反之,二者差异越大,互相关系数越小,多次波干扰越严重,甚至一次波能量被多次波掩盖。因此,二者匹配不好相关性差的地方说明地震数据中多次波发育。
如图7所示,1.2s~2.0s的时间范围内,一次波合成地震记录Sy与实际地震道S符合很好,说明该段实际地震道以一次有效波为主。2.0s~2.4s的时间范围内,一次波合成地震记录Sy与实际地震道S波组特征有较大差异,而含有多次波的正演记录Sy0与实际地震道S符合更好。具体表现在强反射振幅的个数不同,实际地震道S、含有多次波的正演记录Sy0明显有多个强能量振幅。说明2.0s~2.4s内,除一次波反射之外,还包含有较强的多个多次波反射,实际地震资料在该段时间范围内可能残留了较强的多次波干扰。
在本实例1.2s到2.4s的地震记录中,以0.2s为间隔创建6个时间窗口,并求每个窗口中实际地震道S与一次波合成地震记录Sy、包含多次波的正演记录Sy0的互相关系数。
本实例中的相关运算借用了Matlab中的corrcoef函数。如表2所示,2.0~2.2s,2.2~2.4s的时间窗口,一次波正演记录Sy和实际地震道S的互相关系数明显偏小,反观包含多次波的正演记录Sy0和实际地震道S的互相关系数更大,说明在2.0~2.2s,2.2~2.4s的时间窗口内,多次波干扰越严重。
表2
综上,在本实施例中,实际地震数据中多次波时窗范围为:2.0~2.2s,2.2~2.4s。
2)确定产生多次波的来源地层:
识别多次波特征在于,在同一时窗内,对比合成地震记录和实际地震道的波组特征和相关性。理论上,剥层前包含多次波的正演记录和实际地震道二者波组特征更相近,而剥层后包含多次波的正演记录和一次合成地震道二者波组特征更相近。
理论上,在同一时窗内,剥层前包含多次波的正演记录和实际地震道S二者波组特征更相近,而剥层后包含多次波的正演记录和一次合成地震道Sy0二者波组特征更相近。
如图7所示,图7的a为采用本发明实施例的A井处实际地震道S,图7的b到h为采用本发明实施例的A井的合成地震记录。在2.0~2.2s时窗范围内2.1s附近,含有多次波的合成地震记录Sy1较Sy0,波组特征上看无明显变化,说明该层产生的多次波相对较弱,不是主要的多次波来源层,而Sy2较Sy1变化却比较明显,到Sy3、Sy4、Sy5已无明显变化,且接近一次波合成地震记录S的波组特征。说明T2层位是产生2.1s处多次波的主要地层。
同理,在2.2~2.4s时窗范围内2.3s附近的多次波来源于T3层,2.35s附近的多次波来源于T1层。
依次求取每个时窗中,正演记录Sy0、Sy1、Sy2、Sy3、Sy4、Sy5与实际地震道S、一次合成地震道Sy的互相关系数。如表3所示,随着剥层的进行,正演记录Sy0、Sy1、Sy2、Sy3、Sy4、Sy5与实际地震道S的互相关系数在减小,而与一次合成地震道Sy的互相关系数在增大,剥去T2层后corrcoef(S,Sy2)明显减小,而corrcoef(Sy,Sy2)明显增大,说明2.0~2.2s窗口中的多次波来源于T2层。同理,如表4所示,2.2~2.4s窗口中的多次波主要来源于T1,T3层。
表3 2.0~2.2s时窗内的互相系数表
表4 2.2~2.4s时窗内的互相系数表
corrcoef(,) | S | Sy |
Sy0 | 0.8420 | 0.7178 |
Sy1 | 0.7532 | 0.7832 |
Sy2 | 0.7602 | 0.8375 |
Sy3 | 0.7036 | 0.9110 |
Sy4 | 0.6809 | 0.9544 |
Sy5 | 0.6422 | 0.9886 |
综上,在本实施例中,2.0~2.2s时窗内的多次波主要来源于T2地层,即上三叠系须家河底T3x1。2.2~2.4s时窗内的多次波主要来源于T1和T3地层,即侏罗系底J1和下三叠系飞四段底T1f4。
本发明可适用于地震解释工作中,识别中深层具有强反射界面的地区产生的层间多次波。这些强反射界面通常是基岩面、不整合面、火成岩面等,反射波在返回地面遇到这些界面时,可能再次反射后向下传播,形成多次波,这些多次波与深层一次波能量相当,混在一起难以分离。随着勘探的深入,更深深度地层逐渐进入开发阶段,由于深层多层波的干扰,面临高产井地震响应不明确、确立高产井模式困难、有效储层预测的问题。
本发明可分析产生层间多次波的来源层位,使得多次波的压制针对性更强,一定程度上减少了压制多次波的不确定性,更易分离和消除多次波干扰。
本发明通过反射率法正演模拟和波组对比,剥层分析和识别多次波。目前叠后阶段识别多次波的方法相对单一,主要依靠地震剖面上地震反射的外观、形态判断,主观性较强,多解性强,识别的结果说服力不足。本发明采用和测井、地震解释相结合的方法,能实现半定量的多次波识别,从而大大减少了多次波识别的不确定性,更好地分离和消除多次波干扰波。
本发明实施例中,一次波和多次波合成地震记录正演,通过计算合成地震记录和实际地震道的相关系数,实现半定量的多次波识别;本发明可以直接针对某些复杂层位,识别出产生多次波的地层,极大地减少无关数据的分析和运算。
应当注意,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括””、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的单元、装置或模块等,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种多次波识别方法,其特征在于,所述的方法包括:
获取目标区域的叠后地震数据和测井数据;
根据所述测井数据确定各反射界面的一次反射系数;
将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面;
根据确定的强反射界面和一次反射系数确定多次反射系数;
根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录;所述合成地震记录包括:一次波的正演记录和包含多次波的正演记录;
在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波;其中,
所述的根据确定的强反射界面和一次反射系数确定多次反射系数包括:
以等时时间间隔dt划分地层;
利用反射率法确定划分的各地层的上行波能量和下行波能量,将各地层的上行波能量和下行波能量的比值作为全井段包含多次波的反射系数;
利用所述强反射界面筛选所述全井段包含多次波的反射系数确定多次反射系数;
所述的根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录包括:
从所述叠后地震数据中提取统计子波;
统计子波与反射系数序列褶积确定正演记录;其中,所述反射系数序列为:一次反射系数、全井段包含多次波的反射系数以及多次反射系数组成的反射系数序列;
将所述正演记录进行归一化处理生成合成地震记录,所述合成地震记录包括:一次波正演记录Sy,包含多次波的正演记录Syn,其中n=0,1...n;
所述的在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波包括:
在合成地震记录中创建时间窗口;
确定各时间窗口中,实际地震道的正演记录分别与一次波正演记录Sy、包含多次波的正演记录Sy0的相关系数;
实际地震道的正演记录与正演记录Sy0的相关系数大于实际地震道的正演记录与正演记录Sy时,根据多次波的正演记录Syn中波组特征的变化确定产生多次波的地层。
2.如权利要求1所述的多次波识别方法,其特征在于,所述的测井数据包括:声波时差曲线、密度测井曲线以及深度数据。
3.如权利要求2所述的多次波识别方法,其特征在于,所述的根据所述测井数据确定目标区域的一次反射系数包括:
根据声波时差曲线确定速度数据;
根据速度数据和深度数据确定双程旅行时间;
根据速度和密度曲线确定波阻抗;
对所述波阻抗按等时时间间隔dt方波化处理确定方波化的波阻抗;
根据确定的方波化阻抗和下式确定一次反射系数;
其中,Rj为一次反射系数;Z1为方波化的波阻抗;j为采样点序号,Z1j+1为第j+1个采样点的方波化阻抗,Z1j第j个采样点的方波化阻抗。
4.一种多次波识别装置,其特征在于,所述的装置包括:
数据获取模块,用于获取目标区域的叠后地震数据和测井数据;
一次反射系数确定模块,用于根据所述测井数据确定各反射界面的一次反射系数;
强反射界面确定模块,用于将所述一次反射系数中大于预设的反射系数阈值的反射界面作为强反射界面;
多次反射系数确定模块,用于根据确定的强反射界面和一次反射系数确定多次反射系数;
合成地震记录生成模块,用于根据所述一次反射系数、多次反射系数及叠后地震数据确定合成地震记录;所述合成地震记录包括:一次波的正演记录和包含多次波的正演记录;
识别模块,用于在所述合成地震记录中创建时窗窗口,根据各时窗中实际地震道与合成地震记录中各正演记录的相关系数识别多次波;其中,
所述的多次反射系数确定模块包括:
地层划分单元,用于以等时时间间隔dt划分地层;
全井段反射系数确定单元,用于利用反射率法确定划分的各地层的上行波能量和下行波能量,将各地层的上行波能量和下行波能量的比值作为全井段包含多次波的反射系数;
筛选单元,用于利用所述强反射界面筛选所述全井段包含多次波的反射系数确定多次反射系数;
所述的合成地震记录生成模块包括:
子波提取单元,用于从所述叠后地震数据中提取统计子波;
褶积单元,统计子波与反射系数序列褶积确定正演记录;其中,所述反射系数序列为:一次反射系数、全井段包含多次波的反射系数以及多次反射系数组成的反射系数序列;
归一单元,用于将所述正演记录进行归一化处理生成合成地震记录,所述合成地震记录包括:一次波正演记录Sy,包含多次波的正演记录Syn,其中n=0,1...n;
所述的识别模块包括:
时窗创建单元,在合成地震记录中创建时间窗口;
相关系数确定单元,用于确定各时间窗口中,实际地震道的正演记录分别与一次波正演记录Sy、包含多次波的正演记录Sy0的相关系数;
地层识别单元,实际地震道的正演记录与正演记录Sy0的相关系数大于实际地震道的正演记录与正演记录Sy时,根据多次波的正演记录Syn中波组特征的变化确定产生多次波的地层。
5.如权利要求4所述的多次波识别装置,其特征在于,所述的测井数据包括:声波时差曲线、密度测井曲线以及深度数据。
6.如权利要求5所述的多次波识别装置,其特征在于,所述的一次反射系数确定模块包括:
速度确定单元,用于根据声波时差曲线确定速度数据;
时间确定单元,用于根据速度数据和深度数据确定双程旅行时间;
波阻抗确定单元,根据速度和密度曲线确定波阻抗;
方波化单元,用于对所述波阻抗按等时时间间隔dt方波化处理确定方波化的波阻抗;
一次反射系数确定单元,用于根据确定的方波化阻抗和下式确定一次反射系数;
其中,Rj为一次反射系数;Z1为方波化的波阻抗;j为采样点序号,Z1j+1为第j+1个采样点的方波化阻抗,Z1j第j个采样点的方波化阻抗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810579839.3A CN108828664B (zh) | 2018-06-07 | 2018-06-07 | 一种多次波识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810579839.3A CN108828664B (zh) | 2018-06-07 | 2018-06-07 | 一种多次波识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108828664A CN108828664A (zh) | 2018-11-16 |
CN108828664B true CN108828664B (zh) | 2019-12-10 |
Family
ID=64144344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810579839.3A Active CN108828664B (zh) | 2018-06-07 | 2018-06-07 | 一种多次波识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108828664B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111610563B (zh) * | 2019-02-26 | 2023-02-28 | 中国石油天然气股份有限公司 | 一种识别多次波的方法和装置 |
CN111308554B (zh) * | 2020-03-10 | 2021-03-16 | 清华大学 | 一种基于多次波生成层自适应提取的层间多次波预测方法 |
CN118011472A (zh) * | 2022-11-08 | 2024-05-10 | 中国石油天然气集团有限公司 | 一种超深层地震多次波确定方法、装置及计算机设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000041005A1 (en) * | 1999-01-08 | 2000-07-13 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
CN103376466A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | 一种多次波压制方法 |
CN103460074A (zh) * | 2011-03-31 | 2013-12-18 | 埃克森美孚上游研究公司 | 全波场反演中小波估计和多次波预测的方法 |
CN104502974A (zh) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | 一种压制多次波的组合方法及装置 |
GB2521598A (en) * | 2013-12-02 | 2015-07-01 | Statoil Petroleum As | Multi-dimensional deconvolution using exact boundary conditions |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
MXPA06011579A (es) * | 2004-04-07 | 2008-03-11 | Westerngeco Seismic Holdings Ltd | Prediccion generalizada de multiplos de superfie 3-d. |
-
2018
- 2018-06-07 CN CN201810579839.3A patent/CN108828664B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000041005A1 (en) * | 1999-01-08 | 2000-07-13 | Exxonmobil Upstream Research Company | Method for identifying and removing multiples from seismic reflection data |
CN101598809A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种自适应消除线性规则噪声以及多次波干扰的方法 |
CN101630016A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种提高垂直地震剖面成像质量的方法 |
CN103460074A (zh) * | 2011-03-31 | 2013-12-18 | 埃克森美孚上游研究公司 | 全波场反演中小波估计和多次波预测的方法 |
CN103376466A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | 一种多次波压制方法 |
GB2521598A (en) * | 2013-12-02 | 2015-07-01 | Statoil Petroleum As | Multi-dimensional deconvolution using exact boundary conditions |
CN104502974A (zh) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | 一种压制多次波的组合方法及装置 |
Non-Patent Citations (1)
Title |
---|
"地震多次波的几种常见压制方法";魏晨成;《内蒙古石油化工》;20171231(第11、12期);第80-81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108828664A (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8706420B2 (en) | Seismic fluid prediction via expanded AVO anomalies | |
EP3631523B1 (en) | System and method for predicting fault seal from seismic data | |
US7433265B2 (en) | Converted wave energy removal from seismic data | |
US6789018B1 (en) | Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave | |
CN108828664B (zh) | 一种多次波识别方法及装置 | |
CN109115987B (zh) | 一种基于岩石物理模型的流体因子的评价方法及装置 | |
EP0186569B1 (en) | Entropy guided deconvolution of seismic signal | |
Roots et al. | Interferometric seismic imaging around the active Lalor mine in the Flin Flon greenstone belt, Canada | |
Girard et al. | Automated processing strategies for ambient seismic data | |
Maurya et al. | Reservoir characterization using model based inversion and probabilistic neural network | |
CN113687440A (zh) | 基于米兰科维奇旋回的古水深定量恢复方法和存储介质 | |
EP1554607B1 (en) | A method of reconstructing seismic records to obtain high resolution signals | |
AU2018272313A1 (en) | System and method for predicting fault seal from seismic data | |
CN108254785B (zh) | 储层确定方法和装置 | |
MX2014001334A (es) | Sistemas y metodos para detectar ruido del oleaje en un recolector sismico. | |
Pitarka et al. | Using dense array waveform correlations to build a velocity model with stochastic variability | |
CN112130209B (zh) | 岩溶储层预测方法及装置 | |
CN111694055B (zh) | 一种强反射屏蔽下avo属性提取和分析的方法和装置 | |
CN114462703A (zh) | 声波参数曲线预测、测井曲线预测方法及电子设备 | |
EP0397313A2 (en) | Comprehensive system for evaluation of velocities of acoustical energy of sedimentary rocks | |
CN111312272A (zh) | 一种用于减少近井声数据集中噪声信号的产品、方法及系统 | |
CN114966828A (zh) | 一种断溶体储层地震反射特征增强方法 | |
Conticini | Seismic facies quantitative analysis: New tool in stratigraphic interpretation | |
CN113640873B (zh) | 砂岩型铀矿叠前地震预测方法及装置 | |
Orji et al. | A Cepstral Technique for Unmasking Vital Sub-Seismic Stratigraphic Details Embedded in Dense 3D Seismic Data from Niger Delta |
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 |