Zero-offset vertical seismic profile multiple attribute determination method
阅读说明：本技术 零偏垂直地震剖面多次波属性确定方法 (Zero-offset vertical seismic profile multiple attribute determination method ) 是由 李建国 严海滔 李清锋 赵继龙 于 2021-08-30 设计创作，主要内容包括：本发明提供一种零偏垂直地震剖面多次波属性确定方法,(1)输入零偏垂直地震剖面纵波动校正剖面数据Data,时窗参数Time1、Time2；其中,时窗参数Time2＞Time1；(2)用步骤(1)的数据Data、时窗参数Time1、时窗参数Time2确定分时窗叠加属性stack-(Time1)、stack-(Time2)；(3)用步骤(2)的分时窗叠加属性stack-(Time1)、stack-(Time2),确定多次波属性MultiWAtt；(4)解释步骤(3)的多次波属性MultiWAtt,以0.3为阈值识别多次波MultiW。本发明利用输入的垂直地震剖面纵波动校正剖面数据,分时窗叠加构建多次波属性,用来识别多次波,方法简洁、高效。(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 Time2 Time1 、stack Time2 (ii) a (3) Superposing attribute stack by using the time-sharing window in the step (2) Time1 、stack Time2 Determining 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.)
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.
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.