Geophysical Prospecting最新文献

{"title":"Regularization with differential evolution-based common-offset common-reflection surface: A case study for field onshore seismic data","authors":"Tainá Souza,&nbsp;Tiago Barros,&nbsp;Renato Lopes","doi":"10.1111/1365-2478.70007","DOIUrl":"https://doi.org/10.1111/1365-2478.70007","url":null,"abstract":"<p>We introduce a novel approach for seismic pre-stack data regularization using the common-offset common-reflection-surface method with a global attribute search strategy, employing a bio-inspired differential evolution optimization algorithm. We compare the global common-offset common-reflection-surface approach with the conventional sequential attribute search, focusing on their performance in pre-stack data regularization. Tests on synthetic and onshore (field) seismic datasets demonstrate that the global search approach significantly improves performance, enhancing signal-to-noise ratio and coherently filling missing traces. We show that the global common-offset common-reflection-surface method effectively addresses spatial irregularities, reconstructing reflections without artefacts, filling data gaps and highlighting geological details even in complex areas. In contrast, the sequential common-offset common-reflection-surface method, while capable of reconstructing missing traces, shows lower interpolation quality and fails to adequately highlight complex geological features such as high-dip reflections.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1411-1430"},"PeriodicalIF":1.8,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950278","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Three-dimensional forward modelling and inversion of magnetic anomaly using finite-infinite element coupling method with Comsol: A case study of Pingchuan iron-ore deposit, Southwest China","authors":"Yinzhong Ma,&nbsp;Xiangyun Hu,&nbsp;Shuang Liu","doi":"10.1111/1365-2478.70010","DOIUrl":"https://doi.org/10.1111/1365-2478.70010","url":null,"abstract":"<p>We implemented the finite-infinite element coupling method in Comsol for three-dimensional magnetic forward modelling in geophysics exploration. The finite-infinite element coupling method optimizes boundary conditions by combining finite and infinite elements, and its boundary and mesh convergence were validated against analytical solutions. For non-ellipsoidal magnetic bodies without analytical solutions, we compared the forward magnetic anomaly solutions of a cube using the finite-infinite element coupling method and finite volume method with hexahedral meshes. The finite-infinite element coupling method demonstrated higher accuracy compared to the finite volume method under the presence of remanent magnetization, highlighting the impact of spatial discretization on solution quality. Simulation results from the rectangular prism and composite models showed that the finite-infinite element coupling method, through unstructured meshes, enables precise modelling of complex geometries and non-uniform magnetic field distributions. Using Comsol Multiphysics, we successfully applied the finite-infinite element coupling method for the three-dimensional inversion of ore body models and analysed data from the Pingchuan iron ore deposit in Southwest China. This study illustrates the applicability of the finite-infinite element coupling method in advancing magnetic modelling and the effectiveness of Comsol as a tool for characterizing complex geological structures.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1396-1410"},"PeriodicalIF":1.8,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950323","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Group velocity image of singularities for converted reflected and transmitted waves in orthorhombic anisotropic media","authors":"Alexey Stovas,&nbsp;Yuriy Roganov,&nbsp;Vyacheslav Roganov","doi":"10.1111/1365-2478.70005","DOIUrl":"https://doi.org/10.1111/1365-2478.70005","url":null,"abstract":"<p>We define double singularity points, lines and surfaces and their group velocity images for converted reflected and transmitted PS1 and PS2 waves in layered orthorhombic media. This is performed by considering the dominant term in Taylor series for Christoffel equation defined for converted waves. We show that the group velocity image for converted wave has the same functional form as the one for pure wave modes (ellipse in 3D). The position and area of the ellipse in the group velocity domain and offset domain are also analysed. We also expand this method for two-layer model with S1S2 wave singularity in both layers. It results in quartic (S1S1, S2S2, S1S2 and S2S1 waves) characteristic equation. In the group velocity domain, it gives two close curves: quasi-elliptical (for pure S1S1 and S2S2 waves) and very complex curve with deflection points (for converted S1S2 and S2S1 waves).</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1364-1378"},"PeriodicalIF":1.8,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950257","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Slope assisted Physics-informed neural networks for seismic signal separation with applications on ground roll removal and interpolation","authors":"Francesco Brandolin,&nbsp;Matteo Ravasi,&nbsp;Tariq Alkhalifah","doi":"10.1111/1365-2478.70004","DOIUrl":"https://doi.org/10.1111/1365-2478.70004","url":null,"abstract":"<p>The knowledge of the local slope field of prestack seismic data is essential in several seismic signal processing tasks. Building on our previous slope-assisted, physics-informed seismic interpolation framework, dubbed PINNslope, we introduce a series of enhancements that elevate the framework's versatility. This ultimately enables its application to different signal separation problems, with a specific focus on ground roll removal. To begin with, the local slope estimated using our physics-informed neural networks framework is compared against the analytical local slope and those obtained from several conventional slope estimation algorithms. This comparison showcases that our prediction better approximates the analytical one. Second, we use the derived relation between the analytical slope and the physics-informed neural networks estimated slope to constrain the slope estimation network in the ground roll removal problem, predicting only the clean reflections and avoiding the prediction of the ground roll. To address the large difference in the frequency content of the field data, we utilize a time derivative term in the loss function to emphasize the amplitude of the comparatively higher frequency reflection arrivals. Furthermore, we modify the framework loss function and architecture to demonstrate the possibility of predicting two separate components of the seismic data according to the estimate of two local slopes that can span opposite or different ranges of values between each other. The effectiveness of the double slope framework is demonstrated on a proof of concept of the deblending problem and for the interpolation of complex aliased data characterized by conflicting dips, two tasks that were not achievable using our single slope prediction network implementation.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1337-1363"},"PeriodicalIF":1.8,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70004","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950256","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Ore-controlling structures revealed by magnetic and resistivity anomalies in the Linares–La Carolina mining district (Spain)","authors":"J. Rey,&nbsp;R. Mendoza,&nbsp;J. Nuñez,&nbsp;F. J. Martínez-Moreno,&nbsp;M. C. Hidalgo","doi":"10.1111/1365-2478.70019","DOIUrl":"https://doi.org/10.1111/1365-2478.70019","url":null,"abstract":"<p>The effectiveness of magnetic and electromagnetic geophysical methods in delineating geological structures with mining implications was analysed in this study, specifically for structures that affect the possible extension of the lead-based metal phylonian deposits concealed by young sedimentary cover. This study focuses on the ancient mining district of Linares–La Carolina (south-eastern Spain), characterized by the presence of galena (PbS) deposits. In this sector, two of its main mining concessions (Arrayanes and San Miguel–El Mimbre) were selected, which lose continuity under the Plio–Quaternary basin cover of La Garza. The magnetic surveys reveal the presence of anomalies aligned in two directions, NW‒SE and SW‒NE. Calculated magnetic models indicate that these anomalies could be associated with faults. The continuity of the veins is disrupted or displaced when these anomalies are reached. Therefore, the study of these linear magnetic trends can help to locate these faults that control the continuity and displacement of the veins. Joint interpretation of new time-domain electromagnetic stations and electrical resistivity tomography profiles previously obtained reveal linear features that align with the faults interpreted from the magnetic results. This consistency between methods allowed us to corroborate the presence of these structures. The synergistic application of magnetic techniques alongside electrical or electromagnetic methods proves effective in identifying faults in the Arrayanes and San Miguel–El Mimbre mining concessions. The identified faults are significant for prospecting for mineral deposits as they are spatially associated with both vein disappearances and lateral changes in facies and thickness in the Quaternary fill of the La Garza sedimentary basin. This is evidence that the geophysically mapped faults have meaningful control on ore vein distribution.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 5","pages":"1379-1395"},"PeriodicalIF":1.8,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950258","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Discretization of small-scale, stratigraphic heterogeneities and its impact on the seismic response: Lessons from the application of process-based modelling","authors":"Andrea Cuesta Cano,&nbsp;Azin Karimzadanzabi,&nbsp;Joep Elisabeth Anton Storms,&nbsp;Guillaume Rongier,&nbsp;Dirk Jacob Verschuur,&nbsp;Allard Willem Martinius","doi":"10.1111/1365-2478.70015","DOIUrl":"https://doi.org/10.1111/1365-2478.70015","url":null,"abstract":"&lt;p&gt;Reducing the uncertainty of reservoir characterization requires to better identify the small-scale structures of the subsurface from the available data. Studying the seismic response of meter-scale, stratigraphic heterogeneities typically relies on the generation of reservoir models based on outcrop examples and their forward seismic modelling. To bridge geological information and seismic modelling, these methods allocate values of acoustic properties, such as mass-density and P-wave velocity, according to discretized properties like layer-type lithology or facies units. This strategy matches the current workflow in seismic data inversion in industry, where modelling workflows are based on lithofacies distributions. However, from stratigraphic modelling, we know that meter-scale heterogeneities occur within certain facies and lithologies. Here, we evaluate the difference on the seismic response between allocating acoustic properties in a grain size–based, semi-continuous manner versus discretized manners based on lithology and facies classifications. To do so, we generate a reference geological simulation that we populate with acoustic properties, mass-density and P-wave velocity, using three different strategies: (1) based on grain size distribution; (2) based on facies distribution; and (3) based on lithology. The method we propose includes the generation of realistic geological simulations based on stratigraphic modelling and the transformation of its output into acoustic properties, honouring the intra-lithology and intra-facies, small-scale structures. We, then, generate seismic data by applying a forward seismic modelling workflow. The synthetic data show that the grain size–based simulation allows the identification of small-scale, stratigraphic heterogeneities, such as beds with strong density and velocity contrasts. These stratigraphic structures are smoothened or may completely disappear in the facies and lithology discretized simulations and, therefore, are not (well) represented in the synthetic seismic data. Recognizing meter-scale, stratigraphic heterogeneities is relevant for the characterization of the fluid flow in the reservoir. However, current discrete and lithology-based strategies in seismic inversion are not able to resolve such heterogeneities because real subsurface properties are not discrete properties but continuous, unless there are stratigraphic discontinuities such as erosional surfaces or faults. This research works towards a better understanding of the relationship between changes in these continuous properties and the observed seismic data by introducing greater complexity into the discretized geological simulations. Here, we use synthetic seismic images with the goal of eventually aiding in fine-tuning seismic inversion methodologies applied to real seismic data. One pathway is to foster the development of inversion approaches that can leverage stratigraphic modelling to get stronger geological priors and rep","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 4","pages":"1280-1300"},"PeriodicalIF":1.8,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70015","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143845835","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Three-dimensional gravity inversion based on the Hamiltonian Monte Carlo method","authors":"Ya Xu,&nbsp;Wei Chu,&nbsp;Song Huang,&nbsp;Rui Guo,&nbsp;Shupeng Lu,&nbsp;Fangzhou Nan","doi":"10.1111/1365-2478.70016","DOIUrl":"https://doi.org/10.1111/1365-2478.70016","url":null,"abstract":"<p>In geophysics, Bayesian inversion methods are of significant prominence. Here, we present a novel approach utilizing the Hamiltonian Monte Carlo (HMC) method in gravity inversion for elucidating three-dimensional (3D) density structures. HMC provides a multi-dimensional sampling method that demonstrates enhanced optimization efficiency, facilitating the attainment of distant proposals with elevated acceptance probabilities. Its applicability also extends to resolving linear inverse problems. Three synthetic models of cubic bodies, dipping dykes and a combined model were designed for tests. The testes underscore the promising potential of HMC in recovering subsurface density source bodies and giving the uncertainty of the inversion model. Furthermore, an inversion test conducted on the Vinton salt dome yields a reasonable 3D distribution of cap rock, consistent with prior studies in this area. The modelling and field experiments showed that the proposed HMC gravity inversion method had higher accuracy and application potential.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 4","pages":"1301-1314"},"PeriodicalIF":1.8,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143845836","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Sequential multi-dimensional parameter inversion of induction logging data","authors":"Durra H. Saputera,&nbsp;Morten Jakobsen,&nbsp;K. W. A. van Dongen,&nbsp;Nazanin Jahani","doi":"10.1111/1365-2478.70017","DOIUrl":"https://doi.org/10.1111/1365-2478.70017","url":null,"abstract":"<p>Structural information about the subsurface near the borehole can be obtained from reconstructed conductivity distributions. These distributions may be reconstructed via the inversion of deep-sensing electromagnetic induction log data. Unfortunately, these complex media often display anisotropy and structural variations in both horizontal and vertical directions, making the three-dimensional inversion computationally demanding and ill-posed. To address these challenges, we introduce a sequential inversion strategy of deep-sensing electromagnetic induction logging data that is measured while drilling. For the inversion at each logging position, we employ a matrix-free implementation of the adjoint integral equation method and a quasi-Newton algorithm. To tackle the ill-posed nature of the problem, we regularize the inverse problem by employing a multi-dimensional inversion parameter technique that shifts from zero- to three-dimensional parameterization. The model derived from the inversion of the data at multiple positions is incrementally integrated by utilizing the sensitivity data at each logging position. To validate our approach, we tested our method on simulated data using an anisotropic model. These experiments show that this approach produces a good reconstruction of the true conductivity for the whole track while only doing the inversion at a single position at a time.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 4","pages":"1315-1332"},"PeriodicalIF":1.8,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143845834","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multi-synchrosqueezing polynomial chirplet transform and its application to seismic thin interbeds analysis","authors":"Fang Li,&nbsp;Hui Chen,&nbsp;Yuanwei Song,&nbsp;Rui Wang,&nbsp;Bairong Ding","doi":"10.1111/1365-2478.70014","DOIUrl":"https://doi.org/10.1111/1365-2478.70014","url":null,"abstract":"<p>Thin interbeds are important geological structures in seismic exploration, and the analysis of their thickness variations is a key step in seismic interpretation. As a typical signal-processing technology, the time–frequency analysis method maps one-dimensional signals to the two-dimensional time–frequency domain, which can capture the changes of the instantaneous frequency. On this basis, the time–frequency analysis method can be used to analyse the thickness changes of thin interbeds. To analyse the thickness changes more accurately, the adopted time–frequency analysis method needs to have high time–frequency resolution and robustness. This paper proposes a novel method named the multi-synchrosqueezing polynomial chirplet transform. First, the second-order instantaneous frequency estimation operator is obtained through the corrected polynomial chirplet transform. Then, through squeezing and rearranging, the fuzzy time–frequency energy is gradually concentrated on the corresponding second-order multiple instantaneous frequency estimation operator to obtain the multi-synchrosqueezing polynomial chirplet transform. Simulated signals are used to demonstrate that multi-synchrosqueezing polynomial chirplet transform has robustness while improving time–frequency resolution. Simultaneously, simulated and real seismic signals are used to verify that the proposed method can analyse the thickness variation of thin interbeds.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 4","pages":"1268-1279"},"PeriodicalIF":1.8,"publicationDate":"2025-03-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143845986","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Evaluating the applicability of a 2D weighted compact gravity inversion method for determining crustal undulations and thicknesses in the Aegean region","authors":"M. B. Doğan","doi":"10.1111/1365-2478.70001","DOIUrl":"https://doi.org/10.1111/1365-2478.70001","url":null,"abstract":"<p>This study introduces the application of the 2D weighted compact gravity inversion technique to model crustal thickness and intracrustal discontinuities in the Aegean region, encompassing both marine and terrestrial areas over a significant area of 430 km × 333 km. The method utilized advanced spectral analysis and upward continuation techniques to enhance the quality of Bouguer data, effectively mitigating surface noise. The findings reveal a remarkable correlation (over 99.9%) between observed and theoretical data, demonstrating the algorithm's robustness in accurately delineating crustal features. Depth estimates for the Conrad (2.7 g/cm³) and Moho (3.3 g/cm³) discontinuities were obtained, highlighting distinct density variations across discontinuity zones. Furthermore, the relationship between intracrustal discontinuities and seismicity was examined, revealing that earthquakes predominantly follow the Conrad boundary. Notably, this study uniquely produces 2D depth contour maps of Conrad and Moho discontinuities from specific density derived from gravity inversion sections. The results indicate that this method is a valuable tool for understanding crustal dynamics, suggesting potential applications for future tectonic assessments, especially in regional studies. The successful application of this novel technique emphasizes its significance in advancing geophysical modelling and enhances our understanding of isostatic evaluation in the Aegean region for further studies.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 4","pages":"1243-1254"},"PeriodicalIF":1.8,"publicationDate":"2025-03-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.70001","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143846038","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}