{"title":"Enhancing the seismic response of faults by using a deep learning-based method","authors":"Hao Yan, Zhe Yan, Jiankun Jing, Zheng Zhang, Haiying Li, Hanming Gu, Shaoyong Liu","doi":"10.1111/1365-2478.13549","DOIUrl":"10.1111/1365-2478.13549","url":null,"abstract":"<p>The accuracy of fault interpretation is generally influenced by the quality of seismic images. Because of the blurring effect of the migration process, faults with small throws may not be clearly imaged in seismic images, which will impose limitations on the fault detection. To address this issue, we propose a deep learning-based method to enhance faults in poststack seismic images. We generate abundant training samples by convolving the three-dimensional point-spread functions with the noisy reflectivity models. The corresponding labels are synthesized using the one-dimensional seismic wavelet convolution method, simulating conditions with perfect illumination. To train the network for optimal performance, we investigate the impact of different loss functions. Ultimately, we employ a mixed loss function combining structural similarity index measure and gradient difference loss, since the gradient difference loss focuses more on geological edge information, and the structural similarity index measure possesses excellent image perceptual capability and optimization property. Results from one synthetic seismic image and three real seismic data demonstrate that our proposed method can effectively restore the sharpness of fault surfaces, particularly for faults with small displacements. Compared to the structural smoothing method, the network we trained achieves optimal fault enhancement. Furthermore, coherence-based fault images indicate that seismic images enhanced using our method can improve the accuracy of fault interpretation and yield more continuous fault maps.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2615-2633"},"PeriodicalIF":1.8,"publicationDate":"2024-06-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141353600","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 migration of water-bottom-related multiples accelerated by random phase-encoding strategy","authors":"Yanbao Zhang, Feng Hu, Yanzhi Hu","doi":"10.1111/1365-2478.13552","DOIUrl":"10.1111/1365-2478.13552","url":null,"abstract":"<p>Marine seismic multiples contain more structure information than primaries and should be considered in migrations. However, multiple migrations suffer from severe crosstalks generated by interferences among undesirable multiples. It has been proven that the water-bottom-related multiple migration can suppress crosstalks greatly. However, if all associated consecutive-order multiples are considered, the computation cost is extremely high. To settle this issue, a phase-encoding-based multiple migration is proposed. Supergathers are first created by randomly phase-encoding consecutive-order multiples and stacking-encoded multiples. By migrating supergathers, the proposed method can fulfil migrations of all order multiples simultaneously, thereby reducing the computation cost significantly. We use a three-layer model and the Pluto 1.5 model for numerical comparisons. The results reveal that the method can retrieve high-quality images and increase computation efficiency considerably.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2634-2645"},"PeriodicalIF":1.8,"publicationDate":"2024-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141359901","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}
Guangui Zou, Meijiao Wang, Tailang Zhao, Jiasheng She
{"title":"Assessing the impact of pore-fracture structures on permeability sensitivity in tectonic coal using computerized tomography scanning","authors":"Guangui Zou, Meijiao Wang, Tailang Zhao, Jiasheng She","doi":"10.1111/1365-2478.13551","DOIUrl":"10.1111/1365-2478.13551","url":null,"abstract":"<p>The heterogeneity in permeability of coal reservoirs is primarily attributed to the considerable variation in the morphologies and structures of microscopic pore-fractures, shaped by complex geological processes. This study emphasizes the necessity of understanding the impact and governance of these morphological and structural variations in pore-fractures across different types of coal bodies on their permeability. Utilizing computerized tomography scanning and three-dimensional imaging, we examined coal samples from the Datong coalfield in the southeastern Qinshui Basin, Shanxi Province, to characterize the pore-fracture morphologies and structures distinct to various coal-body types based on tomographic data. This introduces a methodology for assessing the influence of microscopic pore-fracture parameters, such as porosity, specific surface area, tortuosity and fractal dimension, on permeability sensitivity. This is achieved through the application of the modified Kozeny–Carman equation and a fractal permeability model. Findings indicate a predominance of slab fractures in raw coal, whereas fragmented coal under weak brittle deformation exhibits small, isolated pore-fractures with minimal diameter and volume and poor connectivity. In contrast, granular coal subjected to strong brittle deformation features extensive clusters of large pore-fractures with significant diameter and volume, enhancing connectivity. Moreover, permeability predictions are refined by integrating the modified Kozeny–Carman equation with tomographic data, offering a more precise understanding of the permeability across different coal bodies.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2772-2783"},"PeriodicalIF":1.8,"publicationDate":"2024-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141382750","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}
Xinpeng Pan, Zhishun Liu, Pu Wang, Lei Huang, Jianxin Liu
{"title":"Bayesian linearized amplitude variation with offset and azimuth inversion and uncertainty analysis in horizontal transversely isotropic media","authors":"Xinpeng Pan, Zhishun Liu, Pu Wang, Lei Huang, Jianxin Liu","doi":"10.1111/1365-2478.13548","DOIUrl":"10.1111/1365-2478.13548","url":null,"abstract":"<p>The stratum can be modelled as a horizontal transversely isotropic medium when a single set of vertically parallel fractures embedded in an isotropic background medium, which facilitates efficient study for fractured reservoirs. Elastic parameters and fracture weaknesses are important parameters to describe the characteristics of fractured reservoirs, and seismic inversion plays a significant role in parameters estimation. The commonly used deterministic inversion methods do not fully utilize the prior information and fails to present the uncertainty analysis of inversion results. To address these shortcomings, we propose a Bayesian linearized amplitude variation with offset and azimuth inversion method tailored for horizontal transversely isotropic media, enabling a more robust analysis of uncertainty. Within the framework of Bayesian inversion, the proposed method successfully derives analytical expressions for the posterior mean and covariance of both elastic parameters and fracture weaknesses. The response characteristics of the anisotropic reflection coefficient are analysed, and it is found that the perturbations of elastic parameters have a greater effect on reflection coefficient compared to fracture weaknesses. Synthetic data examples confirm that the accuracy of estimated P- and S-wave velocities and density surpasses that of fracture weaknesses, and the proposed method still performs well for the case of moderate noise. A field data example demonstrates that the inverted profiles agree well with the logging curve, and the estimated fracture weaknesses display significantly high values in the reservoir area. The estimated reservoir parameters not only contribute to a more accurate representation of the fractured gas-bearing reservoir but also provide insights into the target gas reservoir through its posterior distribution. Both synthetic and field data examples demonstrate the stability and reliability of the proposed method in characterizing fractured reservoirs. We determine that the proposed method provides an available tool for nuanced evaluation of uncertainty for the inversion results, and it is helpful for the fine description of fractured hydrocarbon-bearing reservoirs.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2666-2684"},"PeriodicalIF":1.8,"publicationDate":"2024-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141270500","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 shortest-path-aided fast-sweeping method to improve the accuracy of traveltime calculation in vertically transverse isotropic media","authors":"Jianming Zhang, Liangguo Dong, Chao Huang","doi":"10.1111/1365-2478.13537","DOIUrl":"10.1111/1365-2478.13537","url":null,"abstract":"<p>The high accuracy and efficiency of traveltime calculation are critical in seismic tomography, migration, static corrections, source locations and anisotropic parameter estimation. The fast-sweeping method is an efficient upwind finite-difference approach for solving the eikonal equation. However, the fast-sweeping method is accurate only along the axis directions. In two-dimensional or higher dimensional cases, the accuracy is severely decreased in the diagonal directions due to the numerical errors in these directions. These similar numerical errors also arose in higher order fast-sweeping method and anisotropic fast-sweeping method. To improve the accuracy of traveltime calculation in two-dimensional or higher dimensional space, a shortest-path-aided fast-sweeping method is proposed. The shortest-path-aided solution is embedded into the sweeping process of the standard fast-sweeping method to improve the traveltime accuracy in the diagonal directions. Shortest-path-aided fast-sweeping method is very easy to implement nearly without additional computational cost and memory consumption. Furthermore, this method is easy to extend from two-dimensional to higher dimensional, from low-order to higher-order and from isotropic to anisotropic cases.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2761-2771"},"PeriodicalIF":1.8,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141107908","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":"Diving waves in acoustic factorized orthorhombic media","authors":"Kristoffer Tesdal Galtung, Alexey Stovas","doi":"10.1111/1365-2478.13532","DOIUrl":"10.1111/1365-2478.13532","url":null,"abstract":"<p>Diving waves propagating in the subsurface are massive sources of low-frequency information that can be used to constrain the kinematic component of the velocity model. Compared to reflected waves, less is known about the behaviour of diving waves, especially in the presence of azimuthal anisotropy. Anisotropy is needed to place the events to the correct depths and match travel times in synthetics with recorded data. Obtaining more insights into the influence of anisotropy on diving wave propagation can help to find parameters with a low trade-off for inversion. Here, we derive equations for diving qP-waves in an acoustic factorized anisotropic model with orthorhombic anisotropy. The effects of the anisotropic parameters in the acoustic factorized orthorhombic model are tested by perturbing <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>ε</mi>\u0000 <mn>1</mn>\u0000 </msub>\u0000 <annotation>$epsilon _1$</annotation>\u0000 </semantics></math>, <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>ε</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <annotation>$epsilon _2$</annotation>\u0000 </semantics></math>, <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>η</mi>\u0000 <mn>1</mn>\u0000 </msub>\u0000 <annotation>$eta _1$</annotation>\u0000 </semantics></math>, <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>η</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <annotation>$eta _2$</annotation>\u0000 </semantics></math> and <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>η</mi>\u0000 <mn>3</mn>\u0000 </msub>\u0000 <annotation>$eta _3$</annotation>\u0000 </semantics></math> and observing differences in the ray paths, the effective vertical slowness and the relative geometrical spreading. The properties of diving waves in this model are also compared with those in an acoustic isotropic model and acoustic factorized anisotropic models with elliptical- and vertical transverse isotropic anisotropy. From our analysis, we found that perturbing <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>ε</mi>\u0000 <mn>1</mn>\u0000 </msub>\u0000 <annotation>$epsilon _1$</annotation>\u0000 </semantics></math> and <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>ε</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 <annotation>$epsilon _2$</annotation>\u0000 </semantics></math> has the most significant influence on these characteristics. The <span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>η</mi>\u0000 <mn>1</mn>\u0000 ","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 6","pages":"2300-2316"},"PeriodicalIF":2.6,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140941913","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}
Saulo Pomponet Oliveira, Luan Thanh Pham, Roman Pašteka
{"title":"Regularization of vertical derivatives of potential field data using Morozov's discrepancy principle","authors":"Saulo Pomponet Oliveira, Luan Thanh Pham, Roman Pašteka","doi":"10.1111/1365-2478.13534","DOIUrl":"10.1111/1365-2478.13534","url":null,"abstract":"<p>The calculation of the vertical derivatives of potential field methods can be carried out in a stable manner by Tikhonov regularization, but this procedure requires the appropriate selection of a regularization parameter. For this purpose, we introduce a criterion based on Morozov's discrepancy principle that uses a preliminary approximation given by the vertical derivative of the smoothed data. The smoothing may be performed by a physical or a mathematically based low-pass filter. The filtered data are computed only for estimating the regularization parameter; once it is found, we evaluate the regularized vertical derivative from the original data (not from the smoothed one) in the frequency domain. We verified from experiments with noise-corrupted synthetic data, as well as gravity and magnetic field data, that the regularized vertical derivative has about the same smoothness as the one obtained from filtered data, but true anomalies are more easily distinguished from noise and the shapes of the anomalies are better preserved.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 8","pages":"2880-2892"},"PeriodicalIF":1.8,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140886208","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}
Muhammad Izzatullah, Abdullah Alali, Matteo Ravasi, Tariq Alkhalifah
{"title":"Physics-reliable frugal local uncertainty analysis for full waveform inversion","authors":"Muhammad Izzatullah, Abdullah Alali, Matteo Ravasi, Tariq Alkhalifah","doi":"10.1111/1365-2478.13528","DOIUrl":"10.1111/1365-2478.13528","url":null,"abstract":"<p>Full waveform inversion stands at the forefront of seismic imaging technologies, pivotal in retrieving high-resolution subsurface velocity models. Its application is especially profound when imaging complex geologies such as salt bodies, which are regions notoriously challenging, yet essential given their hydrocarbon potential. However, with the power of full waveform inversion comes the intrinsic challenge of estimating the associated uncertainties. Such uncertainties are crucial in understanding the reliability of subsurface models, particularly in terrains like subsalt regions. Addressing this, we advocate for a nuanced approach employing the Stein variational gradient descent algorithm. Through a judicious use of a limited number of velocity model particles and the integration of random field-based perturbations, our methodology provides a local representation of the uncertainties inherent in full waveform inversion. Our evaluations, based on the Marmousi model, showcase the robustness of the proposed technique. Yet, it is our exploration into salt-intensive terrains, leveraging data from the Sigsbee 2A synthetic model and the Gulf of Mexico, that emphasizes the method's versatility. Findings indicate pronounced uncertainties along salt boundaries and in the deeper subsalt sediments, contrasting the minimal uncertainties in non-salt terrains. However, anomalies like salt canyons present unique challenges, potentially due to the interplay of multi-scattering effects. Emphasizing the scalability and cost-effectiveness of this approach, we highlight its potential for large-scale industrial applications in full waveform inversion, while also underscoring the necessity for prudence when integrating these uncertainty insights into subsequent seismic-driven geological and reservoir modelling.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2718-2738"},"PeriodicalIF":1.8,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140886360","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":"Introducing the area under stress–velocity curve: Theory, measurement and association with rock properties","authors":"Javad Sharifi","doi":"10.1111/1365-2478.13525","DOIUrl":"10.1111/1365-2478.13525","url":null,"abstract":"<p>Since many years ago, ultrasonic velocity has been used to investigate the physical and mechanical behaviour of rocks, thereby playing an important role in reservoir characterization and seismic interpretation. In order to develop the knowledge of ultrasonic tools, I performed a noble analysis on the ultrasonic behaviour of rocks under confining stress and evaluated a distinctive property of porous media that is measured as the area under the stress–velocity curve (here defined as <i>S<sup>*</sup></i>). I further investigated its relationship with elastic and mechanical behaviours of rock. To validate the theoretical framework developed in this work, 20 core plugs from various rock units with complex microstructures were subjected to triaxial compressional tests to calculate their area under the curve. Calculations were made for crack-closing, elastic and post-elastic stages (e.g. pore collapse) along the ultrasonic velocity–stress curve. Moreover, the selected samples had their microstructure investigated by thin-section studies to quantify their porosity and pore type. The results were analysed to check for the effect of pore type on <i>S<sup>*</sup></i> in different stages of the stress–velocity curve. Based on the outputs of the analysis of variance and Pearson's correlation coefficient analysis, the curve had its shape and underlying area closely related to the porosity and pore geometry. Indeed, the results showed that the shale and sandstone with micro cracks and carbonate with stiff pores correspond to smaller and larger areas under the curve in crack-closing and inelastic stages, respectively. Cross-correlating the results to compressibility (inverse of bulk modulus), it was figured out that the calculated area under curve was well consistent with the compressibility. In addition, <i>S<sup>*</sup></i> represents both static and dynamic behaviours of the rock, and the results revealed that the shape and curvature of the stress–velocity curve give valuable information about the rock microstructure. Another finding was the fact that the type of fluid and wave velocity seemingly affect the <i>S<sup>*</sup></i>. Our findings can help interpret wave velocity behaviour in reservoir rocks and other stressful porous media.</p>","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"72 7","pages":"2685-2700"},"PeriodicalIF":1.8,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140886442","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 new highly accurate and efficient pure visco‐acoustic wave equation for tilted transversely isotropic attenuating media","authors":"Lei Xiang, Jianping Huang, Qiang Mao, Xinru Mu","doi":"10.1111/1365-2478.13519","DOIUrl":"https://doi.org/10.1111/1365-2478.13519","url":null,"abstract":"The propagation of seismic waves in attenuating anisotropic media exhibits amplitude dissipation and phase dispersion. To describe its effects, the fractional Laplacian pure visco‐acoustic wave equations capable of producing stable and noise‐free wavefields have been derived. However, except for acoustic approximation, previous wave equations utilize the approximations with lower accuracy in simplifying the denominator of the approximate complex‐valued dispersion relation, resulting in reduced accuracy. To address this concern, we use a combination of complex stiffness coefficients to replace the denominator term of the approximate complex‐valued dispersion relation. This approximation effectively reduces the loss of accuracy caused by ignoring the influence of the velocity anisotropy parameter <jats:italic>ε</jats:italic> and the attenuation anisotropy parameter <jats:italic>ε<jats:sub>Q</jats:sub></jats:italic> in the denominator term, leading to a wave equation with high accuracy in media with large anisotropic parameters <jats:italic>ε</jats:italic> and <jats:italic>δ</jats:italic>. In addition, the new wave equation only contains two high‐order spatial partial derivatives and has high computational efficiency. Theoretical analysis and numerical examples demonstrate that the proposed pure visco‐acoustic tilted transversely isotropic wave equation outperforms the previous pure visco‐acoustic wave equation in terms of simulation accuracy. The newly developed wave equation is well suited for the application of <jats:italic>Q</jats:italic>‐compensated reverse time migration and full waveform inversion in attenuating anisotropic media.","PeriodicalId":12793,"journal":{"name":"Geophysical Prospecting","volume":"18 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140886364","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}