Geophysics最新文献

{"title":"Sparse seismic data regularization in both shot and trace domains using a residual block autoencoder based on the fast Fourier transform","authors":"Alexandre L. Campi, R. Misságia","doi":"10.1190/geo2023-0097.1","DOIUrl":"https://doi.org/10.1190/geo2023-0097.1","url":null,"abstract":"The increasing use of sparse acquisitions in seismic data acquisition offers advantages in cost and time savings. However, it results in irregularly sampled seismic data, adversely impacting the quality of the final images. In this paper, we propose the ResFFT-CAE network, a convolutional neural network with residual blocks based on the Fourier transform. Incorporating residual blocks allows the network to extract both high- and low-frequency features from the seismic data. The high-frequency features capture detailed information, while the low-frequency features integrate the overall data structure, facilitating superior recovery of irregularly sampled seismic data in the trace and shot domains. We evaluated the performance of the ResFFT-CAE network on both synthetic and field data. On synthetic data, we compared the ResFFT-CAE network with the compressive sensing (CS) method utilizing the curvelet transform. For field data, we conducted comparisons with other neural networks, including the convolutional autoencoder (CAE) and U-Net. The results demonstrated that the ResFFT-CAE network consistently outperformed other approaches in all scenarios. It produced images of superior quality, characterized by lower residuals and reduced distortions. Furthermore, when evaluating model generalization, tests using models trained on synthetic data also exhibited promising results. In conclusion, the ResFFT-CAE network shows great promise as a highly efficient tool for the regularizing irregularly sampled seismic data. Its excellent performance suggests potential applications in the preconditioning of seismic data analysis and processing flows.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"11 11","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138585851","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Full-waveform inversion of pure quasi-P waves in titled transversely isotropic media based on different parameterization","authors":"Zhiming Ren, Xue Dai, Lei Wang","doi":"10.1190/geo2023-0367.1","DOIUrl":"https://doi.org/10.1190/geo2023-0367.1","url":null,"abstract":"Full-waveform inversion (FWI) builds subsurface parameter models by minimizing the residuals between the modeled and observed data. Accounting for the effects of anisotropy is critical for high-resolution imaging of complex structures. We develop an acoustic anisotropic FWI method based on a pure quasi P-wave (qP-wave) equation. The equation coefficients and their derivatives with respect to the Thomsen’s anisotropy parameters ([Formula: see text] and [Formula: see text]) are estimated by least-squares (LS) optimization. We derive the functional gradients and analyze the radiation patterns for six parameter classes: the velocity along the symmetry axis [Formula: see text], [Formula: see text] and [Formula: see text], the normal moveout velocity [Formula: see text], anisotropy parameters [Formula: see text] and [Formula: see text], the horizontal velocity [Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text], and [Formula: see text], [Formula: see text] and [Formula: see text]. The parameterization [Formula: see text] has significant tradeoff between [Formula: see text] and [Formula: see text] at the intermediate and wide scattering angles. The anisotropy parameters [Formula: see text] and [Formula: see text] are resolvable at the short scattering angles for the parameterizations [Formula: see text] and [Formula: see text], respectively. The parameter crosstalk for the parameterizations [Formula: see text] and [Formula: see text] is more serious than that for other types of parameterizations. We perform FWI of pure qP-waves in vertical and titled transversely isotropic (VTI and TTI) media. Inversion results on the Overthrust VTI model and the modified BP TTI model show that the velocity, anisotropy parameters and tilt angle can be individually reconstructed when other parameters are sufficiently accurate. The multi-parameter FWI cannot obtain reliable tilt angles for each type of parameterization. The inversion with the parameterization [Formula: see text] produces [Formula: see text], [Formula: see text] and [Formula: see text] models with modest accuracy, whereas the parameterization [Formula: see text] helps to improve the accuracy of [Formula: see text] and [Formula: see text] models.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"357 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138983211","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Ordering cross-spread gathers","authors":"Stewart Trickett","doi":"10.1190/geo2023-0161.1","DOIUrl":"https://doi.org/10.1190/geo2023-0161.1","url":null,"abstract":"Seismic processing on cross-spread gathers is a valuable but underexploited strategy. To do it properly, sources from a single source line and receivers from a single receiver line must be ordered in a physically sensible way, so that adjacent sources or receivers on a surface diagram are physically near each other. But determining such an ordering is a challenge on irregularly acquired land data. I propose novel automatic ordering algorithms using tools from graph theory that minimize large gaps and preserve the sequential patterns found even in highly irregular acquisition. Java source code is available.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"547 ","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139011090","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A frequency-Hankel transform method to extract multimodal Rayleigh wave dispersion spectra from active and passive source surface wave data","authors":"Zhentao Yang, Yao-Chong Sun, Dazhou Zhang, Peng Han, Xiaofei Chen","doi":"10.1190/geo2023-0189.1","DOIUrl":"https://doi.org/10.1190/geo2023-0189.1","url":null,"abstract":"Rayleigh wave dispersion energy spectra have been widely used to extract dispersion curves and invert for underground shear-wave velocity structures for engineering geophysics and seismology. We propose a frequency-Hankel (F-H) transform method to extract high-quality multimodal Rayleigh wave dispersion energy spectra from active and passive source Rayleigh wave data. The F-H transform method is inspired by the frequency-Bessel (F-J) transform method and considers the physical meaning of Green’s functions for Rayleigh wave dispersion analysis. The F-H transform method can naturally avoid crossed artefacts caused by converging waves on F-J spectrograms and obtains more multimodal dispersion spectra of the same quality with fewer Rayleigh wave data than the F-J transform method. Both synthetic and field Rayleigh wave data from active and passive sources for near-surface exploration and ambient noise tomography are used to demonstrate the validity, accuracy and applicability of the F-H transform method. The F-H transform method unifies the F-J transform method and its modifications for active and passive sources Rayleigh wave data. The F-H transform method is a robust and efficient multimodal Rayleigh wave dispersion analysis method for active and passive source Rayleigh wave data.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"213 ","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139011080","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Image-domain seismic inversion by deblurring with invertible Recurrent Inference Machines","authors":"Haorui Peng, Ivan Vasconcelos, M. Ravasi","doi":"10.1190/geo2022-0780.1","DOIUrl":"https://doi.org/10.1190/geo2022-0780.1","url":null,"abstract":"In complex geological settings and in the presence of sparse acquisition systems, seismic migration images manifest as non-stationary blurred versions of the unknown subsurface model. Thus, image-domain deblurring is an important step to produce interpretable and high-resolution models of the subsurface. Most deblurring methods focus on inverting seismic images for their underlying reflectivity by iterative least-squares inversion of a local Hessian approximation; this is obtained by either direct modeling of the so-called point spread functions or by a migration-demigration process. In this work, we adopt a novel deep learning framework, based on invertible Recurrent Inference Machines (i-RIMs), which allows approaching any inverse problem as a supervised learning task informed by the known modeling operator (convolution with point-spread functions in our case): the proposed algorithm can directly invert migrated images for impedance perturbation models, assisted with the prior information of a smooth velocity model and the modeling operator. Because i-RIMs are constrained by the forward operator, they implicitly learn to shape/regularise output models in a training-data-driven fashion. As such, the resulting deblurred images show great robustness to noise in the data and spectral deficiencies (e.g., due to limited acquisition). The key role played by the i-RIM network design and the inclusion of the forward operator in the training process is supported by several synthetic examples. Finally, using field data, we show that i-RIM-based deblurring has great potential in yielding robust, high-quality relative impedance estimates from migrated seismic images. Our approach could be of importance towards future Deep-Learning-based quantitative reservoir characterization and monitoring.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"17 23","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138589283","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Physical property characterization of rocks in the Bayan Obo REE-Nb-Fe deposit, China","authors":"Lili Zhang, Hongrui Fan, Jian Wang, Liang Zhao, Kuifeng Yang, Ya Xu, Yonggang Zhao, Xingwang Xu, Meizhen Hao, Zhanfeng Yang, Xianhua Li","doi":"10.1190/geo2023-0439.1","DOIUrl":"https://doi.org/10.1190/geo2023-0439.1","url":null,"abstract":"Bayan Obo ore deposit is the world’s largest rare-earth element (REE) resource, the second largest niobium (Nb) resource, and also a significant iron (Fe) resource in China. Evaluating resource potential for the deposit has become a focus of global interest. Rock physical properties bridge geophysical exploration and geological modeling; variation in these parameters is necessary for successful geophysical application. REE, Nb, iron, and potassium are mainly hosted in dolomite and slate in the Bayan Obo Group, and REE mineralization is genetically associated with carbonatite. Three physical properties (resistivity, polarizability, and magnetic susceptibility (MS)) of iron ore, slate, dolomite, and carbonatite dike outcrop samples at Bayan Obo are measured and statistically analyzed using three-dimensional reconstruction, one-/two-/three-dimensional kernel density estimation, scatterplot matrix, three-dimensional histogram, and Pearson- and maximum-information-coefficient-based correlation analysis. It is evident that iron ore, iron-mineralized fine-grained dolomite, and iron-mineralized slate are mainly of low resistivity, and iron ore and iron-mineralized fine-grained dolomite have high MS. MS favorably distinguishes iron ore from slate; MS and resistivity distinguish between iron-mineralized fine-grained dolomite and carbonatite dikes. The physical properties and whole rock geochemistry (major and trace elements) jointly demonstrate that MS of iron ore, slate, and dolomite is positively correlated with TFe2O3 content, polarizability is correlated with TFe2O3, SiO2 content is correlated with K2O, and resistivity is correlated with MS and polarizability respectively. Resistivity of iron ore and dolomite is negatively correlated with TFe2O3 content. Resistivity of iron ore is negatively correlated with TFe2O3, total rare-earth element (REE), and Nb, respectively, and correlated with thorium. The methods used have intuitive visual expression and reflect the characteristics of the physical properties and their correlation with the mineralogical composition. The results will be beneficial to determining the geometry of ore-hosting rock masses and providing crucial evidence for the resource evaluation.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"34 33","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138984006","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Estimation of the subsurface EM velocity distribution from diffraction hyperbolas by means of a novel automated picking procedure: Theory and application to glaciological GPR data sets","authors":"M. Dossi, E. Forte, B. Cosciotti, S. Lauro, E. Mattei, E. Pettinelli, M. Pipan","doi":"10.1190/geo2023-0042.1","DOIUrl":"https://doi.org/10.1190/geo2023-0042.1","url":null,"abstract":"We developed an auto-picking algorithm that is designed to automatically detect subsurface diffractors within GPR data sets; to accurately track the hyperbolic diffractions originating from the identified scatterers; and to recover the subsurface EM velocity distribution, among other possible analyses. The proposed procedure presents several advantages with respect to other commonly applied diffraction tracking techniques, since it can be applied with minimal signal pre-processing, thus making it more versatile and adaptable to local conditions; it requires only limited input from the interpreter, in the form of a few thresholds for the tracking parameters, thus making the results more objective; and it does not involve pre-training, as opposed to machine learning algorithms, thus removing the need to gather a large and comprehensive image database of all possible subsurface situations, which would not be necessarily limited to just examples of diffractions. The presented algorithm starts by identifying those signals that are likely to belong to diffraction apexes, which are then used as initial seeds by the auto-tracking process. The horizontal search window used during the auto-tracking process is locally adapted through a rough preliminary estimate of the size of each diffraction. In addition, multiple seeds within the same apex can produce several acceptable hyperbolas tracking the same diffraction phase. The algorithm thus selects the best-fitting ones by assessing several signal attributes, while also removing both redundant hyperbolas and the expected false positives. The algorithm was applied to two glaciological GPR profiles, and it was able to accurately track the vast majority of the recorded diffractions, with very few false positives and negatives. This produced a statistically sound EM velocity distribution, which was used to assess the state of the surveyed alpine glacier.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"27 22","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138983853","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Seismic Adaptive Multiple Subtraction Using a Structure-oriented Matched Filter","authors":"Yuhan Sui, Yue Ma, Lu Liu, Dongliang Zhang, Yubing Li","doi":"10.1190/geo2023-0025.1","DOIUrl":"https://doi.org/10.1190/geo2023-0025.1","url":null,"abstract":"Multiple removal is a crucial step in seismic data processing prior to velocity model building and imaging. After the prediction, adaptive multiple subtraction is employed to suppress multiples (considered noise) in seismic data, thereby highlighting primaries (considered signal). In practice, conventional adaptive subtraction methods fit the predicted and recorded multiples in the least-squares sense using a sliding window, formulating a localized adaptive matched filter. Subsequently, the filter is applied to the prediction to remove multiples from the recorded data. However, such a strategy runs the risk of over attenuating the useful primaries under the minimization energy constraint. To avoid damage to valuable signals, we propose a novel approach that replaces the conventional matched filter with a structure-oriented version. From the predicted multiples, we extract the structural information to be used in the derivation of the adaptive matched filter. The proposed structure-oriented matched filter emphasizes the structures of predicted multiples which helps to better preserve primaries during the subtraction. Synthetic and field data examples demonstrate the efficacy of the proposed structure-oriented adaptive subtraction approach, highlighting its superior performance in multiple removal and primary preservation compared to conventional methods on 2D regularly sampled data.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"65 19","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138594624","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Reverse time migration angle gathers in acoustic anisotropic media using direction vectors","authors":"Kai Yang, Jianfeng Zhang","doi":"10.1190/geo2023-0328.1","DOIUrl":"https://doi.org/10.1190/geo2023-0328.1","url":null,"abstract":"The application of direction vectors in the generation of reverse time migration (RTM) angle gathers in complex acoustic anisotropic media often encounters three main challenges: not pointing to the phase-velocity direction (PVD) of the Poynting vector, inaccuracy due to overlapping wavefields, and instability due to zero points of the direction vector. In general anisotropic media, the normally used Poynting vector indicates the group-velocity direction (GVD), whereas reflection and transmission phenomena rely on the PVD. Anisotropy introduces discrepancies between the GVD and the PVD. To overcome this issue, we employ the so-called PVD vector to directly calculate the PVD from anisotropic wavefields, eliminating the need of the approxi- mated conversion from the GVD to the PVD. To mitigate the inaccuracy problem, we apply the Hilbert transform based wavefield decomposition method to separate over- lapping wavefields into their up/down components, and then we calculate the PVDs using the separated wavefields. To tackle the instability problem, we incorporate the additionally simulated quadrature wavefield during the wavefield decomposition procedure. By combining the direction vector of the quadrature wavefield with that of the original wavefield, we can eliminate the zero points and thus obtain a stabi- lized PVD vector. With those problems solved or alleviated, we present a scheme for the generation of anisotropic RTM angle gathers in complex areas. Two numerical examples utilizing synthetic data sets demonstrate our method’s effectiveness.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"95 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138988314","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Modeling and numerical simulation of a type of pure-viscoacoustic-wave equation in attenuated TI media","authors":"Rong Huang, Zhiliang Wang, Guojie Song, Dan Wang, Xinmin Zhang, Fan Min","doi":"10.1190/geo2022-0670.1","DOIUrl":"https://doi.org/10.1190/geo2022-0670.1","url":null,"abstract":"The anisotropy and attenuation features of subsurface media significantly affect seismic data processing. Ignoring anisotropy and attenuation in seismic wave propagation may result in inaccurate reflector positions, dimming amplitudes, and reduced spatial resolution in the imaging results. Therefore, accurate seismic wave modeling of anisotropy and attenuation is essential for understanding wave propagation in the Earth’s interior. This article derives three pure-viscoacoustic-wave equations from characterizing the Earth’s frequency-independent Q behavior in transversely isotropic (TI) media. Firstly, we propose three time-space domain pure-qP-wave equations in TI media based on different approximation methods, whose coefficients can be determined by different approximation methods and Thomsen’s anisotropic parameters ε, δ. Subsequently, we introduce the Kelvin-Voigt attenuation model into our derived three time-space domain pure-qP-wave equations and then obtain three pure-viscoacoustic-wave equations. To further demonstrate the effectiveness and accuracy of our methods, we give some 2D and 3D numerical tests. The numerical results indicate that the wavefield generated by pure-qP-wave equations and pure-viscoacoustic-wave equations have accurate kinematic characteristics of qP-wave in TI media and attenuated TI media and are free of S- wave artifacts while remaining stable under Thomsen’s anisotropic parameters ε< δ, so our methods have broader applicability compared with some existing methods. At the same time, simulation results of pure-viscoacoustic-wave equations also reflect the absorption and attenuation characteristics of qP-waves in attenuated TI media.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":"91 10","pages":""},"PeriodicalIF":3.3,"publicationDate":"2023-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139025271","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}