CN104570105B - 一种多震源地震资料处理中的速度分析方法 - Google Patents
一种多震源地震资料处理中的速度分析方法 Download PDFInfo
- Publication number
- CN104570105B CN104570105B CN201310499331.XA CN201310499331A CN104570105B CN 104570105 B CN104570105 B CN 104570105B CN 201310499331 A CN201310499331 A CN 201310499331A CN 104570105 B CN104570105 B CN 104570105B
- Authority
- CN
- China
- Prior art keywords
- speed
- velocity
- score
- seism
- carried out
- 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
技术领域
本发明属于石油勘探领域,具体涉及一种多震源地震资料处理中的速度分析方法,用于油气勘探中多震源同步激发(可控震源同步扫描)地震资料的处理。
背景技术
多震源同步激发地震技术是一种高效地震采集技术,它在利用常规采集技术采集一炮资料的时间采集多炮地震记录,其对采集效率的提高是非常明显的。由于在一个很短的时间内有多个震源同步激发,同步震源信号之间的相互干涉是难免的,从而导致资料的信噪比非常低,常规的由速度谱分析获取速度的方法在这里由于信噪比过低而无从解释。而常规地震资料处理中,针对低信噪比资料的速度分析可以采用速度扫描的方法,对整个地区主要层位的叠加速度有个基本了解的基础上再进行速度分析。而多震源地震资料也是一种特低信噪比的资料,因此也可以利用这种技术路线进行速度分析,从而为多震源地震资料的处理奠定基础。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种多震源地震资料处理中的速度分析方法,解决多震源地震资料中干涉炮噪音很强,资料的信噪比极低,速度谱中能量团发散的问题。
本发明是通过以下技术方案实现的:
一种多震源地震资料处理中的速度分析方法,包括:
在工区内按间隔选取目标线,输入最小速度、最大速度和速度增量;
分别对目标线进行常速扫描,得到系列对应叠加剖面;
根据所述叠加剖面分析不同速度在目标线的不同空间位置和时间段的成像情况,得到工区的空间速度分布规律;
根据所述速度分布规律,结合能量速度谱进行速度解释,得到多震源资料的初始速度;
在初始速度的基础上,进行噪音分离与速度迭代。
所述按间隔选取目标线是按等间隔选取的。
所述分别对目标线进行常速扫描,得到系列对应叠加剖面是这样实现的:
在所述最小速度和最大速度之间,按照所述速度增量,对目标线进行扫描,得到系列对应叠加剖面;通过对目标线进行常速扫描,实质是对目标线用系列常速进行叠加,得到系列对应叠加剖面。
与现有技术相比,本发明的有益效果是:在对9个同步源的3D模拟ISS资料进行速度分析的过程中,存在着明显的无能量团的现象,导致速度谱分析无法进行,对其进行精细速度扫描,可以白浅至深发现一些列的叠加同相轴,根据这些同相轴的出现时间和扫描速度,可以对速度谱资料进行解释,从而建立3D叠加速度模型。实际资料的处理也证明了这种方法对速度谱解释的有效性。附图说明
图仁1是3D模拟巧S资料的速度谱图,此图中能量团非常发散。(正常速度谱都是彩色的,根据色彩强度和区域来确定能量团的集中与否。横坐标是速度值,纵坐标是时间,每个速度谱代表不同区域位置的。)
图仁2是3D模拟巧S资料的速度谱图,此图中0.5s以下则基本没有能量团。(横坐标是速度值,纵坐标是时间。)
图2-1是用扫描速度为1650m/s进行速度扫描得到的剖面。
图2-2是用扫描速度为1800m/s进行速度扫描得到的剖面。
图2-3是用扫描速度为2100m/s进行速度扫描得到的剖面。
图2-4是用扫描速度为2350m/s进行速度扫描得到的剖面。
图2-5是用扫描速度为2550m/s进行速度扫描得到的剖面。
图2-6是用扫描速度为2750m/s进行速度扫描得到的剖面。
图2-7是用扫描速度为2900m/s进行速度扫描得到的剖面。
图2-8是用扫描速度为3100m/s进行速度扫描得到的剖面。
图2-9是用扫描速度为3500m/s进行速度扫描得到的剖面。
图3是图1-1速度谱的解释结果。
图4是本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
多震源地震资料中,由于强烈的干涉炮噪音(cross talk)的存在,使得正常生成的速度谱能量不集中,无法分析速度,即在做速度谱分析时,往往会出现速度谱图中基本没有能量团的现象,使得速度谱分析无法进行。为此将常规地震资料处理中针对低信噪比资料的速度扫描方法应用于多震源地震资料的处理,根据速度谱分析的空间密度,在相应的测线位置进行精细速度扫描。
多震源资料由于干扰因素强大,得不到好的聚焦能量的速度谱,并对其进行解释(在能量谱上按一定的空间与时间关系,点击能量团,生成速度叫解释速度。是资料处理的关键步筹之一),以了解整个工区主要地质层位的叠加速度的空间分布规律,在此基础上再对速度谱资料进行分析,从而得到整个工区的叠加速度资料。
如图4所示,本发明的具体实现方法是:
1,在工区内间隔选好目标线(即在工区内每隔50或者20条就有一条目标线),从常速1600到5500(根据工区速度范围定)按速度增量为50,分别对目标线进行叠加。最小速度对应图4中的Vmin,最大速度对应图4中的Vmax,速度增量对应图4中的i50。常速扫描一次,得到按线号排列的所有目标线的系列常速叠加剖面。
2,每个常速叠加后都会生成对应的叠加剖面,分析这些不同速度在目标线的不同空间位置和时间段的成像情况,可以得到目标线的速度精细分布情况。具体如下:如果叠加的常速正好和叠加剖面的某个位置的地下速度吻合,则叠加剖面在该位置就能很好的成像,由此判断地下速度;从叠加剖面的成像在不同扫描速度的分布就可得到叠加剖面的空间速度分布规律。
3,根据目标扫描得到的速度分布,结合能量速度谱进行速度解释,得到好的初始速度(正因为多震元资料的速度谱能量团不集中,难以按空间与时间关系,点击能量团,但现扫描得到的大概初始速度可以投射到对应的空间位的置上.速度分析时,鼠标在分析区域移动,都会提示仔该区的时间及对应的速度,非常方便投射已知的速度控制点),在初始速度的基础上,再进行噪音分离与速度迭代(噪音分离与速度迭代是必需的后续处理,没有较为准确的速度就难以进行噪音分离,噪音部分离速度迭代就不可能越来越准确越来越收敛.)。但通过常速扫描分析后,在此指导下可以进行速度分析。
该方法得以实现,是因为多震源资料的强干扰干扰了速度能量团的生成,
导致速度解释无从下手,难以获得准确的初始速度进行后续工作。但常速扫描通过若干常速叠加,利用叠加本身的去噪功能,目标区域和扫描速度契合时,
剖面就得以成像。通过扫描,可以得到目标区域的空间和时间上的速度分布规律,返回到速度能量谱上去分析速度谱,就可以得到比较准确的速度,在这个速度基础上,再做后续的处理迭代去噪与迭代速度分析,就可以有比较好的效果了。
图1-1和图1-2是3D模拟巧IS资料的两个速度谱图(不同位置的,两者为了说明同一个问题),由此可以看出左图中能量团(指图中的深色部分)非常发散,而由图中0.5s以下则基本没有能量团,很显然这种速度谱无法解释。
图2-1至图2-9是进行速度扫描时几个主要扫描速度(速度分别是:1650m/s、1800m/s、2100m/s、2350m/s、2550m/s、2750m/s、2900m/s、3100m/s、3500m/s)所得到的剖面,由此图中可以看到,白浅至深,出现了一些比较明显的叠加同相轴,很显然我们就能得到工区主要层位的叠加速度(在哪个位置,是哪个常速扫描叠加成像的,那么这个空间位置及对应的时间位置的速度就是该速度.图2-1至图2-9是常速扫描中截出来示范的几张剖面,同相轴好的就是成像好的速度),从而从整体上建立了工区的叠加速度的大概分布规律。根据扫描得到的工区大概的速度分布速度,利用本发明方法对图1-1速度谱图进行解释,得到如图3所示的结果(找到准确的速度后,资料道集动校拉平效果很好)。从而得到工区叠加速度模型(每一次得到的速度都可以根据专业流程转换为叠加速度模型)。
多震源地震技术是近几年来发展起来的地震勘探高效野外采集技术,由于同步激发震源信号之间的相互干涉,使得资料的信噪比极低,在对其进行速度分析时往往因为速度谱中出现的能量团发散或无能量团现象而无法进行。为此将速度扫描技术应用于速度分析过程,通过对工区资料进行精细速度扫描,在了解工区主要反射层位叠加速度的分布规律的基础上,对速度谱资料进行解释,将有助于克服多震源地震基地信噪比资料速度分析的难题。模拟资料与实际资料的速度分析过程与结果证明,这是一种有效的多震源地震速度分析方法。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (3)
1.一种多震源地震资料处理中的速度分析方法,其特征在于:所述方法包括:
在工区内按间隔选取目标线,输入最小速度、最大速度和速度增量;
分别对目标线进行常速扫描,得到系列对应叠加剖面;
根据所述叠加剖面分析不同速度在目标线的不同空间位置和时间段的成像情况,得到工区的空间速度分布规律;
根据所述速度分布规律,结合能量速度谱进行速度解释,得到多震源资料的初始速度;
在初始速度的基础上,进行噪音分离与速度迭代。
2.根据权利要求1所述的多震源地震资料处理中的速度分析方法,其特征在于:所述按间隔选取目标线是按等间隔选取的。
3.根据权利要求2所述的多震源地震资料处理中的速度分析方法,其特征在于:所述分别对目标线进行常速扫描,得到系列对应叠加剖面是这样实现的:
在所述最小速度和最大速度之间,按照所述速度增量,对目标线进行扫描,得到系列对应叠加剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310499331.XA CN104570105B (zh) | 2013-10-22 | 2013-10-22 | 一种多震源地震资料处理中的速度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310499331.XA CN104570105B (zh) | 2013-10-22 | 2013-10-22 | 一种多震源地震资料处理中的速度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570105A CN104570105A (zh) | 2015-04-29 |
CN104570105B true CN104570105B (zh) | 2017-06-20 |
Family
ID=53086619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310499331.XA Active CN104570105B (zh) | 2013-10-22 | 2013-10-22 | 一种多震源地震资料处理中的速度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570105B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116446B (zh) * | 2015-08-11 | 2017-06-13 | 中国石油化工股份有限公司 | 一种用反射层位约束进行速度分析的方法 |
CN107643541B (zh) * | 2016-07-21 | 2019-02-01 | 中国石油化工股份有限公司 | 基于速度模型的速度谱解释方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5097452A (en) * | 1991-05-21 | 1992-03-17 | Western Atlas International, Inc. | Analysis of migration velocity by migration of velocity spectra |
CN1797033A (zh) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | 利用均方根速度提高叠前偏移成像精度的方法 |
-
2013
- 2013-10-22 CN CN201310499331.XA patent/CN104570105B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104570105A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Tracking the rupture of the 2008 Wenchuan earthquake by using the relative back-projection method | |
CN105093292B (zh) | 一种地震成像的数据处理方法和装置 | |
CN110058305A (zh) | 一种基于卷积神经网络的das地震数据降噪方法 | |
CN106154319B (zh) | 一种成像道集的分选方法 | |
CN103649780A (zh) | 耦合的时间-距离依赖的扫频源采集设计及数据去噪 | |
Kurzon et al. | Real‐time automatic detectors of P and S waves using singular value decomposition | |
CN102645670A (zh) | 一种基于叠加响应分析的观测系统优化设计方法 | |
CN104570105B (zh) | 一种多震源地震资料处理中的速度分析方法 | |
CN103675903A (zh) | 一种随机去噪多震源地震波场分离方法 | |
CN101937102A (zh) | 三维观测系统聚焦性能分析方法 | |
CN105510975A (zh) | 提高地震数据信噪比的方法及装置 | |
CN104133244B (zh) | 一种震源横波勘探的近地表静校正处理方法 | |
CN106990434B (zh) | 椭圆展开转换波成像方法及系统 | |
CN111665563B (zh) | 基于聚焦分析的叠前偏移垂向分辨率评价方法 | |
CN103454675B (zh) | 一种基于仪器通道号的消除随机干扰的波场分离处理方法 | |
Kuponiyi et al. | Temporal Variation in Cultural Seismic Noise and Noise Correlation Functions during COVID‐19 Lockdown in Canada | |
CN102749646A (zh) | 一种瑞雷面波深度-频率分析方法 | |
Esestime et al. | 3D shallow water seismic survey planning to deliver sub-salt imaging in South Gabon | |
Park | Roadside passive MASW survey-dynamic detection of source location | |
Xianguo et al. | Application of Broadband Seismic Acquisition under the Complex Exploration Condition | |
CN104062682B (zh) | 一种基于实时能量谱的速度分析方法 | |
CN105116446B (zh) | 一种用反射层位约束进行速度分析的方法 | |
van Baaren et al. | Survey design for high-density, high-productivity, broadband vibroseis point-source and receiver acquisition—a case study | |
Xue et al. | Automatic velocity analysis with physics-constrained optimal surface picking | |
Jiaotong et al. | Time-varient slip sweep excitation time-distance optimization technology |
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 |