NL2020152B1 - Kirchhoff Beam Migration Method Under Complex Topography - Google Patents

Kirchhoff Beam Migration Method Under Complex Topography Download PDF

Info

Publication number
NL2020152B1
NL2020152B1 NL2020152A NL2020152A NL2020152B1 NL 2020152 B1 NL2020152 B1 NL 2020152B1 NL 2020152 A NL2020152 A NL 2020152A NL 2020152 A NL2020152 A NL 2020152A NL 2020152 B1 NL2020152 B1 NL 2020152B1
Authority
NL
Netherlands
Prior art keywords
migration
rays
shot
ray
information
Prior art date
Application number
NL2020152A
Other languages
Dutch (nl)
Other versions
NL2020152A (en
Inventor
Sun Hui
Original Assignee
Univ Southwest Jiaotong
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Univ Southwest Jiaotong filed Critical Univ Southwest Jiaotong
Publication of NL2020152A publication Critical patent/NL2020152A/en
Application granted granted Critical
Publication of NL2020152B1 publication Critical patent/NL2020152B1/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/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/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • 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/51—Migration
    • 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/67—Wave propagation modeling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a Kirchhoff beam migration method under complex topography. The Kirchhoff beam migration method comprises the following steps: inputting relevant parameter files, migration velocity model files and seismic records; dividing the seismic records of each shot into a plurality of different time domain data volumes with a window center as a core, and then decomposing the data into plane waves; tracing rays from shot points along different directions, and calculating the corresponding information of the grid nodes within the scope of a sectional radial beam for storing each ray; tracing the rays from the window center along different directions, and calculating the corresponding information of grid nodes within the scope of the sectional radial beam for storing each ray; selecting the beam from the shot points and the window center for imaging calculation; and accumulating the imaging results of all of the pairs of the beams so as to obtain the migration imaging results. According to the Kirchhoff beam migration method, the coverage of beams in shallow model is increased, and the imaging quality of the Kirchhoff beam migration method in a shallow structure of complex topography model.

Description

Technical field
The invention relates to a seismic migration imaging method, in particular to a Kirchhoff beam migration method under complex topography.
Background Art
Land seismic exploration is usually performed under complex topography (such as mountains, hills, gobi and loess tableland), where problems that the topography causes problems that surface elevation is higher, the seismic data coverage is insufficient, and signal noise ratio (SNR) of collected seismic data is low, often occur. The problems can bring certain difficulty to the seismic data imaging, and affect the computation efficiency and the imaging precision of the seismic imaging.
Chinese Journal of Geophysics ( Stage 4, 2012-) discloses amplitude-preserved Gaussian beam migration under complex topography written by Yue Yubo, et al., and introduces the amplitude-preserved Gaussian beam migration method under complex topography. By taking account of the elevation, dip angle, and actual trace intervals of topography into a local plane wave decomposition step, the migration method is improved for the topography. Simple layered topographic model and Canadian Foothills model are used fortesting Gaussian beam migration under complex topography and Gaussian beam migration achieves superior images.
The 2017 Ph.D. thesis of Jilin University discloses Kirchhoff Dynamic Focused Beam Migration, introduces Kirchhoff dynamic focused beam migration method under complex topography, and introduces the dynamic focused beam propagator into the Kirchhoff beam migration under complex topography to control the width of the beam. The method is applied to Canadian Foothills model, and a superior image to the original method is obtained.CN102590857A discloses a two-way wave prestack depth migration under true topography which firstly obtains the real elevation of seismic data to redefine the velocity model, performs the forward calculation under the true topography conditions, then regulates the seismic data for wave field extrapolation, uses topography boundary conditions, and finally uses the relevant conditions for imaging so as to overcome the impact of topography, and obtain high-quality imaging results under the complex topographic conditions.
It can be seen from the examples that a conventional imaging method can improve the migration imaging results under complex topography to a certain extent, but the realization process is complex, and the calculation efficiency is also low.
Summary of Invention
The invention aims to solve the technical problem of providing a Kirchhoff beam migration method under complex topography. Through adoption of a new beam in a sectional manner and using a cosine square window function in local slant stacking, it not only improves the imaging effect of seismic migration under complex topography, but also improves the processing efficiency of seismic data.
In order to solve the technical problem, the invention adopts the technical scheme:
the Kirchhoff beam migration method under complex topography comprises the following steps: step A: inputting relevant parameter files, migration velocity model and seismic records, wherein, the parameter files comprise grid points, grid spacing, reference frequency, maximum frequency, initial beamwidth, seismic record sampling points, traces of each shot and the spacing of the traces;
step B: dividing the seismic records of each shot into a plurality of different time domain data volumes with a window center as a core according to a distance d =----—----, and then
decomposing the data into plane waves, wherein, V is the average value of the migration velocity filed, is the reference frequency, and /lliax is the maximum frequency;
step C: tracing rays from shot points indifferent directions, and calculating the corresponding information of the grid nodes within the scope of a sectional radial beam for storing each ray, wherein the information comprises traveltime and amplitude; calculating rays from the shot points indifferent dxt 2 ~T=V SPi uT directions, and solving the kinematical ray tracing equation j by a dpi dr
Runge-Kutta method to get the ray information, wherein, x, represents space coordinates, p, represents slowness, τ represents seismic wave traveltime, and v represents the velocity value at discrete points; after the information at the discrete points on the ray is obtained, acquiring the information of the grid nodes within the beam by paraxial approximation, wherein the information includes traveltime and amplitude, calculating the width of the beam by a sectional method:
2jrF a\.g! fam, 2ftVaVgl *0 2\ua. 2πνα//ηα <2Δα° ' f\ 'Λ Fo wherein, Δα is the angle difference of the adjacent rays, and σ is the integral of the velocity along the ray path;
step D: tracing the rays from the window center indifferent directions, and calculating the corresponding information of grid nodes within the scope of the sectional radial beam for storing each ray, wherein the information comprises traveltime and amplitude; similarly, calculating the width of the beam by the sectional method in Step C;
step E: selecting the beam from the shot points and the window center for imaging calculation,
wherein, /,(+) represents the imaging value at the point χ, ps represents the slowness value of rays traced from the shot point, pl)C represents the slowness values of rays traced from the window center, A represents a weight function, and Ds represents a local plane wave decomposition result; and step F: accumulating the imaging results of all of the pairs of the beams so as to obtain the final migration imaging results.
Further, in the step C, the emitting angle range of the rays is from -60 degrees to +60 degrees, and the emitting angle interval between every two adjacent rays is Ap =-----
Compared with conventional technologies, the Kirchhoff beam migration method has the following beneficial effects that through adoption of the sectional radial beam propagator, the coverage of the beam in the shallow layer is enlarged, and the regularity of the migration result in the topographic surface is stronger, and the invention can more clearly reflect the fault structure.
Brief Description of Drawings
Figure 1 is a flow chart of the Kirchhoff beam migration method under complex topography.
Figure 2 is the width diagram of the original Kirchhoff beam migration beam.
Figure 3 is the width diagram of Kirchhoff beam migration beam of the invention.
Figure 4 shows Marmousi model under topography, wherein x represents the horizontal distance and z represents the depth.
Figure 5 shows the original Kirchhoff beam migration result of the Marmousi model under topography.
Figure 6 shows the novel Kirchhoff beam migration result of the Marmousi model under topography.

Claims (2)

EmbodimentsEmbodiments The invention will now be further described in detains in connection with the accompanying drawings and embodiments.The invention will now be further described in the connection with the accompanying drawings and figures. Figure 1 is a flow chart of a Kirchhoff beam migration method under complex topography, and shows an implementation process of the method, wherein the method comprises the following specific steps ofFigure 1 is a flow chart of a Kirchhoff beam migration method under complex topography, and shows an implementation process of the method, in which the method comprises the following specific steps of A) inputting relevant parameter files, migration velocity model files and seismic records, wherein, the parameter files comprise grid points, grid spacing, reference frequency, maximum frequency, initial beamwidth, seismic record sampling points, traces of each shot and the spacing of the traces;A) inserting relevant parameter files, migration velocity model files and seismic records, where, the parameter files include grid points, grid spacing, reference frequency, maximum frequency, initial beamwidth, seismic record sampling points, traces of each shot and the spacing of the traces; B) dividing the seismic records of each shot into a plurality of different time domain data volumes with a window center as a core according to a distance d = . —, and then decomposing the y/fmin./max data into plane waves, wherein, K is the average value of the migration velocity filed, is the reference frequency, and /niax is the maximum frequency;B) dividing the seismic records of each shot into a multiple or different time domain data volumes with a window center as a core according to a distance d =. -, and then decomposing the y / fmin / max data into plane waves, where, K is the average value of the migration velocity filed, is the reference frequency, and / niax is the maximum frequency; C) emitting rays indifferent directions from the shot point, wherein the emitting angle range of the rays is from -60 degrees to +60 degrees, and the angle interval between every two adjacent rays is Δ/? = —1— Ahl , anci solving the kinematical ray tracing equation by a Runge-Kutta method to getC) emitting rays indifferent directions from the shot point, the emitting angle range of the rays is from -60 degrees to +60 degrees, and the angle interval between every two adjacent rays is Δ /? = —1— Ahl, anc i solving the kinematical ray tracing equation by a Runge-Kutta method to get 2Vavg ,/raax the center ray information, as shown below:2V avg , / raax the center ray information, as shown below: GXt wherein, represents space coordinates, p, represents slowness, τ represents seismic wave traveltime, and v represents the velocity value at discrete points; after the information at the discrete points on the ray is obtained, acquiring the information of the grid nodes of the beam by paraxial approximation, wherein the information includes traveltime and amplitude;GX t , represents space coordinates, p, represents slowness, τ represents seismic wave traveltime, and v represents the velocity value at discrete points; after the information at the discrete points on the ray is obtained, acquiring the information from the grid nodes of the beam by paraxial approximation, the information includes traveltime and amplitude; an original calculation formula (as shown in figure 2) for beam width is Η' = 2Δα©, C wherein, Δα is the angle difference of adjacent rays, Ρθ is the velocity value at the starting position of the rays, σ = | Vds is the integral of the velocity along a ray path, but such determination method foran original calculation formula (as shown in figure 2) for beam width is Η '= 2Δα ©, C expect, Δα is the angle difference of adjacent rays, Ρθ is the velocity value at the starting position of the rays, σ = | Vds is the integral of the velocity along a ray path, but such a determination method for J ray the beamwidth is difficult to meet the beam covering demands for the shallow part of the model; calculating the width of the beam by a sectional method (as shown in figure 3):J ray the beam width is difficult to meet the beam covering demands for the shallow part of the model; calculating the width of the beam by a sectional method (as shown in figure 3): j 2nVaVg2πναν8//αΐί>2ΑαΤ W= 2^/®Γ<2Δα^. : *0 *0j 2nV aV g2πν αν8 // αΐί > 2ΑαΤ W = 2 ^ / ® Γ <2Δα ^. : * 0 * 0 D) tracing the rays from the window center indifferent directions, and calculating the corresponding information of grid nodes within the scope of the sectional radial beam for storing each ray, wherein the information comprises traveltime and amplitude; similarly, calculating the width of the beam by the sectional method in C);D) tracing the rays from the window center indifferent directions, and calculating the corresponding information of grid nodes within the scope of the sectional radial beam for malfunction each ray, the information comprises traveltime and amplitude; similarly, calculating the width of the beam by the sectional method in C); E) selecting a beam from the shot point and the window center, and performing imaging calculation according to the formula:E) selecting a beam from the shot point and the window center, and performing imaging calculation according to the formula: Is (x) = Σ ƒ dP° ƒ dPbcA ld (L P = P',t = t') wherein, I,(x) represents the imaging value at the point x, ps represents the slowness value of rays traced from the shot point, phc represents the slowness values of rays traced from the window center, A represents a weight function, and Ds represents a local plane wave decomposition result; andIs (x) = Σ ƒ d P ° ƒ d Pbc A ld (LP = P ', t = t') wherein, I, (x) represents the imaging value at the point x, p s represents the slowness value of rays traced from the shot point, p hc represents the slowness values of rays traced from the window center, a represents a weight function, and D s represents a local plane wave decomposition result; and F) accumulating the imaging results of all of the pairs of the beams so as to obtain the final migration imaging results.F) accumulating the imaging results of all of the pairs of the beams so as to obtain the final migration imaging results. Now, we will verify the scheme and beneficial effects of the invention through the Marmousi model under complex topography. Figure 4 shows the Marmousi velocity model with 737 grid nodes in the horizontal direction and a horizontal grid spacing of 12.5m, and with 750 grid nodes in the vertical direction and a vertical grid spacing of 4m. The seismic records consist of 240 shots with a 90 m spacing and 101 receiving traces per shot with a 25 m spacing. The offset is from 0 m to 2,500 m with 800 sampling points per trace with a sampling interval of 4 ms. Figure 5 shows the Kirchhoff (} beam migration result with a conventional calculation formula of beam width, and figure 6 shows the migration result corresponding to the sectional radial beam width calculation formula of the invention. As seen from the migration result graph, the migration result of the method of the invention has a stronger regularity of the migration energy in the topographic surface, a higher SNR, a clearer 5 geological structure, and the fault structure with white circle can be clearly identified.Now, we will verify the scheme and beneficial effects of the invention through the Marmousi model under complex topography. Figure 4 shows the Marmousi velocity model with 737 grid nodes in the horizontal direction and a horizontal grid spacing of 12.5m, and with 750 grid nodes in the vertical direction and a vertical grid spacing of 4m. The seismic records consist of 240 shots with a 90 m spacing and 101 receiving traces per shot with a 25 m spacing. The offset is from 0 m to 2,500 m with 800 sampling points per trace with a sampling interval or 4 ms. Figure 5 shows the Kirchhoff (} beam migration result with a conventional calculation formula or beam width, and figure 6 shows the migration result corresponding to the sectional radial beam width calculation formula of the invention. As seen from the migration result graph, the migration result or the method of the invention has a stronger regularity of the migration energy in the topographic surface, a higher SNR, a clearer 5 geological structure, and the fault structure with white circle can be clearly identified. The method of the invention is an important seismic data imaging method for complex topography. In view of the coverage shortcoming of conventional beam width calculation formula in the shallow model and the difficulty in meeting the requirements of special conditions under complex topography, the introduction of a sectional beam width calculation method into the Kirchhoff beam 10 migration method under complex topography enhances the imaging effect of the migration method in the shallow part of the model.The method of the invention is an important seismic data imaging method for complex topography. In view of the coverage shortcoming or conventional beam width calculation formula in the shallow model and the difficulty in meeting the requirements of special conditions under complex topography, the introduction of a sectional beam width calculation method into the Kirchhoff beam 10 migration method under complex topography enhances the imaging effect of the migration method in the shallow part of the model. Conclusies:Conclusions: 1. Kirchhoffs stralenbundelmigratiewerkwijze met complexe topografie, met het kenmerk, dat zij de volgende stappen omvat:1. Kirchhoff's beam-beam migration method with complex topography, characterized in that it comprises the following steps: stap A : invoering van relevante parameterbestanden, van een migratiesnelheidsmodel en van seismische registraties, waarbij de parameterbestanden rasterpunten, rasterafstanden, een referentiefrequentie, een maximumfrequentie, een beginbreedte van de stralenbundel, bemonsteringspunten van de seismische registraties, traces van elk schot en de afstand tussen de schoten omvatten:step A: introduction of relevant parameter files, of a migration speed model and of seismic recordings, wherein the parameter files are raster points, raster distances, a reference frequency, a maximum frequency, a starting width of the beam, sampling points of the seismic recordings, traces of each shot and the distance between the beams shots include: stap B : verdeling van de seismische registraties van elk schot in een veelheid van gegevensvolumina van verscheidene tijdsdomeinen met een venstercentrum als kern volgens een aftstand d = . , K----, en daarna min fmax.step B: division of the seismic records of each shot into a plurality of data volumes of various time domains with a window center as a core according to a distance d =. , K ----, and then min fmax. ontbinding van de gegevens in vlakke golven, waarbij Vilve de gemiddelde waarde van de ingediende migratiesnelheid is, /nün de referentiefrequentie is, en fmax de maximumfrequentie is;dissolution of the data in planar waves, where V ilve is the average value of the submitted migration speed, / nün is the reference frequency , and f max is the maximum frequency ; stap C : afbeelding van stralen vanuit schotpunten in verschillende richtingen, en berekening van de overeenstemmende informaties van de rasterknooppunten binnen de omvang van een sectiebundel om elke straal op te slaan, waarbij de informaties de voortplantingstijd en de amplitude omvatten; berekening van stralen vanuit de schotpunten in verschillende richtingen, en oplossing van de kinematische straalafbeeldingsvergelijking bij middel van een Runge-Kutta-werkwijze ten einde de straalinformaties te verkrijgen, waarbij, X, de ruimtecoördinaten voorstelt, pt de traagheid voorstelt, T de voortplantingssnelheid van de seismische golf voorstelt, en V de waaide van de snelheid op discrete punten voorstelt: na het bekomen van de informaties op de discrete punten van de straal, verwerving van de informaties van de rasterknooppunten binnen de stralenbundel door middel van een paraxiale benadering, waarbij de informaties de voortplantingstijd en de amplitude omvatten, berekening van de breedte van de stralenbundel bij middel van een sectiewerkwijze:step C: mapping rays from shot points in different directions, and calculating the corresponding information from the grid nodes within the size of a section beam to store each ray, the information including the propagation time and the amplitude; calculation of rays from the shot points in different directions, and resolution of the kinematic ray image comparison by means of a Runge-Kutta method in order to obtain the beam information, where, X, represents the space coordinates, p t represents the slowness, T the seismic wave, and V represents the wind of the velocity at discrete points: after obtaining the information at the discrete points of the beam, acquiring the information of the grid nodes within the beam by means of a paraxial approach, the information includes the propagation time and the amplitude, calculation of the width of the beam by means of a section method: waarbij, Δα het hoekverschil is van aanliggende stralen, en σ de integraal van de snelheid langsheen het straalpad is;where, Δα is the angular difference of adjacent beams, and σ is the integral of the velocity along the beam path; stap D : afbeelding van de stralen vanuit het venstercentrum in verscheidene richtingen, en berekening van de overeenstemmende informaties van de rasterknooppunten binnen de omvang van de sectiebundel om elke straal op te slaan, waarbij de informaties de vooitplantingstijd en de amplitude omvatten; op een gelijkaardige wijze, berekening van de breedte van de stralenbundel door middel van de sectiewerkwijze in stap C;step D: mapping the rays from the window center in various directions, and calculating the corresponding information from the frame nodes within the size of the section beam to store each ray, the information including the pre-planting time and the amplitude; in a similar manner, calculating the beam width by means of the section method in step C; stap E: selectie van de stralenbundel vanuit de schotpunten en het venstercentrum om de beeldvorming te berekenen.step E: selection of the beam from the shot points and the window center to calculate the image. L (χ) = Σ ƒ dP* ƒ dPbeA Ds (L, p = p'. τ = τ') waarbij, /s.(x) de waarde van de beeldvorming op het punt x voorstelt, ps de waarde van de traagheid van de stralen afgebeeld vanuit het schotpunt voorstelt, pbc de waarden van de traagheid van de stralen afgebeeld vanuit het venstercentrum voorstelt, A een wegingfunctie voorstelt, en Ds het resultaat van de ontbinding van een locale vlakke golf voorstelt; en stap F : verzameling van de beeldvormingresultaten van alle stralenbundelparen ten einde de uiteindelijke beeldvormingresultaten van de migratie te bekomen.L ( χ ) = Σ ƒ d P * ƒ d Pbe A D s (L, p = p '. Τ = τ') where, / s . (X) represents the value of the imaging at the point x, p s the value of the inertia of the rays depicted represents from the shot point, p bc represented the values of the inertia of the beam representing from the window center, a represents a weighting function, and D s represents the result of the dissolution of a local plane wave; and step F: collection of the imaging results of all the beam bundles to obtain the final imaging results of the migration. 2. Kirchhoffs stralenbundelmigratiewerkwijze met complexe topografie volgens conclusie 1, met het kenmerk, dat in stap C: het hoekenbereik van de stralenemissie van -60 graden tot +60 graden gaat, en het interval van de emissiehoeken tussen twee aanliggende stralenKirchhoff's complex topography ray beam migration method according to claim 1, characterized in that in step C: the angle range of the beam emission ranges from -60 degrees to +60 degrees, and the interval of the emission angles between two adjacent beams Figure 1Figure 1 Figure 3Figure 3 5000 x/rri5000 x / rri 12501250 25002500 37503750 50005000 62506250 Figure 4Figure 4 1450014500 1 40001 4000 35003500 30003000 25002500 20002000 15001500 10001000 500 x/'m500 x / m 25002500 Figure 6Figure 6 12501250 37503750 50005000 62506250
NL2020152A 2017-11-06 2017-12-21 Kirchhoff Beam Migration Method Under Complex Topography NL2020152B1 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711077893.XA CN107870355B (en) 2017-11-06 2017-11-06 A kind of kirchhoff type beam deviation method under MODEL OVER COMPLEX TOPOGRAPHY

Publications (2)

Publication Number Publication Date
NL2020152A NL2020152A (en) 2018-01-25
NL2020152B1 true NL2020152B1 (en) 2018-09-11

Family

ID=61023307

Family Applications (1)

Application Number Title Priority Date Filing Date
NL2020152A NL2020152B1 (en) 2017-11-06 2017-12-21 Kirchhoff Beam Migration Method Under Complex Topography

Country Status (4)

Country Link
CN (1) CN107870355B (en)
BE (1) BE1025285B1 (en)
LU (1) LU100531B1 (en)
NL (1) NL2020152B1 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108802821B (en) * 2018-05-28 2019-11-08 中国石油天然气股份有限公司 Three-dimensional undulating surface seismic data migration imaging method, device and system
CN109917454B (en) * 2019-02-19 2020-10-09 中国石油天然气集团有限公司 True earth surface prestack depth migration imaging method and device based on double reference surfaces
CN110780341B (en) * 2019-10-15 2021-08-20 西南交通大学 Anisotropic seismic imaging method
US11474267B2 (en) * 2020-06-11 2022-10-18 China Petroleum & Chemical Corporation Computer-implemented method and system employing compress-sensing model for migrating seismic-over-land cross-spreads
CN111859268B (en) * 2020-06-29 2023-03-28 中铁二院工程集团有限责任公司 Magnetic tensor abnormal spatial domain fast forward modeling method based on grid point lattice
CN112904418B (en) * 2021-01-22 2021-08-17 西南交通大学 Self-adaptive ray encryption type kirchhoff type beam migration seismic wave imaging method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2004304279B2 (en) * 2003-12-12 2008-12-18 Exxonmobil Upstream Research Company Method for seismic imaging in geologically complex formations
CN102103216B (en) * 2009-12-16 2013-07-31 中国石油天然气集团公司 Prestack migration method of two-dimensional Gaussian ray bundle
CN102914791B (en) * 2011-08-05 2015-05-13 中国石油天然气集团公司 Kirchhoff prestack time migration method for processing seismic data of undulating surface
CN103995172B (en) * 2014-06-04 2017-01-18 广东电网公司佛山供电局 Method for on-line monitoring of load current of GIS bus of substation

Also Published As

Publication number Publication date
BE1025285A1 (en) 2019-01-04
LU100531B1 (en) 2018-03-19
CN107870355B (en) 2019-10-11
CN107870355A (en) 2018-04-03
BE1025285B1 (en) 2019-01-11
NL2020152A (en) 2018-01-25

Similar Documents

Publication Publication Date Title
NL2020152B1 (en) Kirchhoff Beam Migration Method Under Complex Topography
CN105549081B (en) Anisotropic medium is total to big gun domain Gaussian beam offset imaging method
CN108363101A (en) A kind of inclined shaft crosshole seismic Gaussian beam pre-stack depth migration imaging method
CN104765064A (en) Microseism interference imaging method
CN108872985B (en) Near-field circumference SAR rapid three-dimensional imaging method
CN104730579B (en) A kind of joint static correcting method of ripple in length and breadth based on calculation of near surface shear velocity inverting
CN102841379B (en) Method for analyzing pre-stack time migration and speed based on common scatter point channel set
CN102944894B (en) Earthquake prestack migration imaging method
CN102914796B (en) Control method for acquiring offset speeds of longitudinal and transverse waves based on Gaussian beam
CN106291687A (en) Anisotropy many ripples Gaussian beam pre-stack depth migration imaging method
CN107656308B (en) A kind of common scattering point pre-stack time migration imaging method based on time depth scanning
CN105093292A (en) Data processing method and device for seismic imaging
CN104215937B (en) Airborne radar space time adaptation processing method based on environment dynamic perception
Yan et al. Full-wave seismic illumination and resolution analyses: A Poynting-vector-based method
CN107132578A (en) A kind of microseism ground monitoring velocity model corrections algorithm
CN106154319A (en) A kind of method for separating of imaging road collection
CN109541677A (en) A kind of delay source pattern calculation method of parameters for protecting stimulating frequency
WO2023213052A1 (en) Deep-sea sound source localization method, computer device and storage medium
NL2024231B1 (en) Anisotropic seismic imaging method
CN104808242B (en) A kind of imaging of seismic data processing method
CN106526664A (en) Time-space domain geometric transform horizontal stacking seismic imaging method
CN107843924B (en) Utilize the seismic source location of moving constraint at the beginning of P wave and focal mechanism joint inversion method
CN110967751B (en) Positioning method of micro-seismic event based on ground shallow well monitoring and storage medium
CN115469362B (en) Energy flow density vector calculation method in seismic exploration
CN104502972B (en) Three component seismic wave integral offset method and apparatus

Legal Events

Date Code Title Description
MM Lapsed because of non-payment of the annual fee

Effective date: 20210101