Wyniki wyszukiwania - BazTech

1

Modelowanie pola falowego przy użyciu równania jednostronnego w ośrodkach o zmiennej anizotropii poprzecznej

Żuławiński K.

Nafta-Gaz

|

2018

|

R. 74, nr 12

894--897

PL

Modelowanie sejsmiczne przy użyciu jednostronnego równania w wersji pseudoakustycznej jest stosunkowo efektywną procedurą dla fal podłużnych propagujących w ośrodku izotropowym i anizotropowym opisanym parametrami Thomsena, a więc cechującym się słabą izotropią poprzeczną. W ośrodkach o zmiennych parametrach anizotropii lub różnych jej typach pojawia się problem aplikowania tej metody. Pionowe wektory falowe po przejściu z dziedziny współrzędnych przestrzennych do dziedziny liczb falowych obliczane są dla całego ośrodka. Modyfikacja wynikająca z występującej w modelowanym ośrodku anizotropii obejmuje więc całość modelu i propagujących w nim fal. Jest to istotne ograniczenie metody pseudospektralnej, ponieważ naturalna wydaje się budowa modelu ośrodka z lokalnie występującą anizotropią lub anizotropiami różnych typów, na przykład modelu z wkładkami anizotropowymi symulującymi spękania lub szczeliny w warstwie. Stosowana przez nas metoda modelowania oparta jest na określeniu czasowej częstości własnej równania dyspersyjnego uzyskanego z pełnego równania falowego, część teoretyczna opisuje detalicznie cały proces, szczególnie intensywne zastosowanie transformaty Fouriera, odgrywającej kluczową rolę w obliczeniach. W artykule pokazano rozwinięcie modelowania rozprzestrzeniania się fali w przypadku podziału ośrodka na części z różnymi typami anizotropii lub anizotropii i izotropii. Oparto go na podziale na obszary o różnych typach anizotropii, w których sygnał propaguje zgodnie z przypisanym mu zestawem liczb falowych w ich dziedzinie i jest łączony w pełne pole falowe we współrzędnych przestrzennych. Zaprezentowano przykłady takiego modelowania z falą wzbudzaną punktowo, dzięki czemu wyraźnie widać, jak czoło fali oddaje charakter ośrodka definiowany przez nachylenie osi symetrii izotropii poprzecznej oraz zmianę parametrów Thomsena. Wskazano również możliwości rozszerzenia metody przez podział ośrodka na większą liczbę obszarów lub ciągły rozkład typów i parametrów anizotropii, co wymaga zastosowania splotu dwuwymiarowego.

EN

Seismic modeling using the one-sided wave equation method with the use of second order pseudo-acoustic equation in wavenumber domain is a relatively efficient procedure for compressional wave propagation in isotropic and anisotropic media. By anisotropic media we mean weak polar anisotropy (transverse isotropy) determined by Thomsen anisotropic parameters. In complex media, where different types of anisotropy coexist, there is the problem with application of this seismic modeling process. Wavenumbers after transformation from spatial coordinates to wavenumber domain are computed for the whole medium. The consequent impact of anisotropy parameters, affects the whole medium and all propagating waves. This is significant limitation of the pseudo-spectral method, since a model of a medium with fragmentary occurring anisotropy or anisotropies of various types seems to be more realistic and useful, the models with sequences of generally horizontal, isotropic layers insert or cracked shale layer insert in isotropic medium for instance. The applied method is based on the solution of the full system of elastic equations, from which the eigenvalue frequency is determined. The theoretical part of this paper sets out the whole procedure in full details, particularly intensive application of the Fourier transform, playing a key role in the calculations. The article presents the development of wave propagation modeling in the case of division of the medium into parts with various types of TI anisotropy or anisotropy and isotropy. The modeled area is divided into distinct parts, in which longitudinal waves propagate according to a proprietary set of wavenumbers and are concatenated into a full wavefield in the spatial coordinate domain. In the presented examples of the described modeling the wave is generated by a point source and thus the wavefronts distinctly map the property of media characterized by the tilt of symmetry axis and the alteration of Thomsen parameters. The capabilities for further development of the method by division into more distinct areas and even by the application of the continuous distribution of parameters and anisotropy types are indicated, though two-dimensional convolution would be necessary.

2

Modelowanie i migracja zero-offset stref uskokowych w ośrodku o azymutalnej anizotropii HTI (horizontal transverse isotropy)

Kostecki A., Żuławiński K.

Nafta-Gaz

|

2017

|

R. 73, nr 12

918--927

PL

W artykule zaprezentowano studium modelowania sekcji czasowych zero-offset oraz odwzorowań migracyjnych dla uskoków warstw anizotropowych HTI (Horizontal Transverse Isotropy). Badano uskoki pionowe i nachylonych warstw w aspekcie wpływu kąta azymutalnego na charakter zero-offsetowej sekcji. Stwierdzono, że w tym przypadku pole falowe zero-offset jest niezależne od kąta azymutalnego, ponieważ kształtują je pionowe przebiegi fal. Odwzorowanie uskoków wykonywano, stosując migrację MG(F-K) w dziedzinie liczb falowych (K) i częstotliwości (F). Uzyskano rozdzielczość 5÷6 m dla częstotliwości dominującej jednostkowego sygnału 80 Hz i sygnału Rickera 30 Hz.

EN

In this paper, we present the study of zero-offset section modeling and migration for the vertical and inclined thrust of anisotropic HTI layers. We investigated the influence of the azimuthal angle on zero-offset sections. In the case of vertical and moderate inclined thrusts, as has been tested and found, the zero-offset sections are independent of the azimuthal angle because the vertical rays of waves, shape the zero-offset wavefield. The imaging of thrusts was carried out by using MG(F-K) migration in wave number – frequency domain. For a dominant frequency of 80 Hz spike signal and 30 Hz Ricker’s signal, a resolution of 5÷6 m was obtained.

3

Propagacja fal podłużnych w ośrodkach anizotropowych na przykładzie struktury uskokowej horstu

Kostecki A., Żuławiński K.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

211--221

PL

W artykule przedstawiono podstawy teoretyczne modelowania „w przód” w ośrodkach anizotropowych typu TI (Transverse Isotropy) zilustrowane przykładami propagacji fal podłużnych w strefie uskokowej. Te same przykłady posłużyły do wykonania wstecznej propagacji, tj. migracji zero-offset, umożliwiającej odwzorowanie założonego ośrodka.

EN

In this paper, we present fundamentals of compressional waves forward modeling in anisotropic media of TI (Transverse Isotropy) type on some examples of these waves propagation. Back propagation i.e. zero-offset migration was used to image the trust zone.

4

Forward and back-propagation of compressional waves in horizontal transverse isotropy (HTI) media

Kostecki A., Żuławiński K.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 212

1--218

EN

The present monograph, describing in detail compressional wave propagation in horizontal transverse isotropy (HTI) media as a function of azimuthal angle Ψ (angle between the two vertical planes: isotropy plane and plane of measurement), comprises the comprehensive study of wave spreading in anisotropic environment. The paper contains both the account of algorithmic relations of seismic modeling and migration and the abundant set of correctness verifying computer simulations as well. The starting point for theoretical considerations is basic relation between stress and strain according to Hooke's law and the elastic wave equation, which lay down the full system of elastic equations. The dispersion equation derived from it in wavenumber domain, allows to get its eigenvalues – temporal frequency and vertical wavenumber. Temporal frequency was used to build our one-way wave equation for forward-propagation modeling, while the vertical wavenumber provides the main element of wave propagator for depth extrapolation, i.e. for seismic migration. As the solution of third degree equation, we show the exact version and two approximate versions. Those approximate solutions may be successfully applied in practice as shown by numerous examples of propagation simulation (errors do not exceed 1.6%). The approximate version of the depth extrapolator was also validated for zero-offset migration. We verified these algorithms in the range of ε [–0.3; 0.3], δ [–0.2; 0.2] parameter variability, i.e. typical properties for rock anisotropy in oil exploration. The zero-offset migration algorithm, proposed in this monograph, is the adapted version of MG(F-K) migration in wavenumber and frequency domain, developed by the authors for compressional and converted wave seismic migration as well as applied in VTI (Vertical Transverse Isotropy) and TTI (Tilted Transverse Isotropy) media. The accuracy and high quality of wave propagation has been verified in plentiful zero-offset modeling and relevant seismic migration experiments for the two models, for the multilayer anticline and for the fault zone. The wavefield images, obtained by the one-way equation in wavenumber domain, are deprived of the noise, including inherent for the wave equations in space domain "diamond shape" noise and multiple waves interferences.

PL

Niniejsza monografia, traktująca o propagacji fal podłużnych w ośrodku o poprzecznie poziomej izotropii HTI (Horizontal Transverse Isotropy) jako funkcji azymutalnego kąta Ψ (pomiędzy pionowymi płaszczyznami: izotropii oraz pomiarową), stanowi obszerne studium właściwości rozprzestrzeniania się fal w ośrodku anizotropowym. Praca zawiera zarówno opis relacji algorytmicznych procesów modelowania i migracji, jak i bogaty zestaw przykładów dokumentujących poprawność symulacji komputerowych. Punktem wyjściowym w rozważaniach teoretycznych jest podstawowy związek pomiędzy naprężeniem a odkształceniem według prawa Hooke'a oraz równanie ruchu falowego, które formułują pełny system równań sprężystych. Wyprowadzone stąd równanie dyspersyjne w dziedzinie liczb falowych pozwala uzyskać wartość własną – częstotliwość czasową oraz pionową liczbę falową. Częstość czasowa posłużyła do sformułowania oryginalnego równania falowego jednostronnego, będącego podstawowym narzędziem modelowania „w przód", natomiast pionowa liczba falowa stanowi główny element propagatora falowego w procesie ekstrapolacji głębokościowej – migracji. W zakresie modelowania zaprezentowano ścisłą wersję rozwiązania równania trzeciego stopnia oraz dwie wersje aproksymacyjne, które – jak wykazały liczne przykłady symulacji propagacji – z powodzeniem mogą być również stosowane w praktyce (błędy ok. 1,6%). Także w zakresie migracji zero-offset stwierdzono przydatność aproksymacyjnej wersji ekstrapolatora. Weryfikację tych algorytmów przeprowadzono w zakresie zmienności parametrów: ε [−0,3; 0,3], δ [−0,2; 0,2], a więc w obszarze stanowiącym podstawowy przedmiot poszukiwań naftowych. Algorytm migracji zero-offset zaproponowany w niniejszej monografii jest zaadaptowaną wersją MG(F-K) migracji w dziedzinie liczb falowych i częstotliwości, opracowaną przez autorów dla izotropowej wersji dla fal podłużnych i przemiennych oraz dla anizotropowych ośrodków typu VTI (Vertical Transverse Isotropy) i TTI (Tilted Transverse Isotropy). Szeroki zakres modelowań sekcji czasowych zero-offset i odpowiadających im obrazów odwzorowań migracyjnych, wykonanych dla modelu wielowarstwowej antykliny i strefy uskokowej, potwierdziły wysoką dokładność i jakość propagacji falowej pozbawionej efektów zakłócających typu diamond shape, będących immanentną cechą równań falowych we współrzędnych przestrzennych oraz fal wielokrotnych w konsekwencji stosowania jednostronnego równania we współrzędnych liczb falowych.

5

Modeling and seismic migration in anisotropic media as a function of azimuthal angle – HTI(Ψ)

Kostecki A., Żuławiński K.

Nafta-Gaz

|

2016

|

R. 72, nr 9

679--690

EN

In this paper, we present a new way of modeling in an anisotropic medium based on a pseudo-acoustic one-way equation, derived from the full equations system of elasticity using eigenvalues of dispersion equation – time-frequency. The method was shown on the examples of signals propagation in anisotropic medium HTI(Ψ) (Horizontal Transverse Isotropy) as function of azimuthal angle Ψ and zero-offset time sections for anticline model. The correctness of the modeling results was verified by new migration MG(F-K) in wavenumbers (K) – frequency (F) domain with a depth operator of extrapolation, which uses a vertical wavenumber derived from dispersion relation.

PL

W artykule przedstawiono nowy sposób modelowania w ośrodkach anizotropowych, bazujący na jednostronnym równaniu falowym wyprowadzonym z pełnego systemu równań sprężystości i stosującym wartości własne równania dyspersyjnego tj. czasową częstotliwość. Metoda została zaprezentowana na przykładach propagacji w ośrodkach HTI(Ψ) (Horizontal Transverse Isotropy) w funkcji kątów azymutalnych Ψ oraz na sekcjach zero-offset dla modelu antykliny. Poprawność rezultatów modelowania została zweryfikowana za pośrednictwem migracji MG(F-K) w dziedzinie liczb falowych (K) i częstotliwości (F) stosującej operator głębokościowej ekstrapolacji z użyciem pionowej liczby falowej wyprowadzonej z relacji dyspersyjnej.

6

The pseudo-acoustic equations of the scalar wavefield in anisotropic media

Kostecki A., Żuławiński K.

Nafta-Gaz

|

2015

|

R. 71, nr 11

811--815

EN

In this paper we present a new formulation of scalar pseudo-acoustic wave equation in arbitrary anisotropic media TI (Transverse Isotropy) type for 2D and 3D cases. These equations, based on precise dispersion relation, determine eigenvalue – temporal frequency as the function of wavenumber and anisotropy parameters. Here we present a few snapshots obtained for different signals by one-way equation.

PL

W artykule przedstawiono nowe sformułowanie skalarnych pseudospektralnych równań falowych w dowolnych ośrodkach anizotropowych typu TI (Transverse Isotropy) dla przypadków 2D i 3D. Równania bazują na dyspersyjnych relacjach wyznaczających wartości własne – częstotliwości w funkcji liczb falowych i parametrów anizotropii na podstawie pełnego systemu równań akustycznych. Zaprezentowano kilka przypadków propagacji falowej określonej jednostronnym równaniem w formie migawkowych zdjęć dla różnych typów sygnału.

7

Modelowanie i migracja sekcji czasowych zero-offsetowych w ośrodkach TTI metodą pseudospektralną

Kostecki A., Żuławiński K.

Nafta-Gaz

|

2014

|

R. 70, nr 12

855--860

PL

W publikacji przedstawiono pseudospektralną metodę modelowania czasowych sekcji zero-offsetowych w anizotropowych ośrodkach TTI (Tilted Transverse Isotropy), opartą na jednostronnym pseudoakustycznym równaniu falowym. Równanie pseudoakustyczne wyprowadzono z dokładnej formuły dyspersyjnej dla ośrodka TTI. Rozpatrzone zostały obydwa przypadki akwizycji pomiarów, to jest „pod upad” i „z upadem” dla antyklinalnego modelu TTI.

EN

In this paper we present a pseudospectral method of the modeling of zero-offset seismic time-sections in anisot-ropic media of the TTI (Tilted Transverse Isotropy) type, based on one-wave pseudoacoustic equation. This equation was derived from a precise dispersion relation for 2D TTI media. Both cases for the acquisition of data along directions for „up-dip” and „down-dip” for two-dimensional anticlinal model TTI were considered. Obtained results were verified by depth migration MG(F-K) in wave number (k) – frequency (f) domain.

8

Odwzorowanie struktur wgłębnych w ośrodkach anizotropowych metodą migracji sejsmicznej

Kostecki A., Półchłopek A., Żuławiński K.

Prace Naukowe Instytutu Nafty i Gazu

|

2013

|

no. 191

1--134

PL

Na podstawie koncepcji MG(F-K) migracji (autorzy: A. Kostecki, A. Półchłopek) w ośrodkach izotropowych sformułowano algorytmy wstecznej propagacji fal w dziedzinie liczb falowych (K) i częstości (F) w ośrodkach anizotropowych VTI (vertical transverse isotropy), HTI (horizontal transverse isotropy) oraz TTI (tilted transverse isotropy). Pionowe liczby falowe kz dla poszczególnych typów ośrodków zostały wyprowadzone z dyspersyjnego równania dla nachylonej poprzecznej izotropii (TTI) rozpatrywanego dla przypadku, gdy oś symetrii leży w płaszczyźnie akwizycji pomiaru w kierunku „pod upad" i „z upadem" oraz w przypadku, gdy oś symetrii jest prostopadła do kierunku rozciągłości laminowanego ośrodka. Dla małych (do 5°) i dużych (powyżej 85°) kątów upadu przedstawiono aproksymacyjne formuły w formie explicite, natomiast w zakresie 5°-85° kątów upadu zastosowano numeryczną procedurę rozwiązywania równania wielomianowego czwartego stopnia dla określenia pionowej liczby falowej w funkcji poziomej liczby falowej, prędkości, parametrów anizotropii i kąta nachylenia laminowanego ośrodka. Test migracji (prestack) przeprowadzony w ośrodku VTI na silnie niejednorodnym modelu anizotropowym Marmousi potwierdził wysoką zdolność algorytmu migracji do odwzorowywania skomplikowanego modelu strukturalno-prędkościowego, Testy migracji zero-offset wykonane zostały na modelach kompozycyjnych łączących strukturalno-geometryczną formę antykliny z modelami poprzecznie izotropowymi: TTI, HTI, VTI. Sekcje czasowe zero-offset dla tych modeli opracowano oryginalną metodą jednostronnego równania falowego z zastosowaniem metody pseudospektralnej. Relacje dla równań falowych uzyskano z dyspersyjnego równania dla modelu TTI, definiując częstości własne dla poszczególnych typów ośrodka. Testy zero-offsetowych odwzorowań migracji wykazały satysfakcjonującą poprawność działania algorytmów i programów anizotropowej migracji.

EN

Based on the conception of MG(F-K) "isotropic" migration (A. Kostecki, A. Półchłopek) were formulated the algorithms for back propagation wave in domain wavenumber (K) — frequency (F) in anisotropic media: TTI (Tilted Transverse Isotropy), VTI (Vertical Transverse Isotropy), HTI (Horizontal Transverse Isotropy). Vertical wavenumber kz for these types media were derived from dispersion relation for TTI considered for this case when axis of symmetric lay in a plane of acquisition in direction "up dip" and "down dip" and in this case when axis symmetric is perpendicular the direction of the strike. For small angles (to 5°) and large angles (greater then 85°)' of inclination were presented aproximation formulas in explicite form, however for range 5°-85° angles used numerical soubroutine for solution of polinomial equation of fourth degree as function of horizontal wavenumber kz, velocity, parameters of anizotropy and angle of dipping. Prestack migration was tested on strongly inhomogeneous vertical transverse isotropy (VTI) Marmousi model. Obtained results confirmed the ability this method to proper map- ping media with complicated structures. For testing the zero offset migration the composition models i.e. anticlinal structure and TI (Transverse Isotropy) were used. The zero-offset time sections for TTI, HTI and VTI models were described by one-way equation wave with pseudospectral method. Relations for pseudoacoustic equations were obtained from eigenvalues of dispersion relations i.e. time frequency for TTI, HTI and VTI models of anisotropy. Experiments on synthetic wavefield have shown the proper imaging of assumed media.

9

Ewolucja algorytmu inwersji spektralnej oparta na wydzieleniu dubletu i normalizacji mocy

Żuławiński K.

Nafta-Gaz

|

2013

|

R. 69, nr 12

903--908

PL

Artykuł przedstawia rozwinięcie inwersji spektralnej opartej na dopasowaniu adaptacyjnym. Zawiera skrócony opis metody i analizuje jej słabe strony. Na podstawie analizy sygnału zgodnej z kryteriami rozdzielczości sugeruje możliwości podniesienia jej wydajności oraz dokładności.

EN

The article presents the development of spectral inversion based on adaptive matching. A short description of the method’s disadvantages is presented. There are suggested methods of precision and efficiency improvement by signal resolvability analysis.

10

Inwersja spektralna oraz dekompozycja spektralna w interpolacji cienkich warstw

Żuławiński K.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

269--275

PL

Inwersję spektralną i dekompozycję spektralną łączy algorytm oparty na transformacji falkowej o zmiennej długości okna. Interpretacja sejsmiczna, wykorzystująca detekcję i śledzenie cienkich warstw oraz ich wyklinowań, zyskuje narzędzie obiektywnej oceny wyboru możliwych z punktu widzenia interpretatora wersji. Pokazane zostanie zastosowanie obu metod na przykładach modelowych i rzeczywistych.

EN

Both spectral inversion and spectral decomposition are joined by the algorithm of wavelet transform with variable size window. Seismic interpretation with use of detection and tracking of thin-beds and its edges gains universal tool for objective choice among various possible versions. The application of spectral inversion and decomposition for seismic models and real data will be shown.

11

Inwersja spektralna oparta na adaptacyjnej metodzie tworzenia rozwinięć sygnałów falkowych

Żuławiński K.

Nafta-Gaz

|

2012

|

R. 68, nr 12

937--948

PL

Przedstawiono metodę inwersji spektralnej sygnału sejsmicznego, przeprowadzono testy metody oraz zaprezentowano wyniki uzyskane dla wybranego profilu sejsmicznego. Omówiono podstawy matematyczne, począwszy od transformaty Gabora, jej porównanie z transformatą Fouriera oraz metodę słownikową rozwinięcia sygnału. W oparciu o te narzędzia zbudowano równanie macierzowe prowadzące do rozwiązania zadania odwrotnego przez stosowane w programowaniu liniowym metody optymalizacji, czyli dopasowanie adaptacyjne.

EN

The method of spectral inversion of seismic signals was presented and tested on theoretical and field data. The fundamental mathematical concepts, beginning with the Gabor transform, in comparison with the Fourier transform, and overcomplete signal representation based on dictionary of expansion functions were treated. On this ground a matrix equation was built and the method of its solution by linear programming field applied.

12

Inwersja cienkich warstw metodą dekompozycji spektralnej

Żuławiński K.

Prace Naukowe Instytutu Nafty i Gazu

|

2010

|

nr 170 wyd. konferencyjne

577--580

PL

Autor prezentuje teoretyczne rozważania na temat wydzielenia cienkich warstw, podstawowy model klina Widessa i jego ograniczenia. Pokazano wpływ częstotliwościowej zawartości waveletu oraz jego kształtu na rozdzielczość sekcji sejsmicznej w obszarze występowania cienkich warstw. Przedstawiono zastosowanie metody dekompozycji spektralnej do wydzielenia uwarstwienia.

13

Ocena efektywności operatora ekstrapolacji pola falowego w głąb ośrodka anizotropowego typu VTI

Półchłopek A., Żuławiński K.

Nafta-Gaz

|

2010

|

R. 66, nr 12

1093-1101

PL

W artykule przedstawiono analizę rezultatów uzyskanych za pomocą operatora ekstrapolacji pola falowego w głąb ośrodka anizotropowego typu VTI z wykorzystaniem MG(F-K) migracji sejsmicznej, na przykładach syntetycznych pól falowych w wersji zero-offsetowej. Z przeprowadzonych eksperymentow wynika, że prezentowany operator jest skutecznym narzędziem do odwzorowania ośrodka anizotropowego typu VTI. Dokonano analizy błędów odwzorowania ośrodka, wynikających z zastosowania błędnie oszacowanych parametrów anizotropii. Omówiono również przyczyny powstawania artefaktów na zmigrowanych przekrojach głębokościowych oraz przedstawiono możliwości ich usuwania.

EN

In this article an analysis of the results obtained with downward wavefield extrapolation operator in depth domain in VTI anisotropic medium with use of MG(F-K) seismic migration is presented. The analysis was conducted on examples of zero-offset synthetic wavefields. The operator appeared to be effective tool for VTI anisotropic media mapping. The analysis of media mapping errors created by mistakenly estimated anisotropy parameters, the agents generating artefacts on migrated seismic profiles in depth domain and methods of their elimination are included.

14

Detekcja i rozdzielczość czasowa, a zdolność rozdzielcza waveletu sejsmicznego w cienkich warstwach

Żuławiński K.

Nafta-Gaz

|

2010

|

R. 66, nr 11

987-992

PL

W artykule omowiono zdolność rozdzielczą waveletu sejsmicznego w różnych aspektach. Sformułowano kryteria rozdzielczości waveletu sejsmicznego dla zgodnie i przeciwnie spolaryzowanych dubletów. Przedstawiono wpływ wyboru częstotliwości na rozdzielczość dubletów w cienkich warstwach.

EN

In this paper wavelet resolving power is discussed from various aspects. The criteria of wavelet resolution are formulated in regard to equal and opposite polarity doublets. The impact of frequency selection on the wavelet resolving power of the doublets in thin-bed layers is put forward.

15

Wysad solny w obrazie falowym rejestrowanym metodą sejsmiki fal odbitych

Jędrzejowska-Tyczkowska H., Żukowska K., Kostecki A., Półchłopek A., Bartoń R., Leginowicz A., Żuławiński K.

Prace Instytutu Nafty i Gazu

|

2008

|

nr 150 wydanie konferencyjne

129--137

PL

Interpretacja geologiczna danych sejsmicznych w niektórych przypadkach budowy geologicznej napotyka na znaczne trudności związane przede wszystkim z brakiem wyraźnej, a często i jakiejkolwiek, jednoznaczności pomiędzy obrazem zarejestrowanego pola falowego, przedstawionego w postaci tzw. sekcji sumowanej (w przestrzeni X, 2T) a spodziewanym obrazem odwzorowywanych struktur geologicznych. Jednym z możliwych przykładów takiej sytuacji jest obraz sejsmiczny wysadu solnego cechujący się, w centralnej partii, chaotyczną, bezrefleksową strefą zapisu drgań sejsmicznych. W publikacji przedstawiono metodykę rozszyfrowania takiego obrazu w sensie geometrii oraz jakościowej i ilościowej oceny pola prędkości fal podłużnych, opartą na interaktywnej analizie i modelowaniu pojedynczych rekordów sejsmicznych.

EN

The large hydrocarbon deposits are often associated with the presence of salt structures due to favorable conditions for accumulation and trapping of oil and gas. It is important to be able to image subsalt prospects correctly although it is very difficult task because of complexity of salt bodies' geometry and a large contrast of velocities. In publication we propose the use of deterministic seismic inversion paralelly with DMVA to solve the problem of subsult structures' imaging.

16

Sejsmika 4D z perspektywy polskich potrzeb i możliwości

Jędrzejowska-Tyczkowska H., Bartoń R., Misiarz P., Żuławiński K.

Prace Instytutu Nafty i Gazu

|

2004

|

nr 130

167--173

PL

W pracy omówiono możliwości zastosowania sekwencji procedur sejsmiki 4D wraz z podstawami teoretycznymi, zarysem historycznym oraz perspektywami. Zawarto również przykłady zastosowania metody 4D na świecie oraz opisano jej znaczenie dla inżynierii złożowej, projektowania i eksploatacji złóż i podziemnych magazynów.

EN

In this paper the fundamental assumptions of procedures of seismic 4D were discussed. Historical outline, prospects, the experience with 4D in Poland and examples of worldwide application in reservoir management of time-lapse seismic was presented.

17

Rola prędkości propagacji fal sprężystych w zadaniach prognozowania ciśnień złożowych na podstawie danych sejsmicznych

Jędrzejowska-Tyczkowska H., Bartoń R., Żuławiński K., Żukowska K.

Prace Instytutu Nafty i Gazu

|

2004

|

nr 130

181--187

PL

W referacie scharakteryzowano krótko genezę powstawania sfref ciśnień anomalnych (podwyższonych, obniżonych) występujących w górotworze, oraz zaprezentowano metodykę prognozowania podwyższonych ciśnień porowych na podstawie danych sejsmicznych — standardowej sekcji sejsmiki powierzchniowej oraz danych pochodzących z pionowego profilowania sejsmicznego (PPS). Podstawowym parametrem fizycznym ośrodka wykorzystywanym w zadaniach prognozowania ciśnień porowych jest prędkość propagacji fal sprężystych. Poprawność określenia tego parametru z danych sejsmicznych (najczęściej poprzez inwersję sejsmiczną), zależy od wielu czynników, w tym między innymi od zakresu częstotliwości obecnych w zapisie zarejestrowanego pola falowego. W pracy przedstawiono ocenę wpływu zwiększenia przedziału częstotliwości sekcji sejsmicznej na wyniki prognozowania ciśnień porowych. Poszerzenie zakresu częstotliwości uzyskano poprzez wprowadzenie do widma trasy sej smicznej składowej wysokoczęstotliwościowej z rejestracji pionowego profilowania sejsmicznego. Wynikająca stąd zmiana prędkości akustycznych ma swoje konsekwencje w zmianie wartości prognozowanych ciśnień. Przedstawiono algorytm i program obliczeniowy oraz wyniki obliczeń w postaci map zależnościP = f(x, h) wraz z ich analizą i interpretacją.

EN

The publication presents one of the possible way in forecasting pore pressure on the basis of seismic data. The special attention was paid an accuracy of determination of velocity.

18

Górnojurajska bioherma jako klaster obiektów - wyniki modelowania sejsmicznego

Misiarz P., Jędrzejowska-Tyczkowska H., Żuławiński K.

Prace Instytutu Nafty i Gazu

|

2004

|

nr 130

175--180

PL

Sejsmiczny obraz bioherm jest wysoce niekoherentny, mogąbyć więc one z powodzeniem lokalizowane, zarówno przy użyciu analizy atrybutów zespolonej trasy sejsmicznej, jak i analizy korelacji oraz wariancji 3D. Nie daje to jednak podstaw do wysnuwania wniosków na temat ich budowy wewnętrznej. Próba rozwiązania tego problemu polegała w pierwszej kolejności na analizie pola prędkości, które zostało obliczone dla obiektu Prendołówka przy zastosowaniu inwersji sejsmicznej. Na podstawie tego pola oraz przy wykorzystaniu danych interpretacyjnych stworzono prosty model biohermy górnojurajskiej zlokalizowanej w podłożu zapadliska przedkarpackiego. Analiza sekcji syntetycznej wykazała, że dotychczasowa interpretacja obiektu nie znalazła umotywowania w obrazie falowym. Na podstawie tnodeli sedymentologicznych analogów bio- herm pochodzących z odsłonięć zaproponowano architekturę budowli, która w nowym modelu potraktowana została jako kompleks obiektów biohermalnych. W wyniku modelowania otrzymano obraz, którego elementy znajdują potwierdzenie w danych rzeczywistych.

EN

Impedance inversion velocities were used to construct simple velocity model of a selected biohermal object. Then, basing on geological models for Upper Jurassic bioherms detailed architecture was proposed. The basic assumption was that we are not dealing with single object but with a cluster of build-ups. There are three intervals in the model reflecting the development of typical bioherms in this area: (a) firstrepresenting initial stage of separate- ("loose") bioherms, (b) main interval of massive ("rocky") limestones and (c) complex of shallow coral-algal-sponge reefs and bioherms. Several features on modeled image were found and compared with real migrated seismic data confirming correctness of some parts of the mode. Explaining internal heterogeneity of bioherms is crucial for treating them not as a single prospect but as a complex of drilling targets.

19

Praktyczne aspekty rozwiązania równania Zoeppritza w projektowaniu pomiarów sejsmiki 3C

Jędrzejowska-Tyczkowska H., Żuławiński K.

Prace Instytutu Górnictwa Naftowego i Gazownictwa

|

2002

|

nr 116 wyd. konferencyjne

147--150

PL

W pracy przedstawiono zastosowanie wyników algorytmu i programu obliczania współczynników odbicia w odniesieniu do projektowania schematu obserwacji w sejsmice 3C. Wykorzystano uogólniony model geologiczny z rejonu Kościana, gdzie przeprowadzone były prace eksperymentalne metodą odbitych fal przemiennych.

EN

The publication presents the results of application of Zoeppritz's equation solution for programming an optimal schemes of observation in 3C seismic.

20

Nowe techniki przetwarzania i interpretacji danych sejsmicznych w zastosowaniu do określenia szczegółowej budowy złoża ropy Radoszyn

Jędrzejowska-Tyczkowska H., Karpoluk Z., Misiarz P., Żukowska K., Żuławiński K., Bartoń R.

Przegląd Geologiczny

|

2001

|

Vol. 49 nr 10

858-866

PL

Artykuł dokumentuje wyniki bardzo szczegółowych studiów i analiz materiałów sejsmicznych i wybranych profilowań geofizycznych w otworach, ukierunkowanych na wyjaśnienie architektury facjalnej złoża ropy naftowej Radoszyn, a więc w pewnym sensie na udokumentowanie przydatności metody sejsmicznej do rozwiązywania bardzo szczegółowych problemów geologicznych. Pozytywne wyniki badań przypisać należy przede wszystkim poszerzonemu specjalistycznemu przetwarzaniu sekcji sejsmicznych, a szczególnie transformacji fazowej i standaryzacji widma częstotliwości oraz interaktywnie prowadzonemu modelowaniu sejsmicznemu, wspomagającemu optymalizację modelu geometryczno-prędkościowego.

EN

This paper presents the results of very detailed studies and analyses of seismic data and selected borehole loggings. The principal purpose of these studies was to explain the fades distribution pattern of the Radoszyn oil field. In some measure, thev were to confirm the usefulness of seismic methods for solving detailed geological problems. The positive results of these studies should be attributed primarily to the broadened specialist processing of seismic sections, especially to the phasic transformation and frequency spectrum standardization, as well as the interactive seismic modeling aiding the geometrical-velocity model optimization.