The Kula Field is the youngest volcanic center in western Türkiye, and consists of various well-preserved volcanic products. Although many geological studies have been conducted in the region, geophysical anomalies have not been studied in detail. Therefore, we analyzed the aeromagnetic anomalies of these volcanic products by performing inversion studies with a recently proposed global optimizer. This study is the first attempt to use success-history-based adaptive differential evolution algorithm (SHADE) for inverting magnetic anomalies. To reduce the computational cost, we introduced the E-SHADE scheme by incorporating an exponential population reduction strategy into the optimizer. A synthetic anomaly study revealed the mathematical nature of the handled inverse problem. Some pre- and post-inversion analyses showed the efficiency of the proposed algorithm. Additionally, we observed that the E-SHADE algorithm produced better results than a commonly used derivative-based local optimizer. Nine profile data sets including magnetic anomalies of some volcanic cones in the Kula region were inverted. It was determined that the basaltic intrusions that allow the mantle material to uplift rapidly are not very deep in the subsurface. Therefore, it is possible that the three-phased volcanism may become active again and generate new alkaline basaltic lava flows in a new phase through these shallow dykes.
{"title":"Investigation of Kula Volcanic Field (Türkiye) Through the Inversion of Aeromagnetic Anomalies Using Success-History-Based Adaptive Differential Evolution with Exponential Population Reduction Strategy","authors":"Yunus Levent Ekinci, Çağlayan Balkaya, Hanbing Ai, Arka Roy, Şenol Özyalin","doi":"10.1007/s00024-024-03569-y","DOIUrl":"https://doi.org/10.1007/s00024-024-03569-y","url":null,"abstract":"<p>The Kula Field is the youngest volcanic center in western Türkiye, and consists of various well-preserved volcanic products. Although many geological studies have been conducted in the region, geophysical anomalies have not been studied in detail. Therefore, we analyzed the aeromagnetic anomalies of these volcanic products by performing inversion studies with a recently proposed global optimizer. This study is the first attempt to use success-history-based adaptive differential evolution algorithm (SHADE) for inverting magnetic anomalies. To reduce the computational cost, we introduced the E-SHADE scheme by incorporating an exponential population reduction strategy into the optimizer. A synthetic anomaly study revealed the mathematical nature of the handled inverse problem. Some pre- and post-inversion analyses showed the efficiency of the proposed algorithm. Additionally, we observed that the E-SHADE algorithm produced better results than a commonly used derivative-based local optimizer. Nine profile data sets including magnetic anomalies of some volcanic cones in the Kula region were inverted. It was determined that the basaltic intrusions that allow the mantle material to uplift rapidly are not very deep in the subsurface. Therefore, it is possible that the three-phased volcanism may become active again and generate new alkaline basaltic lava flows in a new phase through these shallow dykes.</p>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"30 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142267080","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-12DOI: 10.1007/s00024-024-03570-5
Anna Tymińska, Grzegorz Lizurek
The article investigates the reliability of moment tensor (MT) inversion in time domain with use of first P-wave amplitude, a method used to determine the source mechanisms of earthquakes, across four different seismic networks. The study compares the synthetic tests results of MT inversion for two underground mining and two artificial reservoir monitoring seismic networks. The analysis was performed to assesses how consistency and accuracy of the results depend on different factors like: network configuration, events depth, velocity model, focal mechanism of event and applied noise. The findings highlight the impact of network configuration compared to other variables and data quality on the reliability of moment tensor inversion and provide insights into different factors which have to be considered to enhance MT accuracy. The significance of events depth in P-wave amplitude MT inversion and the necessity to consider velocity model influence, especially presence of high velocity gradient, is highlighted by the presented results.
文章研究了利用第一 P 波振幅进行时域矩张量(MT)反演的可靠性,这种方法用于确定四个不同地震台网的震源机制。研究比较了两个地下采矿和两个人工水库监测地震网络的矩张量反演合成测试结果。分析旨在评估结果的一致性和准确性如何取决于不同的因素,如:网络配置、事件深度、速度模型、事件的焦点机制和应用噪声。研究结果强调了与其他变量和数据质量相比,网络配置对矩形张量反演可靠性的影响,并深入分析了提高矩形张量反演准确性必须考虑的不同因素。研究结果凸显了事件深度在 P 波振幅 MT 反演中的重要性,以及考虑速度模型影响的必要性,尤其是高速梯度的存在。
{"title":"Reliability of Moment Tensor Inversion for Different Seismic Networks","authors":"Anna Tymińska, Grzegorz Lizurek","doi":"10.1007/s00024-024-03570-5","DOIUrl":"10.1007/s00024-024-03570-5","url":null,"abstract":"<div><p>The article investigates the reliability of moment tensor (MT) inversion in time domain with use of first P-wave amplitude, a method used to determine the source mechanisms of earthquakes, across four different seismic networks. The study compares the synthetic tests results of MT inversion for two underground mining and two artificial reservoir monitoring seismic networks. The analysis was performed to assesses how consistency and accuracy of the results depend on different factors like: network configuration, events depth, velocity model, focal mechanism of event and applied noise. The findings highlight the impact of network configuration compared to other variables and data quality on the reliability of moment tensor inversion and provide insights into different factors which have to be considered to enhance MT accuracy. The significance of events depth in P-wave amplitude MT inversion and the necessity to consider velocity model influence, especially presence of high velocity gradient, is highlighted by the presented results.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 9","pages":"2787 - 2800"},"PeriodicalIF":1.9,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00024-024-03570-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209908","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The aim of the present work is to evaluate the performance of several MT-InSAR techniques based on satellite SAR data in monitoring ground deformation phenomena affecting complex scenarios such as small islands of volcanic origin. To such purpose, PS, SBAS and IPTA approaches are applied in the study of Lipari, Salina and Vulcano islands belonging to the Aeolian archipelago, southern Tyrrhenian Sea, Italy. The outcomes retrieved from each technique are then discussed in terms of intrinsic features, spatial coverage, linear trend and coherence. Moreover, the accuracy of PS, SBAS and IPTA results are evaluated by comparison with in-situ measurements from the GNSS network managed by INGV-OE and private operators considering different metrics. Experimental results show that in this case there is no preferred MT-InSAR technique in an absolute way but each of them has strengths and drawbacks that have to be taken into account in the monitoring of complex scenarios.
{"title":"On the Monitoring of Small Islands Belonging to the Aeolian Archipelago by MT-InSAR Data","authors":"Marco Polcari, Mimmo Palano, Silvia Puliero, Francesca Silverii, Claudia Spinetti, Cristiano Tolomei","doi":"10.1007/s00024-024-03568-z","DOIUrl":"https://doi.org/10.1007/s00024-024-03568-z","url":null,"abstract":"<p>The aim of the present work is to evaluate the performance of several MT-InSAR techniques based on satellite SAR data in monitoring ground deformation phenomena affecting complex scenarios such as small islands of volcanic origin. To such purpose, PS, SBAS and IPTA approaches are applied in the study of Lipari, Salina and Vulcano islands belonging to the Aeolian archipelago, southern Tyrrhenian Sea, Italy. The outcomes retrieved from each technique are then discussed in terms of intrinsic features, spatial coverage, linear trend and coherence. Moreover, the accuracy of PS, SBAS and IPTA results are evaluated by comparison with in-situ measurements from the GNSS network managed by INGV-OE and private operators considering different metrics. Experimental results show that in this case there is no preferred MT-InSAR technique in an absolute way but each of them has strengths and drawbacks that have to be taken into account in the monitoring of complex scenarios.</p>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"176 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209900","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-09DOI: 10.1007/s00024-024-03567-0
Serguei G. Dobrovolski, Vladislav P. Yushkov, Irina V. Solomonova
Changes in the integral characteristics of the global water exchange, at climatic time scales, are considered as random functions (processes). “Trajectories” obtained as the results of numerical calculations on various, from 34 to 43, climate models (participating at the CMIP-6 “historical” experiment covering the period from 1850 through 2014) are taken as realizations of these processes. Temporal variations of following annually averaged parameters are studied: (1) average evaporation from the ocean surface, (2) precipitation over the ocean, (3) “effective evaporation” from the ocean (difference “evaporation minus precipitation”, on average equal to the water transport from the ocean to land), (4) precipitation over land, (5) evaporation (evapotranspiration) from the land surface, (6) “effective precipitation” over land (or “climatic runoff”: precipitation minus evaporation), and (7) river runoff. It is shown that precipitation over the ocean and evaporation from land largely suppress the monotonous trends in the mean values of evaporation from the ocean and precipitation over land, respectively, at secular time scales. At the same time, this damping does not extend to the trends of the last few decades, which may be due to a combination of a sharp increase in global temperature with explosive volcanic eruptions that preceded this period. An analysis of the time divergence in the model trajectories of each of the components of the global water exchange, as well as the very existence of such divergences, indicates an increase in the uncertainty of processes that is not associated with anthropogenic impact on the climate system.
{"title":"Stochastic Approach to the Evolution of the Global Water Cycle: Results of Historical Experiments on the CMIP-6 Models","authors":"Serguei G. Dobrovolski, Vladislav P. Yushkov, Irina V. Solomonova","doi":"10.1007/s00024-024-03567-0","DOIUrl":"10.1007/s00024-024-03567-0","url":null,"abstract":"<div><p>Changes in the integral characteristics of the global water exchange, at climatic time scales, are considered as random functions (processes). “Trajectories” obtained as the results of numerical calculations on various, from 34 to 43, climate models (participating at the CMIP-6 “historical” experiment covering the period from 1850 through 2014) are taken as realizations of these processes. Temporal variations of following annually averaged parameters are studied: (1) average evaporation from the ocean surface, (2) precipitation over the ocean, (3) “effective evaporation” from the ocean (difference “evaporation minus precipitation”, on average equal to the water transport from the ocean to land), (4) precipitation over land, (5) evaporation (evapotranspiration) from the land surface, (6) “effective precipitation” over land (or “climatic runoff”: precipitation minus evaporation), and (7) river runoff. It is shown that precipitation over the ocean and evaporation from land largely suppress the monotonous trends in the mean values of evaporation from the ocean and precipitation over land, respectively, at secular time scales. At the same time, this damping does not extend to the trends of the last few decades, which may be due to a combination of a sharp increase in global temperature with explosive volcanic eruptions that preceded this period. An analysis of the time divergence in the model trajectories of each of the components of the global water exchange, as well as the very existence of such divergences, indicates an increase in the uncertainty of processes that is not associated with anthropogenic impact on the climate system.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 9","pages":"2873 - 2893"},"PeriodicalIF":1.9,"publicationDate":"2024-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142226668","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-08DOI: 10.1007/s00024-024-03561-6
Mario González-Escobar, Carlos Simón Reyes-Martínez, Cristian A. Gallegos-Castillo, Sergio M. Arregui-Ojeda, Edgar A. Mastache-Román
The Colorado River Delta is an example of a transtensional plate boundary defined by a series of en-echelon faults along the Gulf of California. Near the gulf’s mouth, the southern segment of the Cerro Prieto transform Fault (CPF) and the listric Pangas Viejas Fault (PVF) outline Montague Basin (MB), which contrasts with the pull-apart basins along the gulf’s axis. In this context, we present a series of seismic reflection profiles and known observations to illustrate the MB subsurface architecture. The imaging findings reveal concentrated deformation within the CPF and PVF regions, which exhibit a high density of faults perpendicular to the major faults. The imaging findings reveal concentrated deformation within the CPF and PVF regions, which exhibit a high density of faults perpendicular to the major faults. At the same time, in the central part of the basin, the seismoreflectors appear largely unaltered and only slightly tilted. The MB covers an area of approximately 250 km2, with a length of around 90 km and a width of 30 km between the CPF to the east and the PVF to the west. Based on the data resolution, the basin's depth is believed to be over 5 km. The acoustic basement appears only on the west side of the basin, around 4 km deep. Other important observations include the Yurimori Fault running parallel to the PVF showing significant deformation at the basement level, as well as the Ometepec and La Salina Faults in the southwest oblique to the CPF, which could be the structural link with the Wagner Basin in the Gulf of California. Historical seismicity and one particular seismic profile indicate that deformation in the Indiviso fault occurred before the El Mayor-Cucapah earthquake (2010, Mw 7.2). Although the connection with the gulf’s basins is uncertain, a slight gravimetric high striking NE is observed; this is interpreted as a strike-slip zone, suggesting that the MB could become an overlap gap between marine and terrestrial basins. Our findings provide various controls for fault linkage in a transtensional setting in which continental and oceanic basins interact.
{"title":"Basin Style Variation Along a Transform Fault: Southern Colorado River Delta, Baja California, México","authors":"Mario González-Escobar, Carlos Simón Reyes-Martínez, Cristian A. Gallegos-Castillo, Sergio M. Arregui-Ojeda, Edgar A. Mastache-Román","doi":"10.1007/s00024-024-03561-6","DOIUrl":"https://doi.org/10.1007/s00024-024-03561-6","url":null,"abstract":"<p>The Colorado River Delta is an example of a transtensional plate boundary defined by a series of en-echelon faults along the Gulf of California. Near the gulf’s mouth, the southern segment of the Cerro Prieto transform Fault (CPF) and the listric Pangas Viejas Fault (PVF) outline Montague Basin (MB), which contrasts with the pull-apart basins along the gulf’s axis. In this context, we present a series of seismic reflection profiles and known observations to illustrate the MB subsurface architecture. The imaging findings reveal concentrated deformation within the CPF and PVF regions, which exhibit a high density of faults perpendicular to the major faults. The imaging findings reveal concentrated deformation within the CPF and PVF regions, which exhibit a high density of faults perpendicular to the major faults. At the same time, in the central part of the basin, the seismoreflectors appear largely unaltered and only slightly tilted. The MB covers an area of approximately 250 km<sup>2</sup>, with a length of around 90 km and a width of 30 km between the CPF to the east and the PVF to the west. Based on the data resolution, the basin's depth is believed to be over 5 km. The acoustic basement appears only on the west side of the basin, around 4 km deep. Other important observations include the Yurimori Fault running parallel to the PVF showing significant deformation at the basement level, as well as the Ometepec and La Salina Faults in the southwest oblique to the CPF, which could be the structural link with the Wagner Basin in the Gulf of California. Historical seismicity and one particular seismic profile indicate that deformation in the Indiviso fault occurred before the El Mayor-Cucapah earthquake (2010, Mw 7.2). Although the connection with the gulf’s basins is uncertain, a slight gravimetric high striking NE is observed; this is interpreted as a strike-slip zone, suggesting that the MB could become an overlap gap between marine and terrestrial basins. Our findings provide various controls for fault linkage in a transtensional setting in which continental and oceanic basins interact.</p>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"48 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209901","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-06DOI: 10.1007/s00024-024-03562-5
Urooj Shakir, Aamir Ali, Muyyassar Hussain, Ahmed E. Radwan, Ahmed Abd El Aal
Seismic inversion has been in use for the last two decades to measure inverted impedances using an integrated data set approach. This research focuses on the application of multi-attribute seismic inversion and the geostatistical probabilistic neural network (PNN) approach for determining rock properties and litho-fluid classification in the Mehar-Mazarani Field of the Lower Indus Basin (LIB), Pakistan. The study compares five different inversion techniques, including model-based inversion (MBI), colored inversion (CI), linear sparse spike inversion (LSSI), band-limited inversion (BLI), and maximum likelihood sparse spike inversion (MLSSI). The inverted outputs, such as acoustic P-impedance (Zp), density (ρ), porosity (φ), and shale volume (Vsh), were analyzed in Paleocene and Cretaceous geological complex reservoirs to identify gas-bearing zones. The results indicated the existence of gas between 1630 and 1700 ms (ms) and corresponding depth ranges from approximately 3200 m up to 4200 m with varying thickness. Amongst the inversion techniques, MBI demonstrated greater accuracy, with inverted density volumes showing a strong correlation coefficient of 0.98 and the lowest root mean square error (RMSE) and relative error of 0.10 m/s * g/cc. A geostatistical PNN approach was employed to estimate variations in Vsh and φ within the sand reservoir. MBI again yielded more reliable results, with a strong correlation between the measured and inverted attributes. High φ and low Vsh were observed in predetermined low-impedance zones. Overall, MBI is proven to be the most accurate and reliable technique, providing clear identification of the gas occurrence. This research highlights the effectiveness of seismic inversion, particularly the application of MBI, in determining rock properties and identifying gas-bearing zones within the Mehar-Mazarani gas field.
{"title":"PNN Enhanced Seismic Inversion for Porosity Modeling and Delineating the Potential Heterogeneous Gas Sands via Comparative Inversion Analysis in the Lower Indus Basin","authors":"Urooj Shakir, Aamir Ali, Muyyassar Hussain, Ahmed E. Radwan, Ahmed Abd El Aal","doi":"10.1007/s00024-024-03562-5","DOIUrl":"10.1007/s00024-024-03562-5","url":null,"abstract":"<div><p>Seismic inversion has been in use for the last two decades to measure inverted impedances using an integrated data set approach. This research focuses on the application of multi-attribute seismic inversion and the geostatistical <i>probabilistic neural network</i> (PNN) approach for determining rock properties and litho-fluid classification in the Mehar-Mazarani Field of the Lower Indus Basin (<i>LIB</i>), Pakistan. The study compares five different inversion techniques, including <i>model-based inversion</i> (MBI), <i>colored inversion</i> (CI), <i>linear sparse spike inversion</i> (LSSI), <i>band-limited inversion</i> (BLI), and <i>maximum likelihood sparse spike inversion</i> (MLSSI). The inverted outputs, such as acoustic <i>P-impedance</i> (Zp), <i>density</i> (ρ), <i>porosity</i> (φ), and <i>shale volume</i> (Vsh), were analyzed in Paleocene and Cretaceous geological complex reservoirs to identify gas-bearing zones. The results indicated the existence of gas between 1630 and 1700 <i>ms</i> (ms) and corresponding depth ranges from approximately 3200 m up to 4200 m with varying thickness. Amongst the inversion techniques, MBI demonstrated greater accuracy, with inverted density volumes showing a strong correlation coefficient of 0.98 and the lowest <i>root mean square error</i> (RMSE) and relative error of 0.10 m/s * g/cc. A geostatistical PNN approach was employed to estimate variations in Vsh and φ within the sand reservoir. MBI again yielded more reliable results, with a strong correlation between the measured and inverted attributes. High φ and low Vsh were observed in predetermined low-impedance zones. Overall, MBI is proven to be the most accurate and reliable technique, providing clear identification of the gas occurrence. This research highlights the effectiveness of seismic inversion, particularly the application of MBI, in determining rock properties and identifying gas-bearing zones within the Mehar-Mazarani gas field.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 9","pages":"2801 - 2821"},"PeriodicalIF":1.9,"publicationDate":"2024-09-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00024-024-03562-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209903","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-05DOI: 10.1007/s00024-024-03566-1
Wei Liu, Zhixing Zhai, Zhou Fang
Time-frequency analysis (TFA) technique has a powerful capacity to characterize non-stationary signals. In this paper, a highly energy concentrated TFA method, called multisynchrosqueezing-based S-transform (MSSBST), is proposed for the analysis of seismic data. Herein, we combine S-transform (ST) and multisynchrosqueezing framework by making full use of an iterative reassignment procedure to concentrate time-frequency energy in a stepwise manner. Furthermore, we derive a series of formulas about MSSBST and its inverse transform, which means that the MSSBST allows for signal reconstruction from its time-frequency representation (TFR). The numerical analysis shows that the proposed method not only can effectively enhance the time-frequency energy concentration but also can offer better performance in characterizing non-stationary signals compared with the short-time Fourier transform (STFT), ST and synchrosqueezing S-transform (SSST). Field examples further demonstrate its potential in depicting spectral anomalies related to hydrocarbon reservoir, thus, facilitating seismic interpretation.
时频分析(TFA)技术具有描述非稳态信号特征的强大能力。本文提出了一种能量高度集中的时频分析方法,即基于多同步阙值的 S 变换(MSSBST),用于分析地震数据。在本文中,我们结合了 S 变换(ST)和多同步挤压框架,充分利用迭代重分配程序,逐步集中时频能量。此外,我们还推导出一系列有关 MSSBST 及其逆变换的公式,这意味着 MSSBST 可以从其时频表示(TFR)中重建信号。数值分析表明,与短时傅里叶变换(STFT)、同步傅里叶变换(ST)和同步 S 变换(SSST)相比,所提出的方法不仅能有效提高时频能量集中度,而且在表征非稳态信号方面具有更好的性能。现场实例进一步证明了它在描述与油气储层相关的频谱异常方面的潜力,从而为地震解释提供了便利。
{"title":"A Multisynchrosqueezing-Based S-Transform for Time-Frequency Analysis of Seismic Data","authors":"Wei Liu, Zhixing Zhai, Zhou Fang","doi":"10.1007/s00024-024-03566-1","DOIUrl":"https://doi.org/10.1007/s00024-024-03566-1","url":null,"abstract":"<p>Time-frequency analysis (TFA) technique has a powerful capacity to characterize non-stationary signals. In this paper, a highly energy concentrated TFA method, called multisynchrosqueezing-based S-transform (MSSBST), is proposed for the analysis of seismic data. Herein, we combine S-transform (ST) and multisynchrosqueezing framework by making full use of an iterative reassignment procedure to concentrate time-frequency energy in a stepwise manner. Furthermore, we derive a series of formulas about MSSBST and its inverse transform, which means that the MSSBST allows for signal reconstruction from its time-frequency representation (TFR). The numerical analysis shows that the proposed method not only can effectively enhance the time-frequency energy concentration but also can offer better performance in characterizing non-stationary signals compared with the short-time Fourier transform (STFT), ST and synchrosqueezing S-transform (SSST). Field examples further demonstrate its potential in depicting spectral anomalies related to hydrocarbon reservoir, thus, facilitating seismic interpretation.</p>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"4 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-09-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209904","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-05DOI: 10.1007/s00024-024-03555-4
U. Riccardi, T. Pivetta, A. Fedele, G. Ricciardi, S. Carlino
We present the results of nearly one year of gravity recording acquired at the active caldera of Campi Flegrei (CFc). CFc is one of the three active volcanoes in the Neapolitan area (southern Italy) and is currently the most active one. In fact, the CFc is undergoing a period of unrest characterised by slow uplift of the ground, a peculiar phenomenon known worldwide as bradyseism, accompanied by seismicity and intense fumarolic emissions. Due to the increased intensity of the volcano dynamics, a permanent gravity station equipped with a gPhoneX spring gravimeter was installed to enhance the geodetic monitoring programmes. The purpose of the continuous recordings is to complement the time-lapse observations carried out periodically on networks of benchmarks, in order to continuously monitor the short-term gravity signals. We report on the various processing steps and analyses performed to obtain reliable parameters of the Earth's tides, non-tidal corrections and gravity residuals. The various methodologies employed to investigate the instrumental drift are also elucidated in depth, because it may masquerade the elusive gravity changes resulting from mass fluctuations within the volcanic and geothermal systems. Residual gravity signals, retrieved from the recordings, after reduction of body and ocean tides, polynomial drift, atmospheric, tilt and change of the Earth Orientation Parameters (EOP) effects, appear to be uncorrelated with hydrology, while they show a clear correlation with the most energetic earthquakes, that strongly characterise the current bradyseismic crisis. The residual gravity signals display peculiar trends characterized by steps or offsets (up to about 600 nm/s2) and transients in coincidence of the most energetic volcano-tectonic events (Magnitude > 2.5) and seismic swarms. The steps in the gravity residuals are likely to be of instrumental origin, while the amplitudes of the observed transients are not consistent with co-seismic or volcanological phenomena, for which there is no evidence from other monitoring techniques. Unfortunately, the lack of repeated absolute gravity measurements severely limits our ability to attribute the observed gravity variations to geological sources. From the analysis of the gravity records, reliable tidal gravity models have been derived, which will improve the accuracy of volcano monitoring by allowing a precise reduction of tidal effects for both relative and absolute gravity measurements taken in these volcanic areas.
{"title":"Continuous Gravity Observations at Campi Flegrei Caldera: An Accurate Assessment of Tidal and Non-Tidal Signals and Implications for Volcano Monitoring","authors":"U. Riccardi, T. Pivetta, A. Fedele, G. Ricciardi, S. Carlino","doi":"10.1007/s00024-024-03555-4","DOIUrl":"https://doi.org/10.1007/s00024-024-03555-4","url":null,"abstract":"<p>We present the results of nearly one year of gravity recording acquired at the active caldera of Campi Flegrei (CFc). CFc is one of the three active volcanoes in the Neapolitan area (southern Italy) and is currently the most active one. In fact, the CFc is undergoing a period of unrest characterised by slow uplift of the ground, a peculiar phenomenon known worldwide as <i>bradyseism</i>, accompanied by seismicity and intense fumarolic emissions. Due to the increased intensity of the volcano dynamics, a permanent gravity station equipped with a gPhoneX spring gravimeter was installed to enhance the geodetic monitoring programmes. The purpose of the continuous recordings is to complement the time-lapse observations carried out periodically on networks of benchmarks, in order to continuously monitor the short-term gravity signals. We report on the various processing steps and analyses performed to obtain reliable parameters of the Earth's tides, non-tidal corrections and gravity residuals. The various methodologies employed to investigate the instrumental drift are also elucidated in depth, because it may masquerade the elusive gravity changes resulting from mass fluctuations within the volcanic and geothermal systems. Residual gravity signals, retrieved from the recordings, after reduction of body and ocean tides, polynomial drift, atmospheric, tilt and change of the Earth Orientation Parameters (EOP) effects, appear to be uncorrelated with hydrology, while they show a clear correlation with the most energetic earthquakes, that strongly characterise the current bradyseismic crisis. The residual gravity signals display peculiar trends characterized by steps or offsets (up to about 600 nm/s<sup>2</sup>) and transients in coincidence of the most energetic volcano-tectonic events (Magnitude > 2.5) and seismic swarms. The steps in the gravity residuals are likely to be of instrumental origin, while the amplitudes of the observed transients are not consistent with co-seismic or volcanological phenomena, for which there is no evidence from other monitoring techniques. Unfortunately, the lack of repeated absolute gravity measurements severely limits our ability to attribute the observed gravity variations to geological sources. From the analysis of the gravity records, reliable tidal gravity models have been derived, which will improve the accuracy of volcano monitoring by allowing a precise reduction of tidal effects for both relative and absolute gravity measurements taken in these volcanic areas.</p>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"21 1","pages":""},"PeriodicalIF":2.0,"publicationDate":"2024-09-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209902","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-04DOI: 10.1007/s00024-024-03564-3
A. Sabarinath, T. Kesavavarthini, Meera M. Nair, A. Naga Rajesh
Surface air temperature (SAT) affects both natural systems and human activities, impacting health, agriculture, energy demand, and so on. To investigate and analyze SAT over the region of interest, it is crucial to choose suitable climate models. The study commenced with the evaluation of 42 Coupled Model Intercomparison Project phase 6 (CMIP6) models’ simulations of SAT over India for annual and all four seasons (summer, southwest monsoon, northeast monsoon, and winter) during the historical period 1985 to 2014 with respect to the gridded SAT datasets obtained from the India Meteorological Department (IMD). Multi Model Mean (MMM) of 42 models was included in the evaluation. The evaluation was performed with various statistical metrics such as root mean squared error (RMSE), mean bias error (MBE), correlation coefficient (R), mean squared error (MAE), Taylor skill score (TSS), Brier skill score (BSS), and Interannual variability skill score (IVSS). By the method of estimating Comprehensive Rating Index (CRI), the top-ranking models were identified to be CMCC-CM2-SR5 for the annual and summer season, MIROC6 for the winter season, ACCESS-ESM-1-5 for the southwest monsoon, and NorESM2-LM for northeast monsoon. The novelty of this study lies in the approach of identifying the best ensemble. For each season, statistical metric-wise top-ranked models were picked to develop the best ensemble. Again, the overall ranking of the models along with the best ensemble for each season is determined by estimating CRI. It was observed that for all seasons, the best ensemble falls within the top 3 models’ category. Future projections of SAT under four shared socio-economic pathways (SSP-2.6, 4.5, 7.0, and 8.5) were also analyzed with the best ensemble obtained for each season. The results convey that, the country will witness, especially during the summer season, there will be a 1.160 °C, 1.288 °C and 2.368 °C rise in the mean SAT between historical (1985–2014) and near future (2021–2040), near and mid future (2041–2060), mid and far future (2081–2100) if the pathway, SSP5-8.5 is followed.
{"title":"Assessment of Annual and Seasonal Surface Air Temperature Simulations in CMIP6 Models over India","authors":"A. Sabarinath, T. Kesavavarthini, Meera M. Nair, A. Naga Rajesh","doi":"10.1007/s00024-024-03564-3","DOIUrl":"10.1007/s00024-024-03564-3","url":null,"abstract":"<div><p>Surface air temperature (SAT) affects both natural systems and human activities, impacting health, agriculture, energy demand, and so on. To investigate and analyze SAT over the region of interest, it is crucial to choose suitable climate models. The study commenced with the evaluation of 42 Coupled Model Intercomparison Project phase 6 (CMIP6) models’ simulations of SAT over India for annual and all four seasons (summer, southwest monsoon, northeast monsoon, and winter) during the historical period 1985 to 2014 with respect to the gridded SAT datasets obtained from the India Meteorological Department (IMD). Multi Model Mean (MMM) of 42 models was included in the evaluation. The evaluation was performed with various statistical metrics such as root mean squared error (RMSE), mean bias error (MBE), correlation coefficient (R), mean squared error (MAE), Taylor skill score (TSS), Brier skill score (BSS), and Interannual variability skill score (IVSS). By the method of estimating Comprehensive Rating Index (CRI), the top-ranking models were identified to be CMCC-CM2-SR5 for the annual and summer season, MIROC6 for the winter season, ACCESS-ESM-1-5 for the southwest monsoon, and NorESM2-LM for northeast monsoon. The novelty of this study lies in the approach of identifying the best ensemble. For each season, statistical metric-wise top-ranked models were picked to develop the best ensemble. Again, the overall ranking of the models along with the best ensemble for each season is determined by estimating CRI. It was observed that for all seasons, the best ensemble falls within the top 3 models’ category. Future projections of SAT under four shared socio-economic pathways (SSP-2.6, 4.5, 7.0, and 8.5) were also analyzed with the best ensemble obtained for each season. The results convey that, the country will witness, especially during the summer season, there will be a 1.160 °C, 1.288 °C and 2.368 °C rise in the mean SAT between historical (1985–2014) and near future (2021–2040), near and mid future (2041–2060), mid and far future (2081–2100) if the pathway, SSP5-8.5 is followed.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 9","pages":"2949 - 2971"},"PeriodicalIF":1.9,"publicationDate":"2024-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209905","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-08-31DOI: 10.1007/s00024-024-03563-4
Bappa Mukherjee, Sohan Kar, Kalachand Sain
In the E&P industry, accurate lithology classification is an essential task for successful exploration and production. Geophysical logs provide high-resolution petrophysical properties, but core logging is expensive and traditional techniques may not accurately classify lithologies. We demonstrated a comparative analysis of six ML algorithms: k-nearest neighbor (kNN), support vector machine (SVM), decision tree (DT), random forest (RF), extreme gradient boosting (XGBoost) and artificial neural network (ANN) for the prediction of lithologies from geophysical logs. Here we analysed the wireline logs of eight wells associated with the petroliferous Lakadong-Therria formation of the Bhogpara oil field of the Assam-Arakan Basin. This formation contains eight typical lithologies, namely clay stone, sand stone, calcareous sandstone, shale, calcareous shale, carbonaceous shale, coal and limestone. Performance of the ML algorithms were evaluated through accuracy, precision, recall, F1-score and receiver operating characteristic (ROC) curve. During the training and test phases, the computed overall accuracy of the predicted ML modes exceeded 82% and 71%, respectively. The model accuracy hierarchy was ANN > XGBoost > RF > SVM > DT > kNN during training, and ANN/XGBoost > kNN > DT/RF > SVM during testing. This approach allows interpreters to select the most accurate ML model based on training phase performance. This study provided a clear insight towards generating a supplement for litholog sequence and improving the accuracy and efficiency of lithology prediction in a geologically complex petroleum reservoir using pre-received core derived litholog information at few wells.
在勘探和开发行业,准确的岩性分类是成功勘探和生产的一项基本任务。地球物理测井可提供高分辨率的岩石物理特性,但岩心测井费用昂贵,而且传统技术可能无法对岩性进行准确分类。我们展示了六种 ML 算法的比较分析:k-近邻(kNN)、支持向量机(SVM)、决策树(DT)、随机森林(RF)、极梯度提升(XGBoost)和人工神经网络(ANN),用于从地球物理测井记录预测岩性。在此,我们分析了与阿萨姆-阿拉干盆地博格帕拉油田含油层拉卡东-特里亚地层相关的八口油井的有线测井记录。该地层包含八种典型岩性,即粘土岩、砂岩、钙质砂岩、页岩、钙质页岩、碳质页岩、煤和石灰岩。通过准确度、精确度、召回率、F1-分数和接收者操作特征曲线(ROC)对 ML 算法的性能进行了评估。在训练和测试阶段,计算得出的 ML 模式预测总体准确率分别超过 82% 和 71%。在训练阶段,模型准确率等级为 ANN > XGBoost > RF > SVM > DT > kNN;在测试阶段,模型准确率等级为 ANN/XGBoost > kNN > DT/RF > SVM。这种方法允许解释人员根据训练阶段的表现选择最准确的 ML 模型。这项研究提供了一个清晰的视角,有助于为岩性序列提供补充,并利用在少数油井中预先接收的岩心衍生岩性信息,提高地质复杂的石油储层中岩性预测的准确性和效率。
{"title":"Machine Learning Assisted State-of-the-Art-of Petrographic Classification From Geophysical Logs","authors":"Bappa Mukherjee, Sohan Kar, Kalachand Sain","doi":"10.1007/s00024-024-03563-4","DOIUrl":"10.1007/s00024-024-03563-4","url":null,"abstract":"<div><p>In the E&P industry, accurate lithology classification is an essential task for successful exploration and production. Geophysical logs provide high-resolution petrophysical properties, but core logging is expensive and traditional techniques may not accurately classify lithologies. We demonstrated a comparative analysis of six ML algorithms: k-nearest neighbor (kNN), support vector machine (SVM), decision tree (DT), random forest (RF), extreme gradient boosting (XGBoost) and artificial neural network (ANN) for the prediction of lithologies from geophysical logs. Here we analysed the wireline logs of eight wells associated with the petroliferous Lakadong-Therria formation of the Bhogpara oil field of the Assam-Arakan Basin. This formation contains eight typical lithologies, namely clay stone, sand stone, calcareous sandstone, shale, calcareous shale, carbonaceous shale, coal and limestone. Performance of the ML algorithms were evaluated through accuracy, precision, recall, F1-score and receiver operating characteristic (ROC) curve. During the training and test phases, the computed overall accuracy of the predicted ML modes exceeded 82% and 71%, respectively. The model accuracy hierarchy was ANN > XGBoost > RF > SVM > DT > kNN during training, and ANN/XGBoost > kNN > DT/RF > SVM during testing. This approach allows interpreters to select the most accurate ML model based on training phase performance. This study provided a clear insight towards generating a supplement for litholog sequence and improving the accuracy and efficiency of lithology prediction in a geologically complex petroleum reservoir using pre-received core derived litholog information at few wells.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 9","pages":"2839 - 2871"},"PeriodicalIF":1.9,"publicationDate":"2024-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142209907","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}