首页 > 最新文献

Geophysics最新文献

英文 中文
Sparse seismic data regularization in both shot and trace domains using a residual block autoencoder based on the fast Fourier transform 利用基于快速傅立叶变换的残差块自动编码器在射域和震迹域实现稀疏地震数据正则化
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-09 DOI: 10.1190/geo2023-0097.1
Alexandre L. Campi, R. Misságia
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.
稀疏采集在地震数据采集中的应用越来越多,在节省成本和时间方面具有优势。然而,它会导致地震数据采样不规则,对最终图像的质量产生不利影响。本文提出了一种基于傅里叶变换的残差块卷积神经网络——ResFFT-CAE网络。结合残余块可以使网络从地震数据中提取高频和低频特征。高频特征捕获详细信息,而低频特征整合整体数据结构,有助于在迹线和射孔域中更好地恢复不规则采样的地震数据。我们在综合数据和现场数据上评估了ResFFT-CAE网络的性能。在综合数据上,我们将ResFFT-CAE网络与利用曲线变换的压缩感知(CS)方法进行了比较。对于现场数据,我们与其他神经网络进行了比较,包括卷积自编码器(CAE)和U-Net。结果表明,在所有场景中,ResFFT-CAE网络始终优于其他方法。它产生了高质量的图像,其特点是低残差和减少畸变。此外,在评估模型泛化时,使用合成数据训练的模型进行的测试也显示出有希望的结果。综上所述,ResFFT-CAE网络作为一种对不规则采样地震数据进行正则化的高效工具,具有广阔的应用前景。其优异的性能在地震数据分析和处理流程的预处理方面具有潜在的应用前景。
{"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}
引用次数: 0
Full-waveform inversion of pure quasi-P waves in titled transversely isotropic media based on different parameterization 基于不同参数化的倾斜横向各向同性介质中纯准 P 波的全波形反演
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-09 DOI: 10.1190/geo2023-0367.1
Zhiming Ren, Xue Dai, Lei Wang
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.
全波形反演(FWI)通过最小化建模数据与观测数据之间的残差来建立地下参数模型。考虑各向异性的影响对于复杂结构的高分辨率成像至关重要。我们开发了一种基于纯准 P 波(qP 波)方程的声学各向异性 FWI 方法。方程系数及其与汤姆森各向异性参数([公式:见正文]和[公式:见正文])的导数是通过最小二乘(LS)优化估算的。我们推导出了函数梯度,并分析了六类参数的辐射模式:沿对称轴的速度[式:见正文]、[式:见正文]和[式:见正文],法向移动速度[式:见正文],各向异性参数[式:见正文]和[式:见正文],水平速度[式:见正文]和[式:见正文]:见正文]、水平速度[式:见正文]、[式:见正文]和[式:见正文]、[式:见正文]、[式:见正文]和[式:见正文]、[式:见正文]和[式:见正文]、[式:见正文]和[式:见正文]、[式:见正文]和[式:见正文]。参数化[式:见正文]在中间角和宽散射角时在[式:见正文]和[式:见正文]之间有明显的折衷。各向异性参数[式:见正文]和[式:见正文]在短散射角时分别可通过参数化[式:见正文]和[式:见正文]解决。公式:见正文]和[公式:见正文]的参数串扰比其他类型的参数串扰更为严重。我们对垂直和倾斜横向各向同性(VTI 和 TTI)介质中的纯 qP 波进行了全波反演。Overthrust VTI 模型和改进的 BP TTI 模型的反演结果表明,当其他参数足够精确时,可以单独重建速度、各向异性参数和倾斜角。多参数 FWI 无法为每种参数化类型获得可靠的倾斜角。用参数化[公式:见正文]反演得到的[公式:见正文]、[公式:见正文]和[公式:见正文]模型精度不高,而参数化[公式:见正文]有助于提高[公式:见正文]和[公式:见正文]模型的精度。
{"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}
引用次数: 0
Ordering cross-spread gathers 订购横幅式集装袋
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-08 DOI: 10.1190/geo2023-0161.1
Stewart Trickett
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.
交叉展布采集的地震处理是一种有价值但未得到充分利用的策略。要正确地进行处理,来自单个源线的源和来自单个接收线的接收器必须以物理上合理的方式排序,以便地表图上相邻的源或接收器在物理上彼此靠近。但是,对于不规则采集的陆地数据来说,确定这样的排序是一项挑战。我提出了新颖的自动排序算法,利用图论工具,最大限度地减少大的间隙,并保留即使在高度不规则采集中也能发现的顺序模式。提供 Java 源代码。
{"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}
引用次数: 0
A frequency-Hankel transform method to extract multimodal Rayleigh wave dispersion spectra from active and passive source surface wave data 从主动源和被动源表面波数据中提取多模态瑞利波频散谱的频率-汉克尔变换方法
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-08 DOI: 10.1190/geo2023-0189.1
Zhentao Yang, Yao-Chong Sun, Dazhou Zhang, Peng Han, Xiaofei Chen
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.
雷利波频散能谱已被广泛用于提取频散曲线和反演地下剪切波速度结构,用于工程地球物理和地震学。我们提出了一种频率-汉克尔(F-H)变换方法,用于从主动源和被动源瑞利波数据中提取高质量的多模态瑞利波频散能谱。F-H 变换方法受到频率-贝塞尔(F-J)变换方法的启发,并考虑了用于瑞利波频散分析的格林斯函数的物理意义。与 F-J 变换法相比,F-H 变换法可以自然地避免 F-J 频谱图上的会聚波造成的交叉伪影,并能以更少的瑞利波数据获得更多相同质量的多模态频散谱。为了证明 F-H 变换方法的有效性、准确性和适用性,我们使用了近地表勘探和环境噪声层析成像中主动源和被动源的合成和现场雷利波数据。F-H 变换方法统一了 F-J 变换方法及其对主动源和被动源雷利波数据的修改。F-H 变换法是一种稳健、高效的多模态瑞利波频散分析方法,适用于主动源和被动源瑞利波数据。
{"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}
引用次数: 0
Image-domain seismic inversion by deblurring with invertible Recurrent Inference Machines 利用可逆循环推理机去模糊进行图像域地震反演
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-08 DOI: 10.1190/geo2022-0780.1
Haorui Peng, Ivan Vasconcelos, M. Ravasi
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.
在复杂的地质环境和存在稀疏采集系统的情况下,地震偏移图像表现为未知地下模型的非平稳模糊版本。因此,图像域去模糊是产生可解释和高分辨率地下模型的重要步骤。大多数去模糊方法侧重于通过局部Hessian近似的迭代最小二乘反演地震图像的底层反射率;这可以通过所谓的点扩散函数的直接建模或通过迁移-反迁移过程来获得。在这项工作中,我们采用了一种新的深度学习框架,基于可逆循环推理机(i-RIMs),它允许将任何逆问题作为由已知建模算子(在我们的情况下是与点扩散函数的卷积)通知的监督学习任务来处理:所提出的算法可以直接反演阻抗扰动模型的迁移图像,辅助平滑速度模型和建模算子的先验信息。由于i- rim受到正向运算符的约束,它们隐式地学习以训练数据驱动的方式塑造/正则化输出模型。因此,所得到的去模糊图像对数据中的噪声和光谱缺陷(例如,由于有限的采集)具有很强的鲁棒性。通过几个综合算例证明了i-RIM网络设计和前向算子在训练过程中所起的关键作用。最后,利用现场数据,我们发现基于i- rim的去模糊技术在从偏移地震图像中获得稳健、高质量的相对阻抗估计方面具有很大的潜力。我们的方法可能对未来基于深度学习的定量油藏表征和监测具有重要意义。
{"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}
引用次数: 0
Physical property characterization of rocks in the Bayan Obo REE-Nb-Fe deposit, China 中国巴彦鄂博 REE-Nb-Fe 矿床岩石的物理性质表征
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-06 DOI: 10.1190/geo2023-0439.1
Lili Zhang, Hongrui Fan, Jian Wang, Liang Zhao, Kuifeng Yang, Ya Xu, Yonggang Zhao, Xingwang Xu, Meizhen Hao, Zhanfeng Yang, Xianhua Li
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.
巴彦奥博矿床是世界上最大的稀土元素(REE)资源、第二大铌(Nb)资源,也是中国重要的铁(Fe)资源。评估该矿床的资源潜力已成为全球关注的焦点。岩石物理性质是地球物理勘探和地质建模的桥梁;这些参数的变化是地球物理应用取得成功的必要条件。REE 、Nb、铁和钾主要赋存于巴彦奥博组的白云岩和板岩中,REE 矿化在遗传上与碳酸盐岩有关。利用三维重建、一/二/三维核密度估计、散点图矩阵、三维直方图以及基于皮尔逊和最大信息系数的相关分析,对巴彦奥博的铁矿、板岩、白云岩和碳酸盐岩堤露头样品的三种物理性质(电阻率、极化率和磁感应强度(MS))进行了测量和统计分析。结果表明,铁矿石、铁矿化细粒白云岩和铁矿化板岩的电阻率主要较低,铁矿石和铁矿化细粒白云岩的 MS 值较高。MS值可将铁矿石与板岩区分开来;MS值和电阻率可将铁矿化细粒白云岩与碳酸盐岩尖晶石区分开来。物理性质和全岩地球化学(主要元素和微量元素)共同表明,铁矿石、板岩和白云岩的MS与TFe2O3含量成正相关,极化率与TFe2O3成正相关,SiO2含量与K2O成正相关,电阻率分别与MS和极化率成正相关。铁矿石和白云石的电阻率与 TFe2O3 含量呈负相关。铁矿石的电阻率分别与 TFe2O3、稀土元素总量和铌呈负相关,与钍呈正相关。所采用的方法具有直观形象的表达效果,反映了物理性质的特点及其与矿物成分的相关性。研究结果将有助于确定矿床岩体的几何形状,为资源评价提供重要依据。
{"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}
引用次数: 0
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 通过新型自动拾取程序从衍射双曲线估算地下电磁速度分布:冰川学 GPR 数据集的理论与应用
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-06 DOI: 10.1190/geo2023-0042.1
M. Dossi, E. Forte, B. Cosciotti, S. Lauro, E. Mattei, E. Pettinelli, M. Pipan
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.
我们开发了一种自动挑选算法,旨在自动检测 GPR 数据集中的地下衍射体;准确跟踪源自已识别散射体的双曲衍射;恢复地下电磁速度分布,以及其他可能的分析。与其他常用的衍射跟踪技术相比,所提出的程序具有以下几个优势:只需对信号进行最少的预处理即可使用,因此用途更广,更能适应当地条件;只需解释器以几个跟踪参数阈值的形式提供有限的输入,因此结果更客观;与机器学习算法相比,它不涉及预训练,因此无需收集所有可能的地下情况的大型综合图像数据库,而这些情况不一定仅限于衍射实例。所介绍的算法首先要识别那些可能属于衍射顶点的信号,然后将其作为自动跟踪过程的初始种子。在自动跟踪过程中使用的水平搜索窗口是通过对每个衍射的大小进行粗略的初步估计而进行局部调整的。此外,同一顶点的多个种子可以产生多个可接受的双曲线,跟踪同一衍射相位。因此,该算法通过评估多个信号属性来选择最合适的双曲线,同时还能去除多余的双曲线和预期的误报。该算法应用于两个冰川学 GPR 剖面,能够准确跟踪绝大多数衍射记录,误报和漏报极少。这产生了一个统计上合理的电磁速度分布,用于评估勘测的高山冰川状况。
{"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}
引用次数: 0
Seismic Adaptive Multiple Subtraction Using a Structure-oriented Matched Filter 使用面向结构的匹配滤波器进行地震自适应多重减法
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-06 DOI: 10.1190/geo2023-0025.1
Yuhan Sui, Yue Ma, Lu Liu, Dongliang Zhang, Yubing Li
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.
在速度模型建立和成像之前,多次去除是地震数据处理的关键步骤。预测后,采用自适应多次减法抑制地震数据中的多次(考虑噪声),从而突出原色(考虑信号)。在实践中,传统的自适应减法方法使用滑动窗口对预测和记录的倍数进行最小二乘拟合,形成局部自适应匹配滤波器。随后,将过滤器应用于预测,以从记录的数据中删除倍数。然而,在最小化能量约束下,这种策略有过度衰减有用初级的风险。为了避免损坏有价值的信号,我们提出了一种新的方法,用面向结构的版本取代传统的匹配滤波器。从预测的倍数中提取结构信息,用于自适应匹配滤波器的推导。提出的面向结构的匹配滤波器强调预测倍数的结构,有助于在减法过程中更好地保留原色。合成和现场数据示例证明了所提出的面向结构的自适应减法方法的有效性,与传统方法相比,它在二维常规采样数据的多次去除和初级保存方面表现优异。
{"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}
引用次数: 0
Reverse time migration angle gathers in acoustic anisotropic media using direction vectors 利用方向矢量在声学各向异性介质中进行反向时间迁移角度采集
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-01 DOI: 10.1190/geo2023-0328.1
Kai Yang, Jianfeng Zhang
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.
在复杂的声学各向异性介质中应用方向矢量生成反向时间迁移(RTM)角度集合时,经常会遇到三个主要挑战:没有指向 Poynting 矢量的相位速度方向(PVD)、波场重叠导致的不准确性以及方向矢量零点导致的不稳定性。在一般各向异性介质中,通常使用的 Poynting 向量表示群速度方向(GVD),而反射和透射现象则依赖于 PVD。各向异性会造成 GVD 和 PVD 之间的差异。为了解决这个问题,我们采用了所谓的 PVD 向量来直接计算各向异性波场的 PVD,省去了从 GVD 到 PVD 的近似转换。为了缓解不准确问题,我们采用基于希尔伯特变换的波场分解方法,将重叠波场分离为上下分量,然后利用分离后的波场计算 PVD。为了解决不稳定性问题,我们在波场分解过程中加入了额外模拟的正交波场。通过将正交波场的方向矢量与原始波场的方向矢量相结合,我们可以消除零点,从而获得稳定的 PVD 矢量。随着这些问题的解决或缓解,我们提出了在复杂区域生成各向异性 RTM 角集的方案。利用合成数据集的两个数值示例证明了我们方法的有效性。
{"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}
引用次数: 0
Modeling and numerical simulation of a type of pure-viscoacoustic-wave equation in attenuated TI media 衰减 TI 介质中的一种纯声波方程的建模和数值模拟
IF 3.3 2区 地球科学 Q1 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2023-12-01 DOI: 10.1190/geo2022-0670.1
Rong Huang, Zhiliang Wang, Guojie Song, Dan Wang, Xinmin Zhang, Fan Min
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.
地下介质的各向异性和衰减特征对地震数据处理有很大影响。忽略地震波传播中的各向异性和衰减可能导致反射体位置不准确、振幅减弱以及成像结果的空间分辨率降低。因此,对各向异性和衰减进行精确的地震波建模对于理解地震波在地球内部的传播至关重要。本文通过描述横向各向同性(TI)介质中与频率无关的地球 Q 行为,推导出三个纯声波方程。首先,我们基于不同的近似方法提出了 TI 介质中的三个时空域纯 qP 波方程,其系数可由不同的近似方法和 Thomsen 各向异性参数 ε、δ 确定;然后,我们在推导出的三个时空域纯 qP 波方程中引入 Kelvin-Voigt 衰减模型,进而得到三个纯 viscoacoustic 波方程。为了进一步证明我们方法的有效性和准确性,我们进行了一些二维和三维数值测试。数值结果表明,纯 qP 波方程和纯 viscoacoustic 波方程产生的波场在 TI 介质和衰减 TI 介质中具有准确的 qP 波运动学特性,并且在 Thomsen 各向异性参数 ε< δ 的条件下保持稳定,没有 S 波的伪影,因此与现有的一些方法相比,我们的方法具有更广泛的适用性。同时,纯声波方程的模拟结果也反映了 qP 波在衰减 TI 介质中的吸收和衰减特性。
{"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}
引用次数: 0
期刊
Geophysics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1