CN109581500B - 一种反射地震记录频变速度分析方法 - Google Patents
一种反射地震记录频变速度分析方法 Download PDFInfo
- Publication number
- CN109581500B CN109581500B CN201811546192.0A CN201811546192A CN109581500B CN 109581500 B CN109581500 B CN 109581500B CN 201811546192 A CN201811546192 A CN 201811546192A CN 109581500 B CN109581500 B CN 109581500B
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- seismic
- velocity
- space
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 claims abstract description 40
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 230000001131 transforming effect Effects 0.000 claims abstract description 6
- 230000001419 dependent effect Effects 0.000 claims description 10
- 238000000034 method Methods 0.000 claims description 7
- 238000003306 harvesting Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
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
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
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
本发明涉及一种反射地震记录频变速度分析方法，包括如下步骤：a、利用广义S变换将地震记录由空间‑时间域变换到空间‑时间‑频率域，在频率切片求取时频谱极大值对应时间，提取不同频率地震波的反射时距曲线；b、进行速度扫描，选取计算的理论时距曲线和地震记录时频谱中提取的时距曲线拟合精度最高的速度，即认为是某频率地震波的传播速度。本发明提出了新的反射地震记录频变速度分析方法，克服了传统速度分析方法无法获得速度频散属性的缺点，频变速度分析技术可以利用速度频散属性进行油气储层预测，基于频变速度分析技术可以提高地震处理速度分析、动校正和叠加的精度。
Description
技术领域
本发明涉及地震记录的速度分析方法，尤其是一种反射地震记录频变速度分析方法。
背景技术
在油气储层、矿产资源、工程地质及环境地质反射地震勘探中，地震速度分析技术是反射地震数据处理的重要和必须环节。速度分析质量的优劣直接影响反射地震成像精度。同时，地震波传播速度作为描述地质体的重要地震属性，也被用来对岩性，尤其时储层性质进行定性和定量解释。
储层和海底介质等是典型的孔隙多相介质，地震波在储层和海底介质中的传播速度从在频散已经得到实践证明(基于储层砂岩微观孔隙结构特征的弹性波频散响应分析[J].地球物理学报,2015,58(9):3389-3400.)，地震波速度频散成为储层预测越来越被重视的属性。但是，目前的反射地震速度分析技术研究尚未有报道频变速度分析技术，学者们都是在不同复杂程度的地质模型基础上或不同噪声条件下探讨将所有频率地震波传播速度认为是同一速度的速度分析方法。不依赖频率的速度分析方法获得的速度参数不随频率变化，会导致反射地震资料处理中的动校正和叠加不能取得更好的效果，会导致无法利用速度分析技术获得速度参数的频散属性，从而预测油气储层。
发明内容
本发明的目的是为了解决现有技术存在的缺陷，提供一种反射地震记录频变速度分析方法。
为了实现上述目的，本发明采用的技术方案是：
一种反射地震记录频变速度分析方法，包括如下步骤：
a、利用广义S变换将地震记录由空间-时间域变换到空间-时间-频率域，在频率切片求取时频谱极大值对应时间，提取不同频率地震波的反射时距曲线；
b、进行速度扫描，选取计算的理论时距曲线和地震记录时频谱中提取的时距曲线拟合精度最高的速度，即认为是某频率地震波的传播速度。
进一步，采用广义S变换时频分析技术，反射地震记录u(x,t)的广义S变换为：
这里u(x,t)为地震记录，x表示地震记录道空间，t表示地震记录时间，U(x,τ,f)为空间-时间-频率域地震记录，τ为地震记录时频谱时间，f为频率，λ、p为广义S变换调节参数，exp表示e指数。
进一步，通过空间-时间-频率谱的分析，提取能谱极大值位置，获得不同频率的反射波时距关系，具体步骤如下：
步骤1：在时间域中输入空间-时间域二维地震数据，然后采用广义S变换，将重数据从空间-时间域变换到空间-时间-频率域；
步骤2：基于空间-时间-频率谱的频率切片，提取能谱极大值对应的不同地震道的该频率的地震波履行时间；
[xi,TpF]＝find(U(xi,t,F)＝max(U(xi,t,F)))
[xi,TpF]表示频率为F的地震波在第i地震道的反射到达时为TpF，也即频率为FHz的地震波的反射时距关系，xi为第i道记录的炮检距，U(xi,t,F)为地震记录空间-时间-频率谱在FHz频率处的切片。
进一步，使用速度扫描，计算理论时距曲线与提取的时距曲线差的2范数，范数的极小值对应的速度即为某频率地震波传播速度；
理论时距曲线表示为：
式中：TtFj(xi,F)表示频率为F的地震波在炮检距为xi的第i道检波器的旅行时，t0Fj表示频率为F的地震波在第j个反射层的自激自收旅行时，V0Fj为频率为F的地震波在第j个反射层处的传播速度；
理论时距曲线与提取的时距曲线差的2范数表示为：||TtFj-TpFj||2，范数的极小值对应的速度即为某频率地震波传播速度：VjF＝find(||TtFj-TpFj||2＝min(||TtFj-TpFj||2))
上式表示通过速度扫描，找出使得||TtFj-TpFj||2取最小值的速度，即为某频率地震波传播速度。
本发明的有益效果为：本发明提出了新的反射地震记录频变速度分析方法，克服了传统速度分析方法无法获得速度频散属性的缺点，频变速度分析技术可以利用速度频散属性进行油气储层预测，基于频变速度分析技术可以提高地震处理速度分析、动校正和叠加的精度。
附图说明
图1为本发明实施例中反射地震记录频变速度分析方法流程图；
图2为地震记录图；
图3为地震记录第一道时频谱图；
图4为子波频谱图；
图5为空间-时间-频率域频率切片能谱极值提取50Hz地震波时距曲线图；
图6a为速度扫描提取的地震波为第一层频变速度图；
图6b为速度扫描提取的地震波为第二层频变速度图；
图7为50Hz地震波理论时距曲线与时频谱提取时距曲线拟合图。
具体实施方式
如图1所示，一种反射地震记录频变速度分析方法，实现该方法的步骤主要包括，地震记录的时频变换，利用时频谱的频率切片谱极值提取频变反射地震波时距曲线，基于速度扫描技术获得理论计算时距曲线，通过理论时距曲线与时频谱提取的频变时距曲线差的2范数极小值获得频变速度等。具体步骤如下：
步骤1：地震记录的时频变换。反射地震记录u(x,t)的广义S变换为：
这里u(x,t)为地震记录，x表示地震记录道空间，t表示地震记录时间，U(x,τ,f)为空间-时间-频率域地震记录，τ为地震记录时频谱时间，f为频率，λ、p为广义S变换调节参数，exp表示e指数。
步骤2：利用时频谱的频率切片谱极值提取频变反射地震波时距曲线。通过空间-时间-频率谱的分析，提取能谱极大值位置，获得不同频率的反射波时距关系：(1)在时间域中输入空间-时间域二维地震数据，然后采用广义S变换，将重数据从空间-时间域变换到空间-时间-频率域；(2)基于空间-时间-频率谱的频率切片，提取能谱极大值对应的不同地震道的该频率的地震波履行时间；
[xi,TpF]＝find(U(xi,t,F)＝max(U(xi,t,F)))
[xi,TpF]表示频率为F的地震波在第i地震道的反射到达时为TpF，也即频率为FHz的地震波的反射时距关系，xi为第i道记录的炮检距，U(xi,t,F)为地震记录空间-时间-频率谱在FHz频率处的切片。
步骤3：基于速度扫描技术获得理论计算时距曲线。使用速度扫描，计算理论时距曲线与提取的时距曲线差的2范数，范数的极小值对应的速度即为某频率地震波传播速度。理论时距曲线可表示为：
式中：TtFj(xi,F)表示频率为F的地震波在炮检距为xi的第i道检波器的旅行时，t0Fj表示频率为F的地震波在第j个反射层的自激自收旅行时，V0Fj为频率为F的地震波在第j个反射层处的传播速度。
步骤4：通过理论时距曲线与时频谱提取的频变时距曲线差的2范数极小值获得频变速度。理论时距曲线与提取的时距曲线差的2范数表示为：||TtFj-TpFj||2，范数的极小值对应的速度即为某频率地震波传播速度：
VjF＝find(||TtFj-TpFj||2＝min(||TtFj-TpFj||2))
上式表示通过速度扫描，找出使得||TtFj-TpFj||2取最小值的速度，即为某频率地震波传播速度。
实现该方法具体操作为：
(1)地震记录广义S变换参数分析与时频变换
为了获得更好地地震记录的广义S变化效果，对具体地震记录(图2所示)广义S变换的调节参数进行试验分析，选取分析效果最佳的参数。本次数据为24道地震记录，包含两个反射界面，第一层地震波传播速度为常数，第二层地震波传播速度为频变速度，数据采样率为1000Hz，采样长度为1s，道间距为10米，偏移距为0米。针对本次数据，λ、p分别取值为3和0.9。从地震记录第一道的时频变换的效果可以看出，时频谱(图3所示)在时间和频率域都有很好的分辨率，地震波频率分布与子波频率(图4所示)分布一致。
(2)根据频率切片谱极值提取每个频率反射地震波时距曲线
通过空间-时间-频率谱的分析，提取频率切片的能谱极大值位置，获得不同频率的反射波时距关系(图5所示)。
(3)速度扫描。分析理论时距曲线与时频谱提取的频变时距曲线差的2范数极小值获得频变速度(图6a为第一层频变速度，图6b为第二层频变速度)。按照获得的正确频变速度会使得理论时距曲线和时频谱提取的时距曲线拟合较好(图7所示，两个反射层的50Hz频率地震波时频谱提取的时距曲线和根据提取的频变速度计算的理论时距曲线拟合效果)。
本发明提出了新的反射地震记录频变速度分析方法，克服了传统速度分析方法无法获得速度频散属性的缺点，频变速度分析技术可以利用速度频散属性进行油气储层预测，基于频变速度分析技术可以提高地震处理速度分析、动校正和叠加的精度。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是本发明的原理，在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明的范围内。本发明要求的保护范围由所附的权利要求书及其等同物界定。
Claims (3)
1.一种反射地震记录频变速度分析方法，其特征在于，包括如下步骤：
a、利用广义S变换将地震记录由空间-时间域变换到空间-时间-频率域，在频率切片求取时频谱极大值对应时间，提取不同频率地震波的反射时距曲线；
b、进行速度扫描，选取计算的理论时距曲线和地震记录时频谱中提取的时距曲线拟合精度最高的速度，即使用速度扫描，计算理论时距曲线与提取的时距曲线差的2范数，范数的极小值对应的速度即为某频率地震波传播速度；
理论时距曲线表示为：
式中：TtFj(xi，F)表示频率为F的地震波在炮检距为xi的第i道检波器的旅行时，t0Fj表示频率为F的地震波在第j个反射层的自激自收旅行时，V0Fj为频率为F的地震波在第j个反射层处的传播速度；
理论时距曲线与提取的时距曲线差的2范数表示为：||TtFj-TpFj||2，范数的极小值对应的速度即为某频率地震波传播速度：
VjF＝find(||TtFj-TpFj||2＝min(||TtFj-TpFj||2))
上式表示通过速度扫描，找出使得||TtFj-TpFj||2取最小值的速度，即为某频率地震波传播速度。
2.根据权利要求1所述的一种反射地震记录频变速度分析方法，其特征在于，采用广义S变换时频分析技术，反射地震记录u(x，t)的广义S变换为：
这里u(x，t)为地震记录，x表示地震记录道空间，t表示地震记录时间，U(x，τ，f)为空间-时间-频率域地震记录，τ为地震记录时频谱时间，f为频率，λ、p 为广义S变换调节参数，exp表示e指数。
3.根据权利要求1或2所述的一种反射地震记录频变速度分析方法，其特征在于，通过空间-时间-频率谱的分析，提取能谱极大值位置，获得不同频率的反射波时距关系，具体步骤如下：
步骤1：在时间域中输入空间-时间域二维地震数据，然后采用广义S变换，将重数据从空间-时间域变换到空间-时间-频率域；
步骤2：基于空间-时间-频率谱的频率切片，提取能谱极大值对应的不同地震道的该频率的地震波履行时间；
[xi，TpF]＝find(U(xi，t，F)＝max(U(xi，t，F)))
[xi，TpF]表示频率为F的地震波在第i地震道的反射到达时为TpF，也即频率为F Hz的地震波的反射时距关系，xi为第i道记录的炮检距，U(xi，t，F)为地震记录空间-时间-频率谱在F Hz频率处的切片。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811546192.0A CN109581500B (zh) | 2018-12-18 | 2018-12-18 | 一种反射地震记录频变速度分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811546192.0A CN109581500B (zh) | 2018-12-18 | 2018-12-18 | 一种反射地震记录频变速度分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109581500A CN109581500A (zh) | 2019-04-05 |
| CN109581500B true CN109581500B (zh) | 2020-06-30 |
Family
ID=65930610
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811546192.0A Expired - Fee Related CN109581500B (zh) | 2018-12-18 | 2018-12-18 | 一种反射地震记录频变速度分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109581500B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113156504B (zh) * | 2021-04-09 | 2021-10-15 | 中国科学院地理科学与资源研究所 | 一种地震波速度的确定方法及装置 |
| CN113189643B (zh) * | 2021-04-30 | 2023-08-25 | 西安石油大学 | 一种水工闸门地震动和水流激励的耦合系统及合成方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018499A (en) * | 1997-11-04 | 2000-01-25 | 3Dgeo Development, Inc. | Three-dimensional seismic imaging of complex velocity structures |
| WO2002006856A1 (en) * | 2000-07-17 | 2002-01-24 | Sintef Petroleumsforskning As | Seismic processing with general non-hyperbolic travel-time corre ctions |
| WO2015042386A1 (en) * | 2013-09-20 | 2015-03-26 | Westerngeco Llc | Eikonal solver for quasi p-waves in anisotropic media |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738637B (zh) * | 2008-11-06 | 2012-01-04 | 北京北方林泰石油科技有限公司 | 一种基于速度随频率变化信息的油气检测方法 |
| CN105445795B (zh) * | 2014-08-11 | 2018-05-08 | 中国石油化工股份有限公司 | 倾斜界面的速度分析方法及其系统 |
| CN105785440B (zh) * | 2016-02-29 | 2017-09-15 | 河南理工大学 | 一种矿井槽波双分量地震信号频散曲线提取方法 |
| CN107346034A (zh) * | 2016-05-04 | 2017-11-14 | 中国石油化工股份有限公司 | 基于广义s变换的谱相关系数的q值估计方法 |
| CN106199711A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | Co2驱油地震频散属性提取方法 |
| CN106019374B (zh) * | 2016-07-04 | 2018-06-26 | 中煤科工集团西安研究院有限公司 | 基于反射槽波频散相似度的断层成像方法 |
| CN106291689A (zh) * | 2016-08-18 | 2017-01-04 | 中国石油天然气集团公司 | 一种提取地震数据频散属性的处理方法、装置及预测系统 |
-
2018
- 2018-12-18 CN CN201811546192.0A patent/CN109581500B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018499A (en) * | 1997-11-04 | 2000-01-25 | 3Dgeo Development, Inc. | Three-dimensional seismic imaging of complex velocity structures |
| WO2002006856A1 (en) * | 2000-07-17 | 2002-01-24 | Sintef Petroleumsforskning As | Seismic processing with general non-hyperbolic travel-time corre ctions |
| WO2015042386A1 (en) * | 2013-09-20 | 2015-03-26 | Westerngeco Llc | Eikonal solver for quasi p-waves in anisotropic media |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109581500A (zh) | 2019-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2707526C (en) | Spectral shaping inversion and migration of seismic data | |
| Srivastava et al. | Stochastic inversion of prestack seismic data using fractal-based initial models | |
| US20160178772A1 (en) | High Resolution Estimation of Attenuation from Vertical Seismic Profiles | |
| CN104570067B (zh) | 一种地球物理勘探中相控地震反演方法 | |
| CN106443768B (zh) | 一种叠前深度域合成地震记录制作方法 | |
| Diallo et al. | Characterization of polarization attributes of seismic waves using continuous wavelet transforms | |
| CN103163554A (zh) | 利用零偏vsp资料估计速度和q值的自适应波形的反演方法 | |
| EP3622327B1 (en) | Processing methodology for full-waveform sonic wavefield separation | |
| CN108020863A (zh) | 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法 | |
| Liner et al. | SPICE: A new general seismic attribute | |
| Liu et al. | A review of variational mode decomposition in seismic data analysis | |
| CN109581500B (zh) | 一种反射地震记录频变速度分析方法 | |
| CN110967749A (zh) | 一种vsp地震资料频变q值估计与反q滤波方法 | |
| US20130144534A1 (en) | Multi-Component Spectral Warping Analysis for Thin Bed Resolution | |
| Xue et al. | Q-factor estimation by compensation of amplitude spectra in synchrosqueezed wavelet domain | |
| CN109283581A (zh) | 一种基于深度域地震波频散分析的储层含气性评价方法 | |
| Zhang et al. | Interval Q inversion based on zero-offset VSP data and applications | |
| CN113031070B (zh) | 一种深度域合成地震记录的制作方法 | |
| CN110673211B (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
| Feng et al. | Dispersion calculation method based on S-transform and coordinate rotation for Love channel waves with two components | |
| CN112180464A (zh) | 一种储层物性的识别方法 | |
| Oren et al. | Image-domain DAS 3D VSP elastic transmission tomography | |
| Yan et al. | Combining the synchrosqueezing generalized S-transform of variational mode decomposition with the Teager–Kaiser energy operator to calculate the attenuation gradient for identifying oil and gas reservoirs | |
| Wang et al. | Research on seismic wave attenuation in gas hydrates layer using vertical cable seismic data | |
| CN104570071A (zh) | 一种粘声介质贝叶斯线性ava和avf反演方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200630 Termination date: 20201218 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |