NL1023296C2 - Werkwijze voor het decimeren van seismische sporen. - Google Patents
Werkwijze voor het decimeren van seismische sporen. Download PDFInfo
- Publication number
- NL1023296C2 NL1023296C2 NL1023296A NL1023296A NL1023296C2 NL 1023296 C2 NL1023296 C2 NL 1023296C2 NL 1023296 A NL1023296 A NL 1023296A NL 1023296 A NL1023296 A NL 1023296A NL 1023296 C2 NL1023296 C2 NL 1023296C2
- Authority
- NL
- Netherlands
- Prior art keywords
- receiver
- track
- source
- mesh
- decimer
- Prior art date
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/003—Seismic data acquisition in general, e.g. survey design
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)
Description
Werkwijze voor het decimeren van seismische sporen.
BESCHRIJVING
De uitvinding betreft een werkwijze voor het decimeren van 5 seismische sporen met het oog op het leveren van een 3D diepte afbeelding van de ondergrond in een acceptabel tijdsbestek en tegen acceptabele kosten.
Men herinnert zich, dat de seismische reflectie bestaat uit in de ondergrond voortplanten van een reeks door een bron afgegeven 10 drukgolven en met behulp van ontvangers (hydrofonen bij seismisch onderzoek op zee en geofonen bij seismisch onderzoek aan land) op het oppervlak de signalen te registreren, die worden gereflecteerd door de tussenvlakken tussen de verschillende doorkruiste geologische lagen.
Iedere ontvanger bewerkstelligt een spoor genoemde registratie, welke 15 afhankelijk van de tijd wordt gevormd door een opeenvolging van signalen, die representatief zijn voor de desbetreffende geologische lagen.
Om juist een 2D doorsnede bij verwerving af te beelden weet men dat men moet registreren: voldoende dicht bij elkaar zijnde CMP sporen 20 (gemeenschappelijk middelpunt) langs een doorsnede om juist de helling van de lagen weer te geven; voldoende sporen met verplaatsing (offset) om de verhouding signaal/ruis optimaal te maken en om de diepe horizontale aardlagen juist af te beelden; 25 voldoende proeven per spoor voor het verkrijgen van een tevredenstellende verticale bepaling welke de fijne detectie van geologische lagen mogelijk maakt.
Het blijkt dus duidelijk, dat de omvang van seismische 2D gegevens gelijk is aan de opbrengst van het aantal CMP door het aantal 30 verplaatste, sporen. Voor een 2D lijn van 20 km moet men bijvoorbeeld normaal 50.000 sporen registreren.
1023296 i 2
Bij seismische 3D moet men het aantal sporen per lijn vermenigvuldigen met het aantal lijnen in de Y richting. Zo benodigt een 3D verkrijging over een zone van 25 x 25 km met de klassieke parameters het registreren van ongeveer 1,2.10® sporen. De met een dergelijke 5 verkrijging samenhangend numerieke omvang is ongeveer 100 Giga bit voor een luistertijd van 4 seconden.
Men begrijpt dat een gebruikelijke seismische verkrijging een zeer groot aantal seismische sporen voortbrengt, waarvan het opslaan een groot aantal magnetische dragers noodzakelijk zal maken. Dit is juist 10 een eerste nadeel van de gebruikelijke 2D en 3D verkrijging.
Een ander nadeel blijkt tijdens de behandeling van de geregistreerde sporen. De behandeling bestaat uit het omzetten van de seismische sporen voor het verkrijgen van een beeld van de ondergrond. Men zal hierna twee gebruikelijke behandelingswerkwijzen in het geval van 15 seismisch 3D onderzoek op zee uiteenzetten. Voor vereenvoudiging zal men veronderstellen, dat de gegevens worden behandeld in hun registreereenheid, welke de tijd is.
Gebruikelijke behandeling door in 3D verzameling brengen van 20 middenpunten.
Deze seismische handeling bestaat uit het bewerkstelligen van een voorbehandeling van de sporen op schietpunten, middenpunt in verzameling brengen of binning, een correctie van de helling en van de horizontalisering van de seismische gebeurtenissen, een optelling en 25 tenslotte een migratie. Tijdens deze behandeling worden de sporen van de schietpunten indien eenmaal voorbehandeld gereorganiseerd om te worden gesorteerd in opslagplaatsen (mazen van een raster waarin het opslaan plaatsvindt). Deze fase gaat vergezeld van een regelen van de offsets door sommeren hetgeen het aantal sporen vermindert.
30 Het opslaan brengt echter twee problemen voort: - het sommeren per klasse van offsets vindt onafhankelijk 1 023296 3 4 * van de azimut plaats, zodat de verschillende.seismische trajecten weer in een enkel spoor kunnen worden samengesteld, - de coördinaten van de bron en van de ontvanger worden tijdens het sommeren weer bijgewerkt. Indien men een Pre Stack Depth 5 Migration PSDM (diepte migratie voor totaal) uitvoert zal deze behandeling deze nieuwe coördinaten gebruiken, welke seismische trajecten kunnen weergeven, die zeer verschillend zijn van die, die tijdens het verkrijgen zijn geregistreerd. Men kan dan komen tot een snel heidsafwijking, dan volledig verborgen door dit type behandeling, 10 zoals blijkt uit de bijgaande figuren 1 en 2.
Deze figuren geven een verwervingsraster G weer, welke het opslaan van de behandelingsvolgorde bepaalt. Het raster is verdeeld in opslagplaatsen B in de richting van de lijnen (inline) en in de richting dwars op de lijnen (crossiine).
15 In de figuur 1 heeft men drie segmenten Lal, Li2, La3 weergegeven respectievelijk samenhangend met drie verworven sporen. Deze segmenten omvatten ieder een bron Sal, S^, Sa3, een ontvanger Ral, RaZ, Ra3 en een middenpunt Mal, Ma2, Ma3. Ieder spoor is loodrecht op het vlak van het veld loodrecht op het overeenkomstig middenpunt.
20 De eerste fase van de behandeling bestaat uit het middelen van de coördinaten van de bronpunten en van de ontvangerpunten van de sporen waarvan de middenpunten zich in eenzelfde opslagplaats bevinden. In de figuur 1 is dit het geval van de segmenten LaJ en LaZ, die samenhangen met de eerste en de tweede sporen. In de figuur 2 heeft men 25 de gemiddelde standen weergegeven van de bron S, van de ontvanger R, van het segment L, dat deze verbindt en van een middenpunt M (deze blijft in dezelfde opslagplaats) van het segment. Men heeft zonder modificatie in deze figuur het segment La2, dat samenhangt met het tweede spoor, weergegeven omdat er geen ander spoor bestaat waarvan het middenpunt zich 30 in dezelfde opslagplaats als Ma2 bevindt. Men zal nu opmerken dat het nieuwe segment S-R niet meer representatief is voor de verwerving omdat 1023296 t 4 j (- hij de snel hei dsafwijking A, die geopenbaard is door de verwerving, volledig verbergt.
Gebruikelijke seismische behandeling door decimeren 5 Een ander seismisch behandelingstype van verkregen gegevens, die ertoe neigt de hoeveelheid van te behandelen sporen te verminderen en dientengevolge de tijd en de kosten· van de berekening te verminderen, bestaat uit het zuiver en eenvoudig bepaalde onder deze sporen te decimeren. De meest eenvoudige en de meest klassieke wijze is 10 de sporen op statistische wijze te decimeren. Onder alle verworven sporen behoudt men slechts een monstercollectie, welke afhankelijk van een criterium van aantal is gekozen. In het algemeen wordt dit criterium niet op goed geluk gekozen, want de geofysicus moet het effect, dat de gekozen decimering zal hebben, kunnen voorspellen.
15 Onder de meest gebruikelijk uitgevoerde decimeringen kan men noemen: het elimineren van een spoor op twee in een gemeenschappelijke schietpunt verzameling: deze decimering heeft tot effect een middenpunt op twee te elimineren onder het behouden van de 20 nominale dekking. Hij wordt zeer vaak toegepast want de klassieke tussenruimte tussen middenpunten tijdens de verwerving is 6,25 m. Indien men een spoor op twee elimineert zal de tussenruimte 12,5 m zijn, hetgeen ruim voldoende is met het oog op de Fresnel zone ter hoogte van de doelen; 25 - de eliminatie van een spoor op n in en gemeenschappelijke middenpuntverzameling: dit decimeren wordt bewerkstelligd indien de dekking zeer belangrijk is en indien de verhouding signaal/ruis tevredenstellend is. Deze decimering kan echter een niet verwaarloosbaar effect hebben op het niveau van ruis. Bovendien kan de aanwezigheid van 30 veelvouden tengevolge hebben af te zien van dit genre decimeren, want de algoritmes "anti-veelvouden" maken een belangrijke dekking nodig om 1023296 5 * doeltreffend te zijn; - de eliminatie van een opslagplaats op n: de grootte van de opslagplaats wordt dan verdubbeld in de richting van de lijn. De meest belangrijke hellingen lopen echter gevaar door dit decimeren te worden 5 beïnvloed; - het elimineren van een lijn op n: dit genre decimeren is eenvoudig te bewerkstelligen, maar op grond van opèrationele redenen is het moeilijk een afstand kleiner dan 25 m tussen lijnen te verkrijgen; - het elimineren van een reeks offsets: indien de af te 10 beelden zone niet zeer diep is is het mogelijk de lange offsets die geen groot nut bezitten, af te wijzen. De verhouding signaal/ruis zal minder goed zijn maar de bepaling van de snelheden zal niet bestraft worden.
Al deze types van decimeren hebben het voordeel gemakkelijk te worden toegepast, maar zij bezitten meerdere nadelen: 15 - de waarde van decimeren is in de meeste gevallen van de vorm 2n (2, 4, 8, ...) hetgeen geen grote plooibaarheid biedt, - een enkele parameter wordt gekozen om het decimeren te sturen (spoornummer, waarde van de offset, CMP nummer), - een decimeerrichting wordt bevoorrecht, 20 - in het geval van een gedeeltelijke sommering worden de positioneerinformaties verloren, hetgeen voor een PSDM behandeling schadelijk is, - het is onmogelijk te weten of het decimeren/sommeren een ruimtelijk voldoende monstercollectie heeft behouden. De afgewezen sporen 25 kunnen onmisbaar voor de afbeelding zijn.
SAMENVATTING VAN DE UITVINDING
De huidige uitvinding beoogt tegemoet te komen aan de nadelen van de bekende werkwijze van decimeren en stelt hiertoe een 30 nieuwe werkwijze van decimeren voor bestaande uit: - voor ieder spoor de coördinaten van het overeenkomstige 1023290 · - ——i^ t 6 t koppel bron-ontvanger op te nemen, l - het koppel bron-ontvanger van ieder spoor op te stellen i op een decimeerraster, dat in mazen is verdeeld, - vervolgens alle koppels bron-ontvanger te hergroeperen in 5 verzamelingen van gelijk traject waarbij iedere verzameling de sporen hergroepeert waarvan de seismische trajecten gelijk zijn of nabij gelegen, - de sporen in iedere gelijk traject verzameling te ordenen, 10 - en het eerste spoor in iedere gelijk traject verzameling uit te kiezen.
De oorspronkelijkheid van de werkwijze volgens de uitvinding is gelegen in het feit, dat het decimeren niet meer op statische wijze wordt bewerkstelligd, maar dat deze is gebaseerd op 15 geofysische criteria. Alle sporen van de 3D verwerving zijn gesorteerd volgens het seismische traject, dat bepaald is door de coördinaten van de bron-ontvanger koppels. Deze voorbehandeling maakt het zo mogelijk de hoeveelheid van seismische 3D gegevens op een eenvoudige en doelmatige wijze te verminderen. Eenvoudig, want de geofysicus hoeft niet meer, 20 zoals voordien, moeilijke keuzes te maken (keuze van de richting van decimeren, keuze van de verzameling waarin de sporen zullen worden geëlimineerd). Nu wordt het decimeren op automatische wijze gestuurd op de seismische trajecten en de overbodige sporen worden verwijderd. Doelmatig, enerzijds omdat deze werkwijze het mogelijk maakt de 25 overtollige sporen snel te kiezen en anderzijds omdat deze werkwijze het in het bijzonder mogelijk maakt de tijd voor berekening van de PSDM op de 3D verwervingen aanzienlijk te verminderen zonder overigens de kwaliteit van het resultaat achteruit te doen gaan.
30 KORTE BESCHRIJVING VAN DE FIGUREN
De uitvinding zal beter worden begrepen bij het lezen van 1023296 7 de gedetailleerde beschrijving die gaat volgen en die gemaakt is aan de hand van bijgaande tekeningen, waarin: de figuur 1 een verwervingsraster weergeeft waarop de coördinaten zijn afgebeeld van de met drie verworven sporen samenhangende 5 bronnen en ontvangers; de figuur 2 de werkwijze van decimeren van de sporen volgens de stand van de techniek weergeeft, bestaande uit het bewerkstelligen van de sommeringen van de offsets; de figuur 3 een decimeerraster weergeeft dat dienst doet 10 voor het toepassen van de werkwijze voor decimeren volgens de uitvinding door verzameling van gelijk traject sporen welke werkwijze wordt toegepast bij het voorbeeld van de figuur 1; de figuur 4 een fase weergeeft van de werkwijze volgens de uitvinding die bestaat uit het respectievelijk op ieder bronpunt plaatsen 15 van de vier knopen van de maas, welke deze omlijst en bij ieder ontvangerpunt de vier knopen van de maas, die deze omlijst; de figuur 5 de zestien segmenten toont die respectievelijk een van de knopen van de bronmaas verbinden met een van de knopen van de ontvangermaas; en 20 de figuur 6 de bepaling toont van het meest representatieve bron-ontvanger segment van het met het verworven spoor samenhangend segment.
GEDETAILLEERDE BESCHRIJVING VAN DE UITVINDING 25 Basisidee van de uitvinding
Daar de figuren 1 en 2 reeds zijn beschreven gaat men direct over op de beschrijving van de figuren 3 tot 6.
Wanneer de gecompliceerdheid van de ondergrond zodanig is, dat een seismische behandeling in tijd onaangepast wordt houdt dit in, 30 dat het idee van het middenpunt geen zin meer heeft. Bij seismische dieptebehandeling en meer nauwkeurigheid bij PSDM behandeling zijn de 1023296 t 8 enige indicaties, die in rekening moeten worden genomen, de coördinaten van de bronnen en van de ontvangers van de sporen. Uitgaande van deze coördinaten kan men afhankelijk van een snel heidsmodel het parcours van de seismische golf reconstitueren en zijn voortplantingstijd berekenen.
5 Uitgaande van het idee, dat er voor een gegeven positie van een bron-ontvanger koppel een uniek samenhangend seismisch traject bestaat heeft aanvraagster bedacht dat het oordeelkundig zal zijn een voor de seismische trajecten gestuurde decimering te bewerkstelligen. Deze decimering zal bestaan uit het weglaten van alle bron-zender 10 koppels, die seismische trajecten hebben, die identiek of sterk overeenkomend zijn, behalve een.
Deze fase van de werkwijze volgens de uitvinding is weergegeven door de figuur 3 waarin men het decimeren heeft bewerkstelligd op de bron-ontvanger trajecten op het oppervlak in het 15 geval van dezelfde drie sporen van de figuur 1.
Volgens de werkwijze van de uitvinding is het met het segment La2 samenhangend spoor overtollig want deze is sterk soortgelijk met het met het segment Lai samenhangend spoor. Het segment is dus niet weergegeven op de figuur 3 ofschoon zijn middenpunt Ma2 zich bevindt 20 in een opslagplaats, die afwijkend is van die waar het middenpunt Mal van het overeenkomstige spoor zich bevindt. Na deze eliminatie blijft er op het raster van de figuur 3 het met het eerste spoor samenhangend segment Lal evenals het met het derde spoor samenhangend segment La3, dat behouden is. Men zal opmerken dat dankzij deze werkwijze van decimeren de zone van 25 snel heidsafwijking niet verborgen is.
Men zal op dezelfde wijze met alle sporen te werk gaan. Men sorteert vervolgens de sporen in gelijk traject verzamelingen waarbij iedere verzameling de sporen hergroepeert, die eenzelfde seismisch traject of een nabij gelegen seismisch traject hebben. In iedere 30 verzameling elimineert men alle sporen behalve een.
De overtolligheid van de sporen doet zich veelvuldig voor 1023296 9 * in seismische verwerving op zee. Dit kan zich bijvoorbeeld voordien indien twee opvolgend gevaren lijnen in term van afdekking gedeeltelijk op elkaar liggen of ook in het geval van met stroom samenhangende afdrijfwijziging, van wijziging van de verwervinginrichting, van het gaan 5 om een obstakel, enz.
Toeoassen van de werkwijze
De gecompliceerdheid van het probleem is gelegen in het feit, dat geen enkel bron-ontvanger koppel identiek is tijdens een 10 verwerving. Het is dus noodzakelijk alle trajecten, die nabij gelegen zijn te inventariseren.
Om dit te doen bepaalt men een nieuw raster, decimeerraster D (figuur 4) genaamd, waarvan de mazen afmetingen hebben, die aan de beoordeling van de geofysicus zijn overgelaten, zoals men hieronder zal 15 uiteenzetten. De maat van de maas, uitgedrukt in meters, varieert tussen een minimale maat, welke de maat is van een verwervingsopslagplaats en een maximale maat, die vijfmaal groter kan zijn.
Het decimeerraster wordt bepaald op hetzelfde gebied als het verwervingsraster. De maas van het decimeerraster kan overeenkomend 20 of verschillend zijn van die van het verwervingsraster.
De sporen worden behandeld in de volgorde waarin zij aankomen. Hiertoe plaatst men het decimeerraster bovenop het verwervingsraster Sa-Ra. De verwervingsbron Sa bevindt zich in een maas van het raster, dat bepaald wordt door vier koppels van waarden (x, y), 25 die de coördinaten zijn van de vier knopen of hoeken Slt S2, S3, St van het raster. Op dezelfde wijze behandelt men het verkregen ontvangstpunt Ra en men verkrijgt vier andere koppels (x, y), die de coördinaten zijn van de vier knopen of hoeken R1( Rz, R3, R< van de maas waar Ra zich bevindt. De acht koppels van zo verkregen waarden maken het mogelijk 30 ieder spoor in de term van traject te bepalen. Men kan zo voor alle sporen de offset SaRa en de azimut θβ (bijvoorbeeld ten opzichte van de as 1023296 ♦ 10
Oost-West) berekenen uitgaande van de verwervingscoördinaten.
Vervolgens wordt iedere knoop van de bronmaas in verhouding gebracht met iedere knoop van de ontvangermaas. Zoals de figuur 5 toont zijn zestien combinaties mogelijk (SjRj, S^, SiR3, SjR,, S2Rlt... S4R3, 5 S4R4) die ieder een fictief bron-ontvanger koppel dicht bij het verwervingskoppel SaRa weergeven. Voor ieder van deze zestien combinaties of trajecten berekent men de offset SnR„ en de azimut 0Snto op dezelfde wijze als men heeft gedaan met de werkelijke coördinaten van Sa en Ra.
Men beschikt zo voor ieder van de zestien koppels Snl^ over 10 de volgende informaties:
Afstand tussen Sn en R,,,
Verschil ten opzichte van de werkelijke offset SaRa Azimut van het segment ten opzichte van de as Oost-
West (N270°) 15 Verschil ten opzichte van het verwervingsazimut
Afstand tussen Sa en Sn Afstand tussen Ra en R,,
Zoals de figuur 6 toont maken deze waarden het mogelijk onder de zestien koppels die te verkrijgen, die het met segment SaRa 20 samenhangend werkelijk verkregen spoor het béste weergeeft, d.w.z. het segment, dat de offset en azimut het meest dichtbij de verwerving heeft. Men ziet in de figuur 6, dat in het onderhavige geval dit koppel SjR3 is. Zijn offset SrR3 en zijn azimut 0S1R3 zijn respectievelijk zeer dicht nabij Sa-Ra en 0a van het werkelijke traject.
25 De berekende afstanden Sa-Sn en Ra-R„ kunnen dienst doen voor het beslissen tussen de koppels, die dezelfde afstand hebben in offset en azimut ten opzichte van het werkelijke traject.
Duidelijk zijn hoe groter de mazen van het raster hoe groter de afstanden tussen het berekende spoor en het werkelijke spoor 30 zijn. Uit figuur 5 volgt, dat deze afstanden het moeilijker klein te maken zijn voor korte offsets dan voor lange offsets.
1023296 11
Men kan zo een decimeerraster kiezen met variabele met maasmaat voor het herstellen van een evenwicht in de decimering tussen de lange offsets en de korte offsets.
De werkelijke sporen worden voortaan opgespoord door een 5 tweede stel bron en ontvangercoördinaten, die overeenkomen met de bepaalde knopen van het decimeerraster. De aanvangscoördinaten, die ongewijzigd blijven worden bewaard in labels van de sporen teneinde te worden gebruikt in de PSDM behandeling. De waarden van van de offsets, van de azimuts en van de verplaatsingsafstanden van de bron en 10 ontvangerpunten worden eveneens opgeslagen teneinde de kwaliteit van het decimeren te kunnen controleren.
Teneinde de overtollige sporen te isoleren is het doelmatig alle sporen te sorteren naar toenemende orde volgens drie parameters in de volgende orde: 15 1. de coördinaten x en y van het bronpunt geplaatst op de knoop van het berekende raster, Sn, 2. de coördinaten x en y van de ontvangerplaats geplaatst op de knoop van het berekende raster, R,,, 3. de waarde van verplaatsing (S#-S„)+(Ra-Rj· 20 Men zal opmerken, dat in de term van seismisch traject en dus van afbeelding men op identieke wijze de koppels Sn-R„ en R,„-Sn zal beschouwen.
Dit sorteren maakt het mogelijk de sporen te rangschikken in de gelijk traject verzamelingen, die worden geacht een samenstel van 25 sporen weer te geven, die nagenoeg dezelfde geologische formaties afbeelden. Dit betekent ook, dat de door de seismische golf afgelegde afstand weinig verschillend is en dat de doorlooptijd vergelijkbaar is en dit is des te exacter naar mate de decimeermaas kleiner is.
Het decimeren volgens de uitvinding bezit de volgende 30 voordelen: het houdt rekening met seismische trajecten (vier 1023296 12 ι parameters) hetgeen uitstekend aangepast is voor PSDM behandeling, t zijn toepassen is zeer eenvoudig aangezien slechts de maat van de maas van het decimeerraster te bepalen is, het decimeerpercentage kan voldoende nauwkeurig worden 5 geregeld, het uitvoeren van het decimeren is zeer snel, het decimeren kan plaatsvinden wat ook de volgorde van binnenkomst van de sporen is, de labels van de sporen worden niet gewijzigd zodat de PSDM 10 behandeling kan worden uitgevoerd met de werkelijke verwervings- coördinaten.
Resultaten
Om de werkwijze volgens de uitvinding te testen heeft men 15 een PSDM studie uitgevoerd op een iri een zout houdende structuur gekozen testzone. De test is bewerkstelligd op 461 753 sporen.
Men heeft allereerst de invloed van de afmeting dx x dy van de mazen van het decimeerrooster op het decimeerpercentage bestudeerd. De onderstaande tabel toont dat hoe groter de maas is hoe groter de 20 decimeerwaarde belangrijk is.
Afmeting van dé maas Gemiddeld percentage van dx x dy (meters) decimeren 50 x 50 1,6 60 x 60 2,2 25 70 x 70 3,2 80 x 80 4,6 90 x 90 6,2 _100 x 100__8^2_
Vervolgens is de testzone verdeeld door reeksen offsets die 30 ieder een grootte van 500 meter hebben, teneinde de verdeling van de gedecimeerde sporen te meten. Het verslag is in overeenstemming met de 102329b 13 verwachtingen namelijk dat het decimeren belangrijker is op de korte offsets dan op de lange offsets.
Om de test te vervolgen is een maasafmeting van 70 x 70 m behouden teneinde een decimeerpercentage van 3,2 te verkrijgen. De 5 gedecimeerde gegevens zijn in diepte gemigreerd en het resultaat is op statistische wijze vergeleken met de gemigreerde beelden met een niet gedecimeerd stel gegevens en met hetzelfde gedecimeerde stel (een lijn op twee geëlimineerd). De vergelijking van de drie verkregen beelden toont dat het volgens de uitvinding gestuurd decimeren een minder ruis tonend 10 resultaat geeft dat dichter is bij het beeld, dat men heeft verkregen met alle sporen dan die dat verkregen is met de op statistische wijze gedecimeerde gegevens.
1023296 1
Claims (5)
1. Werkwijze voor het decimeren van seismische sporen in een verwerving van 3D gegevens, bestaande uit: 5. voor ieder spoor de coördinaten van het overeenkomend koppel bron-ontvanger (Sal, RaX) opnemen, - vervolgens alle koppels bron-ontvanger verenigen in verzamelingen van gelijk traject sporen waarbij iedere verzameling de sporen waarvan de seismische trajecten identiek of soortgelijk zijn 10 hergroepeert, gekenmerkt doordat het bestaat uit het op een in mazen verdeeld decimeerraster (D) opstellen van het koppel bron-ontvanger dat overeenkomt met ieder spoor, waarbij de fase van opstelling van de koppels bron-ontvanger van ieder spoor voor ieder spoor bestaan uit: 15. de bronmaas van het decimeerraster (D) waarin zich de bron (Sa) van het verkregen spoor bevindt op te sporen, - de coördinaten (x, y) van de vier knopen (Sla S2, S3, S4) van de genoemde maas op te nemen, - de ontvangermaas van het decimeerrooster waarin zich de 20 ontvanger (Ra) van het verkregen spoor bevindt, op te sporen, - de coördinaten (x, y) van de vier knopen (Rxa R2, R3, R4) van de genoemde maas op te nemen, en - de sporen in iedere gelijk traject verzameling te ordenen en het eerste spoor in ieder gelijk traject verzameling te kiezen waarbij 25 de fase van hergroeperen van de koppels bron-ontvanger in verzamelingen van gelijk traject sporen bestaat uit: - uitgaande van de verwervingscoördinaten voor ieder spoor de corresponderende statische afwijking (Sa-Ra) en de azimut (Θ.) te berekenen, 30. iedere knoop van de bronmaas in verband brengen met iedere knoop van de ontvangermaas waarbij de verkregen zestien 1023296 combinaties (S^,. S^, SjR3, SaR4, S^, ..... S4R3, S4R4) alle fictieve koppels bron-ontvanger dicht bij het verworven koppel (Sa, Ra) weergeven, i - voor ieder van de genoemde zestien combinaties de corresponderende statische afwijking (Sn-Rj en de azimut (0Sn> ,J te 5 berekenen, - onder de zo verkregen zestien koppels, (S^, 0SnRj die af te scheiden, die het dichtst bij het werkelijke traject ($aRa) is, en - dezelfde handelingen als de voorgaande te bewerkstelligen met alle verworven sporen waarvan de bron en de ontvanger zich 10 respectievelijk bevinden in dezelfde bronmaas en ontvangermaas teneinde gelijk traject verzamelingen te vormen.
2. Decimeerwerkwijze volgens conclusie 1, met het kenmerk, dat de fasen van ordenen en van keuze in een gelijk traject verzameling bestaan uit: 15. voor ieder spoor de verplaatsingswaarde (S4-Sn)+(Ra-Rj te bepalen, - de sporen naar volgorde van toenemende verplaatsing te ordenen, en - het met de minimale verplaatsing samenhangende spoor te 20 kiezen.
3. Decimeerwerkwijze volgens conclusie 1, met het kenmerk, dat de afmeting van de maas van het decimeerraster des te groter wordt gekozen naar het decimeerpercentage, dat men wenst te bereiken hoog is.
4. Decimeerwerkwijze volgens conclusie l.met het kenmerk, dat 25 de mazen van het decimeerraster van variabele grootte zijn.
5. Decimeerwerkwijze volgens conclusie 1, met het kenmerk, dat de sporen worden behandeld in de volgorde waar zij aankomen. 1023296
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0205625A FR2839368B1 (fr) | 2002-05-06 | 2002-05-06 | Methode de decimation de traces sismiques pilotee par le trajet sismique |
FR0205625 | 2002-05-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
NL1023296A1 NL1023296A1 (nl) | 2003-11-10 |
NL1023296C2 true NL1023296C2 (nl) | 2005-06-16 |
Family
ID=8871533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NL1023296A NL1023296C2 (nl) | 2002-05-06 | 2003-04-29 | Werkwijze voor het decimeren van seismische sporen. |
Country Status (6)
Country | Link |
---|---|
US (1) | US6865487B2 (nl) |
FR (1) | FR2839368B1 (nl) |
GB (1) | GB2388662B (nl) |
IT (1) | ITMI20030804A1 (nl) |
NL (1) | NL1023296C2 (nl) |
NO (1) | NO20031946L (nl) |
Families Citing this family (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7400552B2 (en) * | 2006-01-19 | 2008-07-15 | Westerngeco L.L.C. | Methods and systems for efficiently acquiring towed streamer seismic surveys |
US8488409B2 (en) | 2007-05-17 | 2013-07-16 | Westerngeco L.L.C. | Acquiring azimuth rich seismic data in the marine environment using a regular sparse pattern of continuously curved sail lines |
US8559265B2 (en) | 2007-05-17 | 2013-10-15 | Westerngeco L.L.C. | Methods for efficiently acquiring wide-azimuth towed streamer seismic data |
US9857491B2 (en) | 2008-05-15 | 2018-01-02 | Westerngeco L.L.C. | Multi-vessel coil shooting acquisition |
US8681580B2 (en) | 2008-05-15 | 2014-03-25 | Westerngeco L.L.C. | Multi-vessel coil shooting acquisition |
US8724426B2 (en) | 2008-06-03 | 2014-05-13 | Westerngeco L.L.C. | Marine seismic streamer system configurations, systems, and methods for non-linear seismic survey navigation |
US9594181B2 (en) | 2008-06-13 | 2017-03-14 | Westerngeco L.L.C. | Filtering and presentation of heading observations for coil shooting |
US9052411B2 (en) | 2008-06-13 | 2015-06-09 | Westerngeco L.L.C. | Method to determine the deviation of seismic equipment from a planned curved path |
SG193172A1 (en) * | 2008-08-11 | 2013-09-30 | Exxonmobil Upstream Res Co | Removal of surface-wave noise in seismic data |
CN102112894B (zh) * | 2008-08-11 | 2015-03-25 | 埃克森美孚上游研究公司 | 用地震表面波的波形评估土壤性质 |
US8483008B2 (en) * | 2008-11-08 | 2013-07-09 | Westerngeco L.L.C. | Coil shooting mode |
US8717846B2 (en) * | 2008-11-10 | 2014-05-06 | Conocophillips Company | 4D seismic signal analysis |
US8681581B2 (en) * | 2009-12-30 | 2014-03-25 | Westerngeco L.L.C. | Randomization of data acquisition in marine seismic and electromagnetic acquisition |
US8711654B2 (en) | 2009-12-30 | 2014-04-29 | Westerngeco L.L.C. | Random sampling for geophysical acquisitions |
US8537638B2 (en) * | 2010-02-10 | 2013-09-17 | Exxonmobil Upstream Research Company | Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration |
US8223587B2 (en) * | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US8756042B2 (en) | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
US8767508B2 (en) | 2010-08-18 | 2014-07-01 | Exxonmobil Upstream Research Company | Using seismic P and S arrivals to determine shallow velocity structure |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
CN103119552B (zh) | 2010-09-27 | 2016-06-08 | 埃克森美孚上游研究公司 | 同时源编码和源分离作为全波场反演的实际解决方案 |
AU2011337143B2 (en) | 2010-12-01 | 2016-09-29 | Exxonmobil Upstream Research Company | Simultaneous source inversion for marine streamer data with cross-correlation objective function |
EP2691795A4 (en) | 2011-03-30 | 2015-12-09 | CONVERGENCE SPEED OF COMPLETE WAVELENGTH INVERSION USING SPECTRAL SHAPING | |
EP2691794A4 (en) | 2011-03-31 | 2015-12-23 | Exxonmobil Upstream Res Co | METHOD FOR WAVELET DETERMINATION AND MULTIPLE PRECODION IN A COMPLETE WAVY-FIELD INVERT |
US9110177B1 (en) * | 2011-06-07 | 2015-08-18 | Russell R. Opfer | Generation and use of seismic survey designs violating Freznel zone constraints |
ES2640824T3 (es) | 2011-09-02 | 2017-11-06 | Exxonmobil Upstream Research Company | Utilización de la proyección sobre conjuntos convexos para limitar la inversión del campo de onda completa |
US8862408B2 (en) | 2011-09-28 | 2014-10-14 | Westerngeco L.L.C. | Determining one or more target positions in an acquisition domain for processing survey data |
US9103942B2 (en) | 2011-10-28 | 2015-08-11 | Westerngeco L.L.C. | Methods and systems for survey designs |
US9176930B2 (en) | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
RU2488145C1 (ru) * | 2012-01-10 | 2013-07-20 | Министерство образования и науки РФ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Уральский государственный горный университет" | Способ построения сейсмических изображений геологической среды |
MY170622A (en) | 2012-03-08 | 2019-08-21 | Exxonmobil Upstream Res Co | Orthogonal source and receiver encoding |
EP2926170A4 (en) | 2012-11-28 | 2016-07-13 | Exxonmobil Upstream Res Co | Q-Tomography with Seismic Reflection Data |
US9229122B2 (en) * | 2013-03-04 | 2016-01-05 | Cggveritas Services Sa | Image-domain 4D-binning method and system |
US10310122B2 (en) | 2013-03-04 | 2019-06-04 | Cgg Services Sas | Increasing similarity between seismic datasets |
RU2615591C1 (ru) | 2013-05-24 | 2017-04-05 | Эксонмобил Апстрим Рисерч Компани | Многопараметрическая инверсия через зависящую от сдвига упругую полноволновую инверсию (fwi) |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
AU2014309376B2 (en) | 2013-08-23 | 2016-11-17 | Exxonmobil Upstream Research Company | Simultaneous sourcing during both seismic acquisition and seismic inversion |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
WO2015171215A1 (en) | 2014-05-09 | 2015-11-12 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
CA2947410A1 (en) | 2014-06-17 | 2015-12-30 | Exxonmobil Upstream Research Company | Fast viscoacoustic and viscoelastic full-wavefield inversion |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
MY182815A (en) | 2014-10-20 | 2021-02-05 | Exxonmobil Upstream Res Co | Velocity tomography using property scans |
EP3234659A1 (en) | 2014-12-18 | 2017-10-25 | Exxonmobil Upstream Research Company | Scalable scheduling of parallel iterative seismic jobs |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
AU2015382333B2 (en) | 2015-02-13 | 2018-01-04 | Exxonmobil Upstream Research Company | Efficient and stable absorbing boundary condition in finite-difference calculations |
MX2017007988A (es) | 2015-02-17 | 2017-09-29 | Exxonmobil Upstream Res Co | Proceso de inversion de campo ondulatorio completo de multifase que genera un conjunto de datos libres de multiples. |
SG11201708665VA (en) | 2015-06-04 | 2017-12-28 | Exxonmobil Upstream Res Co | Method for generating multiple free seismic images |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
US10265872B2 (en) | 2015-09-09 | 2019-04-23 | Berkshire Grey, Inc. | Systems and methods for providing dynamic communicative lighting in a robotic environment |
CA2998522A1 (en) | 2015-10-02 | 2017-04-06 | Exxonmobil Upstream Research Company | Q-compensated full wavefield inversion |
CN108139498B (zh) | 2015-10-15 | 2019-12-03 | 埃克森美孚上游研究公司 | 具有振幅保持的fwi模型域角度叠加 |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
US10712460B2 (en) * | 2018-03-08 | 2020-07-14 | Chevron U.S.A. Inc. | System and method for improving resolution of digital seismic images |
US11733416B2 (en) * | 2020-12-22 | 2023-08-22 | Landmark Graphics Corporation | Automated horizon layer extraction from seismic data for wellbore operation control |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784335A (en) * | 1995-02-16 | 1998-07-21 | Elf Aquitaine Production | Method for making a 3-D cube of near traces from data acquired in a marine reflection seismic survey |
US6026059A (en) * | 1997-11-14 | 2000-02-15 | Petroleum Geo-Services (Us), Inc. | Method of creating common-offset/common-azimuth gathers in 3-D seismic surveys and method of conducting reflection attribute variation analysis |
GB2350428A (en) * | 1999-05-25 | 2000-11-29 | Baker Hughes Inc | Evaluating a seismic survey |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5430689A (en) * | 1991-07-03 | 1995-07-04 | Atlantic Richfield Company | Method for acquiring marine seismic data having statistically distributed azimuths and offsets |
US5511039A (en) * | 1993-06-03 | 1996-04-23 | Western Atlas International, Inc. | Method of performing high resolution crossed-array seismic surveys |
US5402391A (en) * | 1993-10-08 | 1995-03-28 | Geophysical Exploration & Development Corp. | Arrangement of source and receiver lines for three-dimensional seismic data acquisition |
US5924049A (en) * | 1995-04-18 | 1999-07-13 | Western Atlas International, Inc. | Methods for acquiring and processing seismic data |
GB2334334B (en) * | 1996-10-10 | 2000-04-05 | Geco Prakla | Method of generating a fold distribution and of evaluating a seismic survey |
US6018497A (en) * | 1997-02-27 | 2000-01-25 | Geoquest | Method and apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore |
US6106561A (en) * | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
US5963879A (en) * | 1997-11-26 | 1999-10-05 | Schlumberger Technology Corporation | Binning of three dimensional seismic data |
-
2002
- 2002-05-06 FR FR0205625A patent/FR2839368B1/fr not_active Expired - Fee Related
-
2003
- 2003-04-17 IT IT000804A patent/ITMI20030804A1/it unknown
- 2003-04-29 NO NO20031946A patent/NO20031946L/no not_active Application Discontinuation
- 2003-04-29 NL NL1023296A patent/NL1023296C2/nl not_active IP Right Cessation
- 2003-05-06 GB GB0310360A patent/GB2388662B/en not_active Expired - Fee Related
- 2003-05-06 US US10/429,978 patent/US6865487B2/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784335A (en) * | 1995-02-16 | 1998-07-21 | Elf Aquitaine Production | Method for making a 3-D cube of near traces from data acquired in a marine reflection seismic survey |
US6026059A (en) * | 1997-11-14 | 2000-02-15 | Petroleum Geo-Services (Us), Inc. | Method of creating common-offset/common-azimuth gathers in 3-D seismic surveys and method of conducting reflection attribute variation analysis |
GB2350428A (en) * | 1999-05-25 | 2000-11-29 | Baker Hughes Inc | Evaluating a seismic survey |
Also Published As
Publication number | Publication date |
---|---|
GB2388662B (en) | 2004-04-07 |
NO20031946L (no) | 2003-11-07 |
FR2839368A1 (fr) | 2003-11-07 |
NO20031946D0 (no) | 2003-04-29 |
ITMI20030804A1 (it) | 2003-11-07 |
US20040068377A1 (en) | 2004-04-08 |
GB0310360D0 (en) | 2003-06-11 |
FR2839368B1 (fr) | 2004-10-01 |
GB2388662A (en) | 2003-11-19 |
US6865487B2 (en) | 2005-03-08 |
NL1023296A1 (nl) | 2003-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NL1023296C2 (nl) | Werkwijze voor het decimeren van seismische sporen. | |
US5596547A (en) | Prestack seismic migration | |
US7085195B2 (en) | Method for detecting earth formation fractures by seismic imaging of diffractors | |
US9625593B2 (en) | Seismic data processing | |
Chang et al. | High-frequency Rayleigh-wave tomography using traffic noise from Long Beach, California | |
US20130135966A1 (en) | Acquisition and Processing of Multi-Source Broadband Marine Seismic Data | |
Nedimović et al. | Crooked-line 2D seismic reflection imaging in crystalline terrains: Part 1, data processing | |
US7974153B2 (en) | Three-dimensional seismic survey methods using a perturbation pattern to provide bin fractionation | |
WO2009077442A2 (en) | Method of processing data obtained from seismic prospecting | |
CN102590862A (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
GB2451949A (en) | 3D surface related multiple elimination for wide azimuth seismic data | |
RU2255355C2 (ru) | Способ обработки сейсмических данных для повышения пространственного разрешения | |
CN109765615A (zh) | 一种地层品质因子反演方法及装置 | |
AU2006237348B2 (en) | Method of processing seismic data for AVO or AVOA characterisation | |
Zhang et al. | Weighted stacking of ambient noise correlations based on the coherency of 3C data | |
AU2006237350B2 (en) | Seismic data processing method for RMO picking | |
AU721765B2 (en) | Method of generating a fold distribution and of evaluating seismic survey | |
US5629905A (en) | Method of datuming seismic data and method of processing seismic data | |
CA2161048C (en) | Prestack seismic migration | |
Brink et al. | Coverage and binning issues for marine seismic surveys | |
Osinowo | Reprocessing of regional 2D marine seismic data of part of Taranaki basin, New Zealand using Latest processing techniques | |
Xie et al. | 3D wavefront tomography, Part I: NIP wavefront tomography | |
AU2002247318A1 (en) | Seismic processing method to improve spatial resolution | |
Dufour | Refraction static analysis of P-S seismic data using the plus-minus time analysis method. | |
RU2030766C1 (ru) | Способ сейсмической разведки |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AD1A | A request for search or an international type search has been filed | ||
RD2N | Patents in respect of which a decision has been taken or a report has been made (novelty report) |
Effective date: 20050411 |
|
PD2B | A search report has been drawn up | ||
V1 | Lapsed because of non-payment of the annual fee |
Effective date: 20111101 |