CN113703040A - Zero-offset vertical seismic profile multiple attribute determination method - Google Patents
Zero-offset vertical seismic profile multiple attribute determination method Download PDFInfo
- Publication number
- CN113703040A CN113703040A CN202111002080.0A CN202111002080A CN113703040A CN 113703040 A CN113703040 A CN 113703040A CN 202111002080 A CN202111002080 A CN 202111002080A CN 113703040 A CN113703040 A CN 113703040A
- Authority
- CN
- China
- Prior art keywords
- stack
- time2
- time1
- time
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 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
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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
The invention provides a zero-offset vertical seismic profile multiple attribute determination method which comprises the following steps of (1) inputting Data of a longitudinal fluctuation correction profile of a zero-offset vertical seismic profile, and Time window parameters of Time1 and Time 2; wherein, the Time window parameter Time2 is more than Time 1; (2) determining the Time-sharing window superposition attribute stack by using the Data in the step (1), the Time window parameter Time1 and the Time window parameter Time2Time1、stackTime2(ii) a (3) Superposing attribute stack by using the time-sharing window in the step (2)Time1、stackTime2Determining a multiple attribute MultiWAtt; (4) interpreting the multiple attribute MultiWAtt of step (3), and identifying multiple MultiW with 0.3 as a threshold. The method utilizes the input vertical seismic profile longitudinal fluctuation correction profile data and the time-sharing window superposition to construct the multiple attribute for identifying the multiple, and is simple and efficient.
Description
Technical Field
The invention relates to the field of geophysical exploration, which comprises zero-offset vertical seismic profile processing interpretation, multiple identification and attribute determination, in particular to a method for determining the multiple attribute of a zero-offset vertical seismic profile.
Background
Typically, zero-offset vertical seismic sections are used to determine formation velocities, layer position indices, and the like. The zero-offset vertical seismic section contains rich multiple information and can be used for multiple identification and analysis and guidance of ground seismic processing.
The literature and the patent do not find a similar multiple attribute determination method.
Disclosure of Invention
The invention provides a zero-offset vertical seismic profile multiple attribute determination method, which is used for constructing attribute data for identifying multiples on a zero-offset vertical seismic profile. The method utilizes the processing result longitudinal wave motion correction profile data of the zero-offset vertical seismic profile and time-sharing window superposition construction attributes to establish attributes capable of identifying multiples.
The process of the invention comprises the following steps:
the method for determining the multiple attribute of the zero-offset vertical seismic profile comprises the following steps:
(1) inputting Data of a longitudinal fluctuation correction section of a zero-deviation vertical seismic section, and Time window parameters Time1 and Time 2; wherein, the Time window parameter Time2 is more than Time 1;
(2) determining the Time-sharing window superposition attribute stack by using the Data in the step (1), the Time window parameter Time1 and the Time window parameter Time2Time1、stackTime2。
stackTime1＝∑Time1|Data|
stackTime2＝∑Time2|Data|
Wherein, Time1 and Time2 are Time window parameters, Data is vertical fluctuation correction section Data of zero-offset vertical seismic section, stackTime1、stackTime2Is the superposition property of Data in the Time windows Time1 and Time2, Σ is a summation function, and | is an absolute value function.
(3) Superposing attribute stack by using the time-sharing window in the step (2)Time1、stackTime2And determining a multiple attribute MultiWAtt.
MultiWAtt＝|stacktime2-stacktime1|/max(|stacktime2-stacktime1|)
Wherein, the stackTime1、stackTime2Is the Time window Time1,The superposition property of Data in Time2, | | is an absolute value function, and max is a maximum value function.
(4) Interpreting the multiple attribute MultiWAtt of step (3), and identifying multiple MultiW with 0.3 as a threshold.
Where MultiW is 1, 0 is not, and MultiWAtt is a multiple attribute.
The method utilizes the input vertical seismic profile longitudinal fluctuation correction profile data and the time-sharing window superposition to construct the multiple attribute for identifying the multiple, and is simple and efficient.
Drawings
FIG. 1 is a vertical seismic section longitudinal fluctuation correction section of an example input; the abscissa is the depth (unit: meter); the ordinate is time (unit: millisecond).
FIG. 2 is an example of multiple attributes and multiple interpretation; the horizontal axis of the left image is depth (unit: meter), and the vertical axis of the left image is time (unit: millisecond); the abscissa of the right graph is the number of tracks (unit: none) and the ordinate of the right graph is the time (unit: millisecond).
Detailed Description
The specific technical scheme of the invention is described by combining the embodiment.
(1) The Data of the vertical fluctuation correction section of the zero deviation vertical seismic section and the Time window parameters Time1 and Time2 are input.
Wherein, the Time window parameter Time2 > Time 1.
The input vertical seismic profile longitudinal fluctuation correction profile data with zero offset is shown in figure 1.
(2) Determining the Time-sharing window superposition attribute stack by using the Data in the step (1), the Time window parameter Time1 and the Time window parameter Time2Time1、stackTime2。
stackTime1＝∑Time1|Data|
stackTime2＝∑Time2|Data|
Wherein, Time1, Time2 are Time window parameters, Data are zero offset vertical seismic profile longitudinal fluctuation corrected profile Data, stackTime1、stackTime2Is the superposition property of Data in the Time windows Time1 and Time2, Σ is a summation function, and | is an absolute value function.
(3) Superposing attribute stack by using the time-sharing window in the step (2)Time1、stackTime2And determining a multiple attribute MultiWAtt.
MultiWAtt＝|stacktime2-stacktime1|/max(|stacktime2-stacktime1|)
Wherein, the stackTime1、stackTime2Is the superposition property of Data in Time windows Time1 and Time2, | | is an absolute value function, and max is a maximum value function.
(4) Interpreting the multiple attribute MultiWAtt of step (3), and identifying multiple MultiW with 0.3 as a threshold.
Where MultiW is 1, 0 is not, and MultiWAtt is a multiple attribute.
The determined multiple attributes and the multiple interpretation are shown in fig. 2.
Claims (3)
1. The method for determining the multiple attribute of the zero-offset vertical seismic section is characterized by comprising the following steps of:
(1) inputting Data of a longitudinal fluctuation correction section of a zero-deviation vertical seismic section, and Time window parameters Time1 and Time 2; wherein, the Time window parameter Time2 is more than Time 1;
(2) determining the Time-sharing window superposition attribute stack by using the Data in the step (1), the Time window parameter Time1 and the Time window parameter Time2Time1、stackTime2；
(3) Superposing attribute stack by using the time-sharing window in the step (2)Time1、stackTime2Determining a multiple attribute MultiWAtt;
(4) interpreting the multiple attribute MultiWAtt of the step (3), and identifying multiple MultiW by taking 0.3 as a threshold;
where MultiW is 1, 0 is not, and MultiWAtt is a multiple attribute.
2. The method for determining the multiple property of the zero-offset vertical seismic section according to claim 1, wherein the determining in the step (2) is as follows:
stackTime1＝∑Time1|Data|
stackTime2＝∑Time2|Data|
wherein, Time1 and Time2 are Time window parameters, Data is vertical fluctuation correction section Data of zero-offset vertical seismic section, stackTime1、stackTime2Is the superposition property of Data in the Time windows Time1 and Time2, Σ is a summation function, and | is an absolute value function.
3. The method for determining the multiple property of a zero-offset vertical seismic section according to claim 1, wherein the multiple property MultiWAtt of step (3) is;
MultiWAtt＝|stacktime2-stacktime1|/max(|stacktime2-stacktime1|)
wherein, the stackTime1、stackTime2Is the superposition property of Data in Time windows Time1 and Time2, | | is an absolute value function, and max is a maximum value function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111002080.0A CN113703040B (en) | 2021-08-30 | 2021-08-30 | Method for determining multiple wave attribute of zero offset vertical seismic section |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111002080.0A CN113703040B (en) | 2021-08-30 | 2021-08-30 | Method for determining multiple wave attribute of zero offset vertical seismic section |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113703040A true CN113703040A (en) | 2021-11-26 |
| CN113703040B CN113703040B (en) | 2023-12-01 |
Family
ID=78656588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111002080.0A Active CN113703040B (en) | 2021-08-30 | 2021-08-30 | Method for determining multiple wave attribute of zero offset vertical seismic section |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113703040B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| CN105093292A (en) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | Data processing method and device for seismic imaging |
| CN113090251A (en) * | 2021-04-14 | 2021-07-09 | 中油奥博(成都)科技有限公司 | Logging VSP composite data acquisition system based on optical fiber sensing and acquisition processing method |
-
2021
- 2021-08-30 CN CN202111002080.0A patent/CN113703040B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| CN105093292A (en) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | Data processing method and device for seismic imaging |
| CN113090251A (en) * | 2021-04-14 | 2021-07-09 | 中油奥博(成都)科技有限公司 | Logging VSP composite data acquisition system based on optical fiber sensing and acquisition processing method |
Non-Patent Citations (2)
| Title |
|---|
| 李大卫 等: "《地震地质一体化研究中的地震数据处理质量监控方法综述》", 《勘探地球物理进展》, vol. 33, no. 3, pages 160 - 167 * |
| 赵海英 等: "《基于VSP的地震层位综合标定方法》", 《石油地球物理勘探》, vol. 51, pages 84 - 92 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113703040B (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20230333274A1 (en) | Intelligent real-time updating method and system for stratigraphic framework with geosteering-while-drilling | |
| CN107688201B (en) | RBM-based seismic prestack signal clustering method | |
| US10539008B2 (en) | Systems and methods for determining tool center, borehole boundary, and/or mud parameter | |
| US11353616B2 (en) | Methods and systems for processing slowness values from borehole sonic data | |
| Zimmermann et al. | Machine-learning-based automatic well-log depth matching | |
| Tian et al. | Deep learning assisted well log inversion for fracture identification | |
| CN111583407B (en) | Efficient three-dimensional geological modeling intelligent processing method based on paper drilling | |
| Scheidt et al. | A multi-resolution workflow to generate high-resolution models constrained to dynamic data | |
| Gonzalez et al. | Integrated multi-physics workflow for automatic rock classification and formation evaluation using multi-scale image analysis and conventional well logs | |
| CN111914478A (en) | Comprehensive geological drilling well logging lithology identification method | |
| Chang et al. | SegLog: Geophysical logging segmentation network for lithofacies identification | |
| US20150081259A1 (en) | Method of calibrating a geologic model | |
| Pan et al. | Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models | |
| Chen et al. | Integration of cumulative-distribution-function mapping with principal-component analysis for the history matching of channelized reservoirs | |
| CN113703040A (en) | Zero-offset vertical seismic profile multiple attribute determination method | |
| US6859733B1 (en) | Method for determining parameters of earth formations surrounding a well bore | |
| Torres Caceres et al. | Automated well-log depth matching–1D convolutional neural networks vs. classic cross correlation | |
| EP4196825B1 (en) | Machine learning-based differencing tool for hydrocarbon well logs | |
| US20210293132A1 (en) | Process for real time geological localization with greedy monte carlo | |
| Dhara et al. | Seismic image registration using multiscale convolutional neural networks | |
| CN114581556B (en) | Digital map filling method in regional geological survey | |
| US11156740B2 (en) | Electrofacies determination | |
| CN116068644A (en) | Method for improving resolution and noise reduction of seismic data by using generation countermeasure network | |
| CN112987091B (en) | Reservoir detection method, device, electronic equipment and storage medium | |
| Parimontonsakul et al. | A Machine Learning Based Approach to Automate Stratigraphic Correlation through Marker Determination |
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 |