Wyniki wyszukiwania - BazTech

1

Przestrzenna identyfikacja stref o lepszych parametrach zbiornikowych w utworach węglanowych

Bartoń Robert, Urbaniec Andrzej, Filipowska-Jeziorek Kinga

Wiadomości Naftowe i Gazownicze

|

2022

|

R. 25, Nr 11-12 (286)

4-10

PL

W publikacji zaprezentowano metodykę,która pozwoliła wyodrębnić strefy o korzystnych własnościach petrofizycznych umiejscowionych w stropowej partii kompleksu węglanowego górnej jury i dolnej kredy z obszaru środkowej części przedgórza Karpat. W przeprowadzonych badaniach wykorzystano zdjęcie sejsmiczne 3D oraz dane z otworów wiertniczych. Wyodrębnienie obiektów przestrzennych do potencjalnej sekwestracji CO2 było realizowane na podstawie atrybutów sejsmicznych obliczonych z inwersji symultanicznej. Inwersja sejsmiczna jest cennym narzędziem umożliwiającym estymację parametrów fizycznych ośrodka geologicznego z danych sejsmicznych, gdyż pozwala ona na przekształcenie amplitudy refleksów sejsmicznych w fizyczne parametry skał, a w konsekwencji w ilościowy opis złoża. Prędkość propagacji fal sejsmicznych jest jednym z podstawowych parametrów, który najbardziej wiarygodnie charakteryzuje właściwości fizyczne ośrodka geologicznego. Wykonane zostały wykresy krzyżowe impedancji fali podłużnej względem Lambda--Rho (Zp – λρ) oraz Lambda-Rho względem Mu-Rho (λρ – μρ), które w najlepszym stopniu odzwierciedlały zależności pomiędzy parametrami sprężystymi i elastycznymi. Opracowana metodyka może znaleźć zastosowanie zarówno do rozpoznawania stref o korzystniejszych parametrach zbiornikowych, jak również do bardziej zaawansowanych procesów budowy modeli statycznych i dynamicznych analizowanych formacji skalnych.

EN

The paper presents a methodology for distinguishing zones with better petrophysical properties in the uppermost part of the Upper Jurassic and Lower Cretaceous carbonate complex from the Carpathian Foreland area. For these studies 3D seismic survey and well data were used. Identification of spatial objects for potential CO2 sequestration was realized on the basis of seismic attributes calculated from simultaneous inversion. Seismic inversion is a useful tool for the estimation of reservoir properties from seismic data, as it enables the transformation of the amplitude of seismic reflections into physical parameters of rocks and, consequently, into a quantitative description of the reservoir. Propagation of velocity seismic waves is one of the basic parameters that most reliably characterizes the physical properties of a geological formation. Cross plots of longitudinal wave impedance versus Lambda-Rho (Zp – λρ) and Lambda-Rho versus Mu-Rho (λρ – μρ) were made, which best represented the relationships of the elastic parameters. Developed methodology can be applied both for identifying zones with more favorable reservoir parameters, as well as for more advanced processes of construction of static and dynamic models of the analyzed rock formations.

2

Metodyka identyfikacji obiektów perspektywicznych do potencjalnego składowania CO2 na przykładzie utworów węglanowych jury górnej

Bartoń Robert, Urbaniec Andrzej, Filipowska-Jeziorek Kinga

Nafta-Gaz

|

2022

|

R. 78, nr 5

343--357

PL

Celem artykułu jest opracowanie metodyki pozwalającej na wyodrębnienie stref/obiektów geologicznych o korzystniejszych własnościach petrofizycznych na podstawie analizy danych sejsmicznych i otworowych. Do badań tych wykorzystano zdjęcie sejsmiczne 3D z obszaru środkowej części przedgórza Karpat, a szczegółowe analizy prowadzono w obrębie stropowej partii kompleksu węglanowego górnej jury i dolnej kredy. W artykule przedstawiono wyniki analiz przeprowadzonych w obrębie centralnej części wspomnianego wyżej zdjęcia sejsmicznego. Wyodrębnienie obiektów przestrzennych do potencjalnej sekwestracji CO2 było realizowane na podstawie atrybutów sejsmicznych obliczonych z inwersji symultanicznej. Inwersja sejsmiczna jest cennym narzędziem umożliwiającym estymację parametrów fizycznych ośrodka geologicznego z danych sejsmicznych, gdyż pozwala ona na przekształcenie amplitudy refleksów sejsmicznych w fizyczne parametry skał, a w konsekwencji w ilościowy opis złoża. Prędkość propagacji fal sejsmicznych jest jednym z podstawowych parametrów, który najbardziej wiarygodnie charakteryzuje właściwości fizyczne ośrodka geologicznego. Wykonane zostały wykresy krzyżowe atrybutów impedancji fali podłużnej względem Lambda-Rho (Zp – λρ) oraz Lambda-Rho względem Mu-Rho (λρ – µρ), które w najlepszym stopniu odzwierciedlały zależności pomiędzy parametrami sprężystymi i elastycznymi. W obliczeniach wykorzystano opcję horizon probe dostępną w module Geobody Interpretation oprogramowania Petrel. Obliczenia prowadzono dla bramki czasowej obejmującej interwał od wyinterpretowanego horyzontu sejsmicznego odpowiadającego stropowi jury górnej wraz z dolną kredą (J3+K1str) do wartości czasu 120 ms poniżej tego horyzontu. Opracowana metodyka może znaleźć w przyszłości zastosowanie zarówno do rozpoznawania stref o korzystniejszych parametrach zbiornikowych, jak również do bardziej zaawansowanych procesów budowy modeli statycznych i dynamicznych analizowanych formacji skalnych. Wyznaczone obiekty po przeprowadzeniu niezbędnych analiz oraz modelowań mogą zostać wykorzystane do potencjalnego składowania CO2.

EN

The aim of this paper is to develop a methodology to identify geological zones/objects with more favorable petrophysical properties based on analysis of seismic and well data. For these studies 3D seismic image from the middle part of the Carpathian Foreland was used, and detailed analyses were carried out within the top part of the Upper Jurassic and Lower Cretaceous carbonate complex. This paper presents results of performed analysis in the central part of the above-mentioned seismic image. Identification of spatial objects for potential CO2 sequestration was realized on the basis of seismic attributes calculated from simultaneous inversion. Seismic inversion is a useful tool for the estimation of reservoir properties from seismic data, as it enables the transformation of amplitude of seismic reflections into physical parameters of rocks and, consequently, into a quantitative description of the reservoir. Propagation of velocity seismic waves is one of the basic parameters that most reliably characterizes the physical properties of a geological medium. Cross plots of longitudinal wave impedance attributes versus Lambda-Rho (Zp – λρ) and Lambda-Rho versus Mu-Rho (λρ – µρ) were made, which best represented the relationships of the elastic parameters. The horizon probe option available in Geobody Interpretation module of Petrel software was used to perform calculations. A time gate covering the interval from the interpreted seismic horizon corresponding to the Upper Jurassic – Lower Cretaceous top (J3+K1str) to the time value of 120 ms below this horizon was adopted. The developed methodology can be applied in the future both for identification of the zones of better reservoir parameters as well as for more advanced processes of static and dynamic models building for the analyzed rock formations. After conducting necessary analyses and modeling the identified objects can be used for potential CO2 storage.

3

Detekcja stref o lepszych parametrach zbiornikowych na podstawie atrybutów obliczanych z inwersji sejsmicznej

Bartoń Robert, Urbaniec Andrzej, Filipowska-Jeziorek Kinga

Nafta-Gaz

|

2022

|

R. 78, nr 12

845--860

PL

Celem artykułu było wyodrębnienie stref o korzystniejszych parametrach zbiornikowych na podstawie analizy atrybutów sejsmicznych obliczanych z inwersji sejsmicznej w połączeniu z danymi otworowymi. Obliczenia przeprowadzono na wolumenie sejsmicznym ze zdjęcia 3D usytuowanego w środkowej części przedgórza Karpat. Przedmiotem badań była stropowa część kompleksu węglanowego górnej jury. Do identyfikacji stref o preferowanych parametrach wykorzystane zostały wykresy krzyżowe atrybutów impedancji fali podłużnej względem Lambda-Rho (Zp – λρ) oraz Lambda-Rho względem Mu-Rho (λρ – μρ). W obliczeniach wykorzystano opcję horizon probe dostępną w module Geobody Interpretation oprogramowania Petrel. Na podstawie otrzymanych wyników można stwierdzić, że najlepsze dopasowanie obrazu na przekrojach sejsmicznych z profilowaniem porowatości (PHI) w otworze G-4 uzyskano dla wariantu obliczeń Lambda-Rho vs. Mu-Rho, w którym w obrębie wyodrębnionych stref o lepszych parametrach zbiornikowych wydzielono trzy klasy zakresów wartości. Pozwoliło to również w lepszym stopniu zwizualizować zmienność analizowanych parametrów. Wyniki przeprowadzonych badań wskazują, że utwory węglanowe górnej jury w obszarze badań mają zróżnicowane parametry zbiornikowe, przy czym strefy o lepszych własnościach zlokalizowane są głównie w stropowej części tego kompleksu. Generalnie rejon o lepszym rozwoju własności zbiornikowych znajduje się na północny zachód od analizowanego otworu G-4. Można przypuszczać, że jest to związane z obecnością w tym rejonie regionalnych dyslokacji, wokół których doszło do rozwoju szczelinowatości. Najprawdopodobniej dodatkowym czynnikiem odpowiedzialnym za wzrost porowatości w tej strefie jest powierzchnia krasowa, rozwinięta bezpośrednio w stropie utworów jury. Prezentowana metodyka może znaleźć zastosowanie nie tylko do bezpośredniej identyfikacji poziomów skał zbiornikowych o preferowanych parametrach, ale może być też szeroko wykorzystywana w różnego typu analizach i modelowaniach przestrzennych.

EN

The aim of this paper was to distinguish zones with more preferred reservoir parameters based on the analysis of seismic attributes calculated from seismic inversion, combined with well data. The calculations were performed on 3D seismic volumes located in the central part of the Carpathian Foreland. The subject of the study was the upper part of the Upper Jurassic carbonate complex. Cross plots of longitudinal wave impedance attributes against Lambda-Rho (Zp – λρ) and Lambda-Rho against Mu-Rho (λρ – µρ) were used to identify zones with preferred parameters. The calculations used the horizon probe option available in the Geobody Interpretation module in Petrel software. Based on the results, it was concluded that the best match on seismic sections with porosity profiling (PHI) in G-4 well was obtained for the variant of Lambda-Rho vs. Mu-Rho calculations, in which three classes of value ranges were separated within the identified zones with preferable reservoir parameters. It allowed to better visualize the variability of the analyzed parameters. The results of the research indicate that the Upper Jurassic carbonate formations in the study area are characterized by varying reservoir parameters, with the zones with better properties located mainly in the top part of this complex. In general, the region with better development of reservoir properties is located to the northwest of the analyzed G-4 well. It can be assumed that this is due to the presence of regional dislocations in this area, around which fracturing has developed. Most likely, an additional factor responsible for the increase in porosity in this zone is the karst surface, developed directly in the uppermost part of the Jurassic formations. The presented methodology can find application not only for direct identification of reservoir intervals with preferred parameters, but can be widely used in various types of spatial analysis and geophysical modeling.

4

Predicting gas content in coalbed methane reservoirs using seismic waveform indication inversion: a case study from the Upper Carboniferous Benxi Formation, eastern Ordos Basin, China

Wang Jiabei, Li Zhijun, Chen Cen, Meng Cheng

Acta Geophysica

|

2022

|

Vol. 70, no 2

623--638

EN

The identification of gas content is very important for exploration and development of coalbed methane (CBM) reservoirs. As a kind of gas-bearing reservoirs in coal seam, CBM reservoirs usually show strong heterogeneity, which makes the gas content varies greatly in the strata. What’s more, the thin interlayer that is common in coal-bearing formation makes it difficult to predict the favorable gas-bearing distribution based on conventional methods. In this study, a seismic waveform indication inversion method was applied to reveal the gas content of No. 8 coal seam reservoirs in the Upper Carboniferous Benxi Formation of the DJ area in the eastern of the Ordos Basin, China. The first step of this method is to calculate the p-wave impedance inversion volume of No. 8 coal seam. The second step is to build the correlation between elastic parameters and gas content in No. 8 coal seam. Through the statistical analysis based on velocity, density, p-wave impedance and measured gas content data of 16 wells, the fitting formula between p-wave impedance and gas content is obtained, with a highest correlation coefficient up to 96%. The third step is to calculate the gas content data volume from the p-wave impedance inversion volume by the above fitting formula, and then the quantitative plane distribution of gas content in No. 8 coal seam can be predicted. The prediction results indicate that the gas content of No. 8 coal seam can be divided into two Classes. To verify the reliability of the inversion results, the production data of well X11 was applied to verify the gas content which located in the Class I area. The application of seismic waveform indication inversion has provided a precise prediction for the spatial distribution of gas content in CBM reservoirs, serving as a basis for locating and designing wells for CBM development.

5

Seismic inversion based prediction of the elastic parameters of rocks surrounding roadways: a case study from the Shijiazhuang mining area of North China

Shi Suzhen, Feng Jian, Liu Zuiliang, Qi Youchao, Duan Peifei, Han Qi, Zuo Jianping

Acta Geophysica

|

2021

|

Vol. 69, no. 6

2219--2230

EN

Increasing mining depths and intensities have resulted in frequent accidents during roadway excavation, thus increasing the necessity of preemptively determining the geological conditions of rock formations along roadways. Identifying the spatial distributions of lithotypes and the mechanical properties of non-marker bed intervals in coal measure strata are crucial for roadway horizon design, roof management, and efficient excavation. Therefore, based on petrophysical tests, well logging, and seismic data, a comprehensive method is proposed, which uses post-stack wave impedance inversion, petrophysical statistical analyses, and pre-stack simultaneous inversion to predict the lithology and elastic parameters of rocks surrounding roadways. The lithology of the target layer is first divided according to the wave impedances of the inversion body, and the sand lenses within the mudstone interval are described. Based on petrophysical test data, the shear-wave conversion equations of the sand lenses in this target interval and background mudstone are then obtained via marginal limit extrapolation. Finally, pre-stack simultaneous inversion, based on the seismic data and shear-wave conversion curves of different lithologies, is performed to obtain the elastic parameter profile of the target interval in the cross-section of the roadway, laying the foundation for further prediction of the mechanical properties of the surrounding rocks. The results agree with roadway drilling data by>90%.

6

Pore-type identifcation of a heterogeneous carbonate reservoir using rock physics principles: a case study from south-west Iran

Emami Niri Mohammad, Mehmandoost Fariba, Nosrati Hossein

Acta Geophysica

|

2021

|

Vol. 69, no. 4

1241--1256

EN

The characterization of carbonate rocks is a challenging process compared to siliciclastics because of their more intricate pore-space structure. In this study, we applied rock physics methods to a heterogeneous carbonate reservoir located in southwest Iran to identify zones of various pore types, including inter-particle pores, stif (vuggy and moldic) pores, and cracks, from geophysical measurements. We frst constructed two rock physics templates (RPTs) using well logs and used them to quantitatively analyze the efect of various pore types on elastic properties. Using these RPTs, we then implemented an inversion algorithm to estimate the volume fractions of various pores using total porosity and P-wave velocity (Vp) derived from well logs. Next, we have compared the pore-type inversion results and the image log interpretation/core images at the corresponding depth intervals to validate inversion results. This comparison showed that the inversion results are consistent with the measurements. Also, we applied the introduced pore-type inversion algorithm on seismic data to achieve insight into pore-type distribution in the 3D framework of the reservoir under study. The results of these rock physics-based analyses indicate that the inter-particle pores are dominant in the pore-space, while there are stif pores and dispersal cracks in some subzones of the studied depth interval. Additionally, employing the Xu and Payne rock physics modeling procedure, P- and S-wave velocities were estimated using pore-type inversion results at the location of several wells from the studied feld. The calculated mean absolute error and the correlation coefcient indicate a high consistency between the measured and modeled velocities. This research may contribute to permeability prediction and analysis of the diagenetic processes’ impact on reservoir properties.

7

Machine learning based prediction of regularization parameters for seismic inverse problems

Liu Shihuan, Zhang Jiashu

Acta Geophysica

|

2021

|

Vol. 69, no. 3

809--820

EN

Regularization parameter selection (RPS) is one of the most important tasks in solving inverse problems. The most common approaches seek the optimal regularization parameter (ORP) from a sequence of candidate values. However, these methods are often time-consuming because they need to conduct the estimation process on all candidate values, and they are always restricted to solve certain problem types. In this paper, we propose a novel machine learning-based prediction framework (MLBP) for the RPS problem. The MLBP frst generates a large number of synthetic data by varying the inputs with diferent noise conditions. Then, MLBP extracts some pre-defned features to represent the input data and computes the ORP of each synthetic example by using true models. The pairs of ORP and extracted features construct a training set, which is used to train a regression model to describe the relationship between the ORP and input data. Therefore, for newly practical inverse problems, MLBP can predict their ORPs directly with the pre-trained regression model, avoiding wasting computational resources on improper regularization parameters. The numerical results also show that MLBP requires signifcantly less computing time and provides more accurate solutions for diferent tasks than traditional methods. Especially, even though the MLBP trains the regression model on synthetic data, it can also achieve satisfying performance when directly applied to feld data.

8

Global stochastic seismic inversion using turning bands simulation and co-simulation

Sadeghi Mehdi, Amini Navid, Falahat Reza, Sabeti Hamid, Madani Nasser

Acta Geophysica

|

2021

|

Vol. 69, no. 5

1717--1734

EN

In this study, we develop a 3D geostatistical seismic inversion method based on turning bands simulation and co-simulation to estimate acoustic impedance (AI) from seismic and well data. The proposed method uses an iterative approach based on cross-over concept in genetic algorithm optimization to perform global stochastic inversion. The objective function of the optimization algorithm is the measure of correlation coefficients between the modeled and observed seismic data. In the first iteration of the proposed algorithm, the seismic cubes corresponding to the AI realizations generated by the turning bands simulation are compared to the observed seismic data. Subsequently, the local areas of AI models producing the highest correlation coefficients to the observed seismic data are merged to construct a new supplementary cube that is further used as the secondary variable in the turning bands co-simulation for the next iteration. Generation of the seismic cubes and picking the best AI cube via cross-over genetic algorithm approach is performed iteratively until convergence to a constant correlation coefficient between the modeled and observed seismic data. The proposed method is applied to a synthetic dataset that includes 20 wells with known AI logs and 3D stacked seismic data. According to the results, the algorithm converges to a constant correlation coefficient after a few iterations. In addition, it is observed that employing multi-attribute analysis outputs (meta-attributes) during turning bands co-simulation in the initialization step would improve the final global correlation coefficient from 0.774 to 0.815.

9

Rock physics based estimation of quartz content in the Shenhu area, South China Sea

Deng Wei, Liang Jinqiang, He Yulin, Kuang Zenggui, Meng Miaomiao

Acta Geophysica

|

2020

|

Vol. 68, no. 6

1619--1641

EN

Estimation of quartz content (QC) is important for gas hydrate production. However, the existing methods pay more attention to estimate saturations of hydrate or free gas instead of QC. QC is difcult to be estimated because of its limited and unclear infuences on velocities in the hydrate-bearing sediments. A workfow was proposed to estimate QC from core to logging to seismic inversion whose core technologies were an unconsolidated anisotropic model (UAM) and an inverse modeling approach. We used the UAM to construct the quantitative relationships between physical properties including QC and velocity. Then, the velocities of the reservoirs were obtained by logging-constrained seismic inversion. Finally, QC can be scaled by the inverse modeling of the UAM. To build the UAM, we analyzed the physical properties of hydrate reservoirs based on the cores and logging while drilling (LWD) data in the Shenhu area, South China Sea, and characterized unconsolidated sediments with horizontally layered hydrates and gas occurrences. The calculated P-velocities and S-velocities from the UAM agreed with the LWD data when the input variables were QC, porosity, hydrate saturation, and gas saturation. Conversely, for a group of P-velocity and S-velocity from seismic inversion, the corresponding QC can be scaled out as well as the other parameters based on the UAM, which was defned as an inverse modeling. Because the signifcant parameters such as hydrate saturation have been considered as independent variables in the model, we can efectively avoid the correlation between QC and the others. Prediction multiplicity can be reduced. The estimated QC was consistent with the drilling and geological understanding in the feld application, indicating that the method proposed is efective and practical.

10

High-resolution reflectivity inversion based on joint sparse representation

Shi Zhanzhan, Zhou Huailai, Wang Yuanjun, Niu Cong, Huang Rao

Acta Geophysica

|

2019

|

Vol. 67, no. 6

1535--1550

EN

High-resolution reflectivity inversion is termed as a fundamental yet essential step for the prediction of thin-bedded hydrocarbon reservoirs. However, algorithms suffer from two key issues: (1) seismic inversion is an ill-posed problem that has multiple solutions, and the results of trace-by-trace seismic inversion are quite poor in lateral continuity, and (2) algorithm stability is likely to be decreased owing to the noise and distortion associated with the acquisition and processing flows. In the current article, we formulate a new joint sparse representation through the combination with L2,1- norm misfit function, which possesses superior noise robustness, in particular in the presence of outliers. On the basis of the L2,1- norm regularization, this specific approach enforces a common sparsity profile, together with consistently lowering the multiplicity of solution. Subsequent to that, the resultant algorithm is applied to the multi-trace seismic inversion. Besides, the wedge model trial and practical applications suggest that the proposed inversion algorithm is stable, in addition to having good noise robustness and lateral continuity; moreover, the vertical resolution of λ/8 is realized under the noise and outliers interference. The logging data calibration illustrates that the proposed methodology is accurate and credible.

11

Ocena efektywności sposobów zwiększania rozdzielczości danych sejsmicznych w poszukiwaniach węglowodorów w świetle zastosowanej metodyki akwizycji i procedur modyfikacji spektralnej

Jędrzejowska-Tyczkowska H.

Nafta-Gaz

|

2018

|

R. 74, nr 1

3--21

PL

Porównano efektywność sposobów zwiększania rozdzielczości danych sejsmicznych realizowanych na dwóch różnych etapach projektów sejsmicznych 3D, a to w drodze projektowania zaawansowanych schematów akwizycji, umożliwiających uzyskiwanie obrazu falowego o bardzo wysokiej krotności sumowania (rzędu 2260), jak też w drodze stosowania procedur przetwarzania w domenie spektralnej o (autorskiej) nazwie „modyfikacja spektralna”. Artykuł jest czwartą publikacją omawiającą powyższe zagadnienie w optyce innowacyjnej opracowanej przez autorkę w ramach realizacji jednego z zadań programu Blue Gas – Polski Gaz Łupkowy. Przedstawione w artykule porównania wskazują, że zastosowanie określonej metodyki istotnie zależy od jakości materiału sejsmicznego. W niniejszej pracy wykorzystano obraz falowy uzyskany od zleceniodawcy (migracja czasowa na sumie – migration AFTER STACK). Na podstawie wykonanych porównań należy sądzić, że decyzja o wyborze sposobu zwiększania rozdzielczości powinna brać pod uwagę aspekt ekonomiczny końcowej realizacji.

EN

Comparison of the efficiency of the methods for increasing the resolution of seismic data was performed at two different stages of 3D seismic survey design, by creating advanced acquisition schemes to obtain a very high fold-coverage stacking (with the range of 2260), as well as using processing procedures in the spectral domain – “spectral modification”. This article is the fourth publication discussing the issue above, developed by the author within one of the tasks of the Blue Gas – Polski Gaz Łupkowy program. The comparison presented in this paper indicates, that the application of a particular methodology, depends on the quality of the seismic data. In this study a wavefield image obtained from the client (After Stack time migration) was used based on the comparisons made, it should be assumed that the decision on choosing the method of increasing the resolution should take into account the economic aspect of the final implementation.

12

Zastosowanie prędkości interwałowych PPS do modelowania ośrodka geologicznego

Bartoń R.

Nafta-Gaz

|

2016

|

R. 72, nr 3

155--163

PL

Przedmiotem niniejszego artykułu jest prezentacja modelowania strukturalnego ośrodka geologicznego w domenie głębokości przy wykorzystaniu prędkości fal podłużnych pochodzących z azymutalnych pomiarów pionowego profilowania sejsmicznego PPS. W publikacji zostały przedstawione wyniki modelowania pola prędkości na bazie rozkładu prędkości interwałowych, obliczonych na podstawie inwersji sejsmicznej na transformacjach PPS-WPG. Uzyskany model prędkości posłużył do wykonania transformacji czas–głębokość sekcji sejsmicznych i obliczenia map strukturalnych. Przedmiotowa tematyka była podyktowana możliwością szerszego wykorzystania pomiarów PPS w polskim przemyśle naftowym i zwiększenia szczegółowości rozpoznania ośrodka geologicznego (w stosunku do sejsmiki powierzchniowej), co może w istotnym stopniu rzutować na efektywność poszukiwań i eksploatacji węglowodorów oraz na dokładność wierceń otworów poziomych lub kierunkowych.

EN

This article presents results of modeling structural geological medium in depth domain, by using the velocity longitudinal waves, originating from azimuthal Vertical Seismic Profiling (VSP) measurements. The paper presents results of interval velocity modeling based on velocity distribution, calculated from seismic inversion on VSP-CDP transformations. The obtained velocity model was used to perform time-depth transformation of seismic sections and calculations of structural maps. The present topic was dictated by the possibility of wider useage of VSP measurements in the Polish Oil Industry, increasing the detail identification of geological medium (in relation to seismic surface), which may significantly affect the efficiency of exploration, exploitation of hydrocarbons and the accuracy of drilling horizontal or directional holes.

13

Obliczanie inwersji sejsmicznej na azymutalnych transformacjach PPS-WPG fal podłużnych i poprzecznych

Bartoń R.

Nafta-Gaz

|

2015

|

R. 71, nr 10

711--719

PL

W artykule przedstawiono wyniki obliczeń inwersji sejsmicznej wykonywanej na azymutalnych transformacjach PPS-WPG fal podłużnych PP i poprzecznych PS. Realizacja przedmiotowego zadania była prowadzona na zarejestrowanych wieloazymutalnych pomiarach PPS 3C, pochodzących z północnego obszaru Polski. Otrzymane rezultaty użycia inwersji sejsmicznej, obliczanej metodą rekursywną, na transformacjach PPS-WPG wskazują, że metoda ta daje zadowalające wyniki do wyznaczania rozkładu zmienności prędkości akustycznych i gęstości w strefie okołootworowej.

EN

This article presents results of seismic inversion calculation, on VSP-CDP azimuthal transformations for longitudinal PP and transverse PS waves. Implementation of this task was performed on recorded VSP 3C multi-azimuthal measurements from northern Poland. The results of seismic inversion, calculated by recursive method on VSP-CDP transformations, indicate that this method provides satisfactory results for the determination of variability acoustic velocity and density distribution in a wellbore zone.

14

Estymacja parametrów charakteryzujących ośrodek geologiczny za pomocą procedury inwersji genetycznej

Leginowicz A., Pirowska K.

Nafta-Gaz

|

2013

|

R. 69, nr 5

392--400

PL

Estymacja parametrów charakteryzujących ośrodek geologiczny za pomocą inwersji genetycznej przy użyciu danych sejsmicznych i danych geofizyki otworowej pozwala na uzyskanie wolumenu prędkości lub impedancji akustycznej. W publikacji przedstawiono rezultaty inwersji genetycznej opartej na algorytmie, który wykorzystuje połączenie wielowarstwowych sieci neuronowych oraz algorytmów genetycznych. Wykorzystywany w tym rozwiązaniu algorytm wstecznej propagacji błędów jest podstawowym algorytmem uczenia nadzorowanego dla wielowarstwowych jednokierunkowych sieci neuronowych. Istotną zaletą inwersji genetycznej jest możliwość zastosowania jej dla różnych atrybutów petrofizycznych. Danymi wejściowymi mogą być: dane otworowe, atrybuty sejsmiczne, mapy lub inne parametry charakteryzujące ośrodek geologiczny. Procedurę zaaplikowano na danych szkoleniowych firmy Schlumberger oraz na rzeczywistych danych sejsmicznych 3D, pochodzących z przedgórza Karpat.

EN

The estimation of rock properties characterizing reservoirs by genetic inversion using seismic and well data allows us to obtain the volume of velocity or acoustic impedance. This work presents results of genetic inversion within the Petrel system based on a nonlinear multitrace seismic inversion algorithm. In the case of Genetic Inversion, the required inputs are limited to the seismic amplitude, and the Acoustic Impedance well logs used as training data. Indeed no single unique wavelet, neither initial property modeling are needed as inputs prior to the running of the inversion. A genetic algorithm back-propagates the error in order to update the weights for the neural networks. The Petrel inversion module is not restricted to a pure acoustic impedance inversion, but can be extended to any property with some inherent link with respect to the input seismic volume. This analysis was carried out on demo data and on 3D seismic data, from the Carpathians foreland.

15

Detection of salts deposits geometry variation

Kawalec-Latała E.

AGH Journal of Mining and Geoengineering

|

2012

|

Vol. 36, no. 2

161-169

EN

The construction of underground reservoirs is important for energy security and ecology. Salt deposits reveal favourable conditions for the construction of reservoirs. The inhomogeneities within salt deposit as well as changes of lithology, thickness and geometry variation should be especially predicted before the underground reservoir location. The work concerns the possibility of the recognition of the space distribution of anhydrites in the rock salt. Acoustic impedance is one of the basic factors characterising physical features of rocks. Seismic survey is the most economical source of subsurface information. A method of inversion of seismic section to pseudo-acoustic impedance section gives the possibility of the recognition of discontinuities within salt deposits. The paper includes a number of pseudo-acoustic impedance sections constructed for geometrical variations of anhydrites in rock salt. For the modelling of pseudo impedance sections a computer system INVERSE was used.

PL

Cechsztyńskie złoża soli kamiennej rozpatrywane są jako jedno z ważnych, potencjalnych miejsc składowania odpadów radioaktywnych oraz dwutlenku węgla CO2. Korzystne warunki do budowy takich kawern stwarzają złoża pokładowe soli kamiennej o znacznej miąższości (pow. 100 m grubości). Modelowania syntetycznych sekcji pseudoimpedancji akustycznej generowanych z uwzględnieniem różnej konfiguracji przestrzennej anhydrytów zaburzających jednorodność budowy pokładów soli kamiennych, perspektywicznych dla budowy podziemnych zbiorników. Testy wykonano na danych syntetycznych sekcji pseudoimpedancji akustycznej generowanych dla teoretycznego sejsmogeologicznego modelu ośrodka przybliżającego budowę pokładu soli kamiennej. Do konstrukcji modeli sejsmogeologicznych zawierających pokładowe złoża soli wykorzystano parametry sprężyste typowe dla cechsztynu oraz uzyskane na podstawie prędkości propagacji fal sejsmicznych, ustalonych w oparciu o wyniki profilowań geofizycznych w otworach kontrolno-pomiarowych w rejonie projektowanego szybu SW-4. Do modelowania syntetycznych sekcji pseudo-impedancji akustycznej wykorzystano autorski program INWERS.

16

Sejsmiczna inwersja stochastyczna w systemie INPRES

Bartoń R.

Nafta-Gaz

|

2011

|

R. 67, nr 6

385-394

PL

W artykule przedstawiono ogólne informacje oraz podstawowe zasady postępowania przy tworzeniu projektów inwersji sejsmicznej na podstawie danych sejsmicznych i geofizycznych. Zasadnicza część artykułu dotyczy prezentacji wyników obliczania inwersji stochastycznej w Instytucie Nafty i Gazu, realizowanej w systemie INPRES, na danych sejsmicznych i geofizycznych pochodzących z rejonu zapadliska przedkarpackiego oraz Niżu Polskiego.

EN

This article presents the general information and essentials proceedings in creation projects of the seismic inversion based on geophysical and seismic data. The fundamental part of this article concerns the results of creation and calculation of the stochastic inversion realized in the Oil & Gas Institute in the system INPRES, based on seismic and geophysical data, coming from the regions of the Pre-Carpathian depression and Polish Depression.

17

Porównanie porowatości utworów dolomitu głównego wyznaczonych przy pomocy różnych metod inwersji sejsmicznej

Solarski T.

Prace Naukowe Instytutu Nafty i Gazu

|

2010

|

nr 170 wyd. konferencyjne

503--507

PL

Rozkład porowatości w skałach zbiornikowych może być określany na podstawie impedancji akustycznej wyznaczanej z danych sejsmicznych przy pomocy metody inwersji. Poszczególne odmiany tej metody prowadzą jednak do różnych wyników. Głównym celem przedstawianej pracy było porównanie porowatości utworów dolomitu głównego wyznaczonych przy pomocy wybranych metod inwersji sejsmicznej na obszarze zdjęcia sejsmicznego 3D Międzychód — Sieraków.

EN

Porosity distribution in reservoir rocks can be predicted on the basis of acoustic impedance derived from seismic inversion. Different inversion methods can lead to different results. The principal objective of the presented work was comparison of the Main Dolomite porosity distributions determined with use of selected inversion methods.

18

Dynamiczny model prędkości sejsmicznych dla studium wykonalności sejsmiki 4D określony drogą inwersji sejsmicznej i geostatystycznej z danych sejsmicznych i inżynierii złożowej

Jędrzejowska-Tyczkowska H., Malaga M., Żukowska K.

Prace Naukowe Instytutu Nafty i Gazu

|

2010

|

nr 170 wyd. konferencyjne

515--520

PL

W publikacji przedstawiono metodykę dynamicznego prognozowania modelu prędkościowo-głębokościowego, realizowaną z wykorzystaniem inwersji sejsmicznej oraz inwersji geostatycznej danych produkcyjnych (ciśnień). Wyniki stanowią substytut pomiarów sejsmiki 4D.

EN

Dynamic velocities' model requires INPUT data for all reservoir simulation procedures. Solution presented mimics the results of time-lapse seismics.

19

Inwersja sejsmiczna akustyczna i elastyczna fal podłużnych, konwertowanych i poprzecznych w zagadnieniach interpretacji złożowej

Jędrzejowska-Tyczkowska H.

Prace Instytutu Nafty i Gazu

|

2009

|

nr 160

1-85

PL

Celem przedstawionej publikacji jest pokazanie możliwości metody sejsmicznej fal odbitych - podłużnych i konwertowanych - w zadaniach odtworzenia pola prędkości fal poprzecznych. Prędkość propagacji fal sejsmicznych to jeden z parametrów, który najbardziej wiarygodnie charakteryzuje właściwości fizyczne ośrodka geologicznego. W zadaniach określenia właściwości zbiornikowych oraz nasyceń obiektów złożowych szczególnego znaczenia nabiera informacja o prędkości fali poprzecznej, ze względu na jej specyficzną reakcję na obecność fazy ciekłej (i gazowej) w ośrodku. Pomiary geofizyki otworowej oferują bezpośredni pomiar prędkości fali poprzecznej w procedurze "pełne pole falowe" (full wave sonic). Pomiary sejsmiki powierzchniowej - raczej sporadycznie, a często wyłącznie w zadaniach badawczych - wykonywane są w pełnej opcji wieloskładnikowej (multicomponent) - tzw. Sejsmika 3C - przede wszystkim ze względu na odrębny od fali podłużnej P sposób generowania fal poprzecznych S. Znacznie mniej czasochłonne i kosztochłonne jest wykorzystanie do określenia pola prędkości fali poprzecznej pomiarów fal konwertowanych. Podstawową procedurą odtworzenia prędkości sejsmicznych jest inwersja sejsmiczna. W pracy scharakteryzowano różne typy inwersji akustycznej i elastycznej, i na podstawie danych literaturowych pokazano liczne przykłady zastosowania inwersji do charakterystyk własności zbiornikowych. Zaproponowano również odmienne od rozpoznanych w literaturze podejście do oceny prędkości fali poprzecznej, stosując jednocześnie inwersję elastyczną i akustyczną fali konwertowanej w interaktywnej pętli, wraz z modelowaniem sejsmicznym. Ponadto opracowano model testowy, oparty na wynikach interpretacji rzeczywistych danych polowych sejsmiki 2D-3C (Projekt Rajsko-Krzeczów). W fazie finalnej, na podstawie mapy prędkości fal P i S obliczono współczynnik Vp/Vs i oceniono jego skuteczność jako wskaźnika dyskryminującego w ośrodku geologicznym obecność fazy nasyconej węglowodorami.

EN

The goal of this publication is to present the possibilities of the seismic method of longitudinal and converted reflected waves in reconstruction of transverse waves velocity field. The velocity of seismic waves propagation is one of the parameters which describe the physical properties of the geological medium the most authoritatively. In tasks targeted at defining the reservoir properties and saturation of deposit objects the information on transverse wave velocity acquires a particular importance on account of its specific reaction to the liquid (and gaseous) phase in the medium. The measurements of borehole geophysics offer direct measurement of transverse wave velocity in full wave sonic procedure. The measurements of surface seismics are rarely performed in full 3C seismic option, and often exclusively for research purposes, mostly on account of the method of generating transverse S waves independent from longitudinal P wave. Significantly less time-consuming and less expensive is the application of the measurement of converted waves to define the transverse wave velocity field. The basic procedure in reconstruction of seismic velocities is the seismic inversion. Different types of acoustic and elastic inversion are described in the study and, on the basis of available literature specification, numerous applications of inversion are shown to characterize the reservoir properties. Also, a different approach was suggested from those described in literature to the evaluation of transverse wave velocity, using simultaneously elastic and acoustic inversion of converted wave in an interactive loop along with seismic modeling. Moreover, a test model was developed, based on the results of interpretation of the 2D-3C seismic real field information (Rajsko-Krzeczów project). At the final stage, based on the P and S waves velocity map, the coefficient Vp/Vs was calculated and its effectiveness was assessed as an indicator discriminating the hydrocarbon saturated phase in the geological medium.

20

Najnowsze metody inwersji sejsmicznej na przykładzie zdjęcia sejsmicznego Parzęczewo-Czarna Wieś

Dalętka A.

Prace Instytutu Nafty i Gazu

|

2008

|

nr 150 wydanie konferencyjne

153--158

PL

Tematem pracy jest wykorzystanie najnowszych metod inwersji sejsmicznej dla osadów czerwonego spągowca w obszarze zdjęcia sejsmicznego 3D Parzęczewo—Czarna Wieś. Wcześniejsze badania sejsmiczne przeprowadzone na wspomnianym terenie umożliwiły odkrycie złóż gazu typu litologicznego i strukturalno-litologicznego (Ujazd, Paproć—Cicha Góra). Analizy sedymentologiczne i geologiczne oparte na danych z głębokich odwiertów, wskazały na możliwość istnienia złóż związanych z wyklinowaniami serii osadowej czerwonego spągowca na północnym skłonie wału wolsztyńskiego [2]. Wykrycie pułapek strukturalno-litologicznych i określenie linii zasięgu piaskowców saksonu wymaga szczegółowych analiz amplitudowych odbicia od stropu serii osadowej czerwonego spągowca. Autor wykonał kompleksową analizę wspomnianego refleksu w oparciu o zapis sejsmiczny przygotowany z zachowaniem rzeczywistych relacji amplitud. W ramach ww. badań wykonano modelowanie syntetycznych kolekcji CDP oraz inwersję symultaniczną (inwersję sejsmiczną przed składaniem), której parametry oparte zostały na licznych krzywych karotażowych. Przygotowania do projektu oraz rezultaty inwersji zostaną szczegółowo przedstawione w czasie prezentacji.

EN

This paper describes the up-to-date seismic inversion methods that were utilized in the area of 3D seismic project Parzęczewo—Czarna Wieś. The previous 3D seismic investigation in this region enabled to discover the presence of gas reservoirs in the form of structural- stratigraphic and stratigraphic traps (Ujazd, Paproć—Cicha Góra gas field). Numerous geological and sedimentological investigations based on well data indicated that the area in question is still highly prospective. In order to outline stratigraphic traps, a careful amplitude analysis of top reservoir reflection is required. On the basis of pre-stack seismic data, processed with preserved true relative amplitudes, author performed comprehensive study of the Zechstein Base reflection. Within the scope of the project an extensive synthetic gather generation and simultaneous prestack inversion was carried out. The preparations and results will be thoroughly discussed during presentation.