Geophysical Prospecting最新文献

{"title":"Temperature effects on the electrical conductivity of K-feldspar","authors":"Supti Sadhukhan,&nbsp;Tapati Dutta","doi":"10.1111/1365-2478.13605","DOIUrl":"10.1111/1365-2478.13605","url":null,"abstract":"<p>K-feldspar, which constitutes about 60<span></span><math>\u0000 <semantics>\u0000 <mo>%</mo>\u0000 <annotation>$%$</annotation>\u0000 </semantics></math> of the Earth's crust, is crucial for understanding electrical conductivity in porous rocks. Its electrical properties are vital for applications in ceramics, electrical insulation and conductive polymers. In this work, we study the time evolution of electrical conductivity of K-feldspar-rich rocks with varying temperatures, at high and low pH, which has been studied through simulation using time domain random walk. Random walkers, mimicking ions in transport, move in accordance with appropriate hydrodynamic equations, dissolution and precipitation kinetics. Electrical conductivity has been calculated considering variations in the parameters of temperature, fluid pH and the abundance of K-feldspar in rocks. Electrical conductivity is found to increase with temperature up to a critical value, after which it decreases. The sharpness of the rise and fall in electrical conductivity is quantified through a measure defined as the conductivity quality factor <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>Q</mi>\u0000 <mi>σ</mi>\u0000 </msub>\u0000 <annotation>$Q_{sigma }$</annotation>\u0000 </semantics></math>. We find that <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>Q</mi>\u0000 <mi>σ</mi>\u0000 </msub>\u0000 <annotation>$Q_{sigma }$</annotation>\u0000 </semantics></math> increases with a decrease in the availability of K-feldspar mineral. Our simulated results of electrical conductivity show a good match with the experimental trends reported.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3338-3349"},"PeriodicalIF":1.8,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222279","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":"An approach of 2D convolutional neural network–based seismic data fault interpretation with linear annotation and pixel thinking","authors":"Bowen Deng,&nbsp;Guangui Zou,&nbsp;Suping Peng,&nbsp;Jiasheng She,&nbsp;Chengyang Han,&nbsp;Yanhai Liu","doi":"10.1111/1365-2478.13606","DOIUrl":"10.1111/1365-2478.13606","url":null,"abstract":"<p>This article introduces a novel method for geological fault interpretation utilizing a 2D convolutional neural network approach with a focus on coalbed horizons, dealing with the problem as an image classification task. At the beginning, linear annotations, reflecting geological fault features, are applied to multiple seismic sections. By considering texture differences between fault and non-fault areas, we construct samples that represent these distinct zones for training deep neural networks. Initially, fault annotations are transformed into single dots to facilitate pixel-based processing. To depict a specific dot's geological structure, we employ a matrix clipped around the point, determined by a combination of range and step parameters. Convolutional layers generate filters equivalent to seismic data transformation, streamlining the need for analysis and selection of seismic attributes. The article discusses enhancing the efficiency of 2D convolutional neural network–based fault interpretation by optimizing sample selection, data extraction and model construction procedures. Through the incorporation of data from two mining areas (totalling 27.09 km²) in sample creation, the overall accuracy exceeds 0.99. Recognition extends seamlessly to unlabelled sections, showcasing the innovative technical route and methodology of fault interpretation with linear annotation and pixel-based thinking. This study presents a method that integrates planar and raster thinking, transitioning from vision-oriented geological structure annotation to algorithm-oriented pixel location. The proposed 2D convolutional neural network–based matrix-oriented fault/non-fault binary classification demonstrates feasibility and reproducibility, offering a new automated approach for fault detection in coalbeds through convolutional neural network algorithms.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3350-3370"},"PeriodicalIF":1.8,"publicationDate":"2024-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222348","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":"Bayesian time-lapse full waveform inversion using Hamiltonian Monte Carlo","authors":"P. D. S. de Lima,&nbsp;M. S. Ferreira,&nbsp;G. Corso,&nbsp;J. M. de Araújo","doi":"10.1111/1365-2478.13604","DOIUrl":"10.1111/1365-2478.13604","url":null,"abstract":"<p>Time-lapse images carry out important information about dynamic changes in Earth's interior, which can be inferred using different full waveform inversion schemes. The estimation process is performed by manipulating more than one seismic dataset, associated with the baseline and monitors surveys. The time-lapse variations can be so minute and localized that quantifying the uncertainties becomes fundamental to assessing the reliability of the results. The Bayesian formulation of the full waveform inversion problem naturally provides confidence levels in the solution, but evaluating the uncertainty of time-lapse seismic inversion remains a challenge due to the ill-posedness and high dimensionality of the problem. The Hamiltonian Monte Carlo can effectively sample over high-dimensional distributions with affordable computational efforts. In this context, we explore the sequential approach in a Bayesian fashion for time-lapse full waveform inversion using the Hamiltonian Monte Carlo method. The idea relies on integrating the baseline survey information as prior knowledge to the monitor estimation. We compare this methodology with a parallel scheme in perfect and a simple perturbed acquisition geometry scenario considering the Marmousi and a typical Brazilian pre-salt velocity model. We also investigate the correlation effect between baseline and monitor samples on the propagated uncertainties. The results show that samples between different surveys are weakly correlated in the sequential case, while the parallel strategy provides time-lapse images with lower dispersion. Our findings demonstrate that both methodologies are robust in providing uncertainties even in non-repeatable scenarios.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3381-3398"},"PeriodicalIF":1.8,"publicationDate":"2024-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13604","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222291","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":"Imaging of train noise with heavy traffic events recorded by distributed acoustic sensing","authors":"Hanyu Zhang,&nbsp;Lei Xing,&nbsp;Xingpeng Zheng,&nbsp;Tuanwei Xu,&nbsp;Dimin Deng,&nbsp;Mingbo Sun,&nbsp;Huaishan Liu,&nbsp;Shiguo Wu","doi":"10.1111/1365-2478.13597","DOIUrl":"10.1111/1365-2478.13597","url":null,"abstract":"<p>Train noise is a kind of green, non-destructive and strong-energy artificial seismic sources, which is widely used in railway safety monitoring, near-surface imaging and urban underground space exploration. Distributed acoustic sensing is a new seismic acquisition technology, which has the advantages of dense sampling, simple deployment and strong anti-electromagnetic interference ability. In recent years, distributed acoustic sensing has been gradually applied in the fields of urban traffic microseism monitoring, crack detection and underground space imaging. However, previous studies mainly focused on microseism interferometry using train event coda noise, and there is limited research on the workflow of interferometry imaging using distributed acoustic sensing–based heavy train events noise (with short coda windows), which produces an abundant of near-source interference. Aiming at proving the effectiveness of this idea, we investigated a process workflow to get underground shear-velocity structure based on distributed acoustic sensing recorded heavy traffic noise near Qinhuangdao train station. A weighted sliding absolute average method is used to weaken the strong amplitude to the coda wave level and reduce the near-source influence. We demonstrated that the cross-coherence interferometry method, after spectral whitening, has the best effect on sidelobe suppression in the virtual source surface wave shot gathers, through a comparative analysis of cross-correlation and cross-coherence results. For obtaining concentrated energy and strong continuity in phase velocity spectra, we selected the time windows with high spatial coherence and signal-to-noise ratio not less than 1.2 for stacking from 720 time windows in <i>F</i>–<i>K</i> domain. When dividing subarrays to extract pseudo-two-dimensional profile, we set the overlap rate between adjacent time windows to 80% to increase stacking times, enhancing the precision of phase velocity spectra and reducing the errors of picking dispersion curve. Our results show that heavy traffic train events noise (non-pure coda) can be used to detect underground velocity structure with clear dispersion and high inversion reliability. This research provides a new processing flow for distributed acoustic sensing train noise imaging and can be applied in future urban underground space exploration.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3399-3413"},"PeriodicalIF":1.8,"publicationDate":"2024-09-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222280","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":"Envelope normalized reflection waveform inversion","authors":"Yilin Wang, Benxin Chi, Liangguo Dong","doi":"10.1111/1365-2478.13598","DOIUrl":"https://doi.org/10.1111/1365-2478.13598","url":null,"abstract":"The reflection waveform inversion has the capability to reconstruct the background velocity model using only the reflection data by employing a migration/demigration process. Utilizing the waveform discrepancy to update the background velocity model, the conventional reflection waveform inversion method heavily relies on the true‐amplitude migration/demigration technique to reproduce the primary amplitude information from the observed reflections. We can reproduce the amplitude of observed reflections by performing least‐squares reverse time migration to estimate the reflectivity in each iteration. However, this strategy is quite time‐consuming. To avoid the need for the true‐amplitude migration/demigration or least‐squares reverse time migration, we develop an amplitude‐independent reflection waveform inversion method that uses an envelope‐normalized objective function. The envelope‐normalized waveform difference can extract the phase residuals accurately as a function of time. Compared with the global energy–normalized misfit, our proposed envelope‐normalized objective function is essentially a phase‐matched measurement. At the same time, due to the amplitude independence of our proposed objective function, the subsequent weak reflections contribute with a similar weight to the total value of the misfit as the strong early reflections do. This makes it possible to recover the deep subsurface velocity. Synthetic data of the Sigsbee model and marine streamer field data applications validate that our amplitude‐independent reflection waveform inversion method can further improve the resolution and accuracy by aligning the reflection events of synthetic and observed data phase to phase without the need to perform true‐amplitude migration/demigration or least‐squares reverse time migration as in conventional reflection waveform inversion.","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"33 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222315","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":"Improved elastic full‐waveform inversion of ocean bottom node data","authors":"Bo Wu, Gang Yao, Qingqing Zheng, Fenglin Niu, Di Wu","doi":"10.1111/1365-2478.13601","DOIUrl":"https://doi.org/10.1111/1365-2478.13601","url":null,"abstract":"Elastic full‐waveform inversion enables the quantitative inversion of multiple subsurface parameters, significantly enhancing the interpretation of subsurface lithology. Simultaneously, with the ongoing advancements in ocean bottom node technology, the application of elastic full‐waveform inversion to marine ocean bottom node data is receiving increasing attention. This is attributed to the capability of ocean bottom node to acquire high‐quality four‐component data. However, elastic full‐waveform inversion of ocean bottom node data typically encounters two challenges: First, the presence of low S‐wave velocity layers in the seabed leads to weak energy of converted S‐waves, resulting in significantly poorer inversion results for S‐wave velocity compared to those for P‐wave velocity; second, the cross‐talk effect of multiple parameters further exacerbates the difficulty in inverting S‐wave velocity. To effectively recover the S‐wave velocity using ocean bottom node data, we modify the S‐wave velocity gradient in conventional elastic full‐waveform inversion to alleviate the impact of cross‐talk from multiple parameters on the inversion of S‐wave velocity. Furthermore, to invert for density parameters, we adopt a two‐stage inversion strategy. In the first stage, P‐wave and S‐wave velocities are updated simultaneously with a single‐step length. Because the initial density model is far from the true one, density is updated using an empirical relationship derived from well‐log data. In the second stage, velocities and density are updated simultaneously with multi‐step length to further refine the models obtained in the first stage. The high effectiveness of the improved elastic full‐waveform inversion is validated by numerical examples.","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"30 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142227571","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":"Cross‐correlation reflection waveform inversion based on a weighted norm of the time‐shift obtained by dynamic image warping","authors":"Yingming Qu, Shihao Dong, Tianmiao Zhong, Yi Ren, Zizheng Li, Boshen Xing, Yifan Li","doi":"10.1111/1365-2478.13599","DOIUrl":"https://doi.org/10.1111/1365-2478.13599","url":null,"abstract":"The computational efficiency of cross‐correlation reflection waveform inversion can be improved by utilizing the outcomes of reverse time migration instead of the least‐squares reverse time migration results in each iteration. However, the inversion effect of cross‐correlation reflection waveform inversion needs to be optimized as the inversion results may not be optimal. The conventional cross‐correlation operator tends to produce interference values that can compromise the precision of time‐shift estimations. Moreover, the time shift obtained through dynamic image warping can exhibit spiky disturbances, making it difficult to determine accurate time‐shift values. These challenges can cause the inversion process to converge to a local minimum, thereby affecting the quality of the inversion results. To address these limitations, this paper proposes a new approach called cross‐correlation reflection waveform inversion based on dynamic image warping. The proposed method integrates a weighted norm derived from dynamic image warping to effectively regulate the time‐shift values throughout the inversion process. The effectiveness of the proposed cross‐correlation reflection waveform inversion based on the dynamic image warping method is validated through simulations using a simple two‐layer model and a resampled Sigsbee 2A model. A comparative analysis is performed to evaluate the performance of cross‐correlation reflection waveform inversion based on dynamic image warping in mitigating cross‐correlation interference, demonstrating its superior capability compared to the conventional cross‐correlation reflection waveform inversion method.","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"291 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222281","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":"Thomsen, L., 2023, A logical error in Gassmann poroelasticity: Geophysical Prospecting, 71, 649–663. by Leon Thomsen, University of Houston","authors":"","doi":"10.1111/1365-2478.13567","DOIUrl":"https://doi.org/10.1111/1365-2478.13567","url":null,"abstract":"<p>Two figure captions in this paper were in error, confusing compressibility and incompressibility (the figures themselves were correct). The proper figure captions are</p><p>FIGURE 2. Comparison of Berea sandstone data from Hart and Wang (2010) for <i>K</i>_ud − <i>K</i>_fm (as functions of differential pressure, <i>p_d</i> = <i>p</i> − <i>p_F</i>) with predictions from Gassmann theory (Equation 1, using data for <i>K_𝑆</i> (from Equation 14; see also the unnumbered equation from B&amp;K following Equation 17), or from VRH theory), and from B&amp;K theory (Equation 19, using data for <i><span>K</span>_𝑆</i> and for <i>κ_M</i> (from Equation 21)). The Fluid (water) incompressibility <i>K_F</i> is taken as 2.3 GPa.</p><p>FIGURE 4. Comparison of Indiana limestone data from Hart and Wang (2010) for <i>K</i>_ud − <i>K</i>_fm (as functions of differential pressure, <i>p_d</i> = <i>p</i> − <i>p_F</i>) with predictions from Gassmann theory (Equation 1, using data for <i>K_S</i> (from Equation 14; see also the unnumbered equation from B&amp;K following Equation 17), or from VRH theory), and from B&amp;K theory (Equation 19, using data for <i><span>K</span>_𝑆</i> and <i>κ_M</i> (from Equation 21)). The Fluid (water) incompressibility <i>K_F</i> is taken as 2.3 GPa.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2857"},"PeriodicalIF":1.8,"publicationDate":"2024-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13567","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142045334","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":"Centralized feature pyramid-based supervised deep learning for object detection model from GPR data","authors":"Kun Yan,&nbsp;Xianlei Xu,&nbsp;Pengqiao Zhu,&nbsp;Zhaoyang Zhang","doi":"10.1111/1365-2478.13590","DOIUrl":"10.1111/1365-2478.13590","url":null,"abstract":"<p>To address low detection accuracy and speed due to the multisolvability of the ground-penetrating radar signal, we proposed a novel centralized feature pyramid-YOLOv6l–based model to enhance detection precision and speed in road damage and pipeline detection. The centralized feature pyramid was used to obtain rich intra-layer features and improve the network performance. Our proposed model achieves higher accuracy compared with the existing detection models. We also built two new evaluating indexes, relative average precision and relative mean average precision, to fully evaluate the detection accuracy. To verify the applicability of our model, we conducted a road field detection experiment on a ground-penetrating radar dataset we collected and found that the proposed model had good performance in increasing detection precision, achieving the highest mean average precision compared with YOLOv7, YOLOv5 and YOLOx models, with relative mean average precision and frame rate per second at 16.38% and 30.5%, respectively. The detection information for the road damage and pipeline were used to conduct three-dimensional imaging. Our model is suitable for object detection in ground-penetrating radar images, thereby providing technical support for road damage and underground pipeline detection.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3414-3435"},"PeriodicalIF":1.8,"publicationDate":"2024-08-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222299","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":"Blind spectral inversion of seismic data","authors":"Yaoguang Sun,&nbsp;Siyuan Cao,&nbsp;Siyuan Chen,&nbsp;Yuxin Su","doi":"10.1111/1365-2478.13594","DOIUrl":"10.1111/1365-2478.13594","url":null,"abstract":"<p>Reflectivity inversion is a key step in reservoir prediction. Conventional sparse-spike deconvolution assumes that the reflectivity (reflection coefficient series) is sparse and solves for the reflection coefficients by an <i>L</i>1-norm inversion process. Spectral inversion is an alternative to sparse-spike deconvolution, which is based on the odd–even decomposition algorithm and can accurately identify thin layers and reduce the wavelet tuning effect without using constraints from logging data, from horizon interpretations or from an initial model of the reflectivity. In seismic processing, an error exists in wavelet extraction because of complex geological structures, resulting in the low accuracy of deconvolution and inversion. Blind deconvolution is an effective method for solving the problem mentioned above, which comprises seismic wavelet and reflectivity sequence, assuming that the wavelets that affect some subsets of the seismic data are approximately the same. Therefore, we combined blind deconvolution with spectral inversion to propose blind spectral inversion. Given an initial wavelet, we can calculate the reflectivity based on spectral inversion and update the wavelet for the next iteration. During the update processing, we add the smoothness of the wavelet amplitude spectrum as a regularization term, thus reducing the wavelet oscillation in the time domain, increasing the similarity between inverted and initial wavelets, and improving the stability of the solution. The blind spectral inversion method inherits the wavelet robustness of blind deconvolution and high resolution of spectral inversion, which is suitable for reflectivity inversion. Applications to synthetic and field seismic datasets demonstrate that the blind spectral inversion method can accurately calculate the reflectivity even when there is an error in wavelet extraction.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 9","pages":"3436-3447"},"PeriodicalIF":1.8,"publicationDate":"2024-08-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142222289","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}