WO2022239305A1 - 反射法地震探査による受振データの処理方法 - Google Patents
反射法地震探査による受振データの処理方法 Download PDFInfo
- Publication number
- WO2022239305A1 WO2022239305A1 PCT/JP2022/002213 JP2022002213W WO2022239305A1 WO 2022239305 A1 WO2022239305 A1 WO 2022239305A1 JP 2022002213 W JP2022002213 W JP 2022002213W WO 2022239305 A1 WO2022239305 A1 WO 2022239305A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- wave
- water surface
- water
- pseudo
- reflected wave
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013459 approach Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 266
- 238000012545 processing Methods 0.000 claims description 34
- 230000000644 propagated effect Effects 0.000 claims description 17
- 230000003111 delayed effect Effects 0.000 claims description 10
- 230000001902 propagating effect Effects 0.000 claims description 10
- 238000001028 reflection method Methods 0.000 claims description 7
- 230000002250 progressing effect Effects 0.000 abstract 1
- 230000002123 temporal effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 15
- 238000004364 calculation method Methods 0.000 description 10
- 238000003672 processing method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 4
- 238000013213 extrapolation Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 241000283153 Cetacea Species 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001629 suppression 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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- 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
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3843—Deployment of seismic devices, e.g. of streamers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/614—Synthetically generated data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Definitions
- the present disclosure relates to a method of processing received vibration data from reflection method seismic survey using a hypocenter placed in water.
- Reflection seismic survey is known as one of the underwater exploration methods.
- a seismic source towed by a seismic survey vessel and a streamer cable are thrown into the water.
- the seismic source generates sound waves in a predetermined band, and individual geophones (hydrophones) provided on the streamer cable receive the sound waves output from the seismic source.
- the received sound waves include reflected waves of sound waves from strata under the seabed.
- physical properties such as the shape and density of strata can be specified by measuring the arrival time of reflected waves while intermittently changing the positions of the epicenter and the streamer cable (see Patent Document 1).
- the target of underwater exploration is the stratum beneath the waterbed. Therefore, it is possible to increase the intensity of the reflected wave from the stratum as the epicenter is closer to the bottom of the water.
- the geophone since sound waves propagate in all directions from the epicenter, the geophone receives not only the reflected waves of the sound waves that directly reach the stratum from the epicenter (for convenience, they are called direct waves), but also A reflected wave of the arriving sound wave (for convenience, referred to as a water surface reflected wave) is also received. In this case, the same stratum appears at different depths in the cross-sectional view, which is the received vibration data, making appropriate evaluation of the stratum difficult.
- the present disclosure has been made in view of such circumstances, and an object thereof is to provide a processing method capable of reducing the influence of reflected waves on the water surface for seismic data obtained by reflection method seismic survey. .
- a method is a method of processing received vibration data by reflection method seismic survey, wherein the direct wave of the sound wave indicated in the received vibration data is transmitted at a distance twice the water depth of the epicenter in water. is propagated virtually in the advancing direction of the time axis for only the time it advances, and the amplitude is corrected so as to approach the amplitude of the water surface reflected wave of the sound wave shown in the received vibration data, thereby producing a pseudo water surface reflected wave and subtracting the component corresponding to the pseudo water surface reflected wave from the received vibration data.
- the direct wave is a sound wave that is output from the seismic source and reflected by a reflecting surface without being reflected by the water surface
- the water surface reflected wave is output from the seismic source, reflected by the water surface, and then reflected by the water surface. It is a sound wave delayed by a time corresponding to a distance twice the water depth of the epicenter by being reflected on a surface.
- the water surface reflected wave is virtually propagated in the opposite direction of the time axis for a period of time corresponding to a distance twice the water depth of the seismic source in water, and the phase is changed to
- the method may further include calculating a pseudo direct wave by inverting, and adding data indicating the pseudo direct wave to the received vibration data from which the component corresponding to the pseudo water surface reflected wave is subtracted.
- the method virtually propagates the direct wave of the sound wave indicated in the received vibration data in the direction of travel of the time axis for a time required to travel a distance twice the depth of water of a geophone of the sound wave, and , correcting the amplitude so as to approach the amplitude of the second water surface reflected wave of the sound wave indicated in the received vibration data to calculate the second pseudo water surface reflected wave, and calculating the second pseudo water surface reflected wave of the sound wave indicated in the received vibration data
- the direct wave is propagated virtually in the advancing direction of the time axis only for the time required for the sound wave to travel in water a distance twice the water depth of the epicenter and a distance twice the water depth of the geophone of the sound wave, and A third pseudo water surface reflected wave is calculated by correcting the amplitude so as to approach the amplitude of the third water surface reflected wave of the sound wave indicated in the received vibration data, and the second pseudo water surface reflected wave is obtained from the received vibration data.
- the second water-surface-reflected wave is a sound wave that is output from the seismic source, reflected by the reflecting surface, and further reflected by the water surface, thereby being delayed by a time corresponding to a distance twice as large as the water depth of the geophone.
- the third water surface reflected wave is output from the epicenter, reflected by the water surface, reflected by the reflective surface, and further reflected by the water surface for a time corresponding to a distance twice the water depth of the epicenter. and a sound wave delayed by a time corresponding to a distance twice the water depth of the geophone.
- the first water-surface reflected wave is virtually propagated in a direction opposite to the advancing direction of the time axis for a period of time corresponding to a distance twice the water depth of the seismic source in water;
- the first pseudo direct wave is calculated by inverting the phase, and the second water surface reflected wave is shifted from the traveling direction of the time axis for the time corresponding to the distance twice the water depth of the geophone in the water.
- a second pseudo direct wave is calculated by virtually propagating in the opposite direction and further inverting its phase, and the third water surface reflected wave is propagated in water at a distance twice the water depth of the epicenter.
- a third pseudo direct wave is calculated by virtually propagating in a direction opposite to the traveling direction of the time axis for a time corresponding to the corresponding time and a time corresponding to a distance twice the water depth of the geophone, and the first pseudo direct wave is calculated.
- the data representing the wave, the second pseudo direct wave and the third pseudo direct wave are obtained by subtracting the components corresponding to the first pseudo water surface reflected wave, the second pseudo water surface reflected wave and the third pseudo water surface reflected wave. It may be added to the received vibration data.
- FIG. 1 is a block diagram of an underwater exploration system;
- FIG. It is a flow chart which shows a processing method of received vibration data.
- FIG. 4 is a diagram for explaining an example of a cross-sectional view obtained from received vibration data and its processing;
- FIG. 4 is a diagram for explaining the progress of sound waves reflected by reflecting surfaces at the same position;
- It is a graph which shows the relationship between the operating depth (water depth) of an epicenter and signal strength.
- FIG. 4 is a diagram for explaining changes in the reception time of sound waves (converted depth) due to differences in path length; It is a flow chart which shows a processing method of received vibration data. It is a flow chart which shows a processing method of received vibration data. FIG. 4 is a diagram for explaining an example of a cross-sectional view obtained from received vibration data and its processing; It is a flow chart which shows a processing method of received vibration data.
- FIG. 1 is a schematic configuration diagram showing a bottom-of-water exploration system 10 according to this embodiment.
- FIG. 2 is a block diagram of the underwater exploration system 10 according to this embodiment.
- the sea bottom exploration system 10 includes an epicenter 11, a control device 12, and a plurality of geophones 13, and uses a seismic reflection method to explore the structure under the sea bottom 3. .
- the bottom-of-water exploration system 10 of this embodiment assumes operating on the sea.
- the operating environment of the bottom-of-water exploration system 10 is not limited to the ocean, and can be operated in environments such as rivers and lakes where bottom-of-water exploration can be performed.
- the epicenter 11 is set in water with a depth of D1 and outputs sound waves 5.
- the epicenter 11 is, for example, a low-frequency generator having a well-known configuration, and includes a diaphragm (not shown) driven by hydraulic pressure. Vibration of the diaphragm is controlled by a control device 12 to generate sound waves 5 of a desired frequency.
- the diaphragm may be driven by a piezoelectric element or by any other known drive mechanism.
- the epicenter 11 is not limited to the low-frequency generator described above, and may be a sound wave generator such as an air gun.
- a plurality of geophones 13 are installed underwater and arranged in one direction. Each geophone 13 is a so-called hydrophone and receives sound waves 5 in water. A plurality of geophones 13 are connected to each other or held at intervals by a member 15 such as a cable to form a streamer cable 14 .
- the control device 12 is configured as a so-called computer including a central processing unit (CPU), a storage unit, an auxiliary storage device, etc., and is mounted on the seismic survey ship 20, for example.
- the control device 12 controls the generation of the sound wave 5 by the seismic source 11, the frequency, and the like.
- the sound waves generated may be impulse waves, sweep waves, pseudo-random waves, or other well-known waveforms of sound waves. In either case, the control device 12 records the sound wave 5 received by each geophone 13 as signal data.
- the seismic survey ship 20 is equipped with a lifting device 21.
- the hoisting device 21 launches and hoists a towline 23 and a towline 24 .
- the epicenter 11 is connected to the end of the towline 23 and the streamer cable 14 is connected to the end of the towline 24 . Therefore, depending on the operation of the lifting device 21, the seismic center 11 and the streamer cable 14 are thrown into or recovered from the water.
- the seismic survey ship 20 has, for example, a lifting device 21 that supports a towline 23 and a towline 24 connected to the epicenter 11 and the geophone 13 so as to be able to wind up and unwind.
- the equipment on which the underwater exploration system 10 is mounted is not limited to the seismic survey ship 20, and any equipment that can arrange the epicenter 11 and the streamer cable 14 underwater may be used.
- FIG. 5 is a diagram for explaining an example of a cross-sectional view obtained from received vibration data and its processing.
- FIG. 6 is a diagram for explaining the progress of sound waves reflected by the reflecting surface at the same position.
- the signal data recorded in the control device 12 or the like is subjected to basic processing in the reflection method seismic survey such as cross-correlation processing, NMO correction, and common midpoint superimposition. (step S10).
- basic processing in the reflection method seismic survey such as cross-correlation processing, NMO correction, and common midpoint superimposition.
- step S10 it is assumed that strata A, strata B, and strata C (see FIG. 5(c)) are confirmed under the water bottom as the reflecting surface A1 from the received vibration data obtained by the processing of step S10.
- Reflective surface A1 is the interface between the water and the bottom of the water.
- the reflecting surface of the sound wave 5 corresponding to the interface between the stratum A and the stratum B is denoted by B1
- the reflecting surface of the acoustic wave 5 corresponding to the interface between the stratum B and the stratum C is represented by C1.
- the reflecting surface is the interface between water and the stratum (bottom of the water), or the interface between two strata having different geological features.
- FIG. 5(a) is a cross-sectional view 40A showing the structure under the reflective surface A1 (water bottom) obtained by the process of step S10.
- Cross-sectional view 40A shows a reflective surface A1 that is the bottom of the water and two reflective surfaces B1 and C1 below the bottom of the water, as well as a pseudo reflective surface A2 that is a pseudo-bottom and two pseudo reflective surfaces B2 and C2 that are below the pseudo-bottom. ing.
- the epicenter 11 used to acquire the received vibration data showing the cross-sectional view 40A is placed at the water depth D1 (see Fig. 1).
- the water depth D1 is sufficiently deep with respect to the water surface 2 to such an extent that the influence of waves reflected from the water surface, which will be described later, cannot be ignored.
- the geophone 13 receives reflected waves of sound waves that have directly reached strata A (B, C) from the epicenter 11 (referred to as direct waves 6 for convenience).
- the geophone 13 reaches the stratum A (B, C) after being reflected by the water surface 2 from the epicenter 11 (that is, at twice the water depth D1 of the epicenter 11).
- a reflected wave of the sound wave 5 (delayed by a time corresponding to the distance) (for convenience, referred to as a water surface reflected wave 7) is also received.
- a water surface reflected wave 7 A reflected wave of the sound wave 5 (delayed by a time corresponding to the distance) (for convenience, referred to as a water surface reflected wave 7) is also received.
- the received vibration data includes components that appear at different depths and indicate the same stratum, and it is necessary to remove the influence of this reflected wave from the water surface.
- step S20 the direct wave 6 is virtualized in the traveling direction of the time axis for the time q 1 (see FIG. 5B) during which the sound wave 5 travels in water twice as deep as the water depth D1 of the epicenter 11. , and further corrects the amplitude so that it approaches the amplitude of the water surface reflected wave 7 of the sound wave 5 shown in the received vibration data, thereby calculating the pseudo water surface reflected wave.
- a pseudo water surface reflected wave is obtained, for example, by the following calculation.
- the wave field Q of the received vibration data is propagated downward (that is, in the advancing direction of the time axis).
- Equation (1) is the 15° formula for Wave-Extrapolation by Claerbout (Claerbout J. F., Imaging the Earth's Interior, Blackwell Scientific Publications, 1985, p.88), where v is the speed of sound in water and ⁇ is the angular frequency.
- the water surface reflected wave component is removed from the received vibration data (step S30). Specifically, the component corresponding to the pseudo water surface reflected wave obtained in step S20 is subtracted from the received vibration data.
- the component of the water surface reflected wave 7 is separated from the received vibration data by the following equation (2), and the component of the direct wave 6 is Received vibration data P' of only the component is obtained.
- f' denotes a prediction filter for correcting the amplitude of the pseudo water surface reflected wave.
- FISTA Fast Iterative Shrinkage Thresholding Algorithm
- LASSO Least Absolute Shrinkage and Selection Operator
- n denotes the ln norm
- ⁇ denotes a regularization parameter
- f′ is derived by applying the proximity gradient method to equation (3). Substituting the derived f' into the equation (2) and performing calculation, a cross-sectional view (see FIG. 5(c)) of the received vibration data P' of only the component of the direct wave 6 can be obtained.
- FIG. 5(b) is a cross-sectional view 40B based on the component of the pseudo water surface reflected wave.
- a cross-sectional view showing only the component of the pseudo water surface reflected wave corrected by the prediction filter is obtained.
- the components of the cross-sectional view shown in FIG. 5(b) are subtracted from the components of the cross-sectional view 40A shown in FIG. 40C without the component of the water surface reflected wave 7 is obtained. That is, it is possible to obtain received vibration data with an improved S/N ratio by reducing the influence of the water surface reflected wave 7 .
- the calculation method is not limited to the wave field extrapolation method described above, and a well-known calculation method capable of arithmetically simulating the state of the wave field after t seconds can be applied.
- Well-known calculation methods include, for example, the above-mentioned 45° formula by Claerbout and the wave field extrapolation method by Gazdag (J. Gazdag, Modeling of the acoustic wave equation with transform methods, Geophysics, Vol. 46(6), pp. 854-859), time domain finite difference method (FDTD method) or finite element method.
- the water surface reflected wave 7 shown in the original received seismic data is virtualized in the direction opposite to the traveling direction of the time axis for the time q 1 corresponding to the distance twice the water depth D 1 of the hypocenter 11 in the water. and further inverts its phase to calculate a pseudo direct wave (step S40), and converts the data indicating the pseudo direct wave to received vibration data (that is, step (Step S50).
- steps S40 and S50 instead of treating the water surface reflected wave 7 as noise removed from the received vibration data showing the cross section, it is utilized as a component showing the cross section like the direct wave.
- the received vibration data from the water surface reflected wave 7 can be regarded as the received vibration data from a virtual hypocenter 11 (hereinafter referred to as a mirror image hypocenter 11V) located at a mirror image position of the hypocenter 11 with respect to the water surface 2 (Fig. 6(b)). reference). Therefore, the received vibration data from the mirror image source 11V is added to the received vibration data from the direct wave 6 . In other words, it is possible to acquire the received vibration data equivalent to the seismic survey using the actual hypocenter 11 and the mirror image hypocenter 11V with one hypocenter.
- the received vibration data by the water surface reflected wave 7, that is, the received vibration data by the mirror image epicenter 11V is obtained from the equation (4) by using the original vibration data d and the vibration data P′ of only the direct wave 6 component obtained in step S30. be able to.
- the water surface reflected wave 7 is delayed from the direct wave 6 by time q1 (see FIG. 5(b)). Therefore, the wave field R of the received vibration data due to the water surface reflected wave 7 is propagated upward by time q1 (that is, in the direction opposite to the advancing direction of the time axis). In other words, the water surface reflected wave 7 is corrected due to the difference between the actual survey result resulting from the difference in the propagation paths of the two sound waves and the received vibration data obtained from the mirror image source 11V. Equation (1) was used to propagate the wave field Q downward in step S20.
- the wave field R is propagated upward and a pseudo direct wave is obtained.
- the vibration data indicated by the pseudo direct wave (that is, the corrected water surface reflected wave 7) can be accurately added to the vibration data indicating the direct wave.
- the upward propagation process is not limited to the calculation using the formula (1), and various calculations that can be applied to the downward propagation process can be applied.
- FIG. 7 is a graph showing the relationship between the operational depth (water depth) of the epicenter 11 and the signal intensity.
- the horizontal axis indicates the operational depth of the epicenter 11, and the vertical axis indicates the calculated value of the signal strength of the sound waves emitted from the epicenter at a water depth of 750 m.
- a solid line connecting the white circles indicates the signal intensity of the reflected wave obtained by the direct wave 6 .
- a dotted line connecting black circles indicates the signal intensity of the reflected wave obtained by the water surface reflected wave 7 .
- a solid line connecting the black circles indicates the sum of these signal intensities.
- Fig. 7 shows that the signal strength increases as the operating depth of the hypocenter 11 increases, and the effect of using the water surface reflected wave 7 is about 5 dB of the signal strength when the operating depth of the hypocenter 11 is 50 m. At 225 m, the signal strength increases by about 2 dB.
- the epicenter 11 When using a low-frequency generator equipped with a diaphragm or the like as the epicenter 11, the epicenter 11 outputs a sound wave with a relatively long output time, such as a sweep wave. Such sound waves have less sound pressure energy per unit time than impulse waves obtained from air guns.
- a relatively long output time such as a sweep wave.
- Such sound waves have less sound pressure energy per unit time than impulse waves obtained from air guns.
- the present embodiment by increasing the operational depth of the hypocenter 11, sufficient signal strength can be obtained even with a hypocenter with small sound pressure energy such as a low-frequency generator, and the S/N ratio of the received vibration data can be improved. In this case, since the sound pressure energy is reduced, it is possible to reduce the negative impact on the marine environment, such as the ecological environment of marine mammals, while reducing the energy consumption of the underwater exploration system.
- the phase of the water surface reflected wave 7 is inverted with respect to the direct wave 6 . Therefore, the water surface reflected wave 7 interferes with the direct wave 6 propagating to the bottom of the water 3, and a phenomenon (so-called notch) occurs in which a specific frequency component and its integral multiple frequency component are weakened. This tendency becomes more remarkable as the operating depth of the epicenter 11 becomes deeper.
- the received vibration data is separated into the component due to the water surface reflected wave 7 and the component due to the direct wave 6, and then the phases of both are matched and added. can be recovered.
- FIG. 8 is a diagram showing the frequency spectrum (dotted line) of the original received vibration data and the frequency spectrum (solid line) of the received vibration data obtained by performing the processing according to the present embodiment on the original received vibration data.
- the operating depth of hypocenter 11 at this time is 50 m.
- the arrows in the figure indicate the frequencies at which the above notch effect appears. As shown in this figure, according to this embodiment, the influence of the notches can be removed, and the S/N ratio of the received vibration data can be improved.
- FIG. 9 is a schematic configuration diagram showing the bottom-of-water exploration system 10 according to this embodiment.
- FIG. 10 is a diagram for explaining changes in the reception time of sound waves due to differences in path length.
- 11 and 12 are flow charts showing a method for processing received vibration data.
- the configuration of the bottom-of-water exploration system 10 according to the second embodiment is the same as the configuration of the bottom-of-water exploration system 10 according to the first embodiment.
- the direct wave 6 described in the description of the first embodiment is referred to as a direct wave DW
- the water surface reflected wave 7 described in the description of the first embodiment is referred to as a first water surface reflected wave RW1.
- the geophone 13 receives the direct wave DW (direct wave 6) and the first water surface reflected wave RW1 (water surface reflected wave 7).
- the first water surface reflected wave RW1 is received with a delay of time q1 corresponding to a distance twice as large as the water depth D1 of the epicenter 11 .
- the second water surface reflected wave RW2 and the third water surface reflected wave RW3 reaching the geophone 13 are also used for arithmetic processing described later.
- the second water surface reflected wave RW2 is a sound wave that travels as a direct wave DW and is reflected by the water surface 2
- the third water surface reflected wave RW3 travels as a water surface reflected wave RW1 and is further reflected by the water surface 2. It is a reflected sound wave.
- the path length of the second water surface reflected wave RW2 is longer than the path length of the direct wave DW. Also, the path length of the third water surface reflected wave RW3 is longer than the path length of the second water surface reflected wave RW2. That is, each sound wave is shifted from each other on the vertical round-trip time axis.
- the first water-reflected wave RW1 is the water-surface-reflected wave 7 of the first embodiment, which is output from the epicenter 11, reflected by the water surface 2, and delayed by time q1 from the direct wave DW by being reflected by the reflecting surface. sound wave.
- the second water-surface-reflected wave RW2 is output from the epicenter 11, reflected by the reflecting surface, and further reflected by the water surface 2 , thereby delaying the sound wave by a time q2 corresponding to a distance twice the water depth D2 of the geophone 13. is.
- the second water surface reflected wave RW2 is a sound wave that travels as a direct wave DW and is then reflected by the water surface 2 . Therefore, the second water surface reflected wave RW2 appears with a delay of time q2 with respect to the direct wave DW.
- the third water surface reflected wave RW3 is output from the epicenter 11, reflected by the water surface 2, reflected by the reflecting surface, and further reflected by the water surface 2, resulting in a time q corresponding to a distance twice the water depth D1 of the epicenter 11. 1 and the sound wave delayed by the time q 2 corresponding to the distance twice the water depth D 2 of the geophone 13 .
- the third water surface reflected wave RW3 is a sound wave that travels as the first water surface reflected wave RW1 and is then reflected by the water surface 2.
- the size relationship between the time q1 and the time q2 follows the size relationship between the water depth D1 of the epicenter 11 and the water depth D2 of the geophone 13 .
- the water depth D2 of the geophone 13 is a water depth capable of detecting a total of four types of sound waves: the direct wave DW, the first water surface reflected wave RW1, the second water surface reflected wave RW2, and the third water surface reflected wave RW3.
- the received vibration data includes components that appear at four different depths and indicate the same stratum, and it is necessary to remove the influence of each reflected wave from the water surface.
- step S10 the processing from step S10 to step S30 similar to that of the first embodiment is performed using each water surface reflected wave (see FIG. 3).
- received vibration data is created (step S10), and then arithmetic processing in step S20 is performed.
- step S20 a pseudo water surface reflected wave (hereinafter referred to as first pseudo water surface reflected wave) is calculated from the direct wave DW (direct wave 6) and the first water surface reflected wave RW1 (water surface reflected wave 7) (step S21).
- the direct wave DW and the second water surface reflected wave RW2 are used to calculate the second pseudo water surface reflected wave (step S22), and the direct wave DW and the third water surface reflected wave RW3 are used to calculate the third pseudo water surface reflected wave. (step S23).
- the direct wave DW is virtualized in the direction of travel of the time axis for the time q 1 during which the sound wave 5 travels in water twice as deep as the water depth D1 of the epicenter 11.
- the first pseudo water surface reflected wave is calculated by propagating the wave and correcting the amplitude so as to approach the amplitude of the first water surface reflected wave RW1 of the sound wave 5 shown in the received vibration data.
- the direct wave DW is virtualized in the traveling direction of the time axis for the time q2 required for the sound wave 5 to travel a distance twice the water depth D2 of the geophone 13 in the water.
- the direct wave DW travels a distance twice the water depth D1 of the seismic source 11 and a distance twice the water depth D2 of the geophone 13 in water.
- step S30 components corresponding to the first, second, and third pseudo water surface reflected waves are removed from the received vibration data (steps S31 to S33).
- the arithmetic processing in step S30 components corresponding to the first, second, and third pseudo water surface reflected waves are removed from the received vibration data (steps S31 to S33).
- the convolution matrix D′ using the first pseudo water surface reflected wave the first to third water surface reflected waves RW1, RW2, RW1, RW2, Each component of RW3 is separated, and received vibration data P' of only the component of the direct wave DW is obtained.
- FIG. 13 is a diagram for explaining an example of a cross-sectional view obtained from received vibration data and its processing.
- FIG. 13(a) is a cross-sectional view 50A showing the structure under the reflective surface A1 (water bottom) obtained by the process of step S10.
- a cross-sectional view 50A shows a first pseudo-reflecting surface A2, a second pseudo-reflecting surface A3, and a third pseudo-reflecting surface A4 that are pseudo-bottoms in addition to the reflecting surface A1 that is the bottom of the water.
- the cross-sectional view 50A shows a first pseudo-reflecting surface B2, a second pseudo-reflecting surface B3, and a third pseudo-reflecting surface B4 under the pseudo-bottom, in addition to the two reflecting surfaces B1 below the bottom. Further, the cross-sectional view 50A shows the first pseudo-reflecting surface C2, the second pseudo-reflecting surface C3 and the third pseudo-reflecting surface C4 in addition to the two reflecting surfaces C1 under the water bottom.
- Reflective surfaces A1, B1, C1 are indicated by direct waves DW.
- the first pseudo reflective surfaces A2, B2, C2 are indicated by the first water reflected wave RW1.
- the second pseudo reflective surfaces A3, B3, C3 are indicated by the second water surface reflected wave RW2.
- the third pseudo reflecting surface A4, B4, C4 is indicated by the third water surface reflected wave RW3.
- each component of the first to third pseudo water surface reflected waves corrected by the prediction filter is subtracted from the component of the sectional view 50A shown in FIG. 13(a).
- a cross-sectional view 50B without each component of the first to third water surface reflected waves RW1, RW2, and RW3 shown in FIG. 13(b) is obtained. That is, even if the geophone 13 is placed at the water depth D2 where the effects of the second and third water surface reflected waves RW2 and RW3 appear, the effects of the first to third water surface reflected waves RW1, RW2 and RW3 are reduced, and the S/ It is possible to obtain received vibration data with an improved N ratio.
- the geophone 13 (streamer cable 14) can be placed at a depth sufficiently distant from the water surface 2.
- a general streamer cable has a total length of several hundred meters to several kilometers. When such a long streamer cable is towed near the surface of the water, navigation of ships in the water area around the streamer cable is restricted, and there is concern that the streamer cable itself may be damaged.
- the streamer cable 14 can be arranged in sufficiently deep water. Therefore, estimation exploration can be performed while avoiding interference with ships.
- FIG. 14 is a flowchart showing a method for processing received vibration data.
- the processing of steps S41 to S43 is performed as the arithmetic processing of step S40, and the processing of step S51 is performed as the arithmetic processing of step S50.
- the first water-reflected wave RW1 is virtually propagated in the direction opposite to the traveling direction of the time axis for the time q1 corresponding to the distance twice the water depth D1 of the epicenter 11 in the water, and further, By inverting the phase, the first pseudo direct wave is calculated (step S41). This is the same as the arithmetic processing for the water surface reflected wave 7 in the first embodiment.
- the second water surface reflected wave RW2 is virtually propagated in the opposite direction of the time axis for the time q2 corresponding to the distance twice the water depth D2 of the geophone 13 in the water, Furthermore, the second pseudo direct wave is calculated by inverting the phase (step S42). That is, in step S40 of the first embodiment, the water surface reflected wave 7 is replaced with the second water surface reflected wave RW2, and time q1 is replaced with time q2.
- the third water-surface reflected wave RW3 is applied to the sound wave 5 for a time q1 corresponding to a distance twice the water depth D1 of the epicenter 11 in water and a time q2 corresponding to a distance twice the water depth D2 of the geophone 13.
- the third pseudo direct wave is calculated by virtually propagating in the direction opposite to the advancing direction of the time axis (step S43). That is, in step S40 of the first embodiment, the water surface reflected wave 7 is replaced with the third water surface reflected wave RW3, and time q1 is replaced with the sum of time q1 and time q2.
- the third water surface reflected wave RW3 has the same phase as the direct wave DW, so the inversion processing is not performed.
- the data representing the first to third pseudo direct waves, respectively is received vibration data from which the components corresponding to the first to third pseudo water surface reflected waves are subtracted (that is, the processing of steps S30 (S31 to S33) is performed. received vibration data) (step S51).
- the same effect as in the first embodiment can be obtained. That is, by increasing the operating depth of the epicenter 11, even if such sound pressure energy is small, sufficient signal strength can be obtained in seismic surveys. Moreover, the influence of the notch mentioned above can be reduced. Thereby, the S/N ratio of the received vibration data can be improved.
- suppression of sound pressure energy in underwater exploration can reduce adverse effects on marine ecosystems. conserve and sustainably use the oceans and marine resources for sustainable development.”
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
以下、本開示の第1実施形態について説明する。図1は、本実施形態に係る水底探査システム10を示す概略構成図である。図2は、本実施形態に係る水底探査システム10のブロック図である。
水面反射波7による受振データ、つまり鏡像震源11Vによる受振データは、元の受振データdとステップS30で求めた直達波6の成分のみの受振データP´を用いることで、(4)式から求めることができる。
次に本開示の第2実施形態について説明する。なお、本実施形態の説明において、第1実施形態と重複する事項については同一の符号を付し、その説明を省略する。図9は、本実施形態に係る水底探査システム10を示す概略構成図である。図10は、経路長の違いによる音波の受振時間の変化を説明するための図である。図11及び図12は、受振データの処理方法を示すフローチャートである。
Claims (4)
- 水中に音波を出力する震源と前記音波を受振する受振器を用いた反射法地震探査による受振データの処理方法であって、
前記受振データに示された前記音波の直達波を、前記音波が水中において震源の水深の2倍の距離を進行する時間だけ、時間軸の進行方向に仮想的に伝播させ、更に、その振幅を前記受振データに示された前記音波の第1水面反射波の振幅に近付けるように補正することによって第1疑似水面反射波を算出し、
前記受振データから前記第1疑似水面反射波に相当する成分を減じる
ことを含み、
前記直達波は、前記震源から出力され、水面での反射を経ることなく反射面で反射した音波であり、
前記第1水面反射波は、前記震源から出力され、水面で反射し更に前記反射面で反射することによって、前記震源の水深の2倍の距離に相当する時間だけ遅延した音波である、
受振データの処理方法。 - 前記第1水面反射波を、前記音波が水中において前記震源の水深の2倍の距離に相当する時間だけ、時間軸の進行方向と逆方向に仮想的に伝播させ、更に、その位相を反転させることによって第1疑似直達波を算出し、
前記第1疑似直達波を示すデータを、前記第1疑似水面反射波に相当する成分が減じられた前記受振データに加算する
ことを更に含む請求項1に記載の受振データの処理方法。 - 前記受振データに示された前記音波の前記直達波を、前記音波の受振器の水深の2倍の距離を進行する時間だけ、時間軸の進行方向に仮想的に伝播させ、更に、その振幅を前記受振データに示された前記音波の第2水面反射波の振幅に近付けるように補正することによって第2疑似水面反射波を算出し、
前記受振データに示された前記音波の前記直達波を、前記音波が水中において震源の水深の2倍の距離と前記音波の受振器の水深の2倍の距離を進行する時間だけ、時間軸の進行方向に仮想的に伝播させ、更に、その振幅を前記受振データに示された前記音波の第3水面反射波の振幅に近付けるように補正することによって第3疑似水面反射波を算出し、
前記受振データから前記第2疑似水面反射波に相当する成分と前記第3疑似水面反射波に相当する成分とを減じる
ことを更に含み、
前記第2水面反射波は、前記震源から出力され、前記反射面で反射し、更に水面で反射することによって、前記受振器の水深の2倍の距離に相当する時間だけ遅延した音波であり、
前記第3水面反射波は、前記震源から出力され、水面で反射し、前記反射面で反射し、更に水面で反射することによって、前記震源の水深の2倍の距離に相当する時間と前記受振器の水深の2倍の距離に相当する時間だけ遅延した音波である、
請求項1に記載の受振データの処理方法。 - 前記第1水面反射波を、前記音波が水中において前記震源の水深の2倍の距離に相当する時間だけ、時間軸の進行方向と逆方向に仮想的に伝播させ、更に、その位相を反転させることによって第1疑似直達波を算出し、
前記第2水面反射波を、前記音波が水中において前記受振器の水深の2倍の距離に相当する時間だけ、時間軸の進行方向と逆方向に仮想的に伝播させ、更に、その位相を反転させることによって第2疑似直達波を算出し、
前記第3水面反射波を、前記音波が水中において前記震源の水深の2倍の距離に相当する時間と前記受振器の水深の2倍の距離に相当する時間だけ、時間軸の進行方向と逆方向に仮想的に伝播させることによって第3疑似直達波を算出し、
前記第1疑似直達波、前記第2疑似直達波及び前記第3疑似直達波を示すデータを、前記第1疑似水面反射波、第2疑似水面反射波及び第3疑似水面反射波に相当する成分が減じられた前記受振データに加算する
ことを更に含む請求項3に記載の受振データの処理方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2023520768A JPWO2022239305A1 (ja) | 2021-05-12 | 2022-01-21 | |
EP22807012.4A EP4339651A1 (en) | 2021-05-12 | 2022-01-21 | Method for using reflection seismic survey to process received data |
US18/502,427 US20240069231A1 (en) | 2021-05-12 | 2023-11-06 | Method of processing seismic data acquired by reflection seismic survey |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021081267 | 2021-05-12 | ||
JP2021-081267 | 2021-05-12 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US18/502,427 Continuation US20240069231A1 (en) | 2021-05-12 | 2023-11-06 | Method of processing seismic data acquired by reflection seismic survey |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022239305A1 true WO2022239305A1 (ja) | 2022-11-17 |
Family
ID=84028131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2022/002213 WO2022239305A1 (ja) | 2021-05-12 | 2022-01-21 | 反射法地震探査による受振データの処理方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20240069231A1 (ja) |
EP (1) | EP4339651A1 (ja) |
JP (1) | JPWO2022239305A1 (ja) |
WO (1) | WO2022239305A1 (ja) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7616523B1 (en) * | 2008-10-22 | 2009-11-10 | Pgs Geophysical As | Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth |
CN101876715A (zh) * | 2009-10-30 | 2010-11-03 | 中国科学院声学研究所 | 一种拖曳声学阵列的拖船噪声抑制方法 |
US20110058450A1 (en) * | 2005-05-05 | 2011-03-10 | Van Den Berg Peter M | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
US20110166790A1 (en) * | 2004-10-08 | 2011-07-07 | Cggveritas Services Sa | Seismic Processing for the Elimination of Multiple Reflections |
US20130201791A1 (en) * | 2012-02-07 | 2013-08-08 | Gregory Ernest Parkes | Method and system for determining source signatures after source ghost removal |
JP2014137320A (ja) | 2013-01-18 | 2014-07-28 | Ihi Corp | 水域地中探査システム及び水域地中探査方法 |
US20150198730A1 (en) * | 2014-01-10 | 2015-07-16 | Cgg Services Sa | Seismic image acquisition using ghost compensation in kirchhoff migration |
US20160054465A1 (en) * | 2013-06-10 | 2016-02-25 | Downunder Geosolutions Pty Ltd | Seismic data spectrum restoring and broadening |
CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
-
2022
- 2022-01-21 EP EP22807012.4A patent/EP4339651A1/en active Pending
- 2022-01-21 JP JP2023520768A patent/JPWO2022239305A1/ja active Pending
- 2022-01-21 WO PCT/JP2022/002213 patent/WO2022239305A1/ja active Application Filing
-
2023
- 2023-11-06 US US18/502,427 patent/US20240069231A1/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110166790A1 (en) * | 2004-10-08 | 2011-07-07 | Cggveritas Services Sa | Seismic Processing for the Elimination of Multiple Reflections |
US20110058450A1 (en) * | 2005-05-05 | 2011-03-10 | Van Den Berg Peter M | Method for deghosting and water layer multiple reflection attenuation in marine seismic data |
US7616523B1 (en) * | 2008-10-22 | 2009-11-10 | Pgs Geophysical As | Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth |
CN101876715A (zh) * | 2009-10-30 | 2010-11-03 | 中国科学院声学研究所 | 一种拖曳声学阵列的拖船噪声抑制方法 |
US20130201791A1 (en) * | 2012-02-07 | 2013-08-08 | Gregory Ernest Parkes | Method and system for determining source signatures after source ghost removal |
JP2014137320A (ja) | 2013-01-18 | 2014-07-28 | Ihi Corp | 水域地中探査システム及び水域地中探査方法 |
US20160054465A1 (en) * | 2013-06-10 | 2016-02-25 | Downunder Geosolutions Pty Ltd | Seismic data spectrum restoring and broadening |
US20150198730A1 (en) * | 2014-01-10 | 2015-07-16 | Cgg Services Sa | Seismic image acquisition using ghost compensation in kirchhoff migration |
CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
Non-Patent Citations (2)
Title |
---|
CLAERBOUT J.F.: "Imaging the Earth's Interior", 1985, BLACKWELL SCIENTIFIC PUBLICATIONS, pages: 88 |
J. GAZDAG: "Modeling of the acoustic wave equation with transform methods", GEOPHYSICS, vol. 46, no. 6, pages 854 - 859 |
Also Published As
Publication number | Publication date |
---|---|
EP4339651A1 (en) | 2024-03-20 |
US20240069231A1 (en) | 2024-02-29 |
JPWO2022239305A1 (ja) | 2022-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2575274C (en) | System for attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors | |
CA2698020C (en) | Method for operating marine seismic vibrator array to enhance low frequency output | |
US7835225B2 (en) | Method for attenuating particle motion sensor noise in dual sensor towed marine seismic streamers | |
RU2507543C2 (ru) | Совместные интерполяция и подавление волн-спутников в сейсмических данных | |
CN103823235A (zh) | 用于连续数据采集的设备和方法 | |
CN103823243A (zh) | 用于分离在连续数据采集地震勘测期间记录的数据的方法 | |
US8833509B2 (en) | Two component source seismic acquisition and source de-ghosting | |
CA2897141C (en) | Subaqueous underground survey system and subaqueous underground survey method | |
CN1664617A (zh) | 组合海洋地震拖缆中的压力传感器和质点运动传感器的信号的系统 | |
RU2678261C2 (ru) | Сейсмический датчик с датчиками перемещения снижения шума | |
US10996359B2 (en) | Removal of acquisition effects from marine seismic data | |
CA2856335A1 (en) | Separation of simultaneous source data | |
US9453928B2 (en) | Methods and computing systems for processing data | |
CA2838711A1 (en) | A system and method for the removal of shallow water multiples using a hybrid multi-channel prediction method | |
MX2011008294A (es) | Formacion de imagenes con mediciones vectoriales. | |
EA024149B1 (ru) | Способ подавления интерференционных помех в сейсмических данных с помощью системы записи с двумя типами датчиков | |
AU2017378995A1 (en) | System and method for reconstructed wavefield imaging | |
CA3111778A1 (en) | Adaptive receiver deghosting for seismic streamer | |
EP3345020B1 (en) | Method for attenuating multiples | |
GB2481840A (en) | Positioning seismic sources beneath an acoustic ceiling to reduce source ghosts and bubble pulses | |
US10451759B2 (en) | Control system for a marine vibrator | |
US9671513B2 (en) | Seismic image acquisition using ghost compensation in Kirchhoff migration | |
WO2022239305A1 (ja) | 反射法地震探査による受振データの処理方法 | |
CN112424645A (zh) | 利用重合滤波通过回波去混合对地震数据去重影 | |
US20190391290A1 (en) | Dipole Source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22807012 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2023520768 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2022807012 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2022807012 Country of ref document: EP Effective date: 20231212 |