BE1027342B1 - METHOD FOR ANISOTROPIC SEISMIC IMAGING - Google Patents
METHOD FOR ANISOTROPIC SEISMIC IMAGING Download PDFInfo
- Publication number
- BE1027342B1 BE1027342B1 BE20195787A BE201905787A BE1027342B1 BE 1027342 B1 BE1027342 B1 BE 1027342B1 BE 20195787 A BE20195787 A BE 20195787A BE 201905787 A BE201905787 A BE 201905787A BE 1027342 B1 BE1027342 B1 BE 1027342B1
- Authority
- BE
- Belgium
- Prior art keywords
- imaging
- anisotropic
- seismic
- point
- window center
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000013508 migration Methods 0.000 claims abstract description 27
- 230000005012 migration Effects 0.000 claims abstract description 27
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000010304 firing Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002688 persistence Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
-
- 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
- G01V2210/671—Raytracing
-
- 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
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
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)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
De onderhavige uitvinding openbaart een werkwijze voor anisotrope seismische beeldvorming omvattende de stappen in figuur 1. Volgens de werkwijze voor anisotrope seismische beeldvorming, wordt de verhouding van bijdrage van effectieve signalen tot het definitieve beeldvormingsresultaat verhoogd, en wordt de nauwkeurigheid van de berekening van de anisotrope migratiewerkwijze verbeterd.The present invention discloses an anisotropic seismic imaging method comprising the steps in Figure 1. According to the anisotropic seismic imaging method, the ratio of contribution of effective signals to the final imaging result is increased, and the calculation accuracy of the anisotropic migration method is increased. improved.
Description
TECHNISCH DOMEIN De onderhavige uitvinding behoort tot het domein van seismische migratiebeeldvorming en heeft in het bijzonder betrekking op een werkwijze voor anisotrope seismische beeldvorming.TECHNICAL DOMAIN The present invention belongs to the field of seismic migration imaging and relates in particular to a method for anisotropic seismic imaging.
ACHTERGROND In een gebruikelijke werkwijze voor migratiebeeldvorming, wordt een beoogd geologisch lichaam beschouwd als een isotroop medium, terwijl anisotropie gebruikelijkis bij geologische lichamen. Wanneer seismische gegevens met lange offsetafstand en brede azimut verwerkt worden, wordt geen rekening gehouden met de problemen dat migratie-energie niet beter gericht kan worden en de migratieruis verhoogd wordt, wat gemakkelijk veroorzaakt wordt onder invloed van anisotropie. Door deze problemen, kan de nauwkeurigheid van seismische beeldvorming afnemen en kunnen bepaalde problemen ontstaan voor het onderzoek naar olie en gas.BACKGROUND In a common migration imaging method, a target geologic body is considered an isotropic medium, while anisotropy is common in geologic bodies. When processing long offset distance, wide azimuth seismic data, the problems that migration energy cannot be better directed and the migration noise is increased, which is easily caused by anisotropy, is not taken into account. Due to these problems, the accuracy of seismic imaging can decrease and certain problems can arise for oil and gas research.
De Kirchhoff Type Dynamic Focusing Beam Migration werd geopenbaard in een doctoraatsproefschrift van de Jilin University in 2017. In het proefschrift wordt een anisotrope werkwijze van Kirchhoff geïntroduceerd van het type stralenbundelmigratie, en anisotrope straalopsporing wordt in de anisotrope stralenbundelmigratie geïntroduceerd door de werkwijze om anisotrope geologische lichamen te verwerken. Bovendien wordt een beeldvormingstest uitgevoerd op een anisotroop Hessmodel door de anisotrope Kirchhoff stralenbundelmigratiewerkwijze, en worden goede migratieresultaten verkregen.The Kirchhoff Type Dynamic Focusing Beam Migration was disclosed in a doctoral thesis from Jilin University in 2017. In the thesis, an anisotropic Kirchhoff method of ray beam migration type is introduced, and anisotropic beam detection is introduced in the anisotropic ray beam migration by the anisotropic body migration method. to process. In addition, an imaging test is performed on an anisotropic Hess model by the anisotropic Kirchhoff ray beam migration method, and good migration results are obtained.
In ‘Application of Pseudo-Acoustic Prestack Reverse-Time Migration and Imaging Conditions of Anisotropic Medium’ geschreven door Ayizemuguli Ruze en anderen in het vierde nummer van 2017 van Geophysical and Geochemical Exploration Calculation Technology, werd een werkwijze geopenbaard waarbij anisotrope reverse-tijdsmigratiebeeldvorming werd geïntroduceerd, werd een akoestische golfvergelijking van een een VTI-media (variabele timinginjectie) bestudeerd, werden geoptimaliseerde genormaliseerde gezamenlijke correlatie-beeldvormingscondities toegepast, werd de werkwijze gecontroleerd door het Hessmodel voor anisotropie en werd een goed beeldvormingseffect bereikt.In 'Application of Pseudo-Acoustic Prestack Reverse-Time Migration and Imaging Conditions of Anisotropic Medium' written by Ayizemuguli Ruze and others in the fourth issue of 2017 of Geophysical and Geochemical Exploration Calculation Technology, a method was disclosed in which anisotropic reverse time migration imaging was introduced , an acoustic wave equation of a VTI (variable timing injection) media was studied, optimized normalized joint correlation imaging conditions were applied, the method was controlled by the Hess model for anisotropy, and a good imaging effect was achieved.
Uit de bovenstaande voorbeelden, kan worden opgemaakt dat anisotrope gegevenslichamen in zekere mate goed door een gebruikelijke beeldvormingswerkwijze kunnen worden afgebeeld, maar dat de nauwkeurigheid van de beeldvorming nog voor verbetering vatbaar is.From the above examples, it can be seen that anisotropic data bodies can be mapped to some extent by a conventional imaging method, but that imaging accuracy is still capable of improvement.
SAMENVATTING Om de nauwkeurigheid van de berekening van een werkwijze voor anisotrope seismische beeldvorming te verbeteren, verschaft de onderhavige uitvinding een werkwijze voor anisotrope seismische beeldvorming.SUMMARY To improve the computation accuracy of an anisotropic seismic imaging method, the present invention provides an anisotropic seismic imaging method.
De werkwijze voor anisotrope seismische beeldvorming omvat de volgende stappen van stap 1, het lezen in een anisotroop parametermodel, van een P-golf-snelheidsmodel en een parameterbestand; stap 2, het uitvoeren van anisotrope straalopsporing op een schotpunt door een Runge-Kutta-werkwijze in verschillende richtingen, en het berekenen van de informatie over een stralenbundel die met elke straal overeenstemt; stap 3, het verdelen van seismische bestanden van enkelvoudige schoten in een veelheid van gegevenslichamen met een venster als een eenheid; stap 4, het berekenen van de gedeeltelijke afgeleide van de gegevenslichamen in de vensters naar tijd en de gedeeltelijke afgeleide van de gegevenslichamen naar ruimte, en het uitvoeren van lokale vlakgolfontbinding op de seismische bestanden in de vensters; stap 5, het uitvoeren van anisotrope straalopsporing op het venstermiddelpunt in verschillende richtingen, en het berekenen van de informatie over een stralenbundel die met elke straal overeenstemt; stap 6, het uitvoeren van beeldvormingsberekening op alle stralenbundelparen in het schotpunt en het venstermiddelpunt door een nieuwe beeldvormingsformule met een toegevoegde weegfunctie, en stap 7, het samenvoegen van beeldvormingsresultaten van alle stralenbundelparen, om zo een definitief migratiebeeldvormingsresultaat te verkrijgen.The anisotropic seismic imaging method comprises the following steps of step 1, reading into an anisotropic parameter model, a P-wave velocity model and a parameter file; step 2, performing anisotropic beam detection on a shot point by a Runge-Kutta method in different directions, and calculating the information on a beam of beam corresponding to each beam; step 3, dividing single shot seismic files into a plurality of data bodies with one window as a unit; step 4, calculating the partial derivative of the data bodies in the windows by time and the partial derivative of the data bodies by space, and performing local plane wave decomposition on the seismic files in the windows; step 5, performing anisotropic beam search on the window center in different directions, and calculating the information on a beam of beam corresponding to each ray; step 6, performing imaging computation on all ray beam pairs in the shot point and window center by a new imaging formula with an added weighting function, and step 7, merging imaging results of all ray beam pairs, so as to obtain a final migration imaging result.
Het anisotrope parametermodel in stap 1 omvat verder een anisotroop parametermodel en een anisotroop parametermodel; het parameterbestand omvat de afmetingen van een raster, de aanvankelijke stralenbundelbreedte, het aantal seismische kanalen, kanaalafstanden, het aantal monsterpunten in elk kanaal en minimale en maximale frequenties.The anisotropic parameter model in step 1 further includes an anisotropic parameter model and an anisotropic parameter model; the parameter file includes the dimensions of a grid, the initial beam width, the number of seismic channels, channel spacing, the number of sample points in each channel, and minimum and maximum frequencies.
Verder ziin de straalopsporingvergelijkingen in stap 2 de volgende: dx; _ dr Aj P18 781 dp; 1 Od dt TS 2 ÖX; DaDPi8j8xr waarbij x, staat voor de ruimtelijke positie van discrete punten; Pi + Pan en Pp; Staan voor de traagheidscomponent; T staat voor de seismische bewegingstijd; a, wordt berekend door de formule aq, = CGu/P1 €, staat voor de elasticiteitsmodulus, en p staat voor de dichtheid; 8; en 2, staan voor kenmerk-vectorcomponenten, en O een symbool is voor een gedeeltelijk afgeleide; Informatie over overeenstemmende stralenbundels wordt verkregen door een berekeningsformule van stralenbundelbreedte nadat de informatie over centrale stralen gekend is, waarbij de berekeningsformule van de stralenbundelbreedte W als volgt wordt getoond: 5 w=240 7 V, waarbij V, staat voor de snelheidswaarde aan het schotpunt en © staat voor de integraal van een straalpad op basis van snelheid; verder, in stap 3, is de centrale afstand van de vensters gewoonlijk gekozen uit 200 m tot 500 m en bedraagt de persistentielengte van de vensters 1,5 keer de aanvankelijke stralenbundelbreedte. verder, in stap 6, is de beeldvormingsformule met toegevoegde weegfunctie als volgt: 1,0) =F [dp, [dp,A-W,D,(L, P'.T") 7 waarbij /, staat voor de beeldvormingswaarde van een enkelvoudig schot; x staat voor de positie van een beeldvormingspunt; pP, en pP, respectievelijk staan voor traagheidsparameters van stralen afgevuurd vanuit het schotpunt en het venstermiddelpunt; A staat voor de amplitude; D, staat voor het lokale vlakgolfontbindingsresultaat; L staat voor de positie van het venstermiddelpunt r en 7’ staan voor de traagheids- en bewegingsparameters voor lokale hellingsuperpositie; de uitdrukking van de weegcoëfficiënt W, in de beeldvormingsformule wordt getoond als volgt: 2 W.- ape VW) ” Danen WEG) || Een WEG) | waarbij YW staat voor de seismische bestanden, WY, en YW, respectievelijk staan voor de gedeeltelijke afgeleiden van de seismische bestanden naar tijd en ruimte, W staat voor een reeks van punten die voldoen aan de vereisten van traagheid en bewegingstijd, en l,; staat voor seismische bewegingstijd.Furthermore, the beam tracing equations in step 2 are the following: dx; Dr Aj P18 781 dp; 1 Od dt TS 2 ÖX; DaDPi8j8xr where x, represents the spatial position of discrete points; Pi + Pan and Pp; Stand for the inertial component; T represents the seismic movement time; a, is calculated by the formula aq, = CGu / P1 €, represents the elastic modulus, and p represents the density; 8; and 2, represent feature vector components, and O is a symbol for a partially derivative; Information about matching beams is obtained by a beam width calculation formula after knowing the central beam information, the beam beam width calculation formula W is shown as follows: 5 w = 240 7 V, where V, is the velocity value at the shot point and © represents the integral of a velocity-based beam path; further, in step 3, the center distance of the windows is usually selected from 200 m to 500 m and the persistence length of the windows is 1.5 times the initial beam width. further, in step 6, the weighting added imaging formula is as follows: 1.0) = F [dp, [dp, AW, D, (L, P'.T ") 7 where /, represents the imaging value of a single shot; x represents the position of an imaging point; pP, and pP, respectively, represent inertial parameters of beams fired from the firing point and the window center; A represents the amplitude; D represents the local plane wave decomposition result; L represents the position of the window center r and 7 'represent the inertia and motion parameters for local slope superposition, the expression of the weighting coefficient W, in the imaging formula is shown as follows: 2 W.- ape VW) ”Danen WEG) || A WEG) | where YW represents the seismic records, WY, and YW, respectively, represent the partial derivatives of the seismic records by time and space, W represents a series of points meeting the requirements of inertia and movement time, and l ,; vo or seismic movement time.
In vergelijking met de stand van de technniek, heeft de werkwijze voor anisotrope beeldvorming de gunstige effecten dat een nieuwe weegcoëfficiënt aan de beeldvormingsformule toegevoegd wordt, zodanig dat de bijdrageverhouding van effectieve signalen tot een definitief beeldvormingsresultaat verhoogd wordt, en het vermogen tot anti-interferentie en de nauwkeurigheid van de berekening van een anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie, verbeterd worden.Compared to the prior art, the anisotropic imaging method has the beneficial effects of adding a new weighting coefficient to the imaging formula, such that the contribution ratio of effective signals to a final imaging result is increased, and the anti-interference and anti-interference ability. the accuracy of the calculation of a Kirchhoff anisotropic method of the beam migration type can be improved.
KORTE BESCHRIJVING VAN DE TEKENINGEN Fig. 1 is een stroomdiagram van een anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie; Fig 2. is een stroomverdelingsdiagram van een P-golf-snelheidswaarde van een Hessmodel;BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a flow chart of a Kirchhoff anisotropic process of the beam migration type; Fig 2 is a flow distribution diagram of a P wave velocity value of a Hess model;
Fig 3. is een stroomverdelingsdiagram van een anisotroop parameter 6 van het Hessmodel; Fig 4. is een stroomverdelingsdiagram van een anisotroop parameter s van het Hessmodel; Fig. 5 is een vergroot aanzicht van een lokaal beeldvormingsresultaat van een oorspronkelijke anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie van het Hessmodel; Fig. 6 is een vergroot aanzicht van een lokaal beeldvormingsresultaat van een nieuwe anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie van het Hessmodel.Fig 3 is a flow distribution diagram of an anisotropic parameter 6 of the Hess model; Fig 4 is a flow distribution diagram of an anisotropic parameter s of the Hess model; FIG. 5 is an enlarged view of a local imaging result of an original Kirchhoff anisotropic method of the Hess model beam migration type; FIG. 6 is an enlarged view of a local imaging result of a new Kirchhoff anisotropic method of the beam migration type of the Hess model.
GEDETAILLEERDE BESCHRIJVING De onderhavige uitvinding wordt verder in detail beschreven in combinatie met tekeningen en specifieke uitvoeringsvormen.DETAILED DESCRIPTION The present invention is further described in detail in conjunction with drawings and specific embodiments.
Fig. 1 toont een stroomdiagram van een werkwijze voor anisotrope seismische beeldvorming en die specifiek de volgende stappen omvat:FIG. 1 shows a flow diagram of an anisotropic seismic imaging method specifically comprising the following steps:
1. het lezen in een anisotroop parametermodel, van een P-golf-snelheidsmodel en een parameterbestand, waarbij het anisotrope parametermodel een anisotroop T -parametermodel en een anisotroop & -parametermodel omvat; het parameterbestand omvat de afmetingen van een raster, de aanvankelijke stralenbundelbreedte, het aantal seismische kanalen, kanaalafstanden, het aantal monsterpunten in elk kanaal en minimale en maximale frequenties;1. reading in an anisotropic parameter model, a P wave velocity model and a parameter file, wherein the anisotropic parameter model comprises an anisotropic T parameter model and an anisotropic & parameter model; the parameter file includes the dimensions of a grid, the initial beam width, the number of seismic channels, channel spacing, the number of sample points in each channel, and minimum and maximum frequencies;
2. het afvuren van stralen vanuit een schotpunt in verschillende richtingen, waarbij het hoekbereik van de stralen -70 graden tot 70 graden bedraagt, en de hoekafstand AO van de stralen gewoonlijk gekozen wordt als 2 graden tot 4 graden; anisotrope kinematische straalopsporingvergelijkingen worden opgelost door een Runge-Kutta-werkwijze, en de vergelijkingsgroep wordt getoond als volgt: dx; _ dr Aj P18 781 dp; 1 Ca EI UM ; dT 2 0%, DnPiS j8r waarbij x, staat voor de ruimtelijke positie van discrete punten; p; staat voor de traagheidscomponent; 7 staat voor de seismische bewegingstijd; 4, wordt berekend door de formule Gym = Cu/P + ©, Staat voor de elasticiteitsmodulus, en p staat voor de dichtheid; 8; en g, Staat voor kenmerk-vectorcomponenten; het verkrijgen van informatie over overeenstemmende stralenbundels door een berekeningsformule van stralenbundelbreedte nadat de informatie over centrale stralen gekend is, waarbij de berekeningsformule van de stralenbundelbreedte W als volgt getoond wordt: Oo w=2A0— V,2. firing beams from a shot point in different directions, the angular range of the beams being -70 degrees to 70 degrees, and the angular distance A0 of the beams is usually selected as 2 degrees to 4 degrees; anisotropic kinematic beam tracing equations are solved by a Runge-Kutta method, and the comparison group is shown as follows: dx; Dr Aj P18 781 dp; 1 Ca EI UM; dT 2 0%, DnPiS j8r where x, represents the spatial position of discrete points; p; represents the inertial component; 7 represents the seismic movement time; 4, is calculated by the formula Gym = Cu / P + ©, represents the elastic modulus, and p represents the density; 8; and g, represents feature vector components; obtaining information about matching beams by a beam width calculation formula after knowing the central beam information, the beam beam width calculation formula W being shown as follows: Oo w = 2A0— V,
waarbij V, staat voor de snelheidswaarde aan het schotpunt en © staat voor de integraal van een straalpad op basis van snelheid;where V, represents the velocity value at the shot point and © represents the integral of a jet path based on velocity;
3. het verdelen van seismische bestanden van enkelvoudige schoten in een veelheid van gegevenslichamen met een venster als een eenheid, waarbij de centrale afstand van de vensters gewoonlijk gekozen is uit 200 m tot 500 m en de persistentielengte van de vensters 1,5 keer de aanvankelijke stralenbundelbreedte bedraagt;3.dividing single shot seismic files into a plurality of data bodies with a window as a unit, with the center distance of the windows usually selected from 200 m to 500 m and the persistence length of the windows 1.5 times the initial beam width is;
4. het berekenen van de gedeeltelijke afgeleide van de gegevenslichamen in de vensters naar tijd en de gedeeltelijke afgeleide van de gegevenslichamen naar ruimte, en het uitvoeren van lokale vlakgolfontbinding op de seismische bestanden in de vensters;4. calculating the partial derivative of the data bodies in the windows by time and the partial derivative of the data bodies by space, and performing local plane wave decomposition on the seismic files in the windows;
5. het uitvoeren van anisotrope straalopsporing op het venstermiddelpunt in verschillende richtingen, en het berekenen van de informatie over een stralenbundel die met elke straal overeenstemt, waarbij Stap 5 vergelijkbaar is met Stap 2 en Stap 5 en Stap 2 als enige verschil hebben dat de posities van de coördinaten van het schotpunt en het venstermiddelpunt verschillend zijn;5.performing anisotropic beam tracing on the window center point in different directions, and calculating the information of a ray beam corresponding to each ray, Step 5 being similar to Step 2 and Step 5 and Step 2 having the only difference that the positions of the coordinates of the shot point and the window center are different;
6. het uitvoeren van beeldvormingsberekening op alle stralenbundelparen in het schotpunt en het venstermiddelpunt door een nieuwe beeldvormingsformule met een toegevoegde weegfunctie, waarbij een beeldvormingsformule van oorspronkelijke anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie getoond wordt als volgt: 1,.00=> [dp, | dp,A-D,(L,p'.T") waarbij /, staat voor de beeldvormingswaarde van een enkelvoudig schot; x staat voor de positie van een beeldpunt; Ps en P, respectievelijk staan voor traagheidsparameters van stralen afgevuurd vanuit het schotpunt en het venstermiddelpunt; A staat voor de amplitude; D, staat voor het lokale vlakgolf-ontbindingsresultaat; staat voor de positie van het venstermiddelpunt; p' en T' staan voor de traagheids- en bewegingstijdparameters voor lokale hellingsuperpositie; in de oorspronkelijke beeldvormingsformule T-P kunnen gegevenslichamen van het domein verkregen door de lokale vlakgolfontbinding, een effect teweegbrengen op een definitief beeldvormingsresultaat met gelijk gewicht en zonder verschil zolang aan beeldvormingsvoorwaarden voldaan wordt, maar, ongeldige 7—P gegevens van het domein kunnen in de stappen van het lokale vlakgolfontbindingsresultaat geïntroduceerd worden wegens problemen zoals een afvlakkingseffect, en de gegevens hebben een negatief effect op het definitieve beeldvormingsresultaat; In de onderhave uitvinding, wordt een nieuwe weegcoëfficiënt aan de oorspronkelijke beeldvormingsformule toegevoegd, zodanig dat de bijdrageverhouding van effectieve signalen tot een definitief migratieresultaat verhoogd wordt en de nieuwe beeldvormingsformule getoond wordt als volgt: 1,00= > [dp, |dp AW, D(L. p'7") waarbij een uitdrukking van de weegcoëfficiënt W, in de beeldvormingsformule getoond wordt als volgt: 2 W.- ape VW) ” Nae WEE) || ne VG) waarbij YW staat voor de seismische bestanden, WY, en YW, respectievelijk staan voor de gedeeltelijke afgeleiden van de seismische bestanden naar tijd en ruimte, W staat voor een reeks van punten die voldoen aan de vereisten van traagheid en bewegingstijd, en6. performing imaging computation on all beam pairs in the shot point and the window center by a new imaging formula with an added weighting function, showing an imaging formula of Kirchhoff's original anisotropic method of ray beam migration type as follows: 1.00 => [dp , | dp, AD, (L, p'.T ") where /, represents the imaging value of a single shot, x represents the position of a pixel, Ps and P, respectively, represent inertial parameters of beams fired from the firing point and the window center; A stands for amplitude; D stands for local plane wave decomposition result; stands for position of window center; p 'and T' stands for inertial and moving time parameters for local slope superposition; original imaging formula TP allows data bodies from the domain obtained by the local plane wave decomposition, will have an effect on a final imaging result of equal weight and no difference as long as imaging conditions are met, but invalid 7-P data of the domain may be introduced in the steps of the local plane wave decomposition result due to problems such as a flattening effect, and the data has a negative effect on the def initial imaging result; In the present invention, a new weighting coefficient is added to the original imaging formula, such that the contribution ratio of effective signals to a final migration result is increased and the new imaging formula is shown as follows: 1.00 => [dp, | dp AW, D (L. p'7 ") where an expression of the weighting coefficient W, in the imaging formula is shown as follows: 2 W.- ape VW)” Nae WEE) || ne VG) where YW is the seismic files, WY, and YW, respectively, represent the partial derivatives of the seismic files by time and space, W represents a series of points that meet the requirements of inertia and movement time, and
7. het samenvoegen van beeldvormingsresultaten van alle stralenbundelparen, om zo een definitief migratiebeeldvormingsresultaat te verkrijgen.7. merging imaging results from all ray beam pairs to obtain a final migration imaging result.
Controle van de simulatie: Het schema en de gunstige effecten van de onderhavige uitvinding worden gecontroleerd door een anisotroop Hessmodel. Fig. 2, Fig. 3 en Fig. 4 staan respectievelijk voor de verdeling van een P-golf-snelheidswaarde, de verdeling van een paramater ó en de verdeling van een parameter s van het anisotrope Hessmodel. Het model heeft horizontaal 3,617 rasterpunten, en de rasterafstand bedraagt 20 m; het model heeft 1,501 rasterpunten in de lengte, en de rasterafstand bedraagt 20 m.Simulation control: The scheme and beneficial effects of the present invention are controlled by an anisotropic Hess model. FIG. 2, FIG. 3 and FIG. 4 represent the distribution of a P-wave velocity value, the distribution of a parameter δ and the distribution of a parameter s of the anisotropic Hess model. The model has 3,617 grid points horizontally, and the grid spacing is 20 m; the model has 1,501 grid points in the length, and the grid spacing is 20 m.
De gegevensverzameling omvat 720 schoten, de schietmodus is eenzijdig schieten, de schietafstand bedraagt 100 m en de kanaalafstand bedraagt 40 m; en elk kanaal omvat 1,333 monsterpunten en het monsterinterval bedraagt 6 ms.The data collection includes 720 shots, the shooting mode is one-sided shooting, the shooting distance is 100m, and the channel distance is 40m; and each channel includes 1,333 sample points and the sample interval is 6 ms.
FIG. 5 is een oorspronkelijke anisotrope werkwijze van Kirchhoff van het type stralenbundelmigratie, en FIG. 6 is een migratieresultaat van de werkwijze volgens de onderhavige uitvinding.FIG. 5 is an original Kirchhoff anisotropic method of the beam migration type, and FIG. 6 is a migration result of the method of the present invention.
Uit een contrasterend resultaatsdiagram kan gezien worden dat het beeldvormingsresultaat van de onderhavige uitvinding minder migratieruis heeft, hogere signaal-ruisverhouding heeft en duidelijker is in de weergegeven geologische structuur.From a contrast result diagram, it can be seen that the imaging result of the present invention has less migration noise, higher signal-to-noise ratio, and is more apparent in the geological structure shown.
De werkwijze die door de onderhavige uitvinding geopenbaard wordt, is een belangrijke prestack-dieptemigratiewerkwijze van een anisotroop medium; ongeldige gegevens in lokale hellingsuperpositie worden niet speciaal verwerkt volgens de oorspronkelijke beeldvormingsformule; een nieuwe weegcoëfficiënt wordt aan de oorspronkelijke beeldvormingsformule toegevoegd, zodat de verhouding van bijdrage van effectieve signalen tot het definitieve beeldvormingsresultaat verhoogd wordt, en de nauwkeurigheid van de berekening van de anisotrope migratiewerkwijze verbeterd wordt.The method disclosed by the present invention is an important prestack depth migration method of an anisotropic medium; invalid data in local slope superposition is not specially processed according to the original imaging formula; a new weighting coefficient is added to the original imaging formula, so that the ratio of contribution of effective signals to the final imaging result is increased, and the calculation accuracy of the anisotropic migration method is improved.
Claims (2)
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910977970.XA CN110780341B (en) | 2019-10-15 | 2019-10-15 | Anisotropic seismic imaging method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| BE1027342B1 true BE1027342B1 (en) | 2021-01-07 |
Family
ID=68807915
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| BE20195787A BE1027342B1 (en) | 2019-10-15 | 2019-11-13 | METHOD FOR ANISOTROPIC SEISMIC IMAGING |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN110780341B (en) |
| BE (1) | BE1027342B1 (en) |
| NL (1) | NL2024231B1 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114428328B (en) * | 2020-09-21 | 2024-08-06 | 中国石油化工股份有限公司 | Wave packet type depth offset imaging method, device and storage medium |
| CN112630825B (en) * | 2020-12-02 | 2022-08-26 | 中国海洋大学 | Common offset domain Beam prestack time migration imaging method, system, medium and application |
| CN112904418B (en) * | 2021-01-22 | 2021-08-17 | 西南交通大学 | Self-adaptive ray encryption type kirchhoff type beam migration seismic wave imaging method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105549081A (en) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | Anisotropic medium common shot domain Gaussian beam migration imaging method |
| CN106291687A (en) * | 2016-07-21 | 2017-01-04 | 中国地质科学院地质研究所 | Anisotropy many ripples Gaussian beam pre-stack depth migration imaging method |
| NL2020152A (en) * | 2017-11-06 | 2018-01-25 | Univ Southwest Jiaotong | Kirchhoff Beam Migration Method Under Complex Topography |
| CN108363101A (en) * | 2018-02-02 | 2018-08-03 | 西安石油大学 | A kind of inclined shaft crosshole seismic Gaussian beam pre-stack depth migration imaging method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102495426B (en) * | 2011-12-02 | 2014-10-22 | 北京中科联衡科技有限公司 | Kirchhoff integral seismic imaging method |
| US9594176B1 (en) * | 2012-08-01 | 2017-03-14 | Z Terra Inc. | Fast beam migration using plane-wave destructor (PWD) beam forming |
| CN103984019B (en) * | 2014-06-07 | 2017-01-11 | 吉林大学 | Local relevant weighted earthquake beam synthesis method |
| US11041972B2 (en) * | 2017-06-14 | 2021-06-22 | Pgs Geophysical As | Methods and systems to enhance resolution of seismic images |
-
2019
- 2019-10-15 CN CN201910977970.XA patent/CN110780341B/en not_active Expired - Fee Related
- 2019-11-13 NL NL2024231A patent/NL2024231B1/en not_active IP Right Cessation
- 2019-11-13 BE BE20195787A patent/BE1027342B1/en not_active IP Right Cessation
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105549081A (en) * | 2016-01-29 | 2016-05-04 | 中国石油大学(华东) | Anisotropic medium common shot domain Gaussian beam migration imaging method |
| CN106291687A (en) * | 2016-07-21 | 2017-01-04 | 中国地质科学院地质研究所 | Anisotropy many ripples Gaussian beam pre-stack depth migration imaging method |
| NL2020152A (en) * | 2017-11-06 | 2018-01-25 | Univ Southwest Jiaotong | Kirchhoff Beam Migration Method Under Complex Topography |
| CN108363101A (en) * | 2018-02-02 | 2018-08-03 | 西安石油大学 | A kind of inclined shaft crosshole seismic Gaussian beam pre-stack depth migration imaging method |
Non-Patent Citations (3)
| Title |
|---|
| "Advances in Wave Propagation in Heterogenous Earth", vol. 48, 1 January 2007, ELSEVIER, ISBN: 978-0-12-018850-5, ISSN: 0065-2687, article VLASTISLAV CERVENÝ ET AL: "Seismic ray method: Recent developments", pages: 1 - 126, XP055706556, DOI: 10.1016/S0065-2687(06)48001-8 * |
| QIANG LIU ET AL: "The revised ray tracing in anisotropic VTI media and its application of Gaussian beam migration", SPG/SEG 2016 INTERNATIONAL GEOPHYSICAL CONFERENCE, BEIJING, CHINA, 20-22 APRIL 2016, 22 April 2016 (2016-04-22), pages 728 - 730, XP055707480, DOI: 10.1190/IGCBeijing2016-222 * |
| QINGDA LV ET AL: "Efficient Gaussian beam method in time domain for anisotropic media", INTERNATIONAL GEOPHYSICAL CONFERENCE, QINGDAO, CHINA, 17-20 APRIL 2017, 31 May 2017 (2017-05-31), pages 558 - 560, XP055707466, DOI: 10.1190/IGC2017-142 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110780341A (en) | 2020-02-11 |
| NL2024231B1 (en) | 2021-01-27 |
| CN110780341B (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| BE1027342B1 (en) | METHOD FOR ANISOTROPIC SEISMIC IMAGING | |
| CN106353821B (en) | Microseism event positioning method | |
| CN104765064B (en) | A kind of method of microseism interference imaging | |
| CN108241173B (en) | A kind of seismic data offset imaging method and system | |
| CN107132578B (en) | A kind of microseism ground monitoring velocity model corrections algorithm | |
| NL2020152B1 (en) | Kirchhoff Beam Migration Method Under Complex Topography | |
| CN106291687A (en) | Anisotropy many ripples Gaussian beam pre-stack depth migration imaging method | |
| CN112034520B (en) | Anisotropic medium dynamic focusing beam offset imaging method and system | |
| CN110389377B (en) | Microseism offset imaging positioning method based on waveform cross-correlation coefficient multiplication | |
| CN106154319A (en) | A kind of method for separating of imaging road collection | |
| CN107807393A (en) | Separate unit station collection preliminary wave Enhancement Method based on seismic interference method | |
| CN113466933B (en) | Depth weighting-based seismic slope tomography method | |
| CN111352153B (en) | Microseism interference positioning method based on instantaneous phase cross-correlation weighting | |
| CN111665563B (en) | Pre-stack migration vertical resolution evaluation method based on focusing analysis | |
| CN112630830B (en) | Reflection wave full waveform inversion method and system based on Gaussian weighting | |
| CN112748463A (en) | Local offset imaging method based on deep learning illumination analysis | |
| CN110873895A (en) | Variable grid micro-seismic reverse-time interference positioning method | |
| CN112305615B (en) | Seismic data angle domain common imaging point gather extraction method and system | |
| CN114942472A (en) | Offset imaging method and equipment based on uplink ray tracing strategy | |
| CN111505567B (en) | Multi-target tracking method based on single-vector hydrophone direction of arrival estimation | |
| CN108333575B (en) | Gaussian prior and interval constraint based time delay filtering method for mobile sound source | |
| CN115980837A (en) | Angular domain dynamic focusing beam reverse time migration method and device, electronic equipment and medium | |
| CN115993650B (en) | Earthquake interference imaging method based on prismatic waves | |
| Zheng et al. | Application of S-transform threshold filtering in Anhui experiment airgun sounding data de-noising | |
| Peng‐Cheng et al. | NEAR‐FIELD SEISMIC LOCALIZATION USING DELTA T MAPPING |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FG | Patent granted |
Effective date: 20210107 |
|
| MM | Lapsed because of non-payment of the annual fee |
Effective date: 20211130 |