Geophysical Prospecting最新文献

{"title":"Fracture network characterization through fractal dimension and Gutenberg–Richter parameter: Decatur open-source dataset as a study case","authors":"Claudia Pavez-Orrego,&nbsp;Denisse Pastén,&nbsp;Rodrigo Estay","doi":"10.1111/1365-2478.13561","DOIUrl":"10.1111/1365-2478.13561","url":null,"abstract":"<p>The fractal formalisms are well known for providing new understandings regarding the geometrical, spatial, and temporal behaviour of seismicity. Particularly, the fractal dimensions give information about the seismic events self-organization and self-similarity. On the other hand, the Gutenberg–Richter value, known as the <i>b</i>-value, has shown through the years to give handy information regarding the statistical distribution of earthquakes, on-site physical parameters, and geomechanical inputs. The Gutenberg–Richter value (<i>b</i>) and the capacity and correlation fractal dimensions, (<i>D</i>₀ and <i>D</i>₂), of the spatial distribution of earthquake hypocentres interact mathematically for micro- and macro-events. From this interaction, it is possible to obtain new insights into the fracture network development and the microseismicity source characterization in terms of single fractures, fault planes, or densely fractured volumetric spaces. Here we show this interaction for the open-source Decatur CO₂ project seismicity catalogue, comparing it with the results obtained for a natural earthquake catalogue of Illinois, in the United States. The fractal dimension <i>D</i>₀ is calculated using two different methodologies: box-counting and correlation integral partitioning. This last method is also used to calculate <i>D</i>₂. The results presented in this study allow us to describe how the fracture network geometry influences the earthquake complexity. Together with the calculation of the <i>b</i>-value, we present clear indications which show that seismicity recorded in the Illinois tectonic environment partially follows the Aki relationship <i>D</i>₀ ∼ 2<i>b</i>, which is not the case for induced events. In addition, the induced earthquake dataset shows that <i>D</i>₂ &gt; <i>D</i>₀, an anomalous behaviour in terms of the fractal formalisms. All these facts might be used to establish spatial fracture network control techniques and seismicity-type distinctions in CO₂ injection sites located in highly active tectonic areas, respectively.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2844-2856"},"PeriodicalIF":1.8,"publicationDate":"2024-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13561","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141869456","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":"Research on 3D modelling of induced polarization in polarizability anisotropic media","authors":"Jiaxuan Ling,&nbsp;Wei Deng,&nbsp;Shiwei Wei,&nbsp;Qingrui Chen,&nbsp;Lihua He,&nbsp;Siqin Liu,&nbsp;Mengmeng Li","doi":"10.1111/1365-2478.13580","DOIUrl":"10.1111/1365-2478.13580","url":null,"abstract":"<p>To enhance the 3D numerical simulation of the induced polarization method within anisotropic media, our study employs the 2D Fourier transform technique. This technique is utilized to convert the 3D integral of the abnormal potential from the space domain into a 1D integral in the wave number domain. Subsequently, we apply the shape function integration method, which is founded on quadratic interpolation, to resolve the 1D integral equation effectively. This methodology significantly decreases the necessary computational resources and storage while simultaneously harnessing the high efficiency and accuracy of the 1D shape function integration method, as well as the high efficiency of the fast Fourier transform, optimizing the numerical simulation process of the induced polarization method. We validate the accuracy of our algorithmic approach using an equivalent uniform layered model. Furthermore, by employing the sphere model, we conduct a comparison of computation time with the finite element method, thereby demonstrating high efficiency of the proposed algorithm. Utilizing the OpenMP parallel algorithm, we confirm that the proposed algorithm has a high degree of parallelism. We also analyse the differences in the equivalent apparent resistivity and apparent polarizability for various electrical parameters, using a prismatic model as the basis for our analysis. Our results clearly indicate that the anisotropy of the polarizability exerts substantial influence on the observe data. Consequently, the implications of polarizability anisotropy are deemed critical and not be disregarded in the field detection applications.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"3129-3139"},"PeriodicalIF":1.8,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141869458","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":"Numerical simulation of acoustic fields in open boreholes generated by linear phased array acoustic transmitters driven by pulse compression signals","authors":"Shengyue Tao,&nbsp;Xiaohua Che,&nbsp;Wenxiao Qiao,&nbsp;Jiale Wang,&nbsp;Qiqi Zhao","doi":"10.1111/1365-2478.13582","DOIUrl":"10.1111/1365-2478.13582","url":null,"abstract":"<p>Acoustic logging is an important method used to determine formation velocities near boreholes. However, in practice, determining accurate formation velocities from acoustic logging data is challenging because of the presence of various noise interferences. Accordingly, a method to increase the amplitudes of refracted waves in open boreholes is proposed herein on the basis of the directional radiation technology of pulse compression signal–driven linear phased array acoustic transmitters. The waveforms generated by a Ricker monopole acoustic transmitter, linear frequency modulation monopole acoustic transmitter and pulse compression signal–driven linear phased array acoustic transmitter in a fluid-filled open borehole are numerically simulated by employing the finite-difference method. The effects of the pulse compression signal–driven linear phased array parameters on the amplitudes of the refracted compressional and shear waves are studied. Results show that borehole mode waves with the same velocities and dispersion characteristics can be determined using the pulse compression signal–driven linear phased array acoustic and Ricker monopole acoustic transmitters in fluid-filled open boreholes. Pulse compression signal–driven linear phased array acoustic transmitters leverage the advantages of pulse compression and phased array technologies, ensuring that a single element can radiate more acoustic energy, whereas pulse compression signal–driven linear phased array parameters can be modulated to further increase the amplitudes of the refracted compressional and shear waves. Compared with Ricker and linear frequency modulation monopole acoustic transmitters, pulse compression signal–driven linear phased array acoustic transmitters can provide downhole received waveforms of better quality and improved a signal-to-noise ratio of the mode wave dispersion curves obtained using the downhole received waveforms. Because pulse compression signal–driven linear phased array acoustic transmitters use linear frequency modulation drive signals of longer duration, the recording time required for the received waveforms is also longer and the amount of data generated is larger, presenting new challenges for downhole data processing and high-speed data transmission.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2542-2556"},"PeriodicalIF":1.8,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141873305","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":"Elastic full waveform inversion for tilted transverse isotropic media: A multi-step strategy accounting for a symmetry axis tilt angle","authors":"Hengli Song,&nbsp;Yuzhu Liu,&nbsp;Jizhong Yang","doi":"10.1111/1365-2478.13578","DOIUrl":"10.1111/1365-2478.13578","url":null,"abstract":"&lt;p&gt;Transversely isotropic media with a tilted symmetry axis (TTI) exits widely underground due to tectonic movement and mineral orientation. Traditional full waveform inversion (FWI) based on isotropic media or transversely isotropic media with a vertical symmetry axis (VTI) cannot deal with such situations. To address this limitation, TTI–based FWI was developed. However, its practical application faces challenges in estimating the symmetry axis tilt angle &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;. Previous studies have generally assumed that &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; is equal to the strata dip angle, which is incorrect in complex structures such as salt domes and magmatic intrusion zones. Another theoretically robust way to estimate &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; is to treat it as the parameter to be inverted, but there are still some problems unresolved. First, the parameter &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; increases the nonlinearity of the inversion process, and its impact mechanism on inversion is not yet clear. Second, there is severe crosstalk (also known as trade-off or coupling) between parameters, but the current parameter decoupling technique for TTI–based FWI is not mature. To address the first problem, we assess the interaction between &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and other parameters by analysing the radiation patterns in the TTI background. Our analysis reveals that &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;θ&lt;/mi&gt;\u0000 &lt;mi&gt;t&lt;/mi&gt;\u0000 &lt;/msub&gt;\u0000 &lt;annotation&gt;${{theta }_{mathrm{t}}}$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; is most coupled by S-wave vertical velocity &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;msub&gt;\u0000 &lt;mi&gt;v&lt;/mi&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mi&gt;s&lt;/mi&gt;\u0000 &lt;mn&gt;0&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;/msub&gt;\u0000 &lt;","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2486-2503"},"PeriodicalIF":1.8,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141869457","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":"Simultaneous seismic data de-aliasing and denoising with a fast adaptive method based on hybrid wavelet transform","authors":"Peng Zhang,&nbsp;Xiaoying Han,&nbsp;Changle Chen,&nbsp;Xinming Liu","doi":"10.1111/1365-2478.13574","DOIUrl":"10.1111/1365-2478.13574","url":null,"abstract":"<p>Missing data and random noise are prevalent issues encountered during the processing of acquired seismic data. Interpolation and denoising represent economical solutions to address these limitations. Recovering regularly missing traces is challenging because of the spatial aliasing, and the extra difficulty is compounded by the presence of noise. Hence, developing an effective approach to realize denoising and anti-aliasing is important. Projection onto convex sets is an effective method for recovering missing seismic data that is typically used for processing data with a good signal-to-noise ratio. The computational attractiveness of the projection onto convex sets reconstruction approach is compromised by its slow convergence rate. In this study, we aimed to efficiently implement simultaneous seismic data de-aliasing and denoising. We combined a discrete wavelet transform with a seislet transform to construct a hybrid wavelet transform. A new fast adaptive method based on the fast projection onto convex sets method was proposed to recover the missing data and remove random noise. This approach adjusts the projection operator and iterative shrinkage threshold operator. The result is influenced by the threshold value. We enhanced the processing accuracy by adopting an optimal threshold strategy. Synthetic and field data tests indicate the effectiveness of the proposed method.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"3038-3051"},"PeriodicalIF":1.8,"publicationDate":"2024-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141771074","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":"Resolution evaluation of the satellite altimetric gravity anomaly models with shipborne gravity data over the Cosmonaut Sea, East Antarctica (44–52° E)","authors":"Chunguo Yang,&nbsp;Linjiang Qin,&nbsp;Weifeng Ding,&nbsp;Jinyao Gao,&nbsp;Guochao Wu","doi":"10.1111/1365-2478.13577","DOIUrl":"10.1111/1365-2478.13577","url":null,"abstract":"<p>Several satellite gravity anomaly models are freely available to calculate the free-air gravity anomaly in areas where shipborne gravity measurements are scarce. Two models produced by the Technical University of Denmark (DTU17) and the Scripps Institution of Oceanography (SIOv32.1), respectively, were selected to compute the free-air anomalies over the Cosmonaut Sea, East Antarctica. A statistical comparison analysis was performed to evaluate the resolution of satellite gravity anomaly models by comparing them with the shipborne surveying date. The radially averaged energy spectra of free-air anomaly from different sources were calculated and compared over two selected regions to further evaluate the reliability of the data derived from satellite gravity anomaly models. The satellite gravity anomaly models have a better resolution in the ocean basin than in the area near the continental shelf. The comparison analysis revealed that the precision of both DTU17 and SIOv32.1 is close to the shipborne gravity data, but on average, SIOv32.1 is a little bit better than DTU17. The spectral analysis showed that the shipborne measurements may provide higher resolution than the satellite gravity anomaly model at wavelengths shorter than 20 km, and the free-air data derived from SIOv32.1 have better resolution than the one from DTU17. These shipborne datasets will provide contributions for the updates of the Antarctic gravity anomaly and enable new high-resolution combined Earth gravity models to be derived in Antarctica.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"3090-3103"},"PeriodicalIF":1.8,"publicationDate":"2024-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141785126","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":"Two-dimensional controlled source electromagnetic inversion algorithm based on a space domain forward modeling approach","authors":"Iktesh Chauhan,&nbsp;Rahul Dehiya","doi":"10.1111/1365-2478.13575","DOIUrl":"10.1111/1365-2478.13575","url":null,"abstract":"<p>We develop a two-dimensional controlled-source electromagnetic inversion algorithm employing a space domain forward modelling algorithm. The space domain forward modelling algorithm is devised by imposing boundary conditions on the plane perpendicular to the strike direction that passes through the source position. The boundary conditions for various source types are derived using the symmetric/antisymmetric character of the electric and magnetic fields. The benchmarking analysis reveals that roughly eight grids are sufficient for discretizing space in the strike directions for accurate forward response computations. For inverse modelling, the inexact Gauss–Newton optimization technique is utilized. Numerical inversion experiments of synthetic and real-field data clearly demonstrate the versatility and robustness of the developed algorithm. The inversion experimentations also concur with the forward response benchmarking analysis and suggest that only a few grids (around eight) are adequate to discretize space in the strike direction. The developed algorithm is more than one order efficient compared to a wavenumber domain code.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"3052-3066"},"PeriodicalIF":1.8,"publicationDate":"2024-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141771075","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":"Enhancing magnetic source edges using the tilt angle of the analytic-signal amplitudes of the horizontal gradient","authors":"Luan Thanh Pham,&nbsp;Richard S. Smith,&nbsp;Saulo P. Oliveira,&nbsp;Vinicius Theobaldo Jorge","doi":"10.1111/1365-2478.13573","DOIUrl":"10.1111/1365-2478.13573","url":null,"abstract":"<p>Enhancing magnetic data is often complicated due to the non-vertical orientation of the geomagnetic field and the orientation of remanent source magnetization. The complication can be reduced by reducing the data to the pole (mathematically making the geomagnetic field vertical), but this reduction process is problematic. The analytic-signal amplitude can be used to enhance the edges of two-dimensional sources without a reduction to the pole. However, the shape of the analytic-signal amplitude is weakly dependent on the magnetization direction for grid data. This study presents an improved technique, namely the tilt angle of the analytic-signal amplitudes of the horizontal gradient of the vertical integral. This quantity is also only weakly dependent on the magnetization direction and outlines the edges as well or somewhat better than other methods. It also implicitly involves second derivatives of the magnetic field, and we use synthetic data to demonstrate that noise is not amplified as much as it is when using other edge enhancement techniques that implicitly use second derivatives. A dataset of the Apiaí Terrane, Brazil, shows good lateral continuity of features compared with other edge-enhancement methods, and subtle features like faults are easier to identify in the images generated by our new method. Upward continuation of the field, which is normally required, was not necessary to reduce the impact of noise on this field example.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"3026-3037"},"PeriodicalIF":1.8,"publicationDate":"2024-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141771076","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":"Re-visible blind block network: An unsupervised seismic data random noise attenuation method","authors":"Jing Wang,&nbsp;Bangyu Wu,&nbsp;Hui Yang,&nbsp;Bo Li","doi":"10.1111/1365-2478.13559","DOIUrl":"10.1111/1365-2478.13559","url":null,"abstract":"<p>Noise is inevitable when acquiring seismic data, and effective random noise attenuation is crucial for seismic data processing and interpretation. Training and inferencing two-stage deep learning-based denoising methods typically require massive noisy–clean or noisy–noisy pairs to train the network. In this paper, we propose an unsupervised seismic data denoising framework called a re-visible blind block network. It is a training-as-inferencing one-stage method and utilizes only single noisy data for denoising, thereby eliminating the effort to prepare training data pairs. First, we introduce a global masker and a corresponding mask mapper to obtain the denoised result containing all blind block information, enabling simultaneous optimization of all blind blocks via the loss function. The global masker consists of two complementary block-wise masks. It is utilized to mask noisy data to obtain two corrupted data, which are then input into the denoising network for noise attenuation. The mask mapper samples the value of blind blocks in the denoised data and projects it onto the same channel to gather the denoised results of all blind blocks together. Second, the original noisy data are incorporated into the network training process to prevent information loss, and a hybrid loss function is employed for updating the network parameters. Synthetic and field seismic data experiments demonstrate that our proposed method can protect seismic signals while suppressing random noise compared with traditional methods and several state-of-the-art unsupervised deep learning denoising techniques.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2739-2760"},"PeriodicalIF":1.8,"publicationDate":"2024-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141739709","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 method for extracting P-SV-converted wave angle-domain common-image gathers based on elastic-wave reverse-time migration","authors":"Qianqian Ci,&nbsp;Bingshou He","doi":"10.1111/1365-2478.13571","DOIUrl":"10.1111/1365-2478.13571","url":null,"abstract":"<p>Multicomponent seismic technology utilizes the kinematic and dynamic characteristics of reflected P-waves and converted S-waves to reduce ambiguity in seismic exploration. The imaging and inversion accuracy of P-SV-converted waves are important in determining whether multicomponent seismic exploration can achieve higher exploration accuracy than conventional P-wave exploration. Pre-stack inversion of P-SV-converted waves requires precise input of P-SV-converted wave angle-domain common-image gathers. Consequently, the P-SV-converted wave angle-domain common-image gather extraction accuracy will significantly affect the P-SV-converted wave inversion accuracy. However, existing methods for extracting P-SV-converted wave angle-domain common-image gathers are constrained by issues such as the P- and S-wave crosstalk artefacts, low-frequency noises and inaccurate calculation of P-wave incident angles, leading to poor imaging accuracy. We study an angle-domain cross-correlation imaging condition and address three key issues based on this condition: the decoupling of P- and S-waves, the separation of up-going and down-going waves and the precise calculation of P-wave incident angles. Our strategies facilitate high-precision extraction of P-SV-converted wave angle-domain common-image gathers using elastic wave reverse-time migration. In this paper, first, we employ the first-order velocity-dilatation-rotation elastic wave equations to decouple P- and S-waves automatically during source and receiver wavefield extrapolations. Second, we calculate the optical flow vectors of P- and S-waves to ensure stable calculations of wave propagation directions. Based on this, we obtain up-going and down-going waves of P- and S-waves. Meanwhile, we calculate the incident angle of the source P-wave using geometric relations. Lastly, we apply the angle-domain imaging condition to achieve high-precision extraction of P-SV-converted wave angle-domain common-image gathers. Model examples demonstrate the effectiveness and advantages of the proposed method.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2469-2485"},"PeriodicalIF":1.8,"publicationDate":"2024-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141739708","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}