CN102478664B - 一种有效信号无污染的空间采样间隔确定方法 - Google Patents
一种有效信号无污染的空间采样间隔确定方法 Download PDFInfo
- Publication number
- CN102478664B CN102478664B CN 201010560506 CN201010560506A CN102478664B CN 102478664 B CN102478664 B CN 102478664B CN 201010560506 CN201010560506 CN 201010560506 CN 201010560506 A CN201010560506 A CN 201010560506A CN 102478664 B CN102478664 B CN 102478664B
- Authority
- CN
- China
- Prior art keywords
- useful signal
- frequency
- spectra
- interval
- near surface
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及石油地球物理勘探方法,是一种地震勘探观测系统中空间采样间隔(道距)确定方法。根据工区信息建立时间、深度、地球物理模型,根据最小视速度,求取道间距进行模型正演模拟获得地震记录和有效信号记录,利用有效信号的频率及速度得到有效信号在频率-波数谱域分布范围,利用包含近地表噪音与有效信号的正演模拟数据抽取不同道间距的道集,利用有效信号和干扰波的频率和速度计算出包含近地表噪音与有效信号的频率-波数谱,对不同道距的频率-波数谱分布情况判别分析,确定最佳的空间采样间隔。本发明可有效地提高室内地震资料处理噪音压制效果,同时保护有效波,提高地震资料的信噪信噪比,降低勘探费用。
Description
技术领域
本发明涉及石油地球物理勘探方法,是一种地震勘探观测系统中空间采样间隔(道距)确定方法。
背景技术
石油地球物理勘探是基于地球物理学和石油地质学理论,采用相应的地球物理仪器和装备在地球表面(包括陆地和海洋),或者在空中、井中记录地下信息,并通过相应的数据处理和解释获取地下地层的物性(弹性、电性、磁性、密度、放射性)及结构,寻找隐藏在地层中的石油及天然气的方法。
地震勘探是地球物理勘探常用的手段,在实施地震勘探时需要将激发点和接收点按照一定方式组成连续观测的排列片,完成这种组合即观测系统设计。
目前观测系统设计中利用相邻两道的距离小于接收到的最高频率有效信号波长的一半,即空间采样间隔(道距)的设计方法从地震信号采样不产生空间假频、满足横向分辨率要求、防止偏移处理时产生偏移噪音等方面进行道距大小分析。
在进行高密度地震勘探时,在采用单点或小组合接收的情况下,如满足充分采样的要求,使信号不产生空间假频,就要求采用很小的空间采样间隔。这样的空间采样间隔可能会远高于实际地质任务的需求与地震有效信号的限制范围,将会导致勘探投入较大的浪费。另外,过大的采样间隔会造成对有效信号的折叠污染。
发明目的
本发明目的在于提供一种减少勘探成本的投入,避免折叠污染的有效信号无污染的空间采样间隔确定方法。
发明内容
本发明采用以下步骤如实现:
1)根据工区内以往的二维勘探剖面标定的不同反射层的反射时间建立时间模型,根据不同反射层的平均速度及反射时间转化为深度模型,建立工区地球物理模型;
步骤1)所述的工区地球物理模型包括的信息是:速度、密度、岩石的弹性参数、吸收衰减Q因子以及近地表信息。
所述的信息中速度、密度通过测井资料测得数据提供,岩石的弹性参数利用目的层的速度、密度进行计算,吸收衰减Q因采用谱比法计算Q值。
所述的谱比法是根据地震剖面的两个时间上的子波进行频谱分析得到两个振幅谱,对每个频率求两个振幅谱的比值,这个比值对频率的关系表示在这个层段内衰减与频率的关系。
2)根据工区的最小视速度,按下式求取道间距Δx<λmin,N/2,式中λmin,N为噪声的最小波长,利用该道距进行模型正演模拟,获得地震记录;
步骤2)所述的地震记录包括近地表噪音(面波、折射波)和有效信号。
3)通过模型正演模拟获得该区有效信号记录;
步骤3)所述的模型正演模拟是在没有近地表噪音的前提下完成。
4)利用有效信号的频率f及速度v,通过公式K=f/V计算出波数k,得到有效信号在频率-波数谱域中分布范围,或根据实际地震数据,处理得到有效信号的分布范围;
5)利用包含近地表噪音与有效信号的正演模拟数据,抽取不同道间距的道集,利用有效信号和干扰波的频率f和速度v利用公式K=f/V分别计算出包含近地表噪音与有效信号的频率-波数谱;
步骤5)所述的抽取是隔1道提取1道的地震记录信息或隔2道提取1道的地震记录信息或隔3道提取1道的地震记录信息可获得步骤2中2倍、3倍、4倍道距的地震记录或依此类推。
6)对不同道距的频率-波数谱分布情况判别分析,当干扰波的频率-波数谱发生折叠后与有效信号频率-波数谱分布发生重叠,表明干扰波的折叠信号污染有效信号,反之,干扰波的折叠信号没有污染有效信号,将干扰波的折叠信号与有效信号在频率-波数谱即将发生重叠时作为空间采样间隔(道距)的最大上限,确定为最佳的空间采样间隔。
本发明根据线性干扰的折叠频率不污染有效信号的原则,求取的空间采样间隔(道距)。可有效地提高室内地震资料处理噪音压制效果,同时保护有效波,特别是高频反射信息,提高地震资料的信噪信噪比,同时最大限度地降低勘探费用。
附图说明
图1是本发明利用模型正演模拟不同道距的单炮记录;
图2利用本发明资料选择合理道距的示意图。
具体实施方式
以下结合附图具体说明本发明。
本发明采用以下骤如实现:
1)根据工区内以往的二维勘探剖面标定的不同反射层的反射时间建立时间模型,根据不同反射层的平均速度及反射时间转化为深度模型,建立工区地球物理模型;
步骤1)所述的工区地球物理模型包括的信息是:速度、密度、岩石的弹性参数、吸收衰减Q因子以及近地表信息。
所述的信息中速度、密度通过测井资料测得数据提供,岩石的弹性参数利用目的层的速度、密度进行计算,吸收衰减Q因采用谱比法计算Q值。
所述的谱比法是根据地震剖面的两个时间上的子波进行频谱分析得到两个振幅谱,对每个频率求两个振幅谱的比值,这个比值对频率的关系表示在这个层段内衰减与频率的关系。
2)根据工区的最小视速度,按下式求取道间距Δx<λmin,N/2,式中λmin,N为噪声的最小波长,利用该道距进行模型正演模拟,获得地震记录;
步骤2)所述的地震记录包括近地表噪音(面波、折射波)与有效信号。
3)通过模型正演模拟获得该区有效信号记录;
步骤3)所述的模型正演模拟是在没有近地表噪音的前提下完成。
4)利用有效信号的频率f及速度v,通过公式K=f/V计算出波数k,得到有效信号在频率-波数谱域中分布范围(如图1),或根据实际地震数据,处理得到有效信号的分布范围(如图2);
5)利用包含近地表噪音与有效信号的正演模拟数据,抽取不同道间距的道集,利用有效信号和干扰波的频率f和速度v利用公式K=f/V分别计算出包含近地表噪音与有效信号的频率-波数谱(如图1、图2中的a为1个道距频率-波数谱,b为2倍道距频率-波数谱,c为4倍道距频率-波数谱);
步骤5)所述的抽取是隔1道提取1道的地震记录信息或隔2道提取1道的地震记录信息或隔3道提取1道的地震记录信息可获得步骤2中2倍、3倍、4倍道距的地震记录或依此类推。
6)对不同道距的频率-波数谱分布情况判别分析(如图1、图2),当干扰波的频率-波数谱发生折叠后与有效信号频率-波数谱分布发生重叠,表明干扰波的折叠信号污染有效信号,反之,干扰波的折叠信号没有污染有效信号(如图1、图2中a),将干扰波的折叠信号与有效信号在频率-波数谱即将发生重叠时作为空间采样间隔(道距)的最大上限,确定为最佳的空间采样间隔(如图1、图2中b为最佳的空间采样间隔)。
步骤6)主要根据类似图1或图2中不同道距的噪音折叠频率对有效波的影响程度进行道距的选择。图1或图2中a的道距噪音假频的能量较弱,说明噪音波场的保真度较高,采用此道距将增加勘探投入。图1、图2中b的道距干扰波发生折叠,但从图中可以看出其折叠噪音对有效波影响不大,也是就说干扰波对地震有效信号基本没有产生污染,说明采用此道距勘探可保证干扰波压制时不影响效波,采用此道距较为合理。图1、图2中c的道距噪音出现明显的假频,其折叠频率与有效波明显发生严重的混叠,污染了有效波,属典型的采样不足,说明采用此道距勘探在干扰波压制时对有效波产生压制,采用此道距不合适。根据上述干扰波的频率-波数谱发生折叠后与有效信号频率-波数谱分布发生重叠影响程度就可选择道距。
Claims (7)
1.一种有效信号无污染的空间采样间隔确定方法,其特征采用以下步骤实现:
1)根据工区内以往的二维勘探剖面标定的不同反射层的反射时间建立时间模型,根据不同反射层的平均速度及反射时间转化为深度模型,建立工区地球物理模型;
2)根据工区的最小视速度,按下式求取道间距Δx<λmin,N/2,式中λmin,N为噪声的最小波长,利用道间距进行模型正演模拟,获得地震记录;
3)通过模型正演模拟获得该区有效信号记录;
4)利用有效信号的频率f及速度V,通过公式K=f/V计算出波数K,得到有效信号在频率-波数谱域中分布范围,或根据实际地震数据,处理得到有效信号的分布范围;
5)利用包含近地表噪音与有效信号的正演模拟数据,抽取不同道间距的道集,利用有效信号和干扰波的频率f和速度V利用公式K=f/V分别计算出包含近地表噪音与有效信号的频率-波数谱;
6)对不同道间距的频率-波数谱分布情况判别分析,当干扰波的频率-波数谱发生折叠后与有效信号频率-波数谱分布发生重叠,表明干扰波的折叠信号污染有效信号,反之,干扰波的折叠信号没有污染有效信号,将干扰波的折叠信号与有效信号在频率-波数谱即将发生重叠时作为空间采样间隔的最大上限,确定为最佳的空间采样间隔。
2.根据权利要求1的方法,特征是步骤1)所述的工区地球物理模型包括的信息是:速度、密度、岩石的弹性参数、吸收衰减Q因子以及近地表信息。
3.根据权利要求2的方法,特征是所述的信息中速度、密度通过测井资料测得数据提供,岩石的弹性参数利用目的层的速度、密度进行计算,吸收衰减Q因子采用谱比法计算Q值。
4.根据权利要求3的方法,特征是所述的谱比法是根据地震剖面的两个时间上的子波进行频谱分析得到两个振幅谱,对每个频率求两个振幅谱的比值,这个比值对频率的关系表示在这个层段内衰减与频率的关系。
5.根据权利要求1的方法,特征是步骤2)所述的地震记录包括近地表噪音和有效信号。
6.根据权利要求1的方法,特征是步骤3)所述的模型正演模拟是在没有近地表噪音的前提下完成。
7.根据权利要求1的方法,特征是步骤5)所述的抽取是隔1道提取1道的地震记录信息或隔2道提取1道的地震记录信息或隔3道提取1道的地震记录信息,由此获得步骤2)中2倍或3倍或4倍道间距的地震记录,或依此类推。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010560506 CN102478664B (zh) | 2010-11-23 | 2010-11-23 | 一种有效信号无污染的空间采样间隔确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010560506 CN102478664B (zh) | 2010-11-23 | 2010-11-23 | 一种有效信号无污染的空间采样间隔确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102478664A CN102478664A (zh) | 2012-05-30 |
CN102478664B true CN102478664B (zh) | 2013-09-04 |
Family
ID=46091371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010560506 Active CN102478664B (zh) | 2010-11-23 | 2010-11-23 | 一种有效信号无污染的空间采样间隔确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102478664B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375173B (zh) * | 2013-08-15 | 2016-12-07 | 中国石油天然气集团公司 | 一种满足叠前偏移的空间采样确定方法 |
CN111257931B (zh) * | 2020-02-26 | 2020-11-17 | 青岛海洋地质研究所 | 一种去除海洋地震勘探过船干扰噪音的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1282424A (zh) * | 1997-12-19 | 2001-01-31 | 施鲁博格控股有限公司 | 叠加地震痕迹的方法 |
CN1306621A (zh) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | 自适应地震噪声和干扰衰减方法 |
CN1633610A (zh) * | 2002-01-14 | 2005-06-29 | 维斯特恩格科地震控股有限公司 | 处理地震数据 |
CN101382598A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种真三维地震数据线性噪音的压制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7715986B2 (en) * | 2007-05-22 | 2010-05-11 | Chevron U.S.A. Inc. | Method for identifying and removing multiples for imaging with beams |
US20090092003A1 (en) * | 2007-10-08 | 2009-04-09 | Nicolas Goujon | Controlling a seismic survey to reduce the effects of vibration noise |
-
2010
- 2010-11-23 CN CN 201010560506 patent/CN102478664B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1282424A (zh) * | 1997-12-19 | 2001-01-31 | 施鲁博格控股有限公司 | 叠加地震痕迹的方法 |
CN1306621A (zh) * | 1998-05-20 | 2001-08-01 | 施鲁博格控股有限公司 | 自适应地震噪声和干扰衰减方法 |
CN1633610A (zh) * | 2002-01-14 | 2005-06-29 | 维斯特恩格科地震控股有限公司 | 处理地震数据 |
CN101382598A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种真三维地震数据线性噪音的压制方法 |
Non-Patent Citations (5)
Title |
---|
刘志鹏等.径向道变换压制高密度采集资料中的相干噪声.《石油地球物理勘探》.2008,(第03期), |
基于平稳过滤有色噪声地震动模型的MTMD性能评价;李春祥等;《计算力学学报》;20030430(第02期);全文 * |
径向道变换压制高密度采集资料中的相干噪声;刘志鹏等;《石油地球物理勘探》;20080615(第03期);全文 * |
李春祥等.基于平稳过滤有色噪声地震动模型的MTMD性能评价.《计算力学学报》.2003,(第02期), |
毕丽飞,孟宪军等.三维转换波资料处理方法研究及应用.《石油地球物理勘探》.2008,第43卷(第增刊2期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102478664A (zh) | 2012-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
Hartzell et al. | Variability of site response in the Los Angeles urban area | |
CN100349011C (zh) | 地震数据处理中压制与激发源无关的背景噪声的方法 | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
CN102565856B (zh) | 一种基于波动方程正演的近地表噪音压制方法 | |
CN101556338A (zh) | 一种可控震源自适应地表一致性反褶积方法 | |
CN102053263B (zh) | 调查表层结构的方法 | |
CN102645670A (zh) | 一种基于叠加响应分析的观测系统优化设计方法 | |
Shao et al. | Near-surface characterization using high-speed train seismic data recorded by a distributed acoustic sensing array | |
Dean et al. | The coherency of ambient seismic noise recorded during land surveys and the resulting implications for the effectiveness of geophone arrays | |
Deng et al. | The extended range phase shift method for broadband surface wave dispersion measurement from ambient noise and its application in ore deposit characterization | |
CN104422957A (zh) | 一种观测系统优化设计方法 | |
CN102478664B (zh) | 一种有效信号无污染的空间采样间隔确定方法 | |
CN104597497A (zh) | 一种基于叠前瞬时频率属性分析的储层烃类预测方法 | |
Sangster | Dunedin rock and roll: 3D seismic wave velocity modelling for seismic hazard analysis. | |
Lau et al. | Low-frequency source for long-offset, sub-basalt and deep crustal penetration | |
Almholt et al. | High resolution 2D reflection seismic land streamer survey for groundwater mapping: Case study from south east Denmark | |
Hu et al. | Numerical simulation scattered imaging in deep mines | |
Yang et al. | Joint reverse-time imaging condition of seismic towed-streamer and OBN data | |
Burschil et al. | Finite-difference modelling to evaluate seismic P-wave and shear-wave field data | |
Anthymidis et al. | Constraining site response and shallow geophysical structure by ambient noise measurements and 1D numerical simulations: the case of Grevena town (N. Greece) | |
Wu et al. | An Automatic Screening Method for the Passive Surface-Wave Imaging Based on the FK Domain Energy Characteristics | |
Vera | Sesmic Modelling of CO2 in a Sandstone Aquifer, Priddis Alberta | |
Bey | Reservoir Characterization and Seismic Expression of the Clinton Interval over Dominion's Gabor Gas Storage Field in North-East Ohio | |
Urosevic et al. | Acquisition of the Nested 3D seismic survey at Harvey |
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 |