CN103823242B - 一种快速准确选取可控震源滑动时间参数的方法 - Google Patents
一种快速准确选取可控震源滑动时间参数的方法 Download PDFInfo
- Publication number
- CN103823242B CN103823242B CN201210464540.6A CN201210464540A CN103823242B CN 103823242 B CN103823242 B CN 103823242B CN 201210464540 A CN201210464540 A CN 201210464540A CN 103823242 B CN103823242 B CN 103823242B
- Authority
- CN
- China
- Prior art keywords
- time
- controlled source
- sliding time
- record
- sliding
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明是地震勘探中快速准确选取可控震源滑动时间参数的方法,根据已知的可控震源扫描长度和起止扫描频率,计算出避开二次谐波所要求的最小滑动时间,将原始数据在时间域取不同的时间向后延迟,获得不同模拟滑动时间的记录,垂直叠加后与可控震源扫描信号进行互相关,根据对模拟记录谐波干扰的大小和范围,选择干扰能量弱和范围小的模拟滑动时间,将起始时间大于目的层时间的滑动时间确定为滑动时间参数。本发明通过室内模拟方法快速准确地选取滑动时间,实现了可控震源滑动扫描方法施工中滑动时间的快速、准确选择,降低了野外试验周期,确保滑动时间参数选取的合理性。
Description
技术领域
本发明涉及石油勘探工程,是地震勘探中使用可控震源滑动扫描方法生产中快速准确选取可控震源滑动时间参数的方法。
背景技术
目前,在石油地震勘探中,滑动扫描技术是国际上可控震源高效采集技术中最成熟、最常用的技术之一,该技术已获得了良好的勘探效果和经济效益。滑动扫描是指使用多组可控震源同时进行作业,在该技术应用中,谐波干扰是客观存在的,在处理时需要专门的压制技术,在野外通过设计合理的滑动时间也可以有效减弱谐波干扰。滑动时间是可控震源相邻两次震动时间间隔,滑动时间原则上大于地震记录长度即可,与交替扫描相比,相邻的两次振动(扫描),突破了第二次扫描必须等第一次扫描记录结束才能开始的限制,压缩了相邻两次扫描的间隔时间,可提高生产效率。滑动时间参数会影响原始资料品质,目前可控震源滑动扫描施工中,确定滑动时间往往是野外应用2套或多套震源,按照滑动扫描方式分别按照不同的滑动时间进行数据采集,室内对上述数据进行分析对比,以确定合适的滑动时间,一般是通过野外大量的点上试验来确定滑动时间参数,工作量大,资料分析周期长,如果施工区域复杂,生产过程中需要考核试验时,将会进一步影响生产效率,同时,确定的滑动时间参数也具有一定的偶然性。
发明内容
本发明目的在于提供一种可提高生产效率,提高资料品质的快速准确选取可控震源滑动时间参数的方法。
本发明通过以下步骤实现:
1)根据已知的可控震源扫描长度和起止扫描频率,计算出避开二次谐波所要求的最小滑动时间Ti;
步骤1)所述的最小滑动时间Ti采用以下公式计算:
其中:T是扫描长度,fB、fE分别是起止扫描频率;k为常数;
所述的k为常数,k=2。
2)野外现场采集可控震源原始未相关数据;
步骤2)所述的野外现场采集可控震源原始未相关数据是利用可控震源扫描信号S(t),在震源1和震源2位置,利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)。
3)将原始数据Raw2(t)在时间域向后延迟时间i,得到延迟后的记录Raw2(ti),每次取不同的时间i,获得不同模拟滑动时间i的记录Raw2(ti);
所述的时间i大于步骤1)中估算的最小滑动时间Ti,小于或等于扫描长度T,i取整数,根据实际情况i取3-5个值。
4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加,得到对应模拟滑动时间i的数据Raw(ti);
步骤4)所述的垂直叠加采用以下公式:
Raw(ti)=Raw1(t)+Raw2(ti) (2);
5)将Raw(ti)与可控震源扫描信号S(t)进行互相关,得到模拟滑动时间i的地震记录R(ti);
6)根据对模拟记录R(ti)谐波干扰的大小和范围,选择干扰能量弱和范围小的模拟滑动时间,同时,在记录上读取谐波干扰的起始时间,选择干扰起始时间大于目的层时间的滑动时间,确定可控震源滑动时间参数i。
本发明实现了可控震源滑动扫描方法施工中滑动时间的快速、准确选择,减少了野外试验工作量,降低了野外试验费用,缩短了试验周期,确保滑动时间参数选取的合理性。
附图说明
图1是采集可控震源原始未相关记录的观测系统示意图;
图2是用于模拟分析的不同位置的相关记录;
图3是模拟不同滑动时间的相关记录;
图4是模拟不同滑动时间的记录时频图,
A)是模拟滑动时间8s的记录时频图;
B)是模拟滑动时间10s的记录时频图;
C)是模拟滑动时间12s的记录时频图;
D)是模拟滑动时间14s的记录时频图。
具体实施方式
以下结合附图详细说明本发明。
本发明采用以下步骤实现:
1)野外通过一条适当长的二维排列接收平面上不同位置的可控震源常规未相关记录,根据已知的可控震源扫描长度T和起止扫描频率fB、fE,由公式(1)计算出避开二次谐波所要求的最小滑动时间Ti;
式中的k为常数,k=2。
2)野外现场采集可控震源原始未相关数据
步骤2)所述的野外现场采集可控震源原始未相关数据是根据已确定的可控震源参数,利用可控震源扫描信号S(t),在震源1和震源2位置,利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t)。见图1,其中“△”代表接收排列,“●”代表可控震源位置。
3)将原始数据Raw2(t)在时间域向后延迟时间i,该步骤可以在CGG、ProMax、Ω、GeoEast等处理系统上,利用野外静校正应用模块实现。如在GeoEast处理系统上,利用静校正应用模块StApply,在“application potions”选项里面选择“constant correction”,然后在“correction constant”参数里填写常数校正时间-i,单位为ms,生成SEG-Y格式数据,输出文件名Raw2(ti)。
对不同的时间i,重复该步骤,得到不同模拟滑动时间i的记录Raw2(ti)。所述的时间i大于步骤1)中估算的最小滑动时间Ti,小于扫描长度T,i取整数,根据实际情况i取3-5个值。
4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加,得到对应模拟滑动时间i的数据Raw(ti);
步骤4)所述的垂直叠加采用以下公式:
Raw(ti)=Raw1(t)+Raw2(ti) (2);
可在CGG、ProMax、Ω、GeoEast等处理系统上,利用多功能叠加模块将Raw1(t)和Raw2(ti)按同道号进行叠加,形成一个新的炮集。如在GeoEast处理系统上,利用叠加模块VersatiStk,在“input gather type”中选择“source”作为输入道集类型,在“header wordfor stack”参数中选择“trace number”作为叠加道头字。将Raw1(t)和Raw2(ti)按同道号进行叠加,形成一个新的炮集,生成SEG-Y格式数据,输出文件名Raw(ti)。
5)将Raw(ti)与可控震源扫描信号S(t)进行互相关,得到模拟滑动时间i的地震记录R(ti);
可在CGG、ProMax、Ω、GeoEast等处理系统上,利用互相关分析模块,将数据Raw(ti)与扫描信号S(t)互相关,生成互相关后的地震记录R(ti)。如在GeoEast处理系统上,在输入数据Raw(ti)上点击右键选择“Correlation Analysis”启动GeoSeismic Correlation交互相关分析模块,此时数据Raw(ti)已经在交互界面的数据窗口下打开,然后需要在信号窗口里面打开扫描信号S(t),选择交互界面上端的“Correlation”选项里面的“Crosscorrelation”互相关选项,弹出互相关参数定义窗口,选择填写需要计算的地震数据时窗和信号数据时窗等参数,点击“OK”后生成互相关后的地震记录R(ti),以SEG-Y格式数据输出。
分别对每个i值重复步骤3)至步骤5),得到不同模拟滑动时间i的地震记录R(ti);
6)分析不同的模拟记录R(ti)即可准确选取可控震源滑动时间参数i;
对模拟记录R(ti)进行能量和信噪比定量分析统计,确定记录上谐波干扰的大小和范围,选择谐波干扰能量弱和范围小的模拟记录,该模拟记录所对应的模拟滑动时间就是要选取的可控震源滑动时间参数i。
本发明实验实例如下:
可控震源参数是:扫描长度T=14s,fB=6Hz,fE=84Hz,根据公式(1)计算出避开二次谐波所要求的最小滑动时间Ti=7.5s。在该实例中,模拟滑动时间i取8s、10s、12s、14s。
从图2所示的可控震源常规相关记录看,谐波干扰较弱,只有近道震源机械干扰。从图3所示的模拟不同滑动时间的相关记录看,干扰主要来自相邻炮的二阶谐波,随着模拟滑动时间的增加,谐波干扰逐渐减弱。从图4时频图分析,除二阶谐波能量稍强外,其它高阶谐波干扰很弱,对资料品质影响较小。在模拟滑动时间8s的记录和图4A)所示的时频图上,本炮受到后续相邻炮的谐波干扰,干扰最强,而在模拟滑动时间10s、12s、14s的记录上,谐波干扰相对较弱,只是干扰到近道排列,从图4B、图4C、图4D时频图上看,基波已经避开二次谐波干扰。因此,通过分析不同滑动时间模拟记录,当扫描长度为14s时,滑动时间10s的模拟记录上谐波干扰较弱。因此,本实例中滑动时间选择为10s。
Claims (4)
1.一种快速准确选取可控震源滑动时间参数的方法,特点是通过以下步骤实现:
1)根据已知的可控震源扫描长度和起止扫描频率,计算出避开二次谐波所要求的最小滑动时间Ti;
2)野外现场采集可控震源原始未相关数据;
野外现场采集可控震源原始未相关数据是利用可控震源扫描信号S(t),在震源1和震源2位置,利用相同排列记录可控震源原始未相关数据Raw1(t)和Raw2(t);
3)将原始数据Raw2(t)在时间域向后延迟时间i,得到延迟后的记录Raw2(ti),每次取不同的时间i,获得不同模拟滑动时间i的记录Raw2(ti);
4)将由步骤2)得到的Raw1(t)和由步骤3)得到的不同模拟滑动时间i的记录Raw2(ti)分别进行垂直叠加,得到对应模拟滑动时间i的数据Raw(ti);
垂直叠加采用以下公式:
Raw(ti)=Raw1(t)+Raw2(ti) (2);
5)将Raw(ti)与可控震源扫描信号S(t)进行互相关,得到模拟滑动时间i的地震记录R(ti);
6)根据对模拟记录R(ti)谐波干扰的大小和范围,选择干扰能量弱和范围小的模拟滑动时间,同时,在记录上读取谐波干扰的起始时间,选择干扰起始时间大于目的层时间的滑动时间,确定可控震源滑动时间参数i。
2.根据权利要求1所述的方法,特点是步骤1)所述的最小滑动时间Ti采用以下公式计算:
其中:T是扫描长度,fB、fE分别是起止扫描频率;k为常数。
3.根据权利要求2所述的方法,特点是,k=2。
4.根据权利要求1所述的方法,特点是步骤3)所述的时间i大于步骤1)中估算的最小滑动时间Ti,小于扫描长度T,i取3-5个整数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210464540.6A CN103823242B (zh) | 2012-11-16 | 2012-11-16 | 一种快速准确选取可控震源滑动时间参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210464540.6A CN103823242B (zh) | 2012-11-16 | 2012-11-16 | 一种快速准确选取可控震源滑动时间参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823242A CN103823242A (zh) | 2014-05-28 |
CN103823242B true CN103823242B (zh) | 2016-12-21 |
Family
ID=50758390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210464540.6A Active CN103823242B (zh) | 2012-11-16 | 2012-11-16 | 一种快速准确选取可控震源滑动时间参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823242B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104111477B (zh) * | 2014-07-04 | 2017-01-11 | 中国石油天然气集团公司 | 一种用于压制谐波的滑动时间获取方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101627322A (zh) * | 2007-03-09 | 2010-01-13 | 格库技术有限公司 | 可控震源采集方法 |
CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050128874A1 (en) * | 2003-12-15 | 2005-06-16 | Chevron U.S.A. Inc. | Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources |
US20070195644A1 (en) * | 2006-02-21 | 2007-08-23 | Timothy Marples | Methods and Systems for Efficient Compaction Sweep |
-
2012
- 2012-11-16 CN CN201210464540.6A patent/CN103823242B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101627322A (zh) * | 2007-03-09 | 2010-01-13 | 格库技术有限公司 | 可控震源采集方法 |
CN102692649A (zh) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | 一种可控震源地震勘探中消除谐波干扰的方法 |
Non-Patent Citations (2)
Title |
---|
可控震源地震采集技术的进展;倪宇东等;《石油地球物理勘探》;20110630;第46卷(第3期);第249-356页 * |
滑动扫描谐波分析;曹务祥等;《石油地球物理勘探》;20051031;第40卷(第5期);第499-509页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103823242A (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Langston | Spatial gradient analysis for linear seismic arrays | |
CN103582825A (zh) | 同步源地震数据的变频滤波 | |
CN104950326B (zh) | 基于目的层频谱的可控震源非线性扫描信号的设计方法 | |
O’Connell et al. | Interferometric multichannel analysis of surface waves (IMASW) | |
CN106094024A (zh) | 一种可控震源滑动扫描数据采集处理方法及装置 | |
CN105572723B (zh) | 基于自相关子波的可控震源扫描信号的设计方法 | |
CN101923177A (zh) | 多波地震资料的速度分析方法和装置 | |
CN102692649A (zh) | 一种可控震源地震勘探中消除谐波干扰的方法 | |
Rodgers et al. | Broadband waveform modeling of moderate earthquakes in the San Francisco Bay area and preliminary assessment of the USGS 3D seismic velocity model | |
Hoshiba | Real‐time correction of frequency‐dependent site amplification factors for application to earthquake early warning | |
CN109507726A (zh) | 时间域弹性波多参数全波形的反演方法及系统 | |
CN105005075B (zh) | 基于地震频率信息的多波匹配方法 | |
CN104035128A (zh) | 可控震源伪随机扫描信号生成方法 | |
CN106324666A (zh) | 横向各向同性介质转换波角道集抽取方法和装置 | |
CN101839998B (zh) | 一种高精度的叠前深度偏移方法 | |
CN103823242B (zh) | 一种快速准确选取可控震源滑动时间参数的方法 | |
CN108957551B (zh) | 基于重构地面力信号的可控震源谐波压制方法 | |
CN103941289B (zh) | 一种绕射点成像的方法和设备 | |
Yamada et al. | Real-time estimation of fault rupture extent using envelopes of acceleration | |
CN102216810B (zh) | 利用基于射线的原理和有限元原理两者来推导地震波场的系统和方法 | |
Liu et al. | Wave-equation diffraction imaging using pseudo dip-angle gather | |
CN104111477B (zh) | 一种用于压制谐波的滑动时间获取方法及系统 | |
CN105445788B (zh) | 一种基于模型和全局寻优的速度谱自动解释方法 | |
CN104217099B (zh) | 一种用于产生并处理可控震源的扫描信号的方法 | |
CN107607995B (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 |