Papers by Aleksander Wilk
Nafta-Gaz, 2017
Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześ... more Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześniejszych opracowań, który pozwoliłby na bardziej szczegółowe odwzorowanie skomplikowanej budowy geologicznej z przedmiotowego rejonu badań, poprzez opracowanie procedur i sekwencji przetwarzania oraz dobór parametrów, z wykorzystaniem wszelkich dostępnych danych i jednoczesną interpretacją geologiczną. Realizując ten cel, dużą ilość czasu poświęcono na wykonywanie testów oraz dobór odpowiednich procedur przetwarzania, jak również ułożenia ich w odpowiednią sekwencję, przynoszącą możliwie najlepsze efekty. Po każdym etapie prac wykonywano sumę kontrolną, która poddawana była szczegółowej analizie i weryfikacji uzyskanego obrazu, przy zastosowaniu różnych parametrów. Prace przeprowadzono na profilu wytypowanym ze względu na stosunkowo najsłabszy zapis sejsmiczny w porównaniu z innymi dostępnymi profilami z tego rejonu. Słowa kluczowe: przetwarzanie sejsmiczne, poprawki statyczne, prędkość, migracja, interpretacja. Improvement of the accuracy of the geological survey imaging in seismic waveform in the Outer Carpathians area The aim of this study was to obtain a more accurate seismic image in relation to previous studies, which would allow detailed mapping of a complex geological structure of the analyzed area through the development of procedures and a processing sequence, selection of parameters, using all available data and simultaneous geological interpretation. By accomplishing this goal, a great deal of time has been devoted to performing tests and the selection of appropriate processing procedures, as well as arranging them into the appropriate sequences with the best possible results. After each stage of the work, a checksum was performed, which was subjected to detailed analysis and verification of the obtained image, using different parameters. The work was carried out on a profile which was selected due to the relatively weakest seismic data in comparison with other available profiles from this region.
Nafta Gaz, Feb 1, 2016
Do obliczenia gradientu ciśnień porowych wykorzystano moduł Seismic Pore Pressure Modeling firmy ... more Do obliczenia gradientu ciśnień porowych wykorzystano moduł Seismic Pore Pressure Modeling firmy Western-Geco, który został udostępniony autorowi na okres jednego miesiąca. Atutem wspomnianego modułu jest jego interaktywność oraz fakt, iż pracuje on w systemie Petrel firmy Schlumberger, będącym jednym z podstawowych narzędzi interpretacyjnych używanych przez polskie jednostki przemysłowe. Do obliczeń wykorzystano rzeczywiste pomiary ze zdjęcia sejsmicznego 3D z północnego obszaru Polski. Uzyskano zadowalające wyniki potwierdzające użyteczność tego narzędzia w pracach poszukiwawczych. Słowa kluczowe: ciśnienie porowe, gradient ciśnień porowych, naprężenie efektywne, model kompakcji, profilowanie akustyczne, gęstość, prędkość. Predicting pore pressure in the geological medium, based on the module "Seismic Pore Pressure Modeling" of WesternGeco Pore pressure gradient calculations were based on the module "Pore Pressure Seismic Modeling" of WesternGeco, which was made available to the author for a period of one month. The advantage of this module is its interactivity, and the fact that it works in the Schlumberger Petrel system, which is one of the basic tools of interpretation used by our native oil industry. For calculations, actual measurements of the 3D seismic survey of the northern Polish territory were used. Satisfactory results confirming the usefulness of this tool in exploration work was achieved.
Prace Naukowe Instytutu Nafty i Gazu, Jun 20, 2018
Nafta-Gaz, 2021
Głównym celem artykułu jest przedstawienie charakterystyki sejsmicznej utworów o charakterze olis... more Głównym celem artykułu jest przedstawienie charakterystyki sejsmicznej utworów o charakterze olistostromowym z obszaru paleodoliny Szczurowej–Wojnicza na podstawie interpretacji połączonych danych sejsmicznych 3D i 2D. Interpretowany rejon cechuje się skomplikowaną budową geologiczną oraz słabym stopniem rozpoznania wiertniczego, gdyż tylko nieliczne odwierty przewierciły pełen profil utworów miocenu autochtonicznego. Paleodolina Szczurowej–Wojnicza, stanowiąca jeden z najbardziej charakterystycznych elementów budowy strukturalnej obszaru badań, wyerodowana została w utworach węglanowych kredy i jury. W okresie neogenu paleodolina ta wypełniona została miąższą serią osadów klastycznych należących do formacji skawińskiej. W podłożu tej formacji zidentyfikowany został również pakiet charakterystycznych osadów gruboklastycznych o silnie zróżnicowanej miąższości, nazwany serią zlepieńcowo-olistostromową. Punktem wyjścia interpretacji był profil otworu Basowy-1, w którym w materiale rdze...
Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześ... more Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześniejszych opracowań, który pozwoliłby na bardziej szczegółowe odwzorowanie skomplikowanej budowy geologicznej z przedmiotowego rejonu badań, poprzez opracowanie procedur i sekwencji przetwarzania oraz dobór parametrów, z wykorzystaniem wszelkich dostępnych danych i jednoczesną interpretacją geologiczną. Realizując ten cel, dużą ilość czasu poświęcono na wykonywanie testów oraz dobór odpowiednich procedur przetwarzania, jak również ułożenia ich w odpowiednią sekwencję, przynoszącą możliwie najlepsze efekty. Po każdym etapie prac wykonywano sumę kontrolną, która poddawana była szczegółowej analizie i weryfikacji uzyskanego obrazu, przy zastosowaniu różnych parametrów. Prace przeprowadzono na profilu wytypowanym ze względu na stosunkowo najsłabszy zapis sejsmiczny w porównaniu z innymi dostępnymi profilami z tego rejonu.
This article presents the results of 2D seismic reprocessing in the West Pomerania region. The pu... more This article presents the results of 2D seismic reprocessing in the West Pomerania region. The purpose of reprocessing was to improve the imaging of under-Zechstein formations and structures. The current results of seismic processing in this area do not allow for their structural (and more facial) interpretation, despite the efforts undertaken to this end. Reprocessing was carried out in the poststack migration version based on a processing sequence developed at the Seismic Department of the Oil and Gas Institute – National Research Institute. Stacking velocities and a velocity model based on VSP measurements were used for poststack migration. The seismic image on the reprocessed profile has a better projection of the geological structure in relation to the previous study. The seismic image obtained as a result of using the velocity model to poststack migration based on VSP measurements appears to be better than the seismic image obtained as a result of using the stacking velocities...
Do obliczenia gradientu ciśnien porowych wykorzystano modul Seismic Pore Pressure Modeling firmy ... more Do obliczenia gradientu ciśnien porowych wykorzystano modul Seismic Pore Pressure Modeling firmy WesternGeco, ktory zostal udostepniony autorowi na okres jednego miesiąca. Atutem wspomnianego modulu jest jego interaktywnośc oraz fakt, iz pracuje on w systemie Petrel firmy Schlumberger, bedącym jednym z podstawowych na rzedzi interpretacyjnych uzywanych przez polskie jednostki przemyslowe. Do obliczen wykorzystano rzeczywiste pomiary ze zdjecia sejsmicznego 3D z polnocnego obszaru Polski. Uzyskano zadowalające wyniki potwierdzające uzytecznośc tego narzedzia w pracach poszukiwawczych. Slowa kluczowe: ciśnienie porowe, gradient ciśnien porowych, naprezenie efektywne, model kompakcji, profilowanie akustyczne, gestośc, predkośc. Predicting pore pressure in the geological medium, based on the module “Seismic Pore Pressure Modeling” of WesternGeco Pore pressure gradient calculations were based on the module “Pore Pressure Seismic Modeling” of WesternGeco, which was made available to the a...
Nafta-Gaz, 2018
W niniejszym artykule zaprezentowano sposób konstrukcji pola prędkości na potrzeby migracji czaso... more W niniejszym artykule zaprezentowano sposób konstrukcji pola prędkości na potrzeby migracji czasowej 2D po składaniu w trudnych rejonach geologicznych na przykładzie Karpat fliszowych w południowo-wschodniej Polsce. Rejon badań charakteryzuje się dużym stopniem skomplikowania budowy geologicznej, co przekłada się na znaczną trudność w jego odwzorowaniu na sekcjach sejsmicznych. Określenie poprawnego pola prędkości do procedury migracji pozwala prawidłowo odwzorować wgłębną budowę geologiczną na przekroju sejsmicznym. W wyniku migracji opartej na prawidłowym rozpoznaniu rozkładu prędkości uzyskuje się rzeczywiste położenie punktów odbicia od granic nachylonych, usunięcie dyfrakcji, znaczną poprawę rozdzielczości przestrzennej, a zwłaszcza rozdzielczości poziomej analizowanego obrazu sejsmicznego. Nowatorskim rozwiązaniem konstrukcji budowy pola prędkości na potrzeby migracji czasowej 2D po składaniu było wykorzystanie prędkości średnich z pomiarów PPS (pionowe profilowanie sejsmiczne), które ze względu na metodykę pomiaru zawierają informację o anizotropii ośrodka geologicznego. Istotnym elementem w konstrukcji modelu prędkości było zdefiniowanie optymalnego rozkładu prędkości średnich, zarejestrowanych w lokalnych pozycjach otworów wiertniczych z offsetowych pomiarów PPS. Uwzględnienie efektu anizotropii pozwoliło na wiarygodniejszy rozkład pola prędkości i uzyskanie polepszenia obrazu falowego w stosunku do wcześniejszych opracowań. Otrzymany model prędkości stanowił podstawę do odtworzenia skomplikowanej budowy ośrodka geologicznego. W przyjętej przestrzeni obliczeniowej rejonu badań rozpatrywane były dwa modele: model płasko-równoległy bez interpretacji strukturalnej oraz model z interpretacją strukturalną. Dla przyjętych modeli prędkości średnie uzyskane z PPS zostały interpolowane i ekstrapolowane przy użyciu trzech algorytmów: rozkładu Gaussa, krigingu i moving average w systemie Petrel firmy Schlumberger. Na podstawie przetestowanych modeli prędkości dla wybranego profilu sejsmicznego stwierdzono, że optymalny wynik uzyskano w przypadku rozkładu Gaussa z wykorzystaniem modelu z interpretacją strukturalną. Zastosowanie modelu do migracji czasowej 2D po składaniu uwzględniającego anizotropię ośrodka dostarcza bardziej wiarygodnego obrazu ośrodka geologicznego w stosunku do dotychczasowych opracowań, co powinno przekładać się na zwiększenie efektywności w poszukiwaniach węglowodorów oraz ograniczać stopień ryzyka poszukiwawczego.
Nafta-Gaz, 2016
Wyznaczanie rozkładu szczelin w formacjach typu shale gas na podstawie azymutalnych pomiarów sejs... more Wyznaczanie rozkładu szczelin w formacjach typu shale gas na podstawie azymutalnych pomiarów sejsmicznych W artykule została zaprezentowana metodyka wyznaczania rozkładu szczelin i ich azymutów dla dwóch horyzontów sejsmicznych odpowiadających stropowi ogniwa Jantaru i formacji z Sasina, tj. poziomów zbiornikowych typu shale gas, na podstawie azymutalnych pomiarów sejsmicznych 3D. Fala sejsmiczna propagująca przez ośrodek geologiczny uzyskuje wyższą prędkość, przechodząc wzdłuż systemu spękań, a niższą prostopadle do nich. Wykorzystanie tego zjawiska umożliwia przeprowadzenie analiz zmian amplitud względem różnych kierunków (azymutów) zapisu sejsmicznego, pozwalając na określenie stopnia zmienności parametrów sprężystych w układzie HTI (horizontal transverse isotropy). Wyznaczony w ten sposób stopień zmienności umożliwia określenie wielkości oraz kierunków rozkładu szczelin rozpatrywanego ośrodka geologicznego. Stopień rozpoznania intensywności oraz orientacji kierunków spękań (szczelin) jest szczególnie ważny w projektowaniu otworów poziomych oraz w optymalizacji procesów eksploatacji węglowodorów ze złoża. Słowa kluczowe: sejsmika azymutalna, aproksymacja, szczelinowatość, shale gas. Determining distribution of fractures in shale gas formations based on azimuthal seismic survey This article presents results of a methodology of determining the distribution of fractures and their azimuths for two seismic horizons corresponding to the tops of Jantar member and Sasino formation, i.e. the levels of shale gas type reservoir, based on azimuth 3D seismic data. Seismic wave propagating through the geological medium achieved higher velocity going along the system fracture and a lower perpendicular to them. Use of this phenomenon makes it possible to carry out analyzes of changes in amplitude, with respect to different directions (azimuths) of seismic record, allowing to determine the degree of variability of elastic parameters in a HTI (Horizontal Transverse Isotropy). The specified degree of variation makes it possible, to determine the size and direction of fractures distribution in the geological medium. The degree of the intensity and orientation directions of cracks (fractures), is particularly important in the design of horizontal wells and optimizing the exploitation processes of hydrocarbons from the reservoir.
Aplikácia nových technológií interpretácie seizmických dát za účelom prehĺbenia efektívnosti skúš... more Aplikácia nových technológií interpretácie seizmických dát za účelom prehĺbenia efektívnosti skúšok The presentation shows an application of new techniques in the structural and litho-facies interpretation of seismic data. The newest: Halliburton Digital and Consulting Solutions-Landmark and Hampson Russell Company softwares were used.
Nafta-Gaz, 2018
Presentation of the of seismic image analysis results obtained by reprocessing two seismic profil... more Presentation of the of seismic image analysis results obtained by reprocessing two seismic profiles is the main aim of the presented work. The profiles are located in the marginal part of the Outer Carpathians. The mentioned profiles were reprocessed in the Seismic Department of the Oil and Gas Institute-National Research Institute in Krakow, Poland. Proper selection of both the processing sequence and parameters, as well as verification of each stage of processing by simultaneous geological interpretation, resulted in a partly different mapping of the geological structures in comparison with the previous stage. Structural interpretation based on the obtained seismic imagery provides new information that could be used for more thorough interpretation of the Outer Carpathians tectonic units, as well as detailed reconstruction of the fault zones in the analysed area.
Nafta-Gaz
This article presents a construction method of the velocity field for poststack time migration fo... more This article presents a construction method of the velocity field for poststack time migration for 2D seismic calculated on the basis of interval velocities in boreholes and structural interpretation, as well as the results of poststack time migration based on this solution. Three velocity field models have been developed. The models used differ in the way of spatial interpolation and extrapolation in the adopted calculation grid in the depth domain, which was created on the basis of a structural interpretation of 2D seismic profiles. Three methods of interpolation and extrapolation were used: Gaussian distribution, kriging and moving average. The spatial distribution of the interval velocities in the boreholes was made using the Petrel software by Schlumberger. The interval velocities along the analyzed seismic profile were extracted from the computed spatial interval velocity models, and after conversion from the depth to the time domain, they were used for the poststack time migr...
Nafta-Gaz, Jul 1, 2017
Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześ... more Celem przedstawionej pracy było uzyskanie dokładniejszego obrazu sejsmicznego w stosunku do wcześniejszych opracowań, który pozwoliłby na bardziej szczegółowe odwzorowanie skomplikowanej budowy geologicznej z przedmiotowego rejonu badań, poprzez opracowanie procedur i sekwencji przetwarzania oraz dobór parametrów, z wykorzystaniem wszelkich dostępnych danych i jednoczesną interpretacją geologiczną. Realizując ten cel, dużą ilość czasu poświęcono na wykonywanie testów oraz dobór odpowiednich procedur przetwarzania, jak również ułożenia ich w odpowiednią sekwencję, przynoszącą możliwie najlepsze efekty. Po każdym etapie prac wykonywano sumę kontrolną, która poddawana była szczegółowej analizie i weryfikacji uzyskanego obrazu, przy zastosowaniu różnych parametrów. Prace przeprowadzono na profilu wytypowanym ze względu na stosunkowo najsłabszy zapis sejsmiczny w porównaniu z innymi dostępnymi profilami z tego rejonu. Słowa kluczowe: przetwarzanie sejsmiczne, poprawki statyczne, prędkość, migracja, interpretacja. Improvement of the accuracy of the geological survey imaging in seismic waveform in the Outer Carpathians area The aim of this study was to obtain a more accurate seismic image in relation to previous studies, which would allow detailed mapping of a complex geological structure of the analyzed area through the development of procedures and a processing sequence, selection of parameters, using all available data and simultaneous geological interpretation. By accomplishing this goal, a great deal of time has been devoted to performing tests and the selection of appropriate processing procedures, as well as arranging them into the appropriate sequences with the best possible results. After each stage of the work, a checksum was performed, which was subjected to detailed analysis and verification of the obtained image, using different parameters. The work was carried out on a profile which was selected due to the relatively weakest seismic data in comparison with other available profiles from this region.
Nafta-Gaz
Obliczanie objętości materii organicznej (TOC) na podstawie inwersji genetycznej danych sejsmiczn... more Obliczanie objętości materii organicznej (TOC) na podstawie inwersji genetycznej danych sejsmicznych W artykule została zaprezentowana metodyka wyznaczania rozkładu parametru TOC na danych sejsmicznych 3D na podstawie inwersji genetycznej. Zaprezentowany schemat obliczeń zastosowano na danych sejsmicznych pochodzących z obszaru północnej Polski. Głównym celem przeprowadzonych badań było rozpoznanie dwóch, wyinterpretowanych horyzontów sejsmicznych (z których pierwszy odpowiadał powierzchni stropowej ogniwa Jantaru, a drugi stropowi formacji z Sasina), pod kątem wyznaczenia stref potencjalnych sweet spotów. Wyinterpretowany rozkład parametrów fizycznych takich jak: prędkość i gęstość, oraz parametrów petrofizycznych takich jak: porowatość, nasycenie, a także TOC w obrębie kolejnych sekwencji skalnych, ma ścisły związek z parametrami sejsmicznymi, takimi jak czasy przejścia fal odbitych od poszczególnych granic sejsmicznych oraz ich charakterystyk amplitudowych i częstotliwościowych [1, 8]. Zaprezentowana metodyka wyznaczania rozkładu parametru TOC była realizowana na podstawie zmienności amplitud pola falowego, pola prędkości fal podłużnych oraz pomiarów labolatoryjnych TOC w otworach wiertniczych. Uzyskane wyniki w postaci map i rozkładów przestrzennych pozwalają zidentyfikować obszary perspektywiczne o podwyższonych wartościach TOC, a ich dokładność uwarunkowana jest rozdzielczością danych sejsmicznych. Słowa kluczowe: TOC, inwersja genetyczna, impedancja akustyczna, shale gas. Calculation of organic matter volume (TOC) based on genetic inversion of seismic data The article presents the methodology for determining the distribution of TOC parameter on 3D seismic data, based on genetic inversion. The presented calculation scheme was applied on seismic data from the northern area of Poland. The main goal of the study was to recognize two seismic horizons (the first of which corresponded to the top surface of the Jantar Member and the second to the top of the Sasino Formation), in terms of designating the zones of potential sweet spots. The interpreted distribution of physical parameters such as: velocity and density and petrophysical parameters, such as: porosity, saturation and TOC within successive rock sequences, is closely related to seismic parameters, such as transition times of waves reflected from individual seismic boundaries and their amplitude characteristics and frequency [1, 8]. The presented methodology for determining the distribution of the TOC parameter, was based on the variability of the wave field amplitudes, the field of longitudinal wave velocities and laboratory TOC measurements in wellbores. Obtained results in the form of maps and spatial distributions allow to identify perspective areas with increased TOC values, and their accuracy is conditioned by the resolution of seismic data.
Uploads
Papers by Aleksander Wilk