Jie Fei Yang, Xia Luo, Dezhi Liu, Hanming Gu, Ming Sun
{"title":"Seismic noise attenuation method based on low-rank adaptive symplectic geometry decomposition","authors":"Jie Fei Yang, Xia Luo, Dezhi Liu, Hanming Gu, Ming Sun","doi":"10.1111/1365-2478.13504","DOIUrl":"10.1111/1365-2478.13504","url":null,"abstract":"<p>The basic assumption of low-rank methods is that noise-free seismic data can be represented as a low-rank matrix. Effective noise reduction can be achieved through the low-rank approximation of Hankel matrices composed of the data. However, selecting the appropriate rank parameter and avoiding expensive singular value decomposition are two challenges that have limited the practical application of this method. In this paper, we first propose symplectic geometric decomposition that avoids singular value decomposition. The symplectic similarity transformation preserves the essence of the original time sequence as well as the signal's basic characteristics and maintains the approximation of the Hankel matrix. To select an appropriate rank, we construct the symplectic geometric entropy according to the distribution of eigenvalues and search for high-contributing eigenvalues to determine the needed rank parameter. Therefore, we provide an adaptive approach to selecting the rank parameter by the symplectic geometric entropy method. The synthetic examples and field data results show that our method significantly improves the computational efficiency while adaptively retaining more effective signals in complex structures. Therefore, this method has practical application value.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2148-2163"},"PeriodicalIF":2.6,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313662","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":"Comparison of elastic anisotropy in the Middle and Upper Wolfcamp Shale, Midland Basin","authors":"Colin M. Sayers, Sagnik Dasgupta","doi":"10.1111/1365-2478.13503","DOIUrl":"10.1111/1365-2478.13503","url":null,"abstract":"<p>Organic-rich shales contain large amounts of oil and gas and are anisotropic because of fine-scale layering and the partial alignment of organic matter and anisotropic clay minerals with the bedding. An example is the Wolfcamp Shale in the Permian Basin. Elastic anisotropy needs to be accounted for in the characterization of such formations using seismic data and plays a role in hydraulic fracturing and in the evaluation of stress changes and geomechanical effects resulting from production. Using extensive well log data acquired in the Midland Basin, the eastern sub-basin of the Permian Basin, we estimate and compare the elastic anisotropy in the Middle and Upper Wolfcamp Shale by combining data from a vertical pilot well with two lateral wells, one (6SM) drilled in the Middle Wolfcamp and one (6SU) drilled in the Upper Wolfcamp. The data used were acquired at the Hydraulic Fracture Test Site 1, located in the eastern part of the Midland Basin. Thomsen's anisotropy parameter <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 </semantics></math> calculated from the fast and slow shear sonic is higher on average for the 6SM lateral than for 6SU, consistent with there being less carbonate content in 6SM than in 6SU. However, the anisotropy parameter <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 </semantics></math> in some regions with higher carbonate content in well 6SU is higher than in well 6SM. This may indicate the influence of natural fractures. The primary set of steeply dipping fractures observed in the lateral wells at Hydraulic Fracture Test Site 1 acts to increase <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 </semantics></math> if the ratio of the normal-to-shear fracture compliance is less than about 0.5. Sub-horizontal fractures may also increase <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 </semantics></math> and could affect the vertical extent of hydraulic fractures. Relations between elastic moduli <i>C</i><sub>33</sub> and <i>C</i><sub>55</sub> in the Upper and Lower Wolfcamp in a vertical pilot well allow <i>C</i><sub>33</sub> to be predicted in a lateral well using measurements of <i>C</i><sub>55</sub> in that well. Comparison of Thomsen's anisotropy parameters <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 </semantics></math> and <span></span><math>\u0000 <semantics>\u0000 <mi>ε</mi>\u0000 <annotation>$varepsilon $</annotation>\u0000 </semantics></math>, with <span></span><math>\u0000 <semantics>\u0000 <mi>γ</mi>\u0000 <annotation>$gamma $</annotation>\u0000 ","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2317-2328"},"PeriodicalIF":2.6,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313673","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":"Research note: A comparison between normalized controlled-source electromagnetic field components and transfer functions as input data for three-dimensional non-linear conjugate gradient inversion","authors":"Paula Rulff, Thomas Kalscheuer","doi":"10.1111/1365-2478.13488","DOIUrl":"10.1111/1365-2478.13488","url":null,"abstract":"<p>Controlled-source electromagnetic methods are applied to survey the electrical resistivity distribution of the subsurface. This work compares normalized electromagnetic field components and transfer functions such as impedance tensors and vertical magnetic transfer functions generated by two independent source polarizations as input data for three-dimensional inversion. As most other available inversion codes allow for inverting only one of the mentioned input data types, it is unclear which data type is preferable for controlled-source electromagnetic inversion. Our three-dimensional non-linear conjugate gradient inversion code can handle both input data types, facilitating a comparison of normalized electromagnetic field components and transfer functions inversion. Examining inversion results for a three-dimensional synthetic model with two anomalies, we infer that the transfer functions inversion is favourable for recovering the overall resistivity distribution below the receiver sites in fewer iterations. The inversion of normalized electromagnetic field components produces a sharper image of the anomalies and may be capable of detecting the resistivity distribution below the extended sources, which comes at the price of introducing a more heterogeneous background resistivity in the model.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 5","pages":"2005-2012"},"PeriodicalIF":2.6,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13488","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313671","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}
Owen Rohwer Huff, Vemund Stenbekk Thorkildsen, Thomas Larsen Greiner, Jan Erik Lie, Andreas Kjelsrud Evensen, Aina Juell Bugge, Jan Inge Faleide
{"title":"Near offset reconstruction for marine seismic data using a convolutional neural network","authors":"Owen Rohwer Huff, Vemund Stenbekk Thorkildsen, Thomas Larsen Greiner, Jan Erik Lie, Andreas Kjelsrud Evensen, Aina Juell Bugge, Jan Inge Faleide","doi":"10.1111/1365-2478.13505","DOIUrl":"10.1111/1365-2478.13505","url":null,"abstract":"<p>Marine seismic data is often missing near offset information due to separation between the source and receiver cables. To solve this problem, a convolutional neural network is trained on synthetic seismic data to reconstruct the near offset gap. The synthetic data is created using a two-dimensional finite difference method within a heterogeneous velocity model. These synthetics are generated with a source-over-receiver acquisition geometry so that they contain complete near offset data. The convolutional neural network is then trained on input-target synthetic pairs where the inputs are common midpoint gathers with the near offset section removed, and the targets are the same gathers with the near offset section retained. Following training, the robustness of the method is investigated with regards to common midpoint data sorting, normal moveout correction and changes in the velocity model. It is found that training on common midpoint-sorted data results in 2.8 times lower error than training on shot gathers, that normal moveout correction of the training data makes no significant difference in error levels, and that the model can reconstruct realistic near offsets on synthetic data generated 10 km away within the heterogeneous velocity model. In field data testing, first a dataset with source-over-cable acquisition geometry from the Barents Sea is used to compare the reconstructed wavefields to ground truth values. Although the reconstructed amplitudes require minor scaling to match the true values, predictions on this dataset yield 2.5 times lower near offset reconstruction error compared to a simple Radon transform interpolation method. Furthermore, amplitude versus offset gradient and intercept sections from the Barents Sea dataset are estimated with half the error when including the convolutional neural network-predicted near offset data, compared to only using the conventionally-acquirable portion of the data (beyond 112.5 m of offset). In a secondary field data test, a conventional northern North Sea dataset is used to demonstrate how the method may be applied in practice. Here, the convolutional neural network generates more realistic predictions than the Radon method, and the gradient and intercept sections calculated using the convolutional neural network-predicted traces have higher signal-to-noise ratios than the sections calculated using only the original data. The combination of high-quality synthetic training data and interpolation in the common midpoint domain enables near offset reconstruction at significant depth (1 s of two-way traveltime or more), which is demonstrated in both synthetic and field examples.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2164-2185"},"PeriodicalIF":2.6,"publicationDate":"2024-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/1365-2478.13505","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140197629","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}
Ana P. O. Muller, Bernardo Fraga, Matheus Klatt, Jessé C. Costa, Clecio R. Bom, Elisangela L. Faria, Marcelo P. de Albuquerque, Marcio P. de Albuquerque
{"title":"Deep-salt: Complete three-dimensional salt segmentation from inaccurate migrated subsurface offset gathers using deep learning","authors":"Ana P. O. Muller, Bernardo Fraga, Matheus Klatt, Jessé C. Costa, Clecio R. Bom, Elisangela L. Faria, Marcelo P. de Albuquerque, Marcio P. de Albuquerque","doi":"10.1111/1365-2478.13506","DOIUrl":"10.1111/1365-2478.13506","url":null,"abstract":"<p>Delimiting salt inclusions from migrated images during the velocity model building flow is a time-consuming activity that depends on highly human-curated analysis and is subject to interpretation errors or limitations of the images and methods available. We propose a supervised deep learning based method to include three-dimensional salt geometries in the velocity models. We compare two convolutional networks – based on the U-Net architecture – which can process three-dimensional seismic data. One architecture uses three-dimensional convolutional kernels, and the other has convolutional long short-term memory units. Each architecture requires specific preprocessing steps which affects their training and predictive performance. Both proposed architectures use subsurface offset gathers obtained from reverse time migration with an extended imaging condition as input and are trained to predict the salt inclusions. The velocity model used in migration is a reasonable approximation of sediment velocity but without salt inclusions. Thus, the migration model and, consequently, the migrated images are inaccurate due to the absence of the salt inclusion in the model using just the sediment velocity information for the segmentation. A similar salt inclusion methodology was previously validated for two-dimensional approaches; we extend it to the three-dimensional case. Our approach relies on subsurface common image gathers to focus the sediments' reflections around the zero offset and spread salt reflections' energy over large subsurface offsets. The results show that both proposed network models can accurately delineate the salt bodies in the test models, but when evaluating the trained networks for the three-dimensional SEG/EAGE salt model, the architecture with convolutional long short-term memory units has proven to generalize better.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2186-2199"},"PeriodicalIF":2.6,"publicationDate":"2024-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140197628","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 numerical solution of anisotropic poroelastic wave equation in microseismicity: Graphic process unit acceleration and moment tensor implementation","authors":"Jing Zheng, Tiezhu Li, Jingyu Xie, Yuan Sun","doi":"10.1111/1365-2478.13500","DOIUrl":"10.1111/1365-2478.13500","url":null,"abstract":"<p>The accuracy and computational efficiency of full waveform forward modelling in poroelastic media are crucial for microseismic monitoring. It enables intuitive, precise and efficient simulation of subsurface responses, thereby improving the reliability of moment tensor inversion and seismic source mechanism interpretation. Additionally, it reflects the role of fluid effects in waveform evolution. In this paper, based on the Biot mechanism, we derived the first-order velocity–stress equation of poroelastic media and discretized the model using a rotated staggered grid. The rotated staggered grid can well adapt to anisotropic media with high contrast parameters. We provide the finite difference formula based on graphic process unit–acceleration and moment tensor and also provide the graphic process unit workflow for forward modelling of anisotropic poroelastic media. First, two models with different grid sizes were run based on single graphic process unit, 1-thread Central Processing Unit (CPU) and 16-thread CPU. The results show that the speedup factors are approximately 14.3 and 3.7, respectively, compared with the running time of 1-thread CPU and 16-thread CPU. Then, we compare and evaluate the response of three basic source mechanisms (isotropic expansion, double couple and compensated linear vector dipole) in the model. The comparison of analytical and numerical results verifies the effectiveness of the method. Furthermore, wavefield snapshots of two typical anisotropic media (vertical transversely isotropic and horizontal transversely isotropic) are analysed to correspond to different moment tensor sources. The results showed that the source mechanism does not change the isotropic and anisotropic plane and the wave travel time, but it does change the polarization amplitude of the velocity component. The attenuation of slow qP-wave increases along with the increase of the value of viscosity. The effect of permeability on wavefield appears with the opposite effect of viscosity. Finally, the seismic waveform differences between multi-layer elastic media and poroelastic media are compared and analysed. The results showed that the seismic wavefield and waveform of poroelastic media are more complex, the propagation speed of seismic waves is faster, but the attenuation of seismic waves is stronger.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2329-2344"},"PeriodicalIF":2.6,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140170765","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}
Zhuowei Li, Tongtong Mo, Jiawen Song, Benfeng Wang
{"title":"Deblending and interpolation of subsampled blended seismic data based on damped randomized singular spectrum analysis","authors":"Zhuowei Li, Tongtong Mo, Jiawen Song, Benfeng Wang","doi":"10.1111/1365-2478.13507","DOIUrl":"10.1111/1365-2478.13507","url":null,"abstract":"<p>When compared to traditional seismic data acquisition, irregular blended acquisition significantly promotes the acquisition efficiency. Yet, the blending noise of subsampled blended data introduces new obstacles for the subsequent processing of seismic data. Due to the predictability of linear events in the frequency–space domain, the constructed Hankel matrices exhibit low-rank characteristics. However, the blending noise of subsampled blended data increases the rank, so deblending and interpolation can be implemented via rank-reduction algorithms such as the singular spectrum analysis. The significant computing cost of the singular value decomposition, however, makes the traditional singular spectrum analysis inefficient. An alternative algorithm, known as the randomized singular spectrum analysis, employs the randomized singular value decomposition instead of the traditional singular value decomposition for rank-reduction. The randomized singular spectrum analysis significantly enhances the efficiency of the decomposition process, particularly when dealing with large Hankel matrices. There still remains some random noise when using the singular spectrum analysis or randomized singular spectrum analysis for subsampled blended data, because the noise subspace and signal subspace are coupled together. Thus, we incorporate a damping operator into the randomized singular value decomposition and propose a novel damped randomized singular spectrum analysis method. The damped randomized singular spectrum analysis combines the advantages of the randomized singular value decomposition and the damping operator to enhance the computational efficiency and suppress the remaining noise. Moreover, an iterative projected gradient descent strategy is introduced to achieve deblended and interpolated seismic data for subsequent processing. Examples from synthetic data and field data are used to demonstrate the effectiveness and superiority of the proposed damped randomized singular spectrum analysis method, which enhances the accuracy and efficiency during simultaneous deblending and interpolation.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2200-2213"},"PeriodicalIF":2.6,"publicationDate":"2024-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140170860","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}
Hongtao Wang, Jiangshe Zhang, Chunxia Zhang, Li Long, Weifeng Geng
{"title":"Automatic stack velocity picking using a semi-supervised ensemble learning method","authors":"Hongtao Wang, Jiangshe Zhang, Chunxia Zhang, Li Long, Weifeng Geng","doi":"10.1111/1365-2478.13492","DOIUrl":"10.1111/1365-2478.13492","url":null,"abstract":"<p>Picking stack velocity from seismic velocity spectra is a fundamental method in seismic stack velocity analysis. With the increase in the scale of seismic data acquisition, manual picking cannot achieve the required efficiency. Therefore, an automatic picking algorithm is urgently needed now. Despite some supervised deep learning–based picking approaches that have been proposed, they heavily rely on sufficient training samples and lack interpretability. In contrast, utilizing physical knowledge to develop semi-data-driven methods has the potential to efficiently solve this problem. Thus, we propose a semi-supervised ensemble learning method to reduce the reliance on manually labelled data and improve interpretability by incorporating the interval velocity constraint. Semi-supervised ensemble learning fuses the information of the estimated spectrum, nearby velocity spectra and few-shot manual picking to recognize the velocity picking. Test results of both the synthetic and field datasets indicate that semi-supervised ensemble learning achieves more reliable and precise picking than traditional clustering-based techniques and the currently popular convolutional neural network method.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 5","pages":"1816-1830"},"PeriodicalIF":2.6,"publicationDate":"2024-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140156481","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":"Determination of the magnetization direction via correlation between reduced-to-the-pole magnetic anomalies and total gradient of the magnetic potential with vertical magnetization","authors":"Xiange Jian, Shuang Liu, Zuzhi Hu, Yunxiang Liu, Hongzhu Cai, Xiangyun Hu","doi":"10.1111/1365-2478.13499","DOIUrl":"10.1111/1365-2478.13499","url":null,"abstract":"<p>The total magnetization of an underground magnetic source is the vector sum of the induced magnetization and the natural remanent magnetization. The direction of the total magnetization serves as important a priori information in the inversion and processing of magnetic data. We demonstrated that the total gradient of the magnetic potential with vertical magnetization constitutes the envelope of the vertical component of the magnetic field for all directions of the Earth's field and source magnetization. The total gradient of the magnetic potential with vertical magnetization and the reduction-to-the-pole field simultaneously tend to achieve maximum symmetry near the correct total magnetization direction. As a result, the total magnetization direction can be estimated by computing the correlations between the reduction-to-the-pole and the total gradient of the magnetic potential with vertical magnetization. The proposed method yields accurate magnetization directions in synthetic model examples. The total gradient of the magnetic potential with vertical magnetization is less susceptible to data noise than transforms which are derived from the high-order magnetic field derivatives or tensors. The estimation results are slightly affected by changes in the source magnetization direction. In a field example in the Weilasito region (North China), the reduction-to-the-pole fields calculated using the estimated magnetization directions are well centred over the source. The proposed method obtained a more focused magnetization direction than that of a three-dimensional magnetization vector inversion. The total gradient of the magnetic potential with vertical magnetization therefore provides a novel and accurate approach to determine the total magnetization direction from the total field anomaly in a variety of situations.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2377-2402"},"PeriodicalIF":2.6,"publicationDate":"2024-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140107840","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}