CN110907989A - Method and system for reconstructing quasi-ground seismic reflection wave imaging - Google Patents
Method and system for reconstructing quasi-ground seismic reflection wave imaging Download PDFInfo
- Publication number
- CN110907989A CN110907989A CN201811082843.5A CN201811082843A CN110907989A CN 110907989 A CN110907989 A CN 110907989A CN 201811082843 A CN201811082843 A CN 201811082843A CN 110907989 A CN110907989 A CN 110907989A
- Authority
- CN
- China
- Prior art keywords
- virtual source
- seismic
- receiving point
- representing
- record
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005284 excitation Effects 0.000 claims description 19
- 238000000926 separation method Methods 0.000 claims description 9
- 238000005286 illumination Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 11
- 230000008901 benefit Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 240000006829 Ficus sundaica Species 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000004613 tight binding model 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
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 method and system for reconstructing the pseudo-ground seismic reflection wave image are disclosed. The method comprises the following steps: step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions; step 2: calculating a virtual source seismic channel of a virtual source position relative to a virtual source receiving point position; and step 3: repeating the step 2 aiming at each virtual source receiving point position, obtaining a virtual source seismic channel of the virtual source position aiming at each virtual source receiving point position, and obtaining a common virtual source channel set; and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data; and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result. According to the method, the underground illumination range of the conventional VSP data can be effectively expanded by reconstructing the VSP data into the simulated ground seismic data, and the description capability of horizon comparison and tracking and reservoir prediction is greatly enhanced.
Description
Technical Field
The invention relates to the field of seismic exploration and development, in particular to a method and a system for reconstructing simulated ground seismic reflection wave imaging.
Background
The imaging range of conventional VSPs is a limited imaging area of a conical shape, as shown in fig. 1, and thus the description capability for horizon comparison and tracking, as well as in reservoir prediction, will be limited. To expand the imaging range, the multiple information in VSP data can be utilized for direct imaging, but there are many practical difficulties in the implementation of the method technology.
At present, two receivers are often used to record the cross-correlation of wavefields to recover the seismic response at the two receiving points as if received at one of the receiving points as an excitation (virtual source) and the other receiving point, which is collectively referred to as seismic interferometry. The green's function reconstruction based on the correlation-type interchange equation proposed by Wapenaar (2008) can be expressed as:
wherein u isi(xAT) and ui(xBT) represents the seismic records received at surface locations A and B, respectively, representing convolution, i representing the ith excitation source, N representing the number of all sources, G (x)B,xAT) denotes the Green function excited at position A (virtual source), received at position B, G (x)B,xAT) represents the inverse time green's function excited at position a (virtual source), received at position B, and s (t) represents the virtual source wavelet. However, in the method, a false in-phase axis is introduced in the reconstruction process, the amplitude energy of the reconstruction result is weakened, and the signal-to-noise ratio is low. Therefore, there is a need to develop a method and system for reconstructing a pseudo-surface seismic reflection wave.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a method and a system for reconstructing simulated ground seismic reflection wave imaging, which can effectively expand the underground illumination range of conventional VSP data by reconstructing VSP data into simulated ground seismic data, and greatly enhance the description capability of horizon comparison and tracking and reservoir prediction.
According to one aspect of the invention, a method for reconstructing a quasi-surface seismic reflection wave image is provided. The method may include: step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions; step 2: calculating a virtual source seismic trace of the virtual source position relative to a virtual source receiving point position; and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set; and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data; and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
Preferably, the method further comprises the following steps: and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
Preferably, the method further comprises the following steps: and obtaining a plurality of actual seismic sources according to the RVSP seismic data.
Preferably, the step 2 includes: respectively calculating the correlation value of the seismic record of the virtual source position and the seismic record of a virtual source receiving point position under each actual seismic source excitation; and summing the calculated correlation values for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
Preferably, the correlation value between the seismic record of the virtual source position and the seismic record of each virtual source receiving point position is calculated by formula (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
According to another aspect of the invention, a system for reconstructing a pseudo-surface seismic reflection wave imaging is provided, the system comprising: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions; step 2: calculating a virtual source seismic trace of the virtual source position relative to a virtual source receiving point position; and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set; and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data; and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
Preferably, the method further comprises the following steps: and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
Preferably, the method further comprises the following steps: and obtaining a plurality of actual seismic sources according to the RVSP seismic data.
Preferably, the step 2 includes: respectively calculating the correlation value of the seismic record of the virtual source position and the seismic record of a virtual source receiving point position under each actual seismic source excitation; and summing the calculated correlation values for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
Preferably, the correlation value between the seismic record of the virtual source position and the seismic record of each virtual source receiving point position is calculated by formula (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
FIG. 1 shows a schematic diagram of a limited imaging region of a VSP primary reflected wave according to one embodiment of the present invention.
FIG. 2 is a flow chart illustrating the steps of a method of reconstructing a pseudoground seismic reflection imaging in accordance with the present invention.
FIG. 3 shows a schematic diagram of a principle schematic of the reconstruction of pseudo-surface seismic data according to one embodiment of the invention.
FIG. 4 shows a schematic diagram of a VSP forward model according to one embodiment of the invention.
FIG. 5 shows a schematic diagram of VSP shot gather records according to one embodiment of the present invention.
FIG. 6 shows a schematic diagram of a comsumption source gather, according to one embodiment of the present invention.
FIG. 7 shows a schematic diagram of reconstructed seismic imaging results according to an embodiment of the invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
FIG. 1 shows a flow chart of the steps of a method of reconstructing a pseudoground seismic reflection imaging according to the invention.
In this embodiment, the method for reconstructing a pseudo-surface seismic reflection wave imaging according to the present invention may include: step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions; step 2: calculating a virtual source seismic channel of a virtual source position relative to a virtual source receiving point position; and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set; and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data; and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
In one example, further comprising: and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
In one example, further comprising: from the RVSP seismic data, a plurality of actual seismic sources are obtained.
In one example, step 2 comprises: respectively calculating the correlation value of the seismic record at the virtual source position and the seismic record at the virtual source receiving point position under each actual seismic source excitation; and summing the correlation values calculated for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
In one example, the correlation value of the seismic record of the virtual source location and the seismic record of each virtual source receiver location is separately calculated by equation (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
Specifically, the primary reflection wave of the quasi-surface earthquake is obtained by correlating the direct wave and the multiple wave in the VSP data. Based on the seismic interference principle, the reconstruction formula of the quasi-ground seismic reflection wave is shown as formula (1).
FIG. 3 shows a schematic diagram of a principle schematic of the reconstruction of pseudo-surface seismic data according to one embodiment of the invention.
The meaning of equation (1) is intuitively understood from a mathematical point of view, namely: the cross-correlation operation represents the phase subtraction of two operands, i.e.: the arrival time of the multiple wave received at B minus the arrival time of the direct wave received at a results in the excitation (virtual source) of the primary reflected wave received at a at B, as shown in fig. 3, the vertical line represents the actual seismic source, and the horizontal line represents the virtual source position. Due to the received wavefield complexity, the reconstructed seismic response contains not only the primary of interest, but also some artifacts, which may be resolved by wavefield separation or other methods prior to wavefield cross-correlation.
The method for reconstructing the quasi-surface seismic reflection wave imaging according to the invention can comprise the following steps:
step 1: VSP seismic data are converted into RVSP seismic data, wave field separation is carried out to obtain initial seismic data, a plurality of actual seismic sources are obtained according to the RVSP seismic data, a certain receiving point position on the earth surface is set as a virtual source position, and the rest receiving point positions are set as virtual source receiving point positions.
Step 2: respectively calculating the correlation value of the seismic record at the virtual source position and the seismic record at the virtual source receiving point position under each actual seismic source excitation through a formula (1); and summing the correlation values calculated for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
And step 3: and (3) repeating the step (2) aiming at the position of each virtual source receiving point to obtain a virtual source seismic channel of the virtual source position for the position of each virtual source receiving point, and further obtain a common virtual source channel set.
And 4, step 4: and (4) replacing the virtual source position, and repeating the steps 1-3 to obtain the reconstructed simulated ground seismic data.
And 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
According to the method, the VSP data are reconstructed into the simulated ground seismic data, the underground illumination range of the conventional VSP data can be effectively expanded, and the description capability of the method for horizon comparison and tracking and reservoir prediction is greatly enhanced.
Application example
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, a specific application example is given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
FIG. 4 shows a schematic diagram of a VSP forward model according to one embodiment of the invention, with horizontal axis representing horizontal distance in m and vertical axis representing depth in m.
Establishing a VSP forward model as shown in FIG. 4, model size: 2000m 3000m, background velocity 2000m/s, wedge velocity 2400m/s, model size: 2000 x 3000; background speed 2000 m/s; the wedge velocity was 2400m/s and was located in the models (0,500), (1000,500), (1000 ). Surface shot point: 10m, 201 shots in total, receiving in the well: 1000-2000, with a spacing of 10m, for a total of 101. Note that the wedge is located above the receiver in the well, which is outside the VSP reflected wave imaging range.
FIG. 5 shows a schematic diagram of VSP shot gather records with the horizontal axis representing number of tracks and the vertical axis representing time in units of s, according to one embodiment of the invention.
VSP data is generated by utilizing the sound wave finite difference wave field simulation, and the parameters are as follows: the use of Rake wavelets, 40Hz dominant frequency. The sampling rate is 0.5ms and the recording length is 2 s. FIG. 5 shows a VSP shot gather record generated by forward modeling, with the principal co-axial being the direct arrival and the surface-related multiples.
The method for reconstructing the quasi-ground seismic reflection wave imaging comprises the following steps:
step 1: VSP seismic data are converted into RVSP seismic data, wave field separation is carried out to obtain initial seismic data, a plurality of actual seismic sources are obtained according to the RVSP seismic data, a certain receiving point position on the earth surface is set as a virtual source position, and the rest receiving point positions are set as virtual source receiving point positions.
Step 2: respectively calculating the correlation value of the seismic record at the virtual source position and the seismic record at the virtual source receiving point position under each actual seismic source excitation through a formula (1); and summing the correlation values calculated for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
FIG. 6 shows a schematic diagram of a comsumption source gather, according to one embodiment of the present invention.
And step 3: for each virtual source receiving point position, repeating step 2, obtaining a virtual source seismic trace of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source gather, as shown in fig. 6, where the virtual source position of the common virtual source gather is located at the middle position of the survey line, as can be seen from the figure: the reflections from the horizontal and inclined interfaces of the wedge are clearly discernible.
And 4, step 4: and (4) replacing the virtual source position, and repeating the steps 1-3 to obtain the reconstructed simulated ground seismic data.
FIG. 7 shows a schematic diagram of reconstructed seismic imaging results with the horizontal axis representing the number of traces and the vertical axis representing time in units of s, according to one embodiment of the invention.
And 5: imaging is performed on the reconstructed pseudo-surface seismic data to obtain a reconstructed seismic imaging result, as shown in fig. 7, compared with fig. 4, the horizontal interface and the inclined interface of the wedge are well imaged, but the vertical interface cannot be imaged, because only reflected waves from the vertical interface are used as interference waves in the process of reconstructing the pseudo-surface seismic data in VSP observation. This result indicates that: after VSP data are reconstructed into simulated ground seismic data, the imaging range can be expanded, and structures above a receiver in a well can be imaged.
In conclusion, the invention can effectively expand the underground illumination range of the conventional VSP data by reconstructing the VSP data into the simulated ground seismic data, and the description capability for horizon comparison and tracking and in reservoir prediction is greatly enhanced.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
According to an embodiment of the invention, there is provided a system for reconstructing a pseudo-ground seismic reflection wave imaging, the system including: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions; step 2: calculating a virtual source seismic channel of a virtual source position relative to a virtual source receiving point position; and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set; and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data; and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
In one example, further comprising: and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
In one example, further comprising: from the RVSP seismic data, a plurality of actual seismic sources are obtained.
In one example, step 2 comprises: respectively calculating the correlation value of the seismic record at the virtual source position and the seismic record at the virtual source receiving point position under each actual seismic source excitation; and summing the correlation values calculated for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
In one example, the correlation value of the seismic record of the virtual source location and the seismic record of each virtual source receiver location is separately calculated by equation (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
The system can effectively expand the underground illumination range of the conventional VSP data by reconstructing the VSP data into simulated ground seismic data, and greatly enhances the description capability of horizon comparison and tracking and reservoir prediction.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. A method for reconstructing a quasi-surface seismic reflection wave image is characterized by comprising the following steps:
step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions;
step 2: calculating a virtual source seismic trace of the virtual source position relative to a virtual source receiving point position;
and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set;
and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data;
and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
2. The method of reconstructing surface-like seismic reflection imaging as claimed in claim 1, further comprising:
and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
3. The method of reconstructing surface-like seismic reflection imaging as claimed in claim 2, further comprising:
and obtaining a plurality of actual seismic sources according to the RVSP seismic data.
4. A method of reconstructing surface-like seismic reflection imaging as claimed in claim 3, wherein said step 2 includes:
respectively calculating the correlation value of the seismic record of the virtual source position and the seismic record of a virtual source receiving point position under each actual seismic source excitation;
and summing the calculated correlation values for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
5. The method of reconstructing pseudoground seismic reflection imaging according to claim 4, wherein the correlation value of the seismic record of the virtual source location and the seismic record of each virtual source receiver location is calculated separately by equation (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
6. A system for reconstructing a pseudoground seismic reflection imaging, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
step 1: determining initial seismic data, setting a certain receiving point position of the earth surface as a virtual source position, and setting other receiving point positions as virtual source receiving point positions;
step 2: calculating a virtual source seismic trace of the virtual source position relative to a virtual source receiving point position;
and step 3: repeating the step 2 aiming at each virtual source receiving point position to obtain a virtual source seismic channel of the virtual source position for each virtual source receiving point position, and further obtaining a common virtual source channel set;
and 4, step 4: replacing the virtual source position, and repeating the steps 1-3 to obtain reconstructed simulated ground seismic data;
and 5: and imaging the reconstructed simulated ground seismic data to obtain a reconstructed seismic imaging result.
7. The system of claim 6, further comprising:
and converting the VSP seismic data into RVSP seismic data, and performing wave field separation to obtain initial seismic data.
8. The system of claim 7, further comprising:
and obtaining a plurality of actual seismic sources according to the RVSP seismic data.
9. The system for reconstructing pseudoground seismic reflection imaging according to claim 8, wherein said step 2 comprises:
respectively calculating the correlation value of the seismic record of the virtual source position and the seismic record of a virtual source receiving point position under each actual seismic source excitation;
and summing the calculated correlation values for each actual seismic source excitation to obtain a virtual source seismic trace of the virtual source position for a virtual source receiving point position.
10. The system of claim 9, wherein the correlation value between the seismic record of the virtual source location and the seismic record of each virtual source receiver location is calculated by equation (1):
wherein d (x)B,xA,t)ReflectionRepresenting the reflection seismic record excited at surface location a and received at location B, i.e., the correlation of the seismic record at virtual source location a with the seismic record at virtual source receiver point location B,
representing cross-correlation, i representing the ith actual source in the well, N representing the number of all sources, d (x)A,t)directDirect wave, d (x), representing VSP acquisitionB,t)ghostRepresenting multiples of VSP acquisitions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811082843.5A CN110907989A (en) | 2018-09-17 | 2018-09-17 | Method and system for reconstructing quasi-ground seismic reflection wave imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811082843.5A CN110907989A (en) | 2018-09-17 | 2018-09-17 | Method and system for reconstructing quasi-ground seismic reflection wave imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110907989A true CN110907989A (en) | 2020-03-24 |
Family
ID=69813137
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811082843.5A Pending CN110907989A (en) | 2018-09-17 | 2018-09-17 | Method and system for reconstructing quasi-ground seismic reflection wave imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110907989A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111257939A (en) * | 2020-03-26 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source bidirectional wave field reconstruction method and system |
| CN111257938A (en) * | 2020-03-25 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source wave field reconstruction method and system based on wavelet cross-correlation |
| CN112505780A (en) * | 2020-10-27 | 2021-03-16 | 中国石油天然气集团有限公司 | Method and device for correcting formation depth data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
| CN103163554A (en) * | 2013-02-04 | 2013-06-19 | 西安交通大学 | Self-adapting wave form retrieval method through utilization of zero offset vertical seismic profile (VSP) data to estimate speed and Q value |
| CN106291684A (en) * | 2015-06-06 | 2017-01-04 | 中国石油化工股份有限公司 | The seismic response of a kind of blind focus earthquake wave field recovers and virtual source road collection construction method |
| CN106772616A (en) * | 2016-12-20 | 2017-05-31 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of seismic imaging |
| CN107402405A (en) * | 2016-05-18 | 2017-11-28 | 中国石油化工股份有限公司 | Quiet phase virtual source trace gather construction method |
-
2018
- 2018-09-17 CN CN201811082843.5A patent/CN110907989A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
| CN103163554A (en) * | 2013-02-04 | 2013-06-19 | 西安交通大学 | Self-adapting wave form retrieval method through utilization of zero offset vertical seismic profile (VSP) data to estimate speed and Q value |
| CN106291684A (en) * | 2015-06-06 | 2017-01-04 | 中国石油化工股份有限公司 | The seismic response of a kind of blind focus earthquake wave field recovers and virtual source road collection construction method |
| CN107402405A (en) * | 2016-05-18 | 2017-11-28 | 中国石油化工股份有限公司 | Quiet phase virtual source trace gather construction method |
| CN106772616A (en) * | 2016-12-20 | 2017-05-31 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of seismic imaging |
Non-Patent Citations (4)
| Title |
|---|
| KURANG MEHTA ET AL.: "Improving the virtual source method by wavefield separation", 《GEOPHYSICS》 * |
| 倪雪灿: "VSP数据地震干涉成像方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 曹辉等: "虚源法地震技术综述", 《地球物理学进展》 * |
| 陈国金等: "虚源法地震技术及数值模型试验", 《地球物理学进展》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111257938A (en) * | 2020-03-25 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source wave field reconstruction method and system based on wavelet cross-correlation |
| CN111257939A (en) * | 2020-03-26 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source bidirectional wave field reconstruction method and system |
| CN112505780A (en) * | 2020-10-27 | 2021-03-16 | 中国石油天然气集团有限公司 | Method and device for correcting formation depth data |
| CN112505780B (en) * | 2020-10-27 | 2024-05-28 | 中国石油天然气集团有限公司 | Formation depth data correction method and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2018201492B2 (en) | Method and system for separating seismic sources in marine simultaneous shooting acquisition | |
| CN102298156B (en) | For the method and apparatus of deghosting geological data | |
| Toomey et al. | Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 9° 30′ N | |
| US11835672B2 (en) | Non-uniform optimal survey design principles | |
| US9405027B2 (en) | Attentuating noise acquired in an energy measurement | |
| US9541659B2 (en) | Noise removal from 3D seismic representation | |
| US20140078860A1 (en) | Interference noise attenuation method and apparatus | |
| WO2013080128A1 (en) | Separation of simultaneous source data | |
| EP2180345A1 (en) | Removing seismic interference using simultaneous or near simultaneous source separation | |
| EP2823336B1 (en) | Methods and computing systems for processing data | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| CN110907989A (en) | Method and system for reconstructing quasi-ground seismic reflection wave imaging | |
| US20150234066A1 (en) | Method and device for de-blending seismic data using source signature | |
| EP2166379A2 (en) | Interbed seismic multiple prediction | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN113885079A (en) | Elastic wave field decoupling-based high-precision multi-azimuth reverse-time seismic source imaging method | |
| AU2015200555B2 (en) | Correction of sea surface state | |
| CN113740906A (en) | Underwater vertical cable seismic wave interference imaging method and device | |
| CN110967734A (en) | Virtual source reconstruction method and system based on fast Fourier transform | |
| CN108345032A (en) | A kind of weak illumination region high s/n ratio offset imaging method | |
| CN113866821B (en) | Passive source interference offset imaging method and system based on illumination direction constraint | |
| CN110618459B (en) | Seismic data processing method and device | |
| WO2002023222A1 (en) | Illumination weighted imaging condition for migrated seismic data | |
| US11287542B2 (en) | Seismic data processing including surface multiple modeling for ocean bottom cable systems | |
| CN117518251A (en) | Free interface multiple prediction method, device, computing equipment and storage medium |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200324 |
|
| RJ01 | Rejection of invention patent application after publication |