CN106772610B - 一种叠前道集资料增维属性分析方法 - Google Patents
一种叠前道集资料增维属性分析方法 Download PDFInfo
- Publication number
- CN106772610B CN106772610B CN201610998827.5A CN201610998827A CN106772610B CN 106772610 B CN106772610 B CN 106772610B CN 201610998827 A CN201610998827 A CN 201610998827A CN 106772610 B CN106772610 B CN 106772610B
- Authority
- CN
- China
- Prior art keywords
- gather
- dimensional
- cdp
- prestack
- imaging point
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000004141 dimensional analysis Methods 0.000 claims abstract description 9
- 210000001233 cdp Anatomy 0.000 claims description 34
- 238000004637 computerized dynamic posturography Methods 0.000 claims description 34
- 230000008859 change Effects 0.000 claims description 4
- 238000004448 titration Methods 0.000 claims 1
- 230000000007 visual effect Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000008569 process 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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)
- Image Analysis (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种叠前道集资料增维属性分析方法包括：二维叠前共成像点道集参数读取，包括CDP点坐标、间距和道集入射角度间距或道集偏移间距；叠前共成像点道集三维数据体形成方法，读取CDP上的叠前共成像点道集作为三维数据体的Inline线，赋值Inline线间距、CDP间距和坐标，形成三维数据体；叠前共成像点道集三维分析，加载三维数据体到地震解释平台，实现叠前共成像点道集三维分析。该方法与目前常用的叠前共成像点道集属性分析相比，具有快捷、方便、结果更直观，实现了任意CDP点、不同入射角（偏移距）道集的直接对比。
Description
技术领域
本发明属于地球物理勘探地震资料信号分析技术领域，主要是地震资料信号分析中的一种叠前道集资料增维属性分析方法。
背景技术
地震资料信号分析具有相对完善的技术体系，常用的地震解释平台均能够实现数据体的三维分析，而对于叠前共成像点道集资料，受数据记录形式的限制，目前主要针对单个CDP点的道集进行二维信号分析，效率低、道集间的关系分析结果不直观，制约了叠前共成像点道集信息的应用。为实现叠前共成像点道集的高效、全方位属性分析，将二维线的叠前共成像点道集改写成三维形式，再输入地震解释平台，就可进行叠前共成像点道集的全方位属性分析。
发明内容
本发明的目的是针对现存方法的不足之处，发明了一种叠前道集资料增维属性分析方法。该发明在将二维线的叠前共成像点道集改写成三维形式后加载到地震解释平台，快速实现了多维度的道集属性分析、道集间属性对比。
本发明依据地震叠前时间偏移成像道集数据特征，提出将CDP点上的道集数据定义为Inline线，相应的CDP点就组成了Crossline线，从而将一条二维线的地震叠前时间偏移成像道集数据转化为三维数据体，将该三维数据体加载到地震解释平台就可以实现成像道集数据的多方位定量直观属性分析。
本发明的主要流程：
步骤1：二维叠前共成像点道集参数读取，包括CDP点坐标、间距和道集入射角度间距或道集偏移间距；
步骤2：叠前共成像点道集三维数据体形成方法，读取CDP上的叠前共成像点道集作为三维数据体的Inline线，赋值Inline线间距、CDP间距和坐标，形成三维数据体；
步骤3：叠前共成像点道集三维分析，加载三维数据体到地震解释平台，实现叠前共成像点道集三维分析。
上述方案还包括：
步骤2的三维数据体的Inline线为二维线的不同CDP上的叠前共成像点道集，Inline线和CDP间距都赋值二维线的CDP间距，Inline线的首CDP点坐标为对应的原始二维线CDP点的坐标。
对叠前共成像点道集进行体分析包括：刻画所有道集的地震属性随入射角或偏移距变化特征、不同道集和偏移距间的关系。
所述地震解释平台具有叠前共成像点道集的振幅、频率、波形类属性的分析，同时具有叠前共成像点道集间属性关系分析和不同CDP点不同入射角度的属性关系分析。
本发明的主旨在于改写叠前共成像点道集的记录为三维形式，加载到分析手段完善的解释平台，实现多维度的道集属性分析、道集间属性对比。本发明的三维叠前共成像点道集形式能够批量分析道集属性，直观表达道集间的关系，比二维形式的叠前共成像点道集分析具有明显的优越性。
附图说明
图1: 叠前共成像点道集增维分析流程图；
图2: 二维楔形模型叠前共成像点角度域道集常规显示图；
图3: 厚度17米的叠前共成像点角度域道集均方根属性分析图；
图4: 厚度42米的叠前共成像点角度域道集均方根属性分析图；
图5: 入射角0度的叠前共成像点共角度剖面；
图6: 入射角40度的叠前共成像点共角度剖面；
图7: 叠前共成像点道集均方根属性分析平面图；
图8: 沿等均方根值切线图。
具体实施方式
下面结合附图1，以二维楔形地震模型的叠前共成像点角度域道集为例，对叠前道集资料增维属性分析方法进行详细说明。
目前常用的做法，对二维楔形数学模型采用波动方程正演和成像后，得到地震叠前共成像点道集，其记录形式如图2所示，在一条二维线上将叠前共成像点道集按CDP点顺序排放，可以取出任意一个CDP点上的叠前共成像点道集进行二维分析。这样就可以轻松实现对单个CDP点上的叠前共成像点道集分析，但如果要对所有的叠前共成像点道集分析，工作量就非常大且不同叠前共成像点道集间的对比分析的结果也不直观。
为了更为快捷、直观地分析叠前共成像点道集的属性特征，依据叠前共成像点道集的数据特征，构想形成三维格式记录叠前共成像点道集。二维的叠前共成像点道集的数据特征是每个CDP点上的叠前共成像点道集各自独立，反映的是反射点信息随入射角的变化，相互之间没有联系。三维数据格式具有实现数据间相互联系的独特优势，如果将二维的叠前共成像点道集转化成三维形成，就能实现快捷、直观的数据分析。
构造叠前共成像点道集三维数据形式，将每个CDP点上的叠前共成像点道集读入作为Inline线，设定Inline线间距和CDP间距为原始二维线的CDP间距，首CDP点坐标为对应的原始二维线CDP点的坐标，那么Crossline线则是相同入射角的地震数据，另外，任意切线还可以反映不同CDP点、不同入射角的地震信息。这样就将一条二维线的叠前共成像点道集转化为一个三维数据体，三维数据体中线间距和CDP间距相同是为了显示方便，分析时可将CDP间距与实际入射角度间距（或道集偏移间距）换算。
加载三维数据体到地震解释平台，实现叠前共成像点道集三维分析。利用地震解释平台可快速实现所有叠前共成像点道集的振幅、频率、波形类属性的分析，同时实现了叠前共成像点道集间属性关系分析和不同CDP点不同入射角度的属性关系分析。以沿层均方根振幅属性提取说明叠前共成像点道集增维分析的优势。对楔形模型顶面反射提取均方根振幅，可以得到三维数据体所有Inline、Crossline线的均方根振幅信息。这样就对所有叠前共成像点道集进行了均方根振幅属性分析。图3、4是任意抽取2条Inline线，反映不同厚度的模型随入射角增大的均方根振幅减小的特征，但两者之间定量对比关系不明确。这是现有方法能够实现的，但增维分析方法效率明显提高。图5、6是任意抽取的2条Crossline线，反映不同入射角的模型随厚度变化的均方根振幅变化特征，相同厚度不同入射角的模型间的均方根振幅对比关系不明确。这是现有方法能够实现的，但增维分析方法效率明显提高。图7是均方根振幅属性随厚度和入射角变化的平面特征，可以反映不同厚度的模型随入射角增大的均方根振幅变化的定量关系，反映相同厚度不同入射角的模型间的均方根振幅变化的定量关系。厚度小于1/4波长时，均方根振幅随入射角增大变化缓慢；大于1/4波长时，均方根振幅随入射角增大变化快；大于1/2波长时，均方根振幅随入射角增大变化较快，相对平稳。这是增维分析方法的优势，现有方法不能实现的。图8是沿图7平面图中等均方根振幅值切的一条联络线，反映出由不同的厚度、不同的入射角度可以有相同的均方根振幅值，是多对一的映射关系。
通过上述流程实现了叠前共成像点道集增维分析，展现了该方法全面、快捷属性分析特征，实现了现有方法不能实现的叠前共成像点道属性平面分析和任意切线分析，发现了地震属性随厚度、入射角变化的新规律。
Claims (3)
1.一种叠前道集资料增维属性分析方法，其特征在于：
步骤1：二维叠前共成像点道集参数读取，包括CDP点坐标、间距和道集入射角度间距或道集偏移间距；
步骤2：叠前共成像点道集三维数据体形成方法，读取CDP上的叠前共成像点道集作为三维数据体的Inline线，赋值Inline线间距、CDP间距和坐标，形成三维数据体；所述三维数据体的Inline线为二维线的不同CDP上的叠前共成像点道集，Inline线和CDP间距都赋值二维线的CDP间距，Inline线的首CDP点坐标为对应的原始二维线CDP点的坐标；
步骤3：叠前共成像点道集三维分析，加载三维数据体到地震解释平台，实现叠前共成像点道集三维分析。
2.根据权利要求1所述的叠前道集资料增维属性分析方法，其特征在于，对叠前共成像点道集进行体分析包括：刻画所有道集的地震属性随入射角或偏移距变化特征、不同道集和偏移距间的关系。
3.根据权利要求1或2所述的叠前道集资料增维属性分析方法，其特征在于：所述地震解释平台具有叠前共成像点道集的振幅、频率、波形类属性的分析，同时具有叠前共成像点道集间属性关系分析和不同CDP点不同入射角度的属性关系分析。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610998827.5A CN106772610B (zh) | 2016-11-14 | 2016-11-14 | 一种叠前道集资料增维属性分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610998827.5A CN106772610B (zh) | 2016-11-14 | 2016-11-14 | 一种叠前道集资料增维属性分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106772610A CN106772610A (zh) | 2017-05-31 |
| CN106772610B true CN106772610B (zh) | 2020-05-01 |
Family
ID=58973620
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610998827.5A Active CN106772610B (zh) | 2016-11-14 | 2016-11-14 | 一种叠前道集资料增维属性分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106772610B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10557954B2 (en) * | 2017-06-12 | 2020-02-11 | Saudi Arabian Oil Company | Modeling angle domain common image gathers from reverse time migration |
| CN108333625A (zh) * | 2018-02-11 | 2018-07-27 | 中国石油化工股份有限公司 | 时频域叠前共成像点道集属性三维分析方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5587942A (en) * | 1994-03-02 | 1996-12-24 | Exxon Production Research Company | 3D wave equation migration of a 2D grid of seismic data |
| CN104297789A (zh) * | 2014-10-23 | 2015-01-21 | 中国科学院地质与地球物理研究所 | 一种三维倾角域稳相叠前时间偏移方法及系统 |
| CN105372718A (zh) * | 2015-12-01 | 2016-03-02 | 西安石油大学 | 一种2d位场构建3d位场的方法 |
| CN105425306A (zh) * | 2015-12-01 | 2016-03-23 | 西安石油大学 | 基于2d离散小波多尺度分解构建3d位场的方法 |
| CN105445801A (zh) * | 2014-09-01 | 2016-03-30 | 中国石油化工股份有限公司 | 一种消除二维地震资料随机噪音的处理方法 |
-
2016
- 2016-11-14 CN CN201610998827.5A patent/CN106772610B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5587942A (en) * | 1994-03-02 | 1996-12-24 | Exxon Production Research Company | 3D wave equation migration of a 2D grid of seismic data |
| CN105445801A (zh) * | 2014-09-01 | 2016-03-30 | 中国石油化工股份有限公司 | 一种消除二维地震资料随机噪音的处理方法 |
| CN104297789A (zh) * | 2014-10-23 | 2015-01-21 | 中国科学院地质与地球物理研究所 | 一种三维倾角域稳相叠前时间偏移方法及系统 |
| CN105372718A (zh) * | 2015-12-01 | 2016-03-02 | 西安石油大学 | 一种2d位场构建3d位场的方法 |
| CN105425306A (zh) * | 2015-12-01 | 2016-03-23 | 西安石油大学 | 基于2d离散小波多尺度分解构建3d位场的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 三维叠后地震数据可纵向显示的处理方法;苑益军 等;《地球物理学进展》;20081031;第23卷(第5期);第1532-1537页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106772610A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Maggi et al. | An automated time-window selection algorithm for seismic tomography | |
| WO2017024702A1 (zh) | 一种射线弹性参数的反演系统 | |
| EP3158368B1 (en) | Instantaneous isochron attribute-based geobody identification for reservoir modeling | |
| US10073190B2 (en) | Method and system for geophysical modeling of subsurface volumes based on computed vectors | |
| US20130294197A1 (en) | Automatic extraction and characterization of fault and fracture populations | |
| US11231513B2 (en) | Systems and methods for refining estimated parameter values in seismic imaging | |
| AU2019248523B2 (en) | Systems and methods for using probabilities of lithologies in an inversion | |
| CA2847782C (en) | Method and system of subsurface horizon assignment | |
| AU2017366758A1 (en) | Diving wave illumination using migration gathers | |
| CN106772610B (zh) | 一种叠前道集资料增维属性分析方法 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| AU2019246947B2 (en) | Systems and methods for refining estimated effects of parameters on amplitudes | |
| CN106468171A (zh) | 多信息耦合的递进式砂砾岩岩相判别方法 | |
| Zhang et al. | Automatic velocity analysis using interpretable multimode neural networks | |
| Lorentzen et al. | Mapping Cretaceous faults using a convolutional neural network-A field example from the Danish North Sea. | |
| CN108761533B (zh) | 一种确定纵横波速度比的方法、装置及系统 | |
| CN110954956B (zh) | 观测系统的采集痕迹评价方法及计算机可读存储介质 | |
| CN108761535B (zh) | 一种侵入岩分布的识别方法 | |
| US20160274269A1 (en) | Geocellular Modeling | |
| AU2013406187A1 (en) | Geocellular modeling |
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 |