Geophysical Prospecting最新文献

{"title":"Study on the accurate detection method of full-polarimetric ground-penetrating radar faults in mines based on modified Yamaguchi decomposition","authors":"Ran Wang,&nbsp;Fan Cui,&nbsp;Guoqi Dong,&nbsp;Qi Cheng,&nbsp;Guixin Zhang,&nbsp;Xuhao Wang,&nbsp;Xiaoxue Zhao","doi":"10.1111/1365-2478.13636","DOIUrl":"https://doi.org/10.1111/1365-2478.13636","url":null,"abstract":"<p>Obtaining information on tectonic tendencies is a prerequisite for intelligent and accurate mining in mines. In the special mine environment, the co-polarized ground-penetrating radar can only identify the spatial location of faults, and it is difficult to analyse the inclination information of fault structures. This paper proposes a mine full-polarimetric ground-penetrating radar fault tendency detection method based on this. First, based on the stacking characteristics of the coal depositional, this paper analyses the propagation law of the pulse electromagnetic wave in the coal seam and puts forward the assumption of the overlapping echo reflection of the fault structure. The reasonableness of the fault reflection assumption is verified through a numerical simulation study. Second, based on the cutting relationship of the fault to the coal seam, we divided the reflection structure of the fault structure into plane scattering and dihedral angle scattering. We realized the mingled echoes’ decomposition using the improved Yamaguchi decomposition technique. To analyse the applicability of the modified Yamaguchi and Freeman decomposition methods in the identification of fault inclination, we use the upright fault simulation data for the discussion, and we find that the improved Yamaguchi decomposition method is more advantageous in the identification of fault inclination in the mine. The decomposition results based on the simulation data of fault models with different dip angles found that when the dip angle of the fault is less than 90°, the scattering of the fault structure is dominated by planar scattering and dihedral angle scattering; when the dip angle of the fault is greater than or equal to 90°, the scattering of the fault structure is dominated by planar scattering, and the scattering power of the dihedral angle model is zero. By analysing the effect of fault strike on the decomposition results, it is found that the fault strike angle has little effect on the identification of fault tendency. Finally, the application potential of this paper's method is tested by constructing complex numerical models and probing experiments. Therefore, the method proposed in this paper can solve the fault tendency identification under a thick coal seam.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"233-250"},"PeriodicalIF":1.8,"publicationDate":"2024-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143111944","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":"A Ramanujan subspace and dynamic time warping and adaptive singular value decomposition combined denoising method for low signal-to-noise ratio surface microseismic monitoring data in hydraulic fracturing","authors":"Xu-Lin Wang,&nbsp;Jian-Zhong Zhang,&nbsp;Zhong-Lai Huang","doi":"10.1111/1365-2478.13634","DOIUrl":"https://doi.org/10.1111/1365-2478.13634","url":null,"abstract":"<p>Surface microseismic monitoring is widely used in hydraulic fracturing. Real-time monitoring data collected during fracturing can be used to perform surface-microseismic localization, which aids in assessing the effects of fracturing and provides guidance for the process. The accuracy of localization critically depends on the quality of monitoring data. However, the signal-to-noise ratio of the data is often low due to strong coherent and random noise, making denoising essential for processing surface monitoring data. To suppress noise more effectively, this paper introduces a novel denoising method that integrates the Ramanujan subspace with dynamic time warping and adaptive singular value decomposition. The new method consists of two steps: First, a Ramanujan subspace is constructed to suppress periodic noise. Then, dynamic time warping and adaptive singular value decomposition are applied to eliminate remaining coherent and random noise. The method has been evaluated using both synthetic and field data, and its performance is compared with traditional microseismic denoising techniques, including bandpass filtering and empirical mode decomposition.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"251-265"},"PeriodicalIF":1.8,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119310","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":"3D Controlled-source electromagnetic modelling in anisotropic media using secondary potentials and a cascadic multigrid solver","authors":"Kejia Pan,&nbsp;Jinxuan Wang,&nbsp;Xu Han,&nbsp;Zhengyong Ren,&nbsp;Weiwei Ling,&nbsp;Rongwen Guo","doi":"10.1111/1365-2478.13624","DOIUrl":"https://doi.org/10.1111/1365-2478.13624","url":null,"abstract":"<p>Quantitative interpretation of the data from controlled-source electromagnetic methods, whether via forward modelling or inversion, requires solving a considerable number of forward problems, and multigrid methods are often employed to accelerate the solving process. In this study, a new extrapolation cascadic multigrid method is employed to solve the large sparse complex linear system arising from the finite element approximation of Maxwell's equations using secondary potentials. The equations using secondary potentials are discretized by the classic nodal finite element method on nonuniform rectilinear grids. The resulting linear systems are solved by the extrapolation cascadic multigrid method with a new prolongation operator and preconditioned Stabilized bi-conjugate gradient method smoother. High-order interpolation and global extrapolation formulas are utilized to construct the multigrid prolongation operator. The extrapolation cascadic multigrid method with the new prolongation operator is easier to implement and more flexible in application than the original one. Finally, several synthetic examples including layered models, models with anisotropic anomalous bodies or layers, are used to validate the accuracy and efficiency of the proposed method. Numerical results show that the extrapolation cascadic multigrid method improves the efficiency of 3D controlled-source electromagnetic forward modelling a lot, compared with traditional iterative solvers and some state-of-the-art methods or software (e.g., preconditioned flexible generalized minimal residual method, emg3d) in the considered models and grid settings. The efficiency benefit is more evident as the number of unknowns increases, and the proposed method is more efficient at low frequencies. The extrapolation cascadic multigrid method can also be used to solve systems of equations arising from related applications, such as induction logging, airborne electromagnetic, etc.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"284-302"},"PeriodicalIF":1.8,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119312","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 methods of visco-acoustic migration based on the De Wolf approximation and comparison of their migration images","authors":"Huachao Sun,&nbsp;Jianguo Sun,&nbsp;Zhenghui Gao","doi":"10.1111/1365-2478.13637","DOIUrl":"https://doi.org/10.1111/1365-2478.13637","url":null,"abstract":"<p>The viscosity of a medium affects the amplitude attenuation and velocity dispersion of seismic waves. Therefore, it is necessary to consider these factors during migration. First, to eliminate the viscous effect of a medium, we combine the Futterman model with the integral equation of the De Wolf approximation to construct a compensation operator of the De Wolf approximation for a visco-acoustic medium. Next, we use the visco-acoustic screen approximation method to realize the continuation operator then establish a prestack depth migration algorithm. Finally, an error analysis, impulse response test and model test are performed. The results show that three different generalized visco-acoustic screen methods (phase screen method, generalized screen method and extended local Born Fourier method) can satisfactorily compensate for the attenuation of deep interface amplitude. Among these methods, the visco-acoustic extended local Born Fourier method has the highest accuracy and the best compensation effect.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"266-283"},"PeriodicalIF":1.8,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119311","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":"Full-tensor magnetic gradiometry: Comparison with scalar total magnetic intensity, processing and visualization guidelines","authors":"Hernan Ugalde,&nbsp;Bill Morris,&nbsp;Akshay Kamath,&nbsp;Brian Parsons","doi":"10.1111/1365-2478.13629","DOIUrl":"https://doi.org/10.1111/1365-2478.13629","url":null,"abstract":"<p>Full-tensor magnetic gradiometry data have been collected commercially for the last few years. However, to date, there is still no clarity on how to compare these data to scalar total field surveys. Some users display the vertical gradient of the vertical component (<i>B</i>_zz) and compare that to a first vertical derivative of total field with the caveat that ‘they are similar’. Others compute the length of the measured vector and call that total field. We establish the basic formulas to calculate total field from the tensor components and demonstrate this with a real data example from Thompson, Manitoba, Canada. Another key question is whether full-tensor interpolation is required to obtain total field from tensor data. We compare the results from using a commercial full-tensor interpolation algorithm with standard minimum curvature of the tensor components individually and with another open-source code that uses a radial basis function interpolator on the individual tensor components. All three applications produced a total field grid of superior quality to that calculated from a scalar total field survey available for the area of study.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"303-314"},"PeriodicalIF":1.8,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143118404","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":"Inversion of downhole resistivity properties through infrared spectroscopy and whole-rock geochemistry using machine-learning","authors":"Mehdi Serdoun,&nbsp;Frédéric Sur,&nbsp;Gaétan Milesi,&nbsp;Elodie Williard,&nbsp;Pierre Martz,&nbsp;Julien Mercadier","doi":"10.1111/1365-2478.13627","DOIUrl":"https://doi.org/10.1111/1365-2478.13627","url":null,"abstract":"<p>The electrical properties of rocks are widely used in the geophysical exploration of natural resources, such as minerals, hydrocarbons and groundwater. In mining exploration, the primary goal is to map electrically anomalous geological features associated with different mineralization styles, such as clay alteration haloes, metal oxides and sulphides, weathered crystalline rocks or fractured zones. As such, the reconciliation of geophysical data with geological information (geochemistry, mineralogy, texture and lithology) is a critical step and can be performed based on petrophysical properties collected either on core samples or as downhole measurements. Based on data from 189 diamond drill cores collected for uranium exploration in the Athabasca Basin (Saskatchewan, Canada), this paper presents a case study of reconciliation of downhole resistivity probing with core sample geochemistry and short-wave infrared spectroscopy (350–2500 nm) through three successive steps: (i) multivariate analysis of resistivity and other petrophysical properties (porosity, density) against geochemical and infrared spectroscopy information to characterize electrical properties of rocks with respect to other physical parameters, (ii) a machine-learning workflow integrating geochemistry and spectral signatures in order to infer synthetic resistivity logs along with uncertainties. The best model in the basin was Light Gradient-Boosting Machine with pairwise log-ratio, which yielded a coefficient of determination <i>R</i>² = 0.80 (root mean square error = 0.16), and in the basement, support vector regression with data fusion of infrared spectroscopy and pairwise log-ratios on geochemistry yielded <i>R</i>² = 0.82 (root mean square error = 0.35); (iii) the best model was then fitted on an area that was excluded from the original dataset (Getty Russell property) in order to infer synthetic resistivity logs for that zone. Software code is publicly available. This workflow can be re-used for the valorization of legacy datasets.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"355-379"},"PeriodicalIF":1.8,"publicationDate":"2024-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143116436","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":"Underwater unexploded ordnance discrimination based on intrinsic target polarizabilities – A case study","authors":"Erika Gasperikova,&nbsp;Ugo Conti,&nbsp;H. Frank Morrison","doi":"10.1111/1365-2478.13631","DOIUrl":"https://doi.org/10.1111/1365-2478.13631","url":null,"abstract":"<p>Seabed unexploded ordnance that resulted partly from the high failure rate among munitions from more than 80 years ago and from decades of military training and testing of weapons systems poses an increasing concern all around the world. Although existing magnetic systems can detect clusters of debris, they are not able to tell whether a munition is still intact requiring special removal (e.g. in situ detonation) or is harmless scrap metal. The marine environment poses unique challenges, and transferring knowledge and approaches from land to a marine environment has not been easy and straightforward. On land, the background soil conductivity is much lower than the conductivity of the unexploded ordnance and the electromagnetic response of a target is essentially the same as that in free space. For those frequencies required for target characterization in the marine environment, the seawater response must be accounted for and removed from the measurements. The system developed for this study uses fields from three orthogonal transmitters to illuminate the target and four three-component receivers to measure the signal arranged in a configuration that inherently cancels the system's response due to the enclosing seawater, the sea–bottom interface and the air–sea interface for shallow deployments. The system was tested as a cued system on land and underwater in San Francisco Bay – it was mounted on a simple platform on top of a support structure that extended 1 m below and allowed the diver to place metal objects to a specific location even in low-visibility conditions. The measurements were stable and repeatable. Furthermore, target responses estimated from marine measurements matched those from land acquisition, confirming that the seawater and air–sea interface responses were removed successfully. Thirty-six channels of normalized induction responses were used for the classification, which was done by estimating the target principal dipole polarizabilities. Our results demonstrated that the system can resolve the intrinsic polarizabilities of the target, with clear distinctions between those of symmetric intact unexploded ordnance and irregular scrap metal. The prototype system was able to classify an object based on its size, shape and metal content and correctly estimate its location and orientation.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"315-329"},"PeriodicalIF":1.8,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13631","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143116011","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":"Effective elastic properties of fractured rocks considering fracture interactions","authors":"Junxin Guo,&nbsp;Bo-Ye Fu","doi":"10.1111/1365-2478.13632","DOIUrl":"https://doi.org/10.1111/1365-2478.13632","url":null,"abstract":"<p>Fracture interactions are an important factor that affects rock effective elastic properties. We study the fracture interaction effects by combinations of theoretical modelling, numerical simulations and experiments in this work. First, we propose a simplified differential effective medium scheme and self-consistent approximation for effective elastic properties of rocks with aligned fractures, and we compare their results to those of non-interaction approximation. The results show that the predictions by differential effective medium scheme and self-consistent approximation are lower than those by non-interaction approximation. This indicates that the differential effective medium scheme and self-consistent approximation quantify the stress amplification effects but not stress shielding effects. To validate this, we carry out numerical simulations for cases with coplanar cracks and stacked cracks, respectively. We find that the stress shielding effect (stacked cracks) causes a significant increment of effective elastic stiffnesses of fractured rocks. However, the stress amplification (coplanar cracks) has the opposite effect, which induces a slight reduction in rock effective elastic stiffnesses. The differential effective medium scheme quantifies the lower bound for this effect well. Besides numerical simulations, applying theoretical models in experimental measurements also shows a pronounced effect of stress shielding but a small influence of stress amplification on fractured rock elastic properties. This work indicates that the stress shielding is an important effect that affects fractured rock elastic properties. Without considering this effect, the fracture density may be largely underestimated by the seismic or sonic logging inversion. Hence, the models accounting for stress shielding effects need to be developed in the future, which can combine with the above models for the seismic or sonic logging inversion of fracture properties.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"330-344"},"PeriodicalIF":1.8,"publicationDate":"2024-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143116010","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":"Seismic fault detection with sliding windowed differential cepstrum–based coherence analysis","authors":"Qi Ran,&nbsp;Kang Chen,&nbsp;Cong Tang,&nbsp;Long Wen,&nbsp;Ming Zeng,&nbsp;Han Liang,&nbsp;Guang-rong Zhang,&nbsp;Han Xiao,&nbsp;Ya-juan Xue","doi":"10.1111/1365-2478.13633","DOIUrl":"https://doi.org/10.1111/1365-2478.13633","url":null,"abstract":"<p>Cepstral decomposition is beneficial for highlighting certain geological features within the particular quefrency bands which may be deeply buried within the wide quefrency range of the seismic data. Converting seismic traces into the corresponding cepstrum components can better analyse some characteristics of underground strata than the traditional spectral decomposition methods. We propose the sliding windowed differential cepstrum–based coherence analysis approach to delineate the fault features. First, the data are decomposed using a sliding windowed differential cepstrum, which results in multi-cepstrum data of corresponding quefrency of certain bandwidth. These different multi-cepstrum data may highlight the different stratigraphic features in a certain quefrency band. We select the first-order common quefrency volume as the featured attribute. Then, eigenstructure-based coherence is applied on the first-order common quefrency data volume to statistically obtain the fault detection result with a finer and sharper image. Synthetic data and field data examples show that the proposed method has the ability to better visualize all the possible subtle and minor faults present in the data more accurately and discernibly than the traditional coherence method. Compared with the ant-tracking method, the proposed method is more effective in revealing the major faults. It is hoped that this work will complement current fault detection methods with the addition of the cepstral-based method.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 1","pages":"345-354"},"PeriodicalIF":1.8,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115313","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":"Addressing cycle-skipping in full-waveform inversion using acoustic wave energy","authors":"Zhonglei Li,&nbsp;Gang Yao","doi":"10.1111/1365-2478.13630","DOIUrl":"https://doi.org/10.1111/1365-2478.13630","url":null,"abstract":"<p>Limitations in acquisition technologies lead to insufficient low-frequency signals in field seismic data. Local optimization methods are the common approaches for full-waveform inversion. Inaccurate initial velocity models and lack of low-frequency signals in seismic data typically cause the local-gradient-based full-waveform inversion to converge to a local minimum due to cycle-skipping. The existing energy-based objective functions can generate artificial low-frequency signals successfully by squaring the pressure but overlook the law of energy conservation, which may mislead model updates. To overcome this issue, we combine acoustic wave potential energy and kinetic energy to develop a new objective function that fits the acoustic wave energy. The new acoustic-wave-energy-based full-waveform inversion considers the law of energy conservation. The new system creates low-frequency signals to avoid cycle-skipping and produce an accurate smooth background velocity model, which provides a sufficient starting model for conventional full-waveform inversion. Numerical examples demonstrate that the combination of acoustic-wave-energy-based full-waveform inversion and conventional full-waveform inversion can deliver more faithful and accurate final results than conventional full-waveform inversion alone.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"73 3","pages":"960-976"},"PeriodicalIF":1.8,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143513455","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}