CN102879814B - 一种准确的深度域层速度更新方法 - Google Patents
一种准确的深度域层速度更新方法 Download PDFInfo
- Publication number
- CN102879814B CN102879814B CN201110198587.8A CN201110198587A CN102879814B CN 102879814 B CN102879814 B CN 102879814B CN 201110198587 A CN201110198587 A CN 201110198587A CN 102879814 B CN102879814 B CN 102879814B
- Authority
- CN
- China
- Prior art keywords
- velocity
- depth
- root
- layer
- speed
- 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
技术领域
本发明涉及反射波地震数据处理技术,具体的说是一种准确的深度域层速度更新方法。
背景技术
叠前深度偏移技术是现阶段生产中非常重要的地震资料处理技术之一,它能提供深度域中我们能看到的最为准确的地震构造成像资料,因其能直接和井数据建立联系所以在发现规模油田中有着举足轻重的作用。这项技术大规模推广应用的前提是需要建立可靠的深度域层速度场(描述地下介质速度的数据体),叠前深度偏移速度场的建立需要借助偏移(这里的偏移指的是叠前深度偏移,下同)速度分析来完成,因此,偏移速度分析方法的精度和有效性也成为影响偏移结果的重要因素。高精度的偏移速度来源于对初始偏移速度的剩余速度分析,所谓的剩余速度分析是指在已有速度的基础上进一步提高速度精度的速度分析方法。
目前,国内叠前深度偏移层速度分析主要依据层速度扫描的方法实现,该方法因交互操作难以实现、层位定位不准、计算量大、速度更新周期过长等缺点已经被工业界淘汰。国外的深度偏移速度建模软件(例如GeoDepth软件)中已经出现了利用垂直剩余速度谱分析的方式更新层速度,而该项技术在国内还没有出现。
本发明基于中等倾角假设计算剩余速度谱,然后依据剩余速度值更新深度域层速度值。该方法具有计算效率高、更新效果好、实现简单的优点。
发明内容
本发明的目的是提供一种适用于陆相沉积盆地构造成像的深度域层速度更新方法。
本发明通过以下具体步骤实现:
1)将野外地震采集数据经过常规的数据预处理后进行叠前深度偏移处理,得到共反射点道集数据;
2)将共反射点道集数据依次进行切除、能量均衡、道插值和带通滤波处理;
步骤2)中的切除是对浅层拉伸严重的道集和能量异常道进行切除,能量均衡是对浅层与深层、近偏移距与远偏移距道集进行均衡,道插值和带通滤波是带通滤波器的滤波门与共反射点道集的优势频带保持一致。
3)将深度域层速度场利用DIX(C.Hewitt Dix提出的计算公式)公式转化成均方根速度场;
其中,Vint,i为第i层的层速度,n为层数,i为层变量,变化范围从1到n;ti为第i层的垂直双程时间,VR,n为第n层的均方根速度。
4)取均方根速度场每一固定深度间隔共深度点位置的速度值;
5)对每一个共深度点位置的均方根速度场和共反射点道集(CRP)应用以下公式计算偏移剩余速度谱;
其中,ΔZh为深度修正量,v为背景速度,v′为修正速度,z为当前分析点的深度,h为偏移距。该公式直接在深度域中应用,而不需要将其变化到时间域中;
6)对偏移剩余速度谱进行能量均衡处理;
步骤6)所述的能量均衡处理是指浅层与深层的能量进行移动时窗加权处理。
7)利用已有的垂直剩余速度谱确定每一深度位置能量最大的剩余速度值;
8)利用该剩余速度值更新均方根速度场(两者相加),得新的均方根速度场;
9)将更新后的均方根速度场利用如下公式转化成对应位置处的深度域层速度场;
t0,n和t0,n-1分别为第n层和第n-1层的垂直双程传播时间,VR,n和VR,n-1分别为第n层和第n-1层的均方根速度。
步骤9)所述的转化成对应位置处的深度域层速度场时每一层的转化都要结合上覆层的数据,上覆层的速度更新传递到下一层。
比如,已知第三层的均方根速度VR,3和剩余速度ΔV3,则更新的新均方根速度VRn,3=VR,3+ΔV3。第二层的均方根速度VR,2以及对应的时间t0,3、t0,2利用步骤9)中的公式可以计算出第三层的新层速度Vint n,3;此时由于第三层的层速度Vint,3已经更新为Vint n,3,所以必须利用步骤3)中的公式(将Vint n,3代入计算而不是Vint,3)计算第四层的均方根速度VR,4,再加上第四层的剩余速度ΔV4得到第四层新的均方根速度VRn,4=VR,4+ΔV4。再利用第三层的均方根速度VRn,3以及对应的时间t0,4、t0,3和步骤9)中的公式可以计算出第四层的新层速度Vint n,4;如此重复,再进行第五层的速度更新。
10)利用更新后的层速度场进行叠前深度偏移处理得到准确的深度域构造地震成像结果。
本发明同相轴的深距曲线为抛物线加高阶修正量,得到分辨率较高的剩余速度谱,为更准确的叠前深度偏移构造成像提供了可靠依据。
本发明的目的是提供一种适用于陆相沉积盆地构造成像的深度域层速度更新方法,该方法可以为得到准确的地震构造成像方法提供保障。
附图说明
图1为某实际资料某一位置处处理前后的剩余速度谱,其中(a)为未经处理的共反射点道集,(b)为经过处理后的共反射点道集;
图2为初始层速度场;
图3为初始层速度场经过转化后的均方根速度场;
图4利用均方根速度场和共反射点道集得到的高精度剩余速度谱;
图5为利用剩余速度谱得到的剩余速度值更新后的均方根速度场;
图6为将更新后的均方根速度场转化为更新后的层速度场;
图7利用更新后的层速度场进行叠前深度偏移处理后得到的偏移剖面。
具体实施方式
以下结合附图具体说明本发明。
本发明提供了一种准确的深度域层速度更新方法,其具体实施方式为:
1)将野外地震采集数据经过常规的数据预处理后进行叠前深度偏移处理,得到共反射点道集数据,如图1(a)所示,为某一个反射点的共反射点道集数据;
2)将共反射点道集(CRP)数据依次进行切除、能量均衡、道插值和带通滤波处理,如图1(b)所示,为经过处理后的共反射点道集;
步骤2)中的切除是指对浅层拉伸严重的道集和能量异常道进行切除,能量均衡是指对浅层与深层、近偏移距与远偏移距道集数据进行均衡,道插值和带通滤波是指带通滤波器的滤波门与CRP道集的优势频带保持一致。
3)将深度域层速度场(如图2)利用DIX公式转化成均方根速度场(如图3);
4)取均方根速度场每一固定深度间隔共深度点位置的速度值;
5)对每一个共深度点位置的均方根速度场和共反射点道集(CRP)应用以下公式计算偏移剩余速度谱;
其中,ΔZh为深度修正量,v为背景速度,v′为修正速度,z为当前分析点的深度,h为偏移距。该公式直接在深度域中应用,而不需要将其变化到时间域中,如图4中的剩余速度谱纵坐标为深度;
6)对偏移剩余速度谱进行能量均衡处理,如图4所示,为某一位置处的剩余速度谱;
步骤6)所述的能量均衡处理是指浅层与深层的能量进行移动时窗加权处理。
7)利用已有的垂直剩余速度谱确定每一深度位置能量最大的剩余速度值,
8)利用该剩余速度值更新均方根速度场(两者相加),得到更新后的均方根速度场,如图5为更新后的均方根速度场;
9)将更新后的均方根速度场转化成对应位置处的深度域层速度场,如图6为更新后的层速度场;
步骤9)所述的均方根速度转层速度,每一层的转化都需要考虑上覆层的影响,上覆层的速度更新必然会传递到下一层。
10)利用更新后的层速度场进行叠前深度偏移处理得到准确的深度域构造地震成像结果,如图7所示为准确的深度域构造成像结果。
本发明应用实验证实效果良好。
(1)本发明利用计算精度较高的剩余速度计算公式进行剩余速度谱的计算,该公式具有以下特点:①同相轴的深度曲线为抛物线加高阶修正量;②同相轴平直时,成像深度为真深度。
(2)应用本发明能得到精度较高的剩余速度场,同时将更新的均方根速度场转化为更新层速度场时可以同时进行多层速度更新,提高了速度更新的效率。应用更新后的层速度场进行叠前深度偏移即可得到非常准确的构造地震成像结果。
(3)本发明对输入的CRP道集的切除、能量均衡、道插值和带通滤波处理等措施都为得到分辨率较高的剩余速度谱打下良好的基础。
本发明的剩余速度计算公式精度更高,而且计算效率没有明显降低。公式剩余时差的深距曲线为抛物线加高阶修正量,从速度谱的计算得到分辨率较高的剩余速度谱,为复杂偏移速度场的建立提供了可靠方法,进而为偏移构造成像提供了可靠依据。
Claims (1)
1.一种深度域层速度更新方法,其特征在于具体步骤包括:
1)野外地震采集数据,经过常规的数据预处理后进行叠前深度偏移处理,得到共反射点道集数据;
2)共反射点道集数据依次进行切除、能量均衡、道插值和带通滤波处理;
其中所述的的切除是对浅层拉伸严重的道集和能量异常道进行切除;
其中所述的能量均衡是对浅层与深层、近偏移距与远偏移距道集进行均衡;
其中所述的道插值和带通滤波是带通滤波器的滤波门与共反射点道集的优势频带保持一致;
3)将深度域层速度场利用DIX公式转化成均方根速度场;
其中,Vint,i为第i层的层速度,n为层数,i为层变量,变化范围从1到n,ti为第i层的垂直双程时间,VR,n为第n层的均方根速度;
4)取均方根速度场每一固定深度间隔共深度点位置的速度值;
5)对每一个共深度点位置的均方根速度场和共反射点道集用以下公式计算偏移剩余速度谱;
其中,ΔZh为深度修正量,v为背景速度,v'为修正速度,z为当前分析点的深度,h为偏移距;
6)对偏移剩余速度谱进行能量均衡处理;
7)利用已有的垂直剩余速度谱确定每一深度位置能量最大的剩余速度值;
8)利用该剩余速度值与均方根速度场相加,得到该层的新均方根速度;
9)将更新均方根速度场利用以下公式转化成对应位置处的深度域层速度场;
其中,Vint,n为第n层的层速度,
t0,n和t0,n-1分别为第n层和第n-1层的垂直双程传播时间,VR,n和VR,n-1分别为第n层和第n-1层的均方根速度;
10)利用更新后的层速度场进行叠前深度偏移处理得到深度域构造地震成像结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110198587.8A CN102879814B (zh) | 2011-07-15 | 2011-07-15 | 一种准确的深度域层速度更新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110198587.8A CN102879814B (zh) | 2011-07-15 | 2011-07-15 | 一种准确的深度域层速度更新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102879814A CN102879814A (zh) | 2013-01-16 |
CN102879814B true CN102879814B (zh) | 2015-05-13 |
Family
ID=47481207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110198587.8A Active CN102879814B (zh) | 2011-07-15 | 2011-07-15 | 一种准确的深度域层速度更新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102879814B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412329B (zh) * | 2013-08-06 | 2015-09-16 | 中国海洋石油总公司 | 一种提高地震数据分辨率的方法 |
CN103592682B (zh) * | 2013-10-15 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种地震波场能量均衡方法及处理装置 |
CN104570103B (zh) * | 2013-10-16 | 2017-05-03 | 中国石油化工股份有限公司 | 一种低信噪比地震资料的井约束速度谱拾取方法 |
CN105319589B (zh) * | 2014-07-25 | 2018-10-02 | 中国石油化工股份有限公司 | 一种利用局部同相轴斜率的全自动立体层析反演方法 |
CN106483559B (zh) * | 2015-08-24 | 2018-11-23 | 中国石油化工股份有限公司 | 一种地下速度模型的构建方法 |
CN111142159B (zh) * | 2018-11-02 | 2022-04-08 | 中国石油化工股份有限公司 | 中深层异常体地震速度建模方法及系统 |
CN109738945B (zh) * | 2018-11-08 | 2021-01-19 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前深度偏移成果直接生成构造图的方法 |
CN109884693B (zh) * | 2019-01-18 | 2021-10-15 | 潜能恒信能源技术股份有限公司 | 自适应走向速度谱求取方法和系统 |
CN112711064B (zh) * | 2019-10-24 | 2022-11-04 | 中国石油天然气股份有限公司 | 初始层速度建模方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5629904A (en) * | 1994-11-30 | 1997-05-13 | Paradigm Geophysical, Ltd. | Migration process using a model based aperture technique |
CN1797033A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 利用均方根速度提高叠前偏移成像精度的方法 |
US7525874B2 (en) * | 2006-04-21 | 2009-04-28 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
-
2011
- 2011-07-15 CN CN201110198587.8A patent/CN102879814B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5629904A (en) * | 1994-11-30 | 1997-05-13 | Paradigm Geophysical, Ltd. | Migration process using a model based aperture technique |
CN1797033A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 利用均方根速度提高叠前偏移成像精度的方法 |
US7525874B2 (en) * | 2006-04-21 | 2009-04-28 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
Also Published As
Publication number | Publication date |
---|---|
CN102879814A (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879814B (zh) | 一种准确的深度域层速度更新方法 | |
CN103487835B (zh) | 一种基于模型约束的多分辨率波阻抗反演方法 | |
CN101661112B (zh) | 一种准确的地震构造成像方法 | |
CN101251604B (zh) | 二参数转换波速度分析及动校正方法 | |
CN110348135B (zh) | 一种随钻声波测井评价地层渗透率的方法 | |
CN103424777B (zh) | 一种提高地震成像分辨率的方法 | |
CN110318740B (zh) | 一种随钻声波测井评价地层各向异性的方法 | |
CN102540252B (zh) | 基于互相关的高精度中值叠加方法 | |
CN104050147B (zh) | 将时域信号转换成频域信号的方法与系统 | |
CN109444959B (zh) | 全频高精度层速度场建立方法 | |
CN111158049A (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
CN102053267A (zh) | 一种地震剖面资料处理中基于参数反演的vsp波场分离方法 | |
CN104459763A (zh) | 紧支集小波探测地下空洞位置的方法及系统 | |
CN105093278A (zh) | 基于激发主能量优化算法的全波形反演梯度算子的提取方法 | |
CN102798896B (zh) | 一种阵列感应测井仪器的测井信号合成处理方法及其系统 | |
CN101738636A (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
CN108828661B (zh) | 基于地震脉冲响应谱测定场地卓越周期的方法 | |
CN103576197A (zh) | 一种转换波角道集抽取方法 | |
CN105277986A (zh) | 基于自适应匹配滤波算子的可控震源谐波压制方法 | |
CN102043167B (zh) | 一种曲面地表地震资料速度分析的方法 | |
CN103675900A (zh) | 一种确定转换波叠前时间偏移最佳速度剖面的方法 | |
CN107015273B (zh) | 一种基于频率域的煤层Love型槽波压缩方法及装置 | |
CN102914797A (zh) | 一种获得地层各向异性系数的方法及装置 | |
Wu et al. | Analysis and research on blasting network delay of deep-buried diversion tunnel crossing fault zone based on EP-CEEMDAN-INHT | |
CN104678439A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |