BazTech - Yadda

1

Ranking of the utility of selected geostatistical interpolation methods in conditions of highly skewed seismic data distributions: a case study of the Baltic Basin (Poland)

Sowińska-Botor Justyna, Mastej Wojciech, Maćkowski Tomasz

Gospodarka Surowcami Mineralnymi

|

2023

|

T. 39, z. 3

149--172

EN

The suitability of several low-labor geostatistical procedures in the interpolation of highly positively skewed seismic data distributions was tested in the Baltic Basin. These procedures were a combination of various estimators of the model of spatial variation (theoretical variogram) and kriging techniques, together with the initial data transformation to normal distribution or lack thereof. This transformation consisted of logarithmization or normalization using the anamorphosis technique. Two variations of the theoretical variogram estimator were used: the commonly used classical Matheron estimator and the inverse covariance estimator (InvCov), which is robust with regard to non-ergodic data. It was expected that the latter would also be resistant to strongly skewed data distributions. The kriging techniques used included the commonly used ordinary kriging, simple kriging useful for standardized data and the non-linear median indicator kriging technique. It was confirmed that normalization (anamorphosis) is the most useful and less laborious geostatistical procedure of those suitable for such data, which results in a standardized normal distribution. The second, not obvious statement for highly skewed data distributions suggests that the non-ergodic inverted covariance (InvCov) estimator of variogram has an advantage over the Matheron’s estimator. It gives a better assessment of the C0 (nugget effect) and C (sill) parameters of the spatial variability model. Such a conclusion can be drawn from the fact that the higher the estimation of the relative nugget effect L = C0/(C0 + C) using the InvCov estimator, the weaker the correlation between the kriging estimates and the observed values. The values of the coefficient L estimates obtained by using the Matheron’s estimator do not meet this expectation.

PL

W ramach studium przypadku w rejonie basenu bałtyckiego przetestowano przydatność kilku mało pracochłonnych procedur geostatystycznych do interpolacji silnie skośnych rozkładów danych sejsmicznych. Były one kombinacją różnych estymatorów modelu zmienności przestrzennej (wariogramu teoretycznego) i technik krigingu, wraz ze wstępną transformacją danych do rozkładu normalnego lub jej brakiem. Transformacja ta polegała na logarytmowaniu bądź na normalizacji z użyciem techniki anamorfozy. Zastosowano dwie odmiany estymatora wariogramu teoretycznego: powszechnie stosowany klasyczny estymator Matherona oraz estymator odwróconej kowariancji (InvCov) odporny na dane nieergodyczne. Spodziewano się, że ten drugi okaże się również odporny na silnie skośne rozkłady dane. Wśród zastosowanych technik krigingu znalazł się powszechnie stosowany kriging zwyczajny, kriging prosty użyteczny dla danych zestandaryzowanych i nieliniowa technika krigingu wskaźnikowego. Najbardziej użyteczną i mało pracochłonną procedurą geostatystyczną, nadającą się do zastosowania w przypadku takich danych, okazała się normalizacja (anamorfoza), w efekcie której uzyskuje się rozkład normalny standaryzowany. Drugim, nieoczywistym wnioskiem dla silnie skośnych rozkładów danych, jest sugestia, iż estymator InvCov ma przewagę nad estymatorem Matherona, ponieważ pozwala na bardziej realistyczną ocenę parametrów C0 (efektu samorodka) i C (wariancji progowej) modelu zmienności przestrzennej. Taki wniosek można wyciągnąć z faktu, że im wyższa wartość relatywnego efektu samorodków L = C0/(C0 + C) obliczona za pomocą estymatora InvCov, tym słabsza korelacja między wartościami obliczonymi a danymi. Wartości współczynnika L obliczone za pomocą estymatora Matherona nie posiadają tej właściwości.

2

Multichannel antileakage least squares spectral analysis for seismic data regularization beyond aliasing

Ghaderpour Ebrahim

Acta Geophysica

|

2019

|

Vol. 67, no. 5

1349--1363

EN

The antileakage least-squares spectral analysis is a new method of regularizing irregularly spaced data series. This method mitigates the spectral leakages in the least-squares spectrum caused by non-orthogonality of the sinusoidal basis functions on irregularly spaced series, and it is robust when data series are wide-sense stationary. An appropriate windowing technique can be applied to adapt this method to non-stationary data series. When data series present mild aliasing, this method can efectively regularize the data series; however, additional information or assumption is needed when the data series is coarsely sampled. In this paper, we show how to incorporate the spatial gradients of the data series into the method to regularize data series presenting severe aliasing and show its robust performance on synthetic and marine seismic data examples.

3

O nowych rozwiązaniach tektonicznych w „Atlasie geologicznym Polski”

Aleksandrowski P., Mazur S.

Przegląd Geologiczny

|

2017

|

Vol. 65, nr 12

1499--1510

EN

Authorial comprehensive comments and explanations are given to some of the interpretations applied in the tectonic part of the newly published Geological Atlas of Poland (Nawrocki, Becker, 2017) that considerably change the hitherto generally accepted concepts. It should be, however, admitted that most of those "new’" solutions were already proposed in the past by other workers as hypotheses that could not have been tested in the then state of knowledge on Poland’s deep geology and scientific tools at hand. This has now changed with abundant new data obtained with modern seismic techniques and advanced methods of potential field modelling. Using those data, we justify the reasons for, among others, a significant eastward shifting the front of the Variscan Orogen in Poland andfor the accompanying change in position of the division line between the Precambrian and Palaeozoic platforms. We also show the rationale for accepting a far-reaching southwestward extent of the East European Craton’s crystalline basement below the Palaeozoic Platform and for reinterpretation of the Teisseyre-Tornquist Zone’s nature, together with the question of early Palaeozoic terranes in the TESZ and the situation of the Caledonian foredeep at the SW margin of the East-European Craton.

4

Wykorzystanie metody Diving Wave Tomography do obliczenia wartości poprawek statycznych i rzędu dla zdjęcia sejsmicznego 3D z użyciem pakietu DW Tomography (Schlumberger)

Zając J.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

181--186

PL

Podstawowe metody obliczania wartości poprawek statycznych I rzędu dla zdjęć sejsmicznych 2D i 3D bazują na analizie przebiegu fali refrakcyjnej typu „head wave”. Metoda ta zakłada istnienie warstwowanego ośrodka o stałych prędkościach wewnątrzwarstwowych oraz generalnego wzrostu prędkości wraz z głębokością dla kolejnych warstw (refraktorów). Znacznie szerszą klasą ośrodków geologicznych stanowią ośrodki o gradientowym typie wzrostu prędkości z głębokością, dla których metoda „klasycznej refrakcji” okazuje się mniej efektywna. Metoda Diving Wave Tomography rozszerza zakres metod bazujących na analizie pierwszych impulsów (first break), uwzględniając tzw. ciągłą refrakcję promienia sejsmicznego do budowy modelu prędkościowo-głębokościowego ośrodka oraz do obliczenia wartości poprawek statycznych I rzędu. W publikacji przedstawiono wyniki zastosowania powyższej metody na przykładzie zdjęcia sejsmicznego 3D z użyciem pakietu DWTOMO.

EN

Basic methods of computing static corrections for 2D and 3D seismic data are based on analysis of the head wave refraction wave type. This method assumes the system of flat layers (refractors) with constant interval velocity and, in general, increasing velocity with the depth. To handle more commonly met near-surface geology with velocity gradients, Diving Wave Tomography should be used as the classical refraction method is less effective. The Diving Wave Tomography is based on the first arrivals analysis and tomography inversion using continuous refraction seismic rays to build depth-velocity model. In the paper results of the application of this method were presented on the example of seismic 3D data with using the DWTOMO package (Schlumberger).

5

Tłumienie zakłóceń koherentnych w warunkach nieregularnej akwizycji danych sejsmicznych realizowane przy użyciu oprogramowania NUCNS (Schlumberger)

Górowski M.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

231--235

PL

Przedstawiona metoda NUCNS (Non-Uniform Coherent Noise Suppresion, Schlumberger) pozwala na modelowanie zakłóceń koherentnych generowanych w trakcie wzbudzania sygnału sejsmicznego. Zakłócenia te interferują z sygnałem użytecznym i są trudne do usunięcia. Trudność ta może wystąpić w sytuacji, gdy prędkość pozorna fali zakłócającej jest zbliżona do prędkości fali użytecznej (np. fali refrakcyjnej) lub gdy z powodu nieregularnej pozycji punktów odbicia i wzbudzenia, dotychczas stosowane metody są nieefektywne. Zadaniem oprogramowania NUCNS jest skuteczne rozwiązanie wspomnianego problemu, również w przypadku nieregularnej akwizycji danych sejsmicznych. Przykłady zastosowania przedstawione dla zdjęć sejsmicznych z terenu południowej Polski pokazują efektywność metody.

EN

The presented method (NUCNS, Non-Uniform Coherent Noise Suppresion, Schlumberger) allows modeling of coherent noise generated during seismic acquisition. Coherent noise waves interfere with the useful signal and are difficult to remove. This difficulty may occur when the apparent velocity is similar to useful wave velocity (eg. refraction waves) or when due to irregularity of sources and detectors positions previously used methods are not effective. The purpose of the NUCNS program is effective problem solution both for regular and irregular se¬ismic data acquisitions. Examples of applications presented for seismic projects placed in southern Poland show the effectiveness of mentioned method.

6

Porównanie metodyki przetwarzania danych sejsmicznych szerokiego profilu i danych sejsmicznych 2D

Tlałka S.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

175--179

PL

Technologia przetwarzania danych sejsmicznych szerokiego profilu, wykorzystująca narzędzia 3D, pozwala osiągnąć lepsze wyniki w stosunku do konwencjonalnego przetwarzania pojedynczych linii 2D, przy nakładach kosztów akwizycyjnych zwykle niewiele przekraczających koszty standardowej linii 2D. Przedstawiono wybrane procedury przetwarzania szerokiego profilu, podejście w geometrii i schemacie przetwarzania. Pokazano wyniki porównań standardowego przetwarzania linii 2D do wyników z szerokiego profilu oraz podjęto próbę odpowiedzi na pytanie o wartość dodaną technologii sejsmiki szerokiego profilu w stosunku do technologii konwencjonalnej sejsmiki 2D.

XX

At an acquisition costs usually slightly exceeding the 2D budget, the design and shooting a seismic with the geometry of the wide line (WLP – Wide Line Profile), allows to achieve measurable benefits both at the stage of seismic data processing and at the interpretation stage. The wide line seismic data processing technology which uses 3D tools, allows to achieve better results compared to the conventional processing of independent 2D lines. Presented are the wide line processing procedures, and approach to geometry and processing scheme. The comparison between the standard processing of 2D lines and wide line processing is discussed. Additionally attempts to answer the question of added value from wide line technology with comparison to the conventional 2D seismic.

7

Określanie właściwości uszczelniających uskoków na podstawie modelowania z wykorzystaniem danych sejsmicznych i geofizyki otworowej – rozwiązania programowe opracowane w firmie Geofizyka Toruń S.A.

Rudzki M.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

197--200

PL

Referat przedstawia zarys metodyki numerycznego wyznaczania parametrów charakteryzujących własności uszczelniające powierzchni uskokowych, w oparciu o dane sejsmiczne i geofizyki wiertniczej. Przedstawiono ogólny schemat postępowania oraz wyniki zastosowania procedur numerycznych, opracowanych w firmie Geofizyka Toruń S.A.

EN

The paper presents an outline of the method of calculating the parameters characterising the seal properties of faults, basing on seismic and well-logging data. The general workflow, together with the results of numerical procedures developed in Geofizyka Toruń S.A. are presented.

8

Interpolacja i regularyzacja 5D z uwzględnieniem antyaliasingu przestrzennego

Gieras J.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

237--240

PL

Nowoczesne algorytmy regularyzacji i interpolacji danych sejsmicznych wykorzystując idee interpolacji 5D odtwarzają trasy poprzez iteracyjną rekonstrukcję w domenie Fouriera. Na odtworzenie trasy sejsmicznej składa się wygenerowanie parametrów zawartych w nagłówkach tras oraz samego sygnału z rozróżnieniem na kierunki i częstotliwości. W zależności od potrzeb dane mogą zostać odtworzone w żądanej domenie. Przykłady z różnych projektów demonstrują skuteczność metody.

EN

Modern algorithms of interpolation and regularization harnessing idea of 5D reconstruct traces by iterative steps in Fourier domain. To generate new seismic trace headers and seismic signal must be taken into account. This processes can be made for any domain depending on what we need. Below pictures demonstrate how algorithms works.

9

Interakcyjna sektoryzacja i selektywne sumowanie kolekcji pełnego azymutu – nowe narzędzie analizy obrazu strukturalnego

Kobusiński W.

Prace Naukowe Instytutu Nafty i Gazu

|

2016

|

nr 209 wyd. konferencyjne

193--196

PL

Najnowsze technologie obrazowania i analizy danych sejsmicznych w systemie EarthStudy360TM (ES360) odpowiadają na potrzebę lepszego zrozumienia często niejednoznacznych obrazów sejsmicznych. Nowe możliwości związane są z trójwymiarowymi, kątowymi kolekcjami CIG (common image gathers) zachowującymi pełną, rzeczywistą informację azymutalną. Kolekcje te umożliwiają m.in. analizę otrzymanych wyników w oparciu o ich interakcyjną sektoryzację i sumowanie. Prezentacja uwzględnia szczególnie sektoryzację azymutalną w oparciu o kolekcje kierunkowe 3D po migracji głębokościowej ES360 w domenie kątowej. Dyskutowane są aspekty teoretyczne oraz przykłady zastosowania tej technologii do danych rzeczywistych z obszaru Polski.

EN

The latest technologies of imaging and analysis of seismic data in the EarthStudy360TM (ES360) meet the need for better understanding of the often ambiguous seismic images. New opportunities are associated with three-dimensional, angle-domain CIG providing full subsurface azimuth information. These gathers allow, among others, detailed analysis of the results based on their on-the-fly sectorisation and selective stacking. The issue is presented with particular focus on azimuthal sectorisation based on the directional full azimuth angle gathers obtained with ES360 prestack depth migration. Discussed are theoretical aspects and examples of application of this technology to real seismic recorded in Poland.

10

Surface-based Internal Multiple Elimination in the CMP Domain − Theory and Application Strategies on Land Seismic Data

Deng S., Zhao W., Liu Z.

Acta Geophysica

|

2016

|

Vol. 64, no. 6

2114--2135

EN

The data-driven internal multiple elimination (IME) method based on feedback model, which includes CFP-based, surface-based and inversion-based methods, are successfully applied to marine datasets. However, these methods are computationally expensive and not always straightforward on land datasets. In this paper, we first proved that the surface-based IME method, which is the most computationally efficient method among the three methods, can be derived from the CFP theory. Then we extend it to CMP domain under the assumption of locally lateral invariance of the earth, which makes it more computationally efficient. In addition, we proposed applying a time-variant taper based on the first Fresnel zone to predict the multiples more percisely. Besides, the improved S/N ratio and dense offset distribution can be obtained by using the CMP supergather, which makes the CMP-oriented method more suitable for land data. Some practical processing strategies are proposed via case study. The effectiveness of the proposed method is demonstrated with the application to synthetic and field data.

11

Applications of Savitzky-Golay Filter for Seismic Random Noise Reduction

Liu Y., Dang B., Li Y., Lin H., Ma H.

Acta Geophysica

|

2016

|

Vol. 64, no. 1

101--124

EN

This article utilizes Savitzky–Golay (SG) filter to eliminate seismic random noise. This is a novel method for seismic random noise reduction in which SG filter adopts piecewise weighted polynomial via leastsquares estimation. Therefore, effective smoothing is achieved in extracting the original signal from noise environment while retaining the shape of the signal as close as possible to the original one. Although there are lots of classical methods such as Wiener filtering and wavelet denoising applied to eliminate seismic random noise, the SG filter outperforms them in approximating the true signal. SG filter will obtain a good tradeoff in waveform smoothing and valid signal preservation under suitable conditions. These are the appropriate window size and the polynomial degree. Through examples from synthetic seismic signals and field seismic data, we demonstrate the good performance of SG filter by comparing it with the Wiener filtering and wavelet denoising methods.

12

Anizotropia ośrodka geologicznego oraz jej wpływ na zarejestrowane dane sejsmiczne

Leginowicz A.

Nafta-Gaz

|

2014

|

R. 70, nr 9

584--598

PL

W publikacji zawarto obszerne teoretyczne studium dotyczące zagadnień anizotropii ośrodka oraz jej wpływu na zarejestrowane dane sejsmiczne. Przedstawiono teoretyczne podstawy metody wykorzystującej zależność zmienności amplitudy z offsetem oraz azymutem. Omówiono przykłady z literatury opisujące metody analizowania parametrów anizotropowych na podstawie danych sejsmicznych, a w szczególności — metody prędkościowe i amplitudowe (AVOAz) używane w celu scharakteryzowania spękań ośrodka anizotropowego.

EN

Anisotropy connected with fractures is an important issue in the exploration of unconventional reservoirs. This paper presents anisotropic study and its influence on seismic data. Amplitude versus offset and azimuth theoretical bases were presented and bibliography examples, such as velocity and amplitude method for extracting fracture information from seismic data were discussed.

13

Strategie przetwarzania danych sejsmicznych w poszukiwaniu i rozpoznawaniu „złóż" gazu w łupkach

Orzeszko M., Sinoracki A., Gromotka J.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

317--324

PL

Autorzy posteru przedstawią dwie strategie przetwarzania danych sejsmicznych, mających zastosowanie w procesie poszukiwań złóż gazu w łupkach. Pierwsza spośród tych strategii obejmuje przetwarzanie danych w dziedzinie czasu. Jej celem jest przygotowanie materiału do interpretacji strukturalnej, litofacjalnej oraz analiza anizotropii pionowej i poziomej. W przetwarzaniu danych wykorzystano opracowaną i rozwijaną od 2008 roku technologię badania anizotropii; dla aktualnych potrzeb zbudowano pakiet technologiczny do przetwarzania danych w poszukiwaniach gazu w łupkach. Ważniejsze elementy pakietu będą komentowane, skomentowane będą również wymagania sprzętowe (potrzebna moc obliczeniowa). Strategia druga obejmuje przetwarzanie w dziedzinie głębokości i zostanie omówiona skrótowo, ponieważ stanowi przedmiot odrębnego opracowania. Komentowane będą założenia oprogramowania, określone zostaną wymagania adresowane do projektantów i wykonawców akwizycji.

EN

The authors of the poster are about to outline two possible strategies of seismic data processing using in shale gas research. First approach base on data processing in time domain. The objective of this method is to prepare input data for structural and litofacial interpretation, as we1l as vertical and horizontal anisotropy analysis. The technology of data processing with regard to anisotropy, designed and developed since 2008, has been used with several updates based on current knowledge and technical possibilities. A uniform package of technology for seismic data processing has been built. Its basic elements will be presented. The second strategy base on processing in depth domain and will be only briefly covered as it is a topic of separate presentation. Aspects of software design, acquisition and equipment requirements will be outlined.

14

Przetwarzanie danych sejsmiki w poszukiwaniach złóż niekonwencjonalnych : ilustracje kluczowych elementów pakietu łupkowego

Kobusiński W., Podolak M.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

433--437

PL

Poszukiwanie gazu ziemnego w skałach łupkowych spowodowało potrzebę określania anizotropii azymutalnej ośrodka oraz znacznego zwiększenia precyzji informacji geologicznych uzyskiwanych z danych sejsmicznych. To z kolei wymaga zastosowania wielu nowych algorytmów przetwarzania danych sejsmicznych. Kluczowe techniki oraz działanie realizujących je modułów przetwarzania ilustrowane są przykładami z praktycznego ich stosowania. Prezentacja obejmuje: głębokościową migrację przed składaniem według wspólnego kąta odbicia (CRAM), korektę efektu cieni przyuskokowych, metodę konwersji atrybutów efektywnych w atrybuty interwałowe, rozkład sumy WPG na składową refleksyjną i dyfrakcyjną.

EN

Search for the shale gas brought need to define azimuthal anisotropy of rock medium, and demand for essential increase of precision expected from seismic data. That, in turn, created need of many new algorithms for seismic data processing. Key techniques, and ideas behind their implementation, are illustrated with examples of practical application. Presented are: prestack depth migration in common reflection angle domain (CRAM), correction for fault shadows, method of converting effective attributes into interval models, decomposition of classical CMP stack into components of reflectivity and diffractivity.

15

Przetwarzanie danych sejsmicznych z wykorzystaniem nowej technologii przetwarzania, poprawiającej koherencję sygnału sejsmicznego

Masiukiewicz D., Lubomski S.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

289--294

PL

Podstawowym celem przetwarzania danych sejsmicznych jest uzyskanie jak najwyższej jakości sekcji sejsmicznych. Zastosowanie standardowego schematu przetwarzania danych nie zawsze jest wystarczająco skuteczne, szczególnie w rejonach o niskim stosunku sygnału do zakłóceń. W Ośrodku Obliczeniowym Geofizyki Toruń wdrożono własną technologię poprawiającą koherencję sygnału ECP (Enhanced Coherency Processing) możliwą do zastosowania przed składaniem, po składaniu, dla 2D i 3D. W niniejszej pracy zaprezentowane zostały przykłady jej działania na danych sejsmicznych syntetycznych i rzeczywistych.

EN

The basic aim and goal of seismic data processing is to obtain the highest quality seismic section. Using a standard schema of data processing is not always sufficiently effective, especially in regions with low signal to noise ratio. In Seismic Data Processing Department of Geofzyka Toruń proprietary technology was implemented which improves the signal coherence ECP (Enhanced Coherency Processing). ECP can be used before stacking, as well as after stacking, for 2D and 3D. In this paper examples of obtained results of ECP were presented on synthetic and real seismic data.

16

Przestrzenne modelowanie parametrów petrofizycznych w oparciu o atrybuty sejsmiczne, rejon zdjęcia sejsmicznego Góra Ropczycka-Iwierzyce 3D

Zych I., Smółka-Gnutek P.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

233--238

PL

Celem przestrzennego modelowania było dostarczenie szczegółowej informacji o rozkładzie parametrów petrofizycznych w kompleksie utworów miocenu. Obszar objęty analizą znajduje się w brzeżnej strefie nasunięcia karpackiego — w strefie znanych złóż gazu ziemnego — Sędziszów (Zagórzyce), Nosówka oraz Góra Ropczycka. W wyniku badań zweryfikowano dotychczasowe rozpoznanie złóż i wyznaczono kierunki dalszych prac poszukiwawczych. Główną zaletą przedstawionej analizy jest integracja szczegółowych danych geologicznych i petrofizycznych (profilowania geofizyki wiertniczej, wyniki prób złożowych, pomiary laboratoryjne) z danymi sejsmicznymi (m.in. atrybuty sejsmiczne, impedancja akustyczna). W wyniku analizy otrzymano przestrzenne modele rozkładu porowatości PHI, gęstości RHOB oraz wielkości naturalnego promieniowania GR. Na podstawie wyników interpretacji prac sejsmicznych 3D Góra Ropczycka—Iwierzyce 3D wykonano przestrzenny model strukturalny kompleksu miocenu o pionowej rozdzielczości 10 m. Dostępne dane wiertnicze w postaci krzywych geofizycznych przeskalowano do rozdzielczości przestrzennego modelu strukturalnego. Do kierowania rozkładem analizowanych danych petrofizycznych wykorzystano przede wszystkim informacje z szeregu atrybutów sejsmicznych. Estymację wykonano za pomocą zaawansowanych algorytmów dostępnych w programie Petrel (Neural Net, Gaussian Random Function Simulation). W wyniku badań uzyskano informację o rozkładzie własności petrofizycznych w każdym punkcie przestrzennego modelu utworów kompleksu miocenu. Na podstawie dostępnej informacji złożowej z przestrzeni modelu wyselekcjonowano także anomalne obszary korelujące się z objawami węglowodorów w otworach wiertniczych. Pozwoliło to na weryfikację dotychczasowego rozpoznania złóż gazu Sędziszów, Nosówka i Góra Ropczycka oraz wskazanie nowych stref o podobnym zapisie sejsmicznym.

EN

The purpose of spatial modeling was to provide detailed information about the distribution of petrophysical parameters in the Miocene complex. The area of study is located in the marginal zone of the Carpathian Foredeep where many deposits of natural gas such as Sędziszów (Zagórzyce), Nosówka and Góra Ropczycka are known. The tests verified existing deposits and set directions for further development. The main advantage of the presented analysis is the integration of detailed geological and petrophysical data (borehole geophysics, laboratory measurements) with the seismic data (seismic attributes, acoustic impedance, etc.). As the result of the analysis spatial models were obtained such as porosity distribution PHI, density RHOB and gamma ray distribution GR. On the basis of structural interpretation of 3D seismic data, the structure framework of Miocene complex was built with the vertical resolution of 10 m. All well data was scaled to the resolution of structural model. A number of seismic attributes were used to steer the distribution of petrophysical well data. Estimation was performed by using advanced algorithms available in Petrel (Neural Net, Gaussian Random Function Simulation). The research provided information about the distribution of petrophysical properties at every point of the Miocene complex. On the basis of all available reservoir information, anomalous areas were extracted. Thanks to this analysis existing gas fields: Sędziszów, Nosówka, Góra Ropczycka could be verified and new areas with similar seismic characteristic were identified.

17

Praktyka zwiększania pionowej rozdzielczości zapisu sejsmiki powierzchniowej

Kowalski H., Podolak M.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

65--67

PL

Zwiększanie pionowej rozdzielczości zapisu sejsmicznego, pomimo długiej historii doskonalenia metod przetwarzania danych sejsmicznych, pozostaje jednym z najważniejszych zadań do rozwiązania. Prezentowane są analizy porównawcze paru metod dostępnych obecnie w warsztacie geofizyka przetwarzającego przemysłowo dane sejsmiczne. Prezentowana jest również nowa metoda, od nazwiska autora nazwana HJT, stwarzająca nowe możliwości w porównaniu do technik dotychczas stosowanych.

EN

Inspite of long-time record of developments in seismic data processing, attempts to enhance vertical resolution of recorded seismic data remains one of tasks to be solved. Presented are results of comparative analysis of several methods available today in industrial-level workshop of seismic data processor. Moreover, the new method is being presented, called HJT from the name of its developer. That seems to open new opportunities when comparing it to the other available methods.

18

Odtworzenie pożądanej rozdzielczości danych sejsmicznych w aspekcie realizowanych zadań geologicznych

Jędrzejowska-Tyczkowska H.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

57--63

PL

W publikacji przedstawiono wyniki zastosowania opracowanego w INiG sposobu zwiększenia rozdzielczości sejsmicznego pola falowego (patent nr 394446). Efekty realizacji pokazano na dwóch różnych obiektach (z USA i Polski), w odniesieniu do dwóch różnych zadań geologicznych. Celem prezentowanego rozwiązania jest: a) zwiększenie rozdzielczości zapisywanych krzywych pomiarowych, co pozwoli odkrywać znacznie mniejsze obiekty geologiczne; b) eliminacja efektu zaniku energii drgań sprężystych, proporcjonalnego w realiach propagacji drgań do długości przebytej drogi, co bezpośrednio wpływa na zwiększenie głębokości propagacji fali sejsmicznej. Dotychczas stosowane sposoby, wykorzystując zabiegi o nazwach filtracja (w domenie częstotliwości) oraz dekonwolucja (w domenie czasu), bazowały wyłącznie na zapisach drgań sprężystych dokonywanych na powierzchni ziemi. Zaproponowany sposób rozpoznania ośrodka geologicznego wykorzystuje dodatkowe źródło informacji w postaci krzywych pomiarowych zapisanych w otworach wiertniczych, do określenia niezbędnego zakresu częstotliwości umożliwiających dokładną analizę i opis ośrodka. Aktualnie, niezwykle wygórowane wymagania i wyzwania dotyczące efektywności metody sejsmicznej, związane z nowymi geologicznymi celami poszukiwań (np. złoża niekonwencjonalne typu shale gas i tight gas) spowodowały, iż znane, stosowane dotąd metody przetwarzania sygnałów nie spełniają oczekiwanej roli.

EN

In publication theoretical and practical results of seismic data resolution improvement are presented (Patent No. 394446). Various techniques and technologies have been implemented in geological exploration. One technique commonly referred to as the "seismic method", includes generating vibrations in the earth and recording reflected signals. Generally, attempts to improve the seismic methods for geological exploration have been aimed at increasing the resolution of recorded measuring curves. In particular, filtration (in the frequency domain) and deconvolution (in the time domain) were used and based exclusively on the recordings of elastic vibrations made on the surface of the Earth. Presently, strict requirements and challenges related to the efficiency of the seismic method with new geological targets of exploration (for example, unconventional deposits of shale gas and tight gas) have exposed previously used methods of signal processing/conversion as less effective at locating such targets.

19

Możliwości zobrazowania struktur podcechsztyńskich przy wykorzystaniu procedur przetwarzania danych sejsmicznych na obszarach Niżu Środkowoeuropejskiego

Rabiasz M., Hyjek P., Sikorski Ł.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

325--330

PL

Eksploatacja gazu łupkowego w Polsce może stać się zasadniczym źródłem tego surowca w niedalekiej przyszłości. Istnienie bogatych złóż gazu łupkowego w Polsce nie jest sprawą nową, jednak brak technologii wydobywczej, jak i trudności w pozyskaniu dobrej jakości zapisu sejsmicznego głęboko leżących formacji geologicznych oraz ich interpretacja geologiczna spowodowały, że dopiero teraz złoża te mogą stać się interesujące pod względem poszukiwawczym. Na jakość obrazowania sejsmicznego ma wpływ niejednorodność budowy przypowierzchniowej i związane z tym odpowiednie podejście do powierzchniowego rozwiązania poprawek statycznych. Niemniej ważne jest odpowiednie dobranie procedur processingowych, które mają za zadanie jak najlepiej oczyścić materiał z fal zakłócających oraz wyeliminować odbicia wielokrotne skutecznie maskujące wgłębną budowę geologiczną. Celem artykułu jest przedstawienie wyników wybranych sekwencji przetwarzania danych sejsmicznych, poczynając od estymacji powierzchniowo zgodnych poprawek statycznych poprzez odpowiedni dobór procedur processingowych, mających na celu uzyskanie jak najlepszej jakości obrazu sejsmicznego głęboko zalegających, perspektywicznych formacji paleozoicznych.

EN

Drilling for shale gas in Poland may well be a major source of this resource in the near future. Rich layers of gas-bearing shale have been around for ages. Until recently this resource was not of interest because of lack of production technology and considerable difficulties in getting good quality seismic imaging and geological interpretation of that data. This is now changing. The quality of the seismic imaging available is impacted by the heterogeneity of the near-surface structure and the approach chosen to deal with statics. Choosing the proper processing procedures to best reduce interference and eliminate multiple waves is equally important in order to show the underlying geological structure. This article aims to present the results of select seismic data processing sequences. Beginning with estimation of surface consistent long-wave statics, through choosing the right processing procedures aimed at producing the best quality image of deep and prospective Paleozoic formations.

20

Metoda szacowania pola prędkości ośrodka geologicznego na podstawie danych sejsmicznych z wykorzystaniem algorytmu genetycznego

Pirowska K.

Prace Naukowe Instytutu Nafty i Gazu

|

2012

|

nr 182 wyd. konferencyjne

277--281

PL

Zaprezentowana została bazująca na algorytmie genetycznym metoda szacowania prędkości propagacji fal podłużnych w izotropowym ośrodku geologicznym na podstawie danych sejsmicznych. W oparciu o wyznaczone z rekordów polowych przebiegi czasowe poszukiwane były parametry dotyczące modelu prędkościowo-głębokościowego. Omówiono rezultaty przetestowania metody na wybranym profilu rzeczywistych danych sejsmicznych.

EN

The method presented in this paper enables estimating velocity of P-waves propagation in isotropic geological medium from seismic survey data. Estimation is conducting with application of genetic algorithm. Parameters which characterize velocity-depth model are searched based on hodographs picked from field data. The technique was tested on real profile of seismic data. Results and obserwations are shown.