PL237618B1 - Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections - Google Patents
Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections Download PDFInfo
- Publication number
- PL237618B1 PL237618B1 PL426304A PL42630418A PL237618B1 PL 237618 B1 PL237618 B1 PL 237618B1 PL 426304 A PL426304 A PL 426304A PL 42630418 A PL42630418 A PL 42630418A PL 237618 B1 PL237618 B1 PL 237618B1
- Authority
- PL
- Poland
- Prior art keywords
- seismic
- velocity field
- time domain
- domain
- migration
- Prior art date
Links
Description
Opis wynalazkuDescription of the invention
Dziedzina technikiThe field of technology
Przedmiotem wynalazku jest sposób budowy pola prędkości dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych, zarówno 2D jak i 3D, z wykorzystaniem pionowych profilowań sejsmicznych, uzyskanych w ramach badań sejsmicznych.The subject of the invention is a method of constructing a velocity field for the needs of temporal migration after assembling seismic sections, both 2D and 3D, with the use of vertical seismic profiling obtained during seismic surveys.
Badania sejsmiczne stanowią na dzień dzisiejszy najlepsze źródło informacji o budowie geologicznej skorupy ziemskiej. Badania te obejmują trzy zasadnicze etapy:Today, seismic surveys are the best source of information about the geological structure of the earth's crust. This research involves three main steps:
- wzbudzanie drgań sprężystych (sejsmicznych) na powierzchni lub pod powierzchnią ziemi i rejestracja energii w postaci amplitud tych drgań odbitych od powierzchni (granic) struktur geologicznych. Wzbudzanie odbywa się za pomocą ładunków wybuchowych w płytkich otworach, bądź za pomocą urządzeń mechanicznych (wibratory), a rejestracja odbywa się za pomocą odbiorników (czujniki elektryczne tzw. geofony). Wynik rejestracji przedstawiany jest w postaci tras sejsmicznych, będących zapisem następujących po sobie amplitud odbić zależnych od czasu propagacji fal. Rejestracja odbywa się w domenie czasu. Wyróżnia się dwa podstawowe schematy rejestracji (geometrii pomiaru): 2D i 3D. Rejestracja 2D odbywa się wzdłuż linii, przez co uzyskuje się informację o ośrodku geologicznym wzdłuż linii rejestracji, natomiast rejestracja 3D obejmuje pewną objętość ośrodka, przez co uzyskuje się informację o budowie przestrzennej badanego ośrodka geologicznego;- excitation of elastic (seismic) vibrations on the surface or below the ground and registration of energy in the form of amplitudes of these vibrations reflected from the surface (boundaries) of geological structures. Excitation is carried out by means of explosives in shallow holes, or by means of mechanical devices (vibrators), and registration is carried out by receivers (electrical sensors, so-called geophones). The recording result is presented in the form of seismic routes, which are a record of successive reflection amplitudes depending on the wave propagation time. Registration takes place in the time domain. There are two basic registration schemes (measurement geometry): 2D and 3D. 2D registration takes place along the line, which provides information about the geological center along the registration line, while the 3D registration covers a certain volume of the center, thus obtaining information about the spatial structure of the studied geological center;
- przetwarzanie zarejestrowanych tras sejsmicznych z wykorzystaniem specjalistycznych programów i sprzętu komputerowego (maszyn obliczeniowych), w wyniku którego otrzymuje się sekcje sejsmiczne czyli obraz amplitudowy (sejsmiczny) budowy strukturalnej badaneg o ośrodka geologicznego;- processing of registered seismic routes with the use of specialized programs and computer equipment (computing machines), as a result of which seismic sections are obtained, i.e. an amplitude (seismic) image of the structure of the studied geological center;
- interpretację geologiczną otrzymanych w procesie przetwarzania sekcji sejsmicznych, która jest realizowana w celu uzyskania możliwie jak najszerszego spektrum informacji o ośrodku geologicznym (budowa strukturalna, obraz facji, własności petrofizyczne itp.).- geological interpretation of seismic sections obtained in the processing of seismic sections, which is carried out in order to obtain the broadest possible spectrum of information about the geological center (structural structure, facies image, petrophysical properties, etc.).
Z powyższego widać, że proces przetwarzania jest niezwykle istotny i ma bezpośredni wpływ na wyniki interpretacji geologicznej.It can be seen from the above that the processing process is extremely important and has a direct impact on the results of the geological interpretation.
Jednym z głównych etapów mających wpływ na właściwe odtworzenie budowy geologicznej w procesie przetwarzania zarejestrowanych danych sejsmicznych jest proces migracji. Proces migracji polega na umiejscowieniu zarejestrowanych amplitud odbić sejsmicznych od granic strukturalnych w ich rzeczywistej pozycji oraz ich dobre zogniskowanie. Podstawowym parametrem używanym w procesie migracji jest pole prędkości, bądź czasy przebiegu promienia sejsmicznego (od źródła do odbiornika), w zależności od algorytmu jakim proces migracji jest realizowany. Należy nadmienić, iż przejście od czasów przebiegu promienia sejsmicznego do pola prędkości i odwrotnie realizowane jest poprzez podstawową formułę matematyczną, wykorzystywaną w pracy przez specjalistów zajmujących się przetwarzaniem danych sejsmicznych. Niewłaściwe określenie prędkości do procesu migracji skutkuje niewłaściwym umiejscowieniem elementów strukturalnych zarówno w pionie jak i w poziomie, ich zniekształceniem, a często ich eliminacją albo generacją nieistniejących. Taki obraz sejsmiczny nie przedstawia rzeczywistego ośrodka geologicznego, jest mało czytelny, o słabej rozdzielczości i trudny do interpretacji.One of the main stages influencing the proper reconstruction of the geological structure in the process of processing recorded seismic data is the migration process. The migration process consists in locating the registered amplitudes of seismic reflections from the structural boundaries in their actual position and their good focus. The basic parameter used in the migration process is the velocity field or the time of the seismic ray (from the source to the receiver), depending on the algorithm of the migration process. It should be mentioned that the transition from the times of the course of the seismic ray to the velocity field and vice versa is carried out by the basic mathematical formula, used at work by specialists in the processing of seismic data. Incorrect determination of the speed for the migration process results in improper location of structural elements both vertically and horizontally, their distortion, and often their elimination or generation of non-existent ones. Such a seismic image does not represent the real geological medium, it is not readable, with poor resolution and difficult to interpret.
Proces migracji można realizować w domenie czasu (migracja czasowa) albo głębokości (migracja głębokościowa) wykorzystując procedurę przed składaniem (prestack), bądź po składaniu (poststack). W procedurze przed składaniem migruje się trasy sejsmiczne, a potem sumuje (składa) otrzymując sekcję sejsmiczną, a w procedurze po składaniu najpierw sumuje się trasy sejsmiczne otrzymując sekcję sejsmiczną, a następnie się ją migruje. Procedura sumowania polega na dodawaniu wszystkich tras sejsmicznych zarejestrowanych przez odbiorniki od wspólnego punktu środkowego, w wyniku czego dostaje się jedną trasę sumaryczną. Zabieg ten stosuje się w celu wzmocnienia zarejestrowanych amplitud odbić sejsmicznych. Kolekcja sumarycznych tras sejsmicznych tworzy sekcję sejsmiczną.The migration process can be carried out in the time domain (time migration) or in depth (depth migration) using the pre-stack or post-stack procedure. In the pre-fold procedure, the seismic routes are migrated and then summed (folded) to obtain a seismic section, and in the post-folding procedure, the seismic routes are first summed to obtain a seismic section and then migrated. The summation procedure adds all seismic routes recorded by the receivers from a common midpoint, resulting in one summary route. This treatment is used to strengthen the recorded amplitudes of seismic reflections. A collection of summary seismic routes forms a seismic section.
Istnieje wiele rozwiązań budowy pola prędkości na potrzeby migracji sejsmicznej realizowanej przed i po składaniu, które lepiej lub gorzej sprawdzają się w realizacji tego celu, jednakże najczęściej stosowanym rozwiązaniem jest rozwiązanie oparte o analizę prędkości realizowaną w powszechnie stosowanych programach przeznaczonych do przetwarzania sejsmicznego.There are many solutions for the construction of the velocity field for the needs of seismic migration carried out before and after assembly, which are better or worse for this purpose, however, the most commonly used solution is the one based on velocity analysis implemented in commonly used seismic processing programs.
Dane z rejestracji pionowych profilowań sejsmicznych wykorzystywane były dotychczas w interpretacji danych sejsmicznych, interpretacji litostratygraficznej, śledzenia niewidocznych granic sejsmicznych, bądź słabo rejestrujących się w metodzie sejsmiki powierzchniowej, jak również do budoThe data from the registration of vertical seismic profiling have been used so far in the interpretation of seismic data, lithostratigraphic interpretation, tracking of invisible seismic boundaries or poorly recorded in the surface seismic method, as well as for construction
PL 237 618 B1 wy pola prędkości na potrzeby transformacji czas - głębokość. Metodyka rejestracji pionowych profilowań sejsmicznych wykorzystuje propagację drgań sprężystych (fal) podobnie jak do opisanych powyżej powierzchniowych badań sejsmicznych z tą różnicą, że odbiorniki rejestrujące amplitudę tych drgań znajdują się w otworze wiertniczym. Na ogół stosuje się trzy źródła wzbudzania w różnych odległościach od otworu i różnych azymutach. Wynikiem są zarejestrowane amplitudy fal padających i odbitych dla każdego punktu wzbudzania (trasy sejsmiczne).Speed field for the time-depth transformation. The methodology of registering vertical seismic profiling uses the propagation of elastic vibrations (waves) similar to the above-described surface seismic surveys, with the difference that the receivers recording the amplitude of these vibrations are located in the borehole. Generally, three excitation sources are used with different distances from the hole and different azimuths. The result is the recorded amplitudes of the incident and reflected waves for each excitation point (seismic traces).
Z opisu zgłoszenia patentowego US 005696735A i US 6002642A znana jest metoda wykorzystania pomiarów offsetowych prędkości średnich w otworach do procesu migracji. Pomiary te są realizowane na potrzeby określenia czasów przebiegów promieni sejsmicznych, bądź pola prędkości na potrzeby migracji zarówno prestack, jak i postack w domenie czasu i głębokości. Metoda ta określa geometrię pomiaru offsetowych prędkości średnich w zależności od geometrii badanego ośrodka geologicznego, oraz wykorzystuje zarejestrowane czasy i ich interpolację liniową bądź interpolację funkcjami sklejanymi (spline). Parametrem do migracji może być zarówno prędkość jak i czas rejestracji. Metoda ta uwzględnia również anizotropię prędkości ośrodka geologicznego.From the patent application US 005696735A and US 6002642A there is known a method of using offset average velocity measurements in the wells for the migration process. These measurements are carried out to determine the times of seismic rays or the velocity field for the needs of both prestack and postack migration in the time and depth domain. This method determines the geometry of the measurement of offset average velocities depending on the geometry of the geological medium under study, and uses the recorded times and their linear interpolation or interpolation with spline functions. The migration parameter can be both the speed and the recording time. This method also takes into account the velocity anisotropy of the geological medium.
Sposób według wynalazku wykorzystuje już istniejące pomiary, co znacznie redukuje koszty w stosunku do znanych rozwiązań. Dokonuje się interpolacji i ekstrapolacji prędkości. Interpolacja i ekstrapolacja prędkości realizowana jest w oparciu o obraz strukturalny ośrodka geologicznego, co zapewnia poprawne odwzorowanie i umiejscowienie struktur w przestrzeni oraz bardziej dokładne i szczegółowe odwzorowanie ośrodka geologicznego badanego rejonu. Pole prędkości zbudowane sposobem według wynalazku ma zastosowanie dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych.The method according to the invention uses already existing measurements, which significantly reduces costs as compared to known solutions. Speeds are interpolated and extrapolated. Velocity interpolation and extrapolation is carried out on the basis of the structural image of the geological center, which ensures correct mapping and location of structures in space as well as more accurate and detailed mapping of the geological center of the studied area. The velocity field constructed by the method according to the invention is applicable for the purpose of temporary migration after assembling the seismic sections.
Istotą rozwiązania według wynalazku jest to, że dane z rejestracji pionowych profilowań sejsmicznych poddaje się przetwarzaniu przy pomocy programu komputerowego, uzyskując wartości średnich prędkości offsetowych, które poddaje się przestrzennej interpolacji i ekstrapolacji przy użyciu programu komputerowego, w oparciu o model strukturalny w domenie czasu lub w domenie głębokości uzyskując przestrzenne pole prędkości średnich. Następnie za pomocą ogólnie znanych algorytmów matematycznych konwertuje się uzyskane pole prędkości średnich na pole prędkości interwałowych w domenie czasu, które wczytuje się jako parametr do migracji czasowej po składaniu sekcji sejsmicznych.The essence of the solution according to the invention is that the data from the registration of vertical seismic profiling are processed with a computer program, obtaining the values of average offset speeds, which are subjected to spatial interpolation and extrapolation with the use of a computer program, based on a structural model in the time domain or in the time domain. the depth domain, obtaining the spatial field of average velocities. The resulting mean velocity field is then converted into a time domain interval velocity field with the aid of generally known mathematical algorithms, which is read as a parameter for time migration after assembling the seismic sections.
Jeżeli model strukturalny jest w domenie głębokości, to przed dokonaniem konwersji na pole prędkości interwałowych w domenie czasu dokonuje się konwersji pola prędkości średnich w domenie głębokości na pole prędkości średnich w domenie czasu.If the structural model is in the depth domain, then before converting to a time domain interval velocity field, the depth domain average velocity field is converted to a time domain average velocity field.
Wynalazek skutecznie wychodzi naprzeciw potrzebie budowy poprawnego pola prędkości dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych 2D lub 3D, które uwzględnia budowę strukturalną badanego ośrodka geologicznego. W sposobie według wynalazku wykorzystuje się już istniejące rejestracje pionowych profilowań sejsmicznych w otworach wiertniczych, bez konieczności dodatkowych pomiarów prędkości dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych. Istotną zaletą sposobu według wynalazku jest również to, iż w rejestracji pionowych profilowań sejsmicznych zawarta jest informacja o anizotropii (zmiana parametru w zależności od kierunku) prędkości ośrodka geologicznego, a której nie uwzględniają standardowe analizy prędkości stosowane dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych 2D lub 3D. Brak informacji o anizotropii w polu prędkości dla potrzeb migracji czasowej po składaniu sekcji sejsmicznych 2D lub 3D prowadzi do błędów w obrazowaniu strukturalnym ośrodka geologicznego na sekcji sejsmicznej. Przestrzenna interpolacja i ekstrapolacja offsetowych prędkości średnich zawierających informację o anizotropii prędkości ośrodka geologicznego w oparciu o wyinterpretowane i dowiązane do reperów stratygraficznych w otworach granice strukturalne zapewnia poprawne odwzorowanie i umiejscowienie struktur w przestrzeni oraz bardziej dokładne i szczegółowe odwzorowanie ośrodka geologicznego badanego rejonu. W konsekwencji umożliwia to bardziej precyzyjne przedstawienie obrazu geologicznego badanego rejonu w stosunku do procesów migracji, opartych o standardowe procedury analizy prędkości zastosowane w modułach istniejących w powszechnie stosowanych oprogramowaniach przeznaczonych do przetwarzania sejsmicznego.The invention effectively meets the need to build the correct velocity field for the needs of temporary migration after assembling 2D or 3D seismic sections, which takes into account the structural structure of the geological center under study. The method according to the invention uses the already existing recordings of vertical seismic profiling in boreholes, without the need for additional velocity measurements for the purposes of time migration after assembling the seismic sections. An important advantage of the method according to the invention is also the fact that the registration of vertical seismic profiling includes information about the anisotropy (change of a parameter depending on the direction) of the speed of the geological medium, which is not taken into account by standard velocity analyzes used for the needs of temporal migration after assembling 2D seismic sections or 3D. Lack of information on anisotropy in the velocity field for the purpose of temporal migration after assembling 2D or 3D seismic sections leads to errors in the structural imaging of the geological medium on the seismic section. Spatial interpolation and extrapolation of offset average velocities containing information on the anisotropy of the geological medium velocity based on the structural boundaries interpreted and linked to the stratigraphic benchmarks in the boreholes ensures the correct mapping and location of structures in space and a more accurate and detailed mapping of the geological center of the studied area. As a consequence, it enables a more precise presentation of the geological image of the surveyed area in relation to migration processes, based on standard velocity analysis procedures used in modules existing in commonly used seismic processing software.
Wynalazek stosuje się identycznie dla geometrii pomiaru 2D jak i 3D.The invention applies identically to 2D and 3D measurement geometries.
Figury rysunku 1-17 przedstawiają:Figures 1-17 show:
fig. 1 - schemat pomiarowy sejsmiki powierzchniowej, fig. 2 - zapis tras sejsmicznych w postaci amplitud drgań sprężystych odbitych od granic strukturalnych ośrodka geologicznego, zarejestrowanego przez jeden odbiornik,Fig. 1 - surface seismic measurement diagram, Fig. 2 - record of seismic routes in the form of amplitudes of elastic vibrations reflected from the structural boundaries of the geological medium, recorded by one receiver,
PL 237 618 B1 fig. 3 - geometria rejestracji pionowych profilowań sejsmicznych, fig. 4 - schemat powierzchniowy rejestracji pionowych profilowań sejsmicznych; PW1, PW2,Fig. 3 - registration geometry of vertical seismic profiling, Fig. 4 - surface diagram of registration of vertical seismic profiling; PW1, PW2,
PW3 - punkty wzbudzania; W1 - otwór wiertniczy, fig. 5 - zarejestrowany pomiar w postaci tras sejsmicznych z jednego punktu wzbudzania w metodzie pionowych profilowań sejsmicznych, fig. 6 - obliczone offsetowe prędkości średnie w domenie czasu po przetworzeniu zarejestrowanych pól falowych dla trzech punktów wzbudzania w metodzie pionowych profilowań sejsmicznych, fig. 7 - obliczone offsetowe prędkości średnie w domenie głębokości po przetworzeniu zarejestrowanych pól falowych dla trzech punktów wzbudzania w metodzie pionowych profilowań sejsmicznych, fig. 8 - model strukturalny w domenie czasu, fig. 9 - model strukturalny w domenie głębokości, fig. 10 - przekrój sejsmiczny obrazujący pole prędkości interwałowych w domenie czasu dla potrzeb migracji czasowej po składaniu sekcji sejsmicznej (fig. 14), uzyskane w oparciu o analizy prędkości zrealizowane przy pomocy programu SeisSpace, fig. 11 - przekrój sejsmiczny w domenie czasu obrazujący pole prędkości średnich w domenie czasu dla potrzeb migracji czasowej po składaniu sekcji sejsmicznej (fig. 14), uzyskane w oparciu o wynalazek, fig. 12 - przekrój sejsmiczny obrazujący pole prędkości średnich w domenie głębokości dla potrzeb migracji czasowej po składaniu sekcji sejsmicznej (fig. 14), uzyskane w oparciu o wynalazek, fig. 13 - przekrój sejsmiczny obrazujący pole prędkości interwałowych w domenie czasu dla potrzeb migracji czasowej po składaniu sekcji sejsmicznej (fig. 14), uzyskane w oparciu o wynalazek, fig. 14 - sekcja sejsmiczna 2D obrazująca ośrodek geologiczny po procesie sumowania w domenie czasu, fig. 15 - sekcja sejsmiczna 2D obrazująca ośrodek geologiczny po procesie migracji czasowej po składaniu otrzymany po zastosowaniu pola prędkości, uzyskanego w oparciu o analizy prędkości zrealizowane przy pomocy programu SeisSpace, fig. 16 - sekcja sejsmiczna 2D obrazująca ośrodek geologiczny po migracji czasowej po składaniu uzyskany po zastosowaniu pola prędkości uzyskanego w oparciu o wynalazek, fig. 17 - blokowy schemat realizacji sposobu według wynalazku.PW3 - excitation points; W1 - borehole, Fig. 5 - recorded measurement in the form of seismic routes from one excitation point in the vertical seismic profiling method, Fig. 6 - calculated offset average velocities in the time domain after processing the recorded wave fields for three excitation points in the vertical seismic profiling method Fig. 7 - calculated offset mean velocities in the depth domain after processing the recorded wave fields for three excitation points in the vertical seismic profiling method, Fig. 8 - time-domain structural model, Fig. 9 - depth-domain structural model, Fig. 10 - seismic cross-section showing the time-domain interval velocity field for the purpose of temporary migration after assembling the seismic section (Fig. 14), obtained on the basis of velocity analyzes carried out with the SeisSpace program, Fig. 11 - seismic cross-section in the time domain showing the average velocity field in time domain for temporary migration after folding the seismic section (Fig. 14), obtained based on the invention, Fig. 12 - seismic section showing the mean velocity field in the depth domain for the purpose of time migration after assembling the seismic section (Fig. 14), obtained based on the invention, Fig. 13 - seismic section showing the field time domain interval velocities for the purpose of temporal migration after assembling the seismic section (Fig. 14), obtained based on the invention, Fig. 14 - 2D seismic section showing the geological medium after the time domain summation process, Fig. 15 - 2D seismic section showing geological center after the temporal migration process after folding obtained after applying the velocity field obtained on the basis of velocity analyzes carried out with the SeisSpace program, Fig. 16 - 2D seismic section showing the geological center after temporary migration after folding, obtained after applying the velocity field obtained on the basis of the invention, fig. 17 is a block diagram of the method according to the invention.
Realizacja wynalazku odbywa się przy pomocy dostępnych na rynku oprogramowań przeznaczonych do interpretacji sejsmicznej i przetwarzania sejsmicznego z użyciem sprzętu komputerowego.The invention is implemented using commercially available software for seismic interpretation and seismic processing with the use of computer hardware.
Realizację sposobu według wynalazku przedstawia schemat blokowy - fig. 17.The implementation of the method according to the invention is shown in a block diagram - Fig. 17.
P r z y k ł a d 1 - demonstracja sposobu według wynalazku na danych przemysłowych według schematu blokowego - fig. 17Example 1 - demonstration of the method according to the invention on industrial data according to a block diagram - Fig. 17
Z pomierzonych tras sejsmicznych - blok A, (fig. 5) metodą pionowych profilowań sejsmicznych PPS (fig. 3) i ich przetworzeniu przy pomocy programu komputerowego UNIVERS dokonano identyfikacji (pikowania) czasów przejścia fali podłużnej od punktu wzbudzania do odbiornika (czasy pierwszych wstąpień) (fig. 5). Czasy te zostały użyte do obliczenia offsetowych prędkości średnich w domenie głębokości - blok B, w oparciu o relację Vsr=S/T (S - droga przejścia fali podłużnej od punktu wzbudzania do odbiornika, T - czas przejścia fali od punktu wzbudzania do odbiornika) (fig. 7) dla każdego punktu wzbudzenia wybranego otworu (fig. 4). Powyższą czynność wykonano dla wszystkich otworów w rejonie badań, które posiadają rejestrację pionowych profilowań sejsmicznych. W wyniku tej operacji otrzymano zbiór offsetowych prędkości średnich w domenie głębokości. Prędkości te przekonwertowano do domeny czasu - blok C, przy pomocy programu komputerowego UNIVERS (fig. 6). W celu uzyskania informacji o prędkościach poza punktami wzbudzania dokonano przestrzennej interpolacji i ekstrapolacji - blok E, przy pomocy programu komputerowego Petrel w oparciu o model strukturalny w domenie czasu - blok D (fig. 8), którego granice strukturalne są dowiązane do reperów stratygraficznych w otworach występujących w obszarze badań sejsmicznych. Model strukturalny uzyskano w wyniku interpretacji sejsmicznej na bazie wstępnie przetworzonych sekcji sejsmicznych. W wyniku powyższej interpolacji i ekstrapolacji offsetowych prędkości średnich otrzymano przestrzenne pole prędkości średnich w domenie czasu - blok F (fig. 11).From the measured seismic routes - block A, (Fig. 5) using the PPS vertical seismic profiling method (Fig. 3) and their processing using the UNIVERS computer program, identification (quilting) of the longitudinal wave transit times from the excitation point to the receiver (first ascent times) was made (fig. 5). These times were used to calculate the offset mean velocities in the depth domain - block B, based on the relation Vsr = S / T (S - the path of the longitudinal wave from the excitation point to the receiver, T - the wave transit time from the excitation point to the receiver) ( Fig. 7) for each pick-up point of the selected opening (Fig. 4). The above operation was performed for all boreholes in the test area that have vertical seismic profiling records. As a result of this operation, a set of offset average speeds in the depth domain was obtained. These speeds were converted to the time domain - block C, by means of the UNIVERS computer program (Fig. 6). In order to obtain information about velocities beyond the excitation points, spatial interpolation and extrapolation - block E, with the use of the Petrel computer program, based on a structural model in the time domain - block D (Fig. 8), was performed, the structural boundaries of which are related to the stratigraphic benchmarks in the wells occurring in the area of seismic research. The structural model was obtained as a result of a seismic interpretation based on pre-processed seismic sections. As a result of the above interpolation and extrapolation of the offset average velocities, a spatial field of average velocity in the time domain was obtained - block F (Fig. 11).
Następnie przy pomocy programu komputerowego SeisSpace pole prędkości średnich w domenie czasu (fig. 11) przekonwertowano na pole prędkości interwałowych w domenie czasu - blok G (fig. 13), gdyż w tej postaci wczytuje się je jako parametr do migracji czasowej po składaniu w programie komputerowym SeisSpace - blok H. Figury 10 i 13 przedstawiają pola prędkości interwałowych w domenie czasu zastosowane do migracji czasowej po składaniu sekcji sejsmicznej (fig. 14). FiguThen, using the SeisSpace computer program, the time domain average velocity field (Fig. 11) was converted to the time domain interval velocity field - block G (Fig. 13), because in this form it is loaded as a parameter for time migration after folding in the program computerized SeisSpace - block H. Figures 10 and 13 show the time domain interval velocity fields used for temporal migration after assembling the seismic section (Fig. 14). Figu
PL 237 618 B1 ra 10 przedstawia pole prędkości interwałowych w domenie czasu otrzymane w programie SeisSpace, za pomocą którego została zrealizowana przy pomocy programu komputerowego SeisSpace i przedstawiona w przykładzie migracja czasowa po składaniu (fig. 15). Natomiast figura 13 przedstawia pole prędkości interwałowych w domenie czasu otrzymane według wynalazku, które stanowiło parametr do migracji czasowej po składaniu zrealizowanej za pomocą programu SeisSpace (fig. 16). Na figurze 13 widoczne jest zdecydowanie większe zróżnicowanie wartości pola prędkości interwałowych w domenie czasu (wartość wyrażona jest odcieniem szarości), którego zmienność oddaje budowę modelu strukturalnego.10 shows the time domain interval velocity field obtained in the SeisSpace program, with which the SeisSpace computer program was performed and the time migration after folding shown in the example (Fig. 15). In contrast, FIG. 13 shows the time domain interval velocity field obtained according to the invention, which was the parameter for post-folding temporal migration performed with the SeisSpace program (FIG. 16). Figure 13 shows a much greater variation in the value of the interval velocity field in the time domain (the value is expressed in a shade of gray), the variability of which reflects the structure of the structural model.
P r z y k ł a d 2 - demonstracja sposobu według wynalazku na danych przemysłowych według schematu blokowego - fig. 17Example 2 - demonstration of the method according to the invention on industrial data according to a block diagram - Fig. 17
Z pomierzonych tras sejsmicznych - blok A (fig. 5) metodą pionowych profilowań sejsmicznych PPS (fig. 3) i ich przetworzeniu przy pomocy programu komputerowego UNIVERS dokonano identyfikacji (pikowania) czasów przejścia fali podłużnej od punktu wzbudzania do odbiornika (czasy pierwszych wstąpień) (fig. 5). Czasy te zostały użyte do obliczenia offsetowych prędkości średnich w domenie głębokości - blok B, w oparciu o relację Vsr=S/T (S - droga przejścia fali podłużnej od punktu wzbudzania do odbiornika, T - czas przejścia fali od punktu wzbudzania do odbiornika) (fig. 7) dla każdego punktu wzbudzenia wybranego otworu (fig. 4). Powyższą czynność wykonano dla wszystkich otworów w rejonie badań, które posiadają rejestrację pionowych profilowań sejsmicznych. W wyniku tej operacji otrzymano zbiór offsetowych prędkości średnich w domenie głębokości. W celu uzyskania informacji o prędkościach poza punktami wzbudzania dokonano przestrzennej interpolacji i ekstrapolacji - blok E’, przy pomocy programu komputerowego Petrel w oparciu o model strukturalny w domenie głębokości - blok D’ (fig. 9), którego granice strukturalne są dowiązane do reperów stratygraficznych w otworach występujących w obszarze badań sejsmicznych. Model strukturalny uzyskano w wyniku interpretacji sejsmicznej na bazie wstępnie przetworzonych sekcji sejsmicznych. W wyniku powyższej interpolacji i ekstrapolacji offsetowych prędkości średnich otrzymano przestrzenne pole prędkości średnich w domenie głębokości - blok F’ (fig. 12), które przekonwertowano na pole prędkości średnich do domeny czasu - blok F”, za pomocą programu komputerowego Petrel (fig. 11).From the measured seismic routes - block A (Fig. 5) using the PPS vertical seismic profiling method (Fig. 3) and their processing using the UNIVERS computer program, identification (quilting) of the longitudinal wave transit times from the excitation point to the receiver (first ascent times) ( Fig. 5). These times were used to calculate the offset mean velocities in the depth domain - block B, based on the relation Vsr = S / T (S - the path of the longitudinal wave from the excitation point to the receiver, T - the wave transit time from the excitation point to the receiver) ( Fig. 7) for each pick-up point of the selected opening (Fig. 4). The above operation was performed for all boreholes in the test area that have vertical seismic profiling records. As a result of this operation, a set of offset average speeds in the depth domain was obtained. In order to obtain information about velocities beyond the excitation points, spatial interpolation and extrapolation - block E ', with the use of the Petrel computer program, based on the structural model in the depth domain - block D' (Fig. 9), was performed, the structural boundaries of which are related to the stratigraphic benchmarks. in boreholes occurring in the area of seismic surveys. The structural model was obtained as a result of a seismic interpretation based on pre-processed seismic sections. As a result of the above interpolation and extrapolation of the offset average velocities, the spatial velocity field of average velocities in the depth domain - block F '(Fig. 12) was obtained, which was converted into the average velocity field to the time domain - block F ", by means of the Petrel computer program (Fig. 11). ).
Następnie przy pomocy programu komputerowego SeisSpace pole prędkości średnich w domenie czasu (fig. 11) przekonwertowano na pole prędkości interwałowych w domenie czasu - blok G’ (fig. 13), gdyż w tej postaci wczytuje się je jako parametr do migracji czasowej po składaniu w programie komputerowym SeisSpace - blok H. Figury 10 i 13 przedstawiają pola prędkości interwałowych w domenie czasu zastosowane do migracji czasowej po składaniu (fig. 14). Figura 10 przedstawia pole prędkości interwałowych w domenie czasu otrzymane w wyniku analizy prędkości w programie SeisSpace, za pomocą którego została zrealizowana przy pomocy programu komputerowego SeisSpace i przedstawiona w przykładzie migracja czasowa po składaniu (fig. 15). Natomiast figura 13 przedstawia pole prędkości interwałowych w domenie czasu otrzymane według wynalazku, które stanowiło parametr do migracji czasowej po składaniu zrealizowanej za pomocą programu SeisSpace (fig. 16). Na figurze 13 widoczne jest zdecydowanie większe zróżnicowanie wartości pola prędkości interwałowych w domenie czasu (wartość wyrażona jest odcieniem szarości), którego zmienność oddaje budowę modelu strukturalnego.Then, using the SeisSpace computer program, the time domain average velocity field (Fig. 11) was converted to the time domain interval velocity field - block G '(Fig. 13), because in this form it is loaded as a parameter for time migration after folding into SeisSpace computer program - block H. Figures 10 and 13 show the time domain interval velocity fields used for post-composite temporal migration (Figure 14). Figure 10 shows the time domain interval velocity field obtained from the speed analysis in the SeisSpace program, with which the SeisSpace computer program was performed and the time migration after folding shown in the example (Figure 15). In contrast, FIG. 13 shows the time domain interval velocity field obtained according to the invention, which was the parameter for post-folding temporal migration performed with the SeisSpace program (FIG. 16). Figure 13 shows a much greater variation in the value of the interval velocity field in the time domain (the value is expressed in a shade of gray), the variability of which reflects the structure of the structural model.
Wyniki zastosowania pola prędkości do migracji czasowej po składaniu według sposobu będącego przedmiotem wynalazku przedstawiono na figurze 14, 15 i 16. Figura 14 przedstawia sekcję sejsmiczną 2D w postaci sumy przed procesem migracji czasowej po składaniu, figura 15 przedstawia sekcję sejsmiczną 2D po procesie migracji czasowej po składaniu po zastosowaniu pola prędkości opartego o standardowe procedury analiz prędkości, a figura 16 przedstawia sekcję sejsmiczną 2D po procesie migracji czasowej po składaniu po zastosowaniu pola prędkości uzyskanego w oparciu o wynalazek.The results of the application of the velocity field to post-spotting temporal migration according to the inventive method are shown in Figures 14, 15 and 16. Figure 14 shows a 2D seismic section as a sum before the post-splicing temporal migration process, Figure 15 shows a 2D seismic section after a post-splicing time migration process. folding after applying a velocity field based on standard velocity analysis procedures, and figure 16 shows a 2D seismic section after a post-folding temporal migration process after applying a velocity field obtained according to the invention.
Porównanie sekcji sejsmicznej (fig. 15) - przed zastosowaniem sposobu będącego przedmiotem wynalazku, oraz sekcji sejsmicznej (fig. 16) - po zastosowaniu sposobu będącego przedmiotem wynalazku - wyraźnie manifestuje wzrost przejrzystości obrazu sejsmicznego, rozdzielczości oraz detali obrazujących budowę geologiczną. Dynamika amplitud jest wyrównana, co przekłada się na ciągłość w zapisie amplitudowym odwzorowującym granice strukturalne. Widoczny jest również wzrost nieciągłości obrazujących uskoki.The comparison of the seismic section (Fig. 15) - before the application of the method according to the invention, and the seismic section (Fig. 16) - after the application of the method according to the invention - clearly demonstrates the increase in the clarity of the seismic image, resolution and details depicting the geological structure. The dynamics of the amplitudes is even, which translates into continuity in the amplitude record mapping the structural boundaries. An increase in the discontinuities depicting faults is also visible.
Claims (1)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PL426304A PL237618B1 (en) | 2018-07-11 | 2018-07-11 | Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PL426304A PL237618B1 (en) | 2018-07-11 | 2018-07-11 | Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| PL426304A1 PL426304A1 (en) | 2020-01-13 |
| PL237618B1 true PL237618B1 (en) | 2021-05-04 |
Family
ID=69161595
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PL426304A PL237618B1 (en) | 2018-07-11 | 2018-07-11 | Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections |
Country Status (1)
| Country | Link |
|---|---|
| PL (1) | PL237618B1 (en) |
-
2018
- 2018-07-11 PL PL426304A patent/PL237618B1/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| PL426304A1 (en) | 2020-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2457513C2 (en) | Methods and systems for processing microseismic data | |
| US10577926B2 (en) | Detecting sub-terranean structures | |
| Bakulin et al. | Smart DAS upholes for simultaneous land near-surface characterization and subsurface imaging | |
| US10884148B2 (en) | Amplitude compensation of reverse time migration (RTM) gathers for AVO/AVA analysis | |
| JP2974331B2 (en) | How to obtain a model that represents a heterogeneous medium, especially a stratum | |
| Haldorsen et al. | Locating microseismic sources using migration-based deconvolution | |
| Hornby et al. | Interferometric imaging of a salt flank using walkaway VSP data | |
| US5754493A (en) | Analysis of velocity data | |
| CN105182408A (en) | Manufacturing method and device for synthesizing earthquake record | |
| EA032186B1 (en) | Seismic adaptive focusing | |
| US6430508B1 (en) | Transfer function method of seismic signal processing and exploration | |
| CN111856577B (en) | Method for reducing calculation amount of reverse-time migration earth surface offset gather | |
| Dafni et al. | Scattering and dip angle decomposition based on subsurface offset extended wave-equation migration | |
| RU2541091C1 (en) | Method of producing apriori hodograph for performing lithologic and stratigraphic referencing | |
| US8144543B2 (en) | Seismic data processing method for RMO picking | |
| CN115469362B (en) | Energy flow density vector calculation method in seismic exploration | |
| CN111352153A (en) | Microseism interference positioning method based on instantaneous phase cross-correlation weighting | |
| Mildner et al. | Target-oriented velocity analysis using Marchenko-redatumed data | |
| PL237618B1 (en) | Method of building a velocity field for temporary migration after assembling of 2D or 3D seismic sections | |
| CN115061197A (en) | Two-dimensional sea surface ghost wave water body imaging measurement method, system, terminal and flow measurement equipment | |
| Zhang et al. | Detecting and locating microseismic events with stacking velocity analysis for surface monitoring | |
| Sun et al. | Anisotropic parameters estimation for VTI media based on first-arrival stereotomography | |
| Kessler et al. | Depth Imaging–More than PSDM | |
| RU2221262C1 (en) | Process of seismic prospecting for examination of sedimentary mantle in presence of heavily broken-up acoustically stiff boundaries ( variants ) | |
| Ziramov et al. | Kevitsa Ni-Cu-PGE deposit, North Finland? A seismic case study |