首页 > 最新文献

Geophysics最新文献

英文 中文
A depth-variant seismic wavelet extraction method for basis pursuit inversion with impedance trend constraint 基于阻抗趋势约束的基序反演深度变异地震小波提取方法
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-23 DOI: 10.1190/geo2023-0255.1
R. Cai, Chengyu Sun, Zhen’an Yao, Shizhong Li
The seismic images produced by pre-stack depth migration show more accurate subsurface structures than time images, resulting in a growing need for depth-domain inversion. However, due to the strong non-stationarity exhibited by depth-domain seismic data, time-domain inversion methods based on the convolutional model cannot be directly applied in the depth domain. To address this issue, we have developed a method for extracting a depth-variant seismic wavelet, which is then combined with a non-stationary convolutional model to enable direct inversion of the depth-domain acoustic impedance. First, we extend the Morlet wavelet to the depth domain and propose an orthogonal matching pursuit spectral decomposition method using the depth-domain Morlet wavelet. We then investigate the waveforms and wavenumber spectra similarities between the depth-domain Morlet wavelet and depth-domain Ricker wavelet and extract depth-variant Ricker wavelets from the depth-wavenumber spectrum. We add a depth-domain impedance trend constraint to the conventional basis pursuit inversion to enhance the lateral continuity of the inversion results. Then, we attain direct inversion of the depth-domain acoustic impedance. Tests of synthetic and field data demonstrate that the proposed method achieves high-accuracy inversion results while maintaining high computational efficiency, highlighting our approach's effectiveness and strong reservoir characterization potential.
叠前深度迁移产生的地震图像比时间图像能显示更精确的地下结构,因此对深度域反演的需求日益增长。然而,由于深度域地震数据表现出很强的非稳态性,基于卷积模型的时域反演方法无法直接应用于深度域。为解决这一问题,我们开发了一种提取深度变异地震小波的方法,然后将其与非稳态卷积模型相结合,从而实现深度域声阻抗的直接反演。首先,我们将 Morlet 小波扩展到深度域,并提出了一种使用深度域 Morlet 小波的正交匹配追寻频谱分解方法。然后,我们研究了深度域 Morlet 小波和深度域 Ricker 小波之间的波形和波谱相似性,并从深度波谱中提取了深度变异 Ricker 小波。我们在传统的基序反演中加入了深域阻抗趋势约束,以增强反演结果的横向连续性。然后,我们实现了深度域声阻抗的直接反演。对合成数据和现场数据的测试表明,所提出的方法既能获得高精度的反演结果,又能保持较高的计算效率,凸显了我们方法的有效性和强大的储层表征潜力。
{"title":"A depth-variant seismic wavelet extraction method for basis pursuit inversion with impedance trend constraint","authors":"R. Cai, Chengyu Sun, Zhen’an Yao, Shizhong Li","doi":"10.1190/geo2023-0255.1","DOIUrl":"https://doi.org/10.1190/geo2023-0255.1","url":null,"abstract":"The seismic images produced by pre-stack depth migration show more accurate subsurface structures than time images, resulting in a growing need for depth-domain inversion. However, due to the strong non-stationarity exhibited by depth-domain seismic data, time-domain inversion methods based on the convolutional model cannot be directly applied in the depth domain. To address this issue, we have developed a method for extracting a depth-variant seismic wavelet, which is then combined with a non-stationary convolutional model to enable direct inversion of the depth-domain acoustic impedance. First, we extend the Morlet wavelet to the depth domain and propose an orthogonal matching pursuit spectral decomposition method using the depth-domain Morlet wavelet. We then investigate the waveforms and wavenumber spectra similarities between the depth-domain Morlet wavelet and depth-domain Ricker wavelet and extract depth-variant Ricker wavelets from the depth-wavenumber spectrum. We add a depth-domain impedance trend constraint to the conventional basis pursuit inversion to enhance the lateral continuity of the inversion results. Then, we attain direct inversion of the depth-domain acoustic impedance. Tests of synthetic and field data demonstrate that the proposed method achieves high-accuracy inversion results while maintaining high computational efficiency, highlighting our approach's effectiveness and strong reservoir characterization potential.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139957217","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
Exploring the challenges of interpreting near-surface towed TEM data on saline permafrost 探索解读盐碱冻土近地表拖曳 TEM 数据的挑战
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-21 DOI: 10.1190/geo2023-0221.1
Thomas Højland Lorentzen, M. A. Kass, Johanna Scheer, Soňa Tomaškovičová, Anders Vest Christiansen, Pradip Kumar Maury, T. Ingeman‐Nielsen
This paper analyzes a continuous near-surface transient electromagnetic (TEM) survey on permafrost, around Ilulissat in western Greenland, an area characterized by continuous saline permafrost. The TEM data is severely affected by induced polarization (IP) effects, causing a large range of decay shapes. We identify seven unique decay shapes: “oversteepened”, ”sign-change“, ”all-negative”, “double-sign-change”, “no apparent IP”, “flat-spot”, and “positive-nonmonotonic”, of which the last two have not previously been identified in the scientific literature. A clear spatial dependency of the decay shapes is observed. Inversion of the data is carried out using a Cole-Cole model which poses a highly nonunique inversion problem with an extreme starting model dependency. A series of inversion and forward modeling experiments demonstrate these challenges, and show that a low-resistivity and highly chargeable layer with time constant ([Formula: see text]) values between [Formula: see text] and[Formula: see text] and frequency exponent (C) values above 0.74 is needed to fit the data used in the inversion experiment. Forward modeling further shows that low [Formula: see text] and high C values are needed to reproduce the observed flat-spot and positive-nonmonotonic decay shapes. Based on these observations, we attribute the IP effects to the orientational polarization of ice in the soil column. This mechanism allows for low-resistivity, high-chargeability layers due to partially frozen saline sediments; a combination which is difficult to explain by the IP mechanisms traditionally considered. Based on forward modeling of a realistic structural model of a layered sediment pack over dipping bedrock, we are able to reproduce all the observed decay shapes. This model provides a consistent framework for qualitative interpretation of the entire data set, and evidence for the presence of saline deposits in the central parts of the sedimentary basins.
本文分析了在格陵兰西部伊卢利萨特周围对永久冻土进行的连续近地表瞬变电磁(TEM)勘测,该地区的特点是存在连续的盐碱永久冻土。TEM 数据受到诱导极化(IP)效应的严重影响,导致衰减形状范围很大。我们发现了七种独特的衰减形状:"过度陡峭"、"符号变化"、"全负"、"双符号变化"、"无明显 IP"、"平点 "和 "正-非单调",其中后两种以前在科学文献中从未发现过。衰减形状具有明显的空间依赖性。使用科尔-科尔模型对数据进行反演,这是一个高度非唯一的反演问题,具有极端的起始模型依赖性。一系列反演和前向建模实验证明了这些挑战,并表明需要一个时间常数([公式:见正文])值介于[公式:见正文]和[公式:见正文]之间、频率指数(C)值高于 0.74 的低电阻率和高电荷层才能适合反演实验中使用的数据。前向建模进一步表明,需要低[式:见正文]和高 C 值来重现观测到的平斑和正非单调衰减形状。根据这些观测结果,我们将 IP 效应归因于土壤柱中冰的定向极化。这种机制允许部分冻结的盐沉积物产生低电阻率、高电荷率层;传统的 IP 机制很难解释这种组合。基于对倾斜基岩上的层状沉积物群的现实结构模型的前向建模,我们能够再现所有观测到的衰减形状。该模型为整个数据集的定性解释提供了一个一致的框架,并为沉积盆地中部存在盐渍沉积提供了证据。
{"title":"Exploring the challenges of interpreting near-surface towed TEM data on saline permafrost","authors":"Thomas Højland Lorentzen, M. A. Kass, Johanna Scheer, Soňa Tomaškovičová, Anders Vest Christiansen, Pradip Kumar Maury, T. Ingeman‐Nielsen","doi":"10.1190/geo2023-0221.1","DOIUrl":"https://doi.org/10.1190/geo2023-0221.1","url":null,"abstract":"This paper analyzes a continuous near-surface transient electromagnetic (TEM) survey on permafrost, around Ilulissat in western Greenland, an area characterized by continuous saline permafrost. The TEM data is severely affected by induced polarization (IP) effects, causing a large range of decay shapes. We identify seven unique decay shapes: “oversteepened”, ”sign-change“, ”all-negative”, “double-sign-change”, “no apparent IP”, “flat-spot”, and “positive-nonmonotonic”, of which the last two have not previously been identified in the scientific literature. A clear spatial dependency of the decay shapes is observed. Inversion of the data is carried out using a Cole-Cole model which poses a highly nonunique inversion problem with an extreme starting model dependency. A series of inversion and forward modeling experiments demonstrate these challenges, and show that a low-resistivity and highly chargeable layer with time constant ([Formula: see text]) values between [Formula: see text] and[Formula: see text] and frequency exponent (C) values above 0.74 is needed to fit the data used in the inversion experiment. Forward modeling further shows that low [Formula: see text] and high C values are needed to reproduce the observed flat-spot and positive-nonmonotonic decay shapes. Based on these observations, we attribute the IP effects to the orientational polarization of ice in the soil column. This mechanism allows for low-resistivity, high-chargeability layers due to partially frozen saline sediments; a combination which is difficult to explain by the IP mechanisms traditionally considered. Based on forward modeling of a realistic structural model of a layered sediment pack over dipping bedrock, we are able to reproduce all the observed decay shapes. This model provides a consistent framework for qualitative interpretation of the entire data set, and evidence for the presence of saline deposits in the central parts of the sedimentary basins.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140445149","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}
引用次数: 1
Improving Parsimonious Refraction Interferometry through U-Net-based Correction of First-Arrival Traveltimes 通过基于 U-Net 的首到旅行时间校正改进准假折射干涉测量法
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-21 DOI: 10.1190/geo2023-0435.1
Ganghoon Lee, Sukjoon Pyun
Parsimonious Refraction Interferometry (PRI) is a useful technique for generating dense virtual traveltimes when the seismic survey is performed with sparse acquisition geometry. However, PRI has limitations in accurately calculating the traveltimes of direct and diving waves, leading to inaccurate velocity structures when applying the first-arrival traveltime tomography (FATT). To address this issue, we propose an improved approach using a deep learning (DL) network called U-Net. We first examine the feasibility of the proposed algorithm analytically through the traveltime interpretation of a simple model. Then, the U-Net model is trained on various datasets to learn the relationship between PRI results and actual traveltimes. Subsequently, the trained network corrects the traveltime errors of the PRI results. As a result, we can obtain accurate first-arrival traveltimes using only two shot gathers without additional information such as infilled shots. The proposed technique enables virtual traveltime corrections, allowing for improved FATT results, even in cases where dense shots are difficult to deploy. Numerical results demonstrate that the proposed method can achieve comparable accuracy to picked traveltime data, indicating its high effectiveness in increasing the resolution of FATT results. The proposed approach maximizes the cost-saving benefits of PRI and can be advantageous in obtaining high-resolution FATT results when dense shot geometry is unavailable.#xD;
准折射干涉测量法(PRI)是一种有用的技术,可在地震勘探采用稀疏采集几何时生成密集的虚拟行波时间。然而,PRI 在精确计算直达波和潜波的旅行时间方面存在局限性,导致在应用首次到达旅行时间层析(FATT)时速度结构不准确。为解决这一问题,我们提出了一种使用深度学习(DL)网络(U-Net)的改进方法。我们首先通过对一个简单模型的旅行时间解释,分析检验了所提算法的可行性。然后,在各种数据集上训练 U-Net 模型,学习 PRI 结果与实际旅行时间之间的关系。随后,经过训练的网络会修正 PRI 结果的旅行时间误差。因此,我们只需使用两个镜头采集即可获得精确的首次到达旅行时间,而无需填充镜头等附加信息。所提出的技术可实现虚拟旅行时间校正,从而改进 FATT 结果,即使在难以部署高密度镜头的情况下也是如此。数值结果表明,所提出的方法可以达到与采到的旅行时间数据相当的精度,这表明它在提高 FATT 结果的分辨率方面非常有效。所提出的方法最大限度地发挥了 PRI 节省成本的优势,在没有高密度射电几何图形的情况下,也有利于获得高分辨率的 FATT 结果;
{"title":"Improving Parsimonious Refraction Interferometry through U-Net-based Correction of First-Arrival Traveltimes","authors":"Ganghoon Lee, Sukjoon Pyun","doi":"10.1190/geo2023-0435.1","DOIUrl":"https://doi.org/10.1190/geo2023-0435.1","url":null,"abstract":"Parsimonious Refraction Interferometry (PRI) is a useful technique for generating dense virtual traveltimes when the seismic survey is performed with sparse acquisition geometry. However, PRI has limitations in accurately calculating the traveltimes of direct and diving waves, leading to inaccurate velocity structures when applying the first-arrival traveltime tomography (FATT). To address this issue, we propose an improved approach using a deep learning (DL) network called U-Net. We first examine the feasibility of the proposed algorithm analytically through the traveltime interpretation of a simple model. Then, the U-Net model is trained on various datasets to learn the relationship between PRI results and actual traveltimes. Subsequently, the trained network corrects the traveltime errors of the PRI results. As a result, we can obtain accurate first-arrival traveltimes using only two shot gathers without additional information such as infilled shots. The proposed technique enables virtual traveltime corrections, allowing for improved FATT results, even in cases where dense shots are difficult to deploy. Numerical results demonstrate that the proposed method can achieve comparable accuracy to picked traveltime data, indicating its high effectiveness in increasing the resolution of FATT results. The proposed approach maximizes the cost-saving benefits of PRI and can be advantageous in obtaining high-resolution FATT results when dense shot geometry is unavailable.#xD;","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140444053","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
Inversion of Fracture Weakness Parameters Based on the 3CVSP Data?Part I: HTI Media Composed of A Single Fracture Set in the Isotropic Background Media 基于 3CVSP 数据的断裂软弱性参数反演:第一部分:各向同性背景介质中由单一断裂组构成的 HTI 介质
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-19 DOI: 10.1190/geo2023-0539.1
Yuyong Yang, Alexey Stovas, Qiaomu Qi, Huailai Zhou
Natural fractures in oil and gas reservoirs are a crucial factor that cannot be ignored, as they significantly influence the reservoir's petrophysical properties and hydrocarbon development. A horizontal transversely isotropic (HTI) medium composed of a single fracture set in an isotropic background is a typical anisotropic medium. Meanwhile, the shear wave splitting (SWS) is a sensitive response of such anisotropic media, resulting in the generation of fast and slow shear waves. The normal and tangential fracture weaknesses are crucial parameters that characterize the anisotropy of fractured media. We proposed an inversion method for fracture weakness based on three-component vertical seismic profiling (3CVSP) data. Firstly, assuming weak anisotropy and an HTI medium containing single fracture set, we derived a first-order linear approximation of the travel times of the converted fast and slow shear waves (PS1- and PS2-waves) with respect to fracture weakness parameters in the phase velocity domain. By solving for the horizontal projection of the slowness vector, approximate equations of the travel times of the PS1- and PS2-waves were converted from phase velocity domain to the group velocity domain. Furthermore, we devised an inversion workflow consisting of three primary steps: 1. pre-processing the VSP data to derive the travel times and azimuth of the HTI medium; 2. constructing a forward model with undetermined fracture weakness parameters; 3. following the establishment of the objective function, conducting the inversion for the fracture weakness parameters. We demonstrated the reliability of the method through numerical examples and synthetic 3CVSP data. The inversion errors are primarily influenced by the azimuth angle, with minimal influence from the receiver depth. Furthermore, a collective set of inverted results derived from all geophones are more stable and accurate than individual geophones. The application to actual 3CVSP data further confirmed the effectiveness of our approach.
油气藏中的天然裂缝是一个不容忽视的关键因素,因为它们对油气藏的岩石物理特性和油气开发有重大影响。水平横向各向同性(HTI)介质是一种典型的各向异性介质,由各向同性背景中的单一裂缝组成。同时,剪切波分裂(SWS)是此类各向异性介质的敏感响应,会产生快慢剪切波。法向和切向断裂弱度是表征断裂介质各向异性的关键参数。我们提出了一种基于三分量垂直地震剖面(3CVSP)数据的裂缝软弱性反演方法。首先,假设各向异性较弱且 HTI 介质包含单个断裂组,我们推导出转换后的快慢剪切波(PS1 波和 PS2 波)在相速度域中的移动时间与断裂软弱性参数的一阶线性近似值。通过求解慢度矢量的水平投影,PS1 波和 PS2 波的行程时间近似方程从相速度域转换到群速度域。此外,我们还设计了一套反演工作流程,包括三个主要步骤:1. 预处理 VSP 数据,得出 HTI 介质的传播时间和方位角;2.我们通过数值示例和合成的 3CVSP 数据证明了该方法的可靠性。反演误差主要受方位角的影响,而接收器深度的影响很小。此外,与单个检波器相比,由所有检波器得出的集合反演结果更加稳定和准确。对实际 3CVSP 数据的应用进一步证实了我们方法的有效性。
{"title":"Inversion of Fracture Weakness Parameters Based on the 3CVSP Data?Part I: HTI Media Composed of A Single Fracture Set in the Isotropic Background Media","authors":"Yuyong Yang, Alexey Stovas, Qiaomu Qi, Huailai Zhou","doi":"10.1190/geo2023-0539.1","DOIUrl":"https://doi.org/10.1190/geo2023-0539.1","url":null,"abstract":"Natural fractures in oil and gas reservoirs are a crucial factor that cannot be ignored, as they significantly influence the reservoir's petrophysical properties and hydrocarbon development. A horizontal transversely isotropic (HTI) medium composed of a single fracture set in an isotropic background is a typical anisotropic medium. Meanwhile, the shear wave splitting (SWS) is a sensitive response of such anisotropic media, resulting in the generation of fast and slow shear waves. The normal and tangential fracture weaknesses are crucial parameters that characterize the anisotropy of fractured media. We proposed an inversion method for fracture weakness based on three-component vertical seismic profiling (3CVSP) data. Firstly, assuming weak anisotropy and an HTI medium containing single fracture set, we derived a first-order linear approximation of the travel times of the converted fast and slow shear waves (PS1- and PS2-waves) with respect to fracture weakness parameters in the phase velocity domain. By solving for the horizontal projection of the slowness vector, approximate equations of the travel times of the PS1- and PS2-waves were converted from phase velocity domain to the group velocity domain. Furthermore, we devised an inversion workflow consisting of three primary steps: 1. pre-processing the VSP data to derive the travel times and azimuth of the HTI medium; 2. constructing a forward model with undetermined fracture weakness parameters; 3. following the establishment of the objective function, conducting the inversion for the fracture weakness parameters. We demonstrated the reliability of the method through numerical examples and synthetic 3CVSP data. The inversion errors are primarily influenced by the azimuth angle, with minimal influence from the receiver depth. Furthermore, a collective set of inverted results derived from all geophones are more stable and accurate than individual geophones. The application to actual 3CVSP data further confirmed the effectiveness of our approach.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140451076","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
Noise Suppression of ATEM Data using Minimum Curvature Method 使用最小曲率法抑制 ATEM 数据的噪声
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-19 DOI: 10.1190/geo2022-0107.1
Wang Hao, Shengqing Xiong, Li Fei, Wanyin Wang, Feng Bin, Jifeng Zhang, Yang Min
The airborne transient electromagnetic (ATEM) method is a geophysical exploration technique that offers several advantages, including rapid exploration, minimal interference from complex terrains, dense data sampling, and environmental friendliness. However, the full-frequency time-domain secondary field signal acquisition of ATEM is susceptible to various electromagnetic interferences, resulting in noise in the collected data. This noise complicates subsequent processing and interpretation. To address this issue, we have developed a one-dimensional minimum-curvature difference equation with unequal spacing based on the minimum-curvature differential equations. Through research, we propose an ATEM noise-suppression method based on multiple single-step and multiple superposition-step implicit iterative formats. In addition, we have investigated the convergence of explicit and implicit iterative schemes via Fourier spectrum analysis, and proved that the implicit iterative schemes converge. This method effectively suppresses noise in both synthetic and real ATEM data during processing. Our results show that the minimum-curvature method is particularly effective in suppressing sferic and Gaussian noises in ATEM data. As a result, our proposed method provides high-quality and reliable ATEM data for further data processing and interpretation. Moreover, the minimum-curvature method is also capable of suppressing noise in ground and ground–space transient electromagnetic data, as well as magnetotelluric data noise. Therefore, this method exhibits a wide range of potential applications.
机载瞬变电磁(ATEM)方法是一种地球物理勘探技术,具有勘探速度快、受复杂地形干扰小、数据采样密集、环境友好等优点。然而,ATEM 的全频时域二次场信号采集容易受到各种电磁干扰的影响,导致采集的数据出现噪声。这种噪声使后续处理和解释变得复杂。为解决这一问题,我们在最小曲率微分方程的基础上,开发了不等间距的一维最小曲率差分方程。通过研究,我们提出了一种基于多个单步和多个叠加步隐式迭代格式的 ATEM 噪声抑制方法。此外,我们还通过傅立叶频谱分析研究了显式和隐式迭代方案的收敛性,并证明了隐式迭代方案的收敛性。这种方法在处理过程中能有效抑制合成和真实 ATEM 数据中的噪声。结果表明,最小曲率法在抑制 ATEM 数据中的干扰和高斯干扰方面尤为有效。因此,我们提出的方法为进一步的数据处理和解释提供了高质量和可靠的 ATEM 数据。此外,最小曲率法还能抑制地面和地空瞬变电磁数据中的噪声,以及磁法数据噪声。因此,该方法具有广泛的应用潜力。
{"title":"Noise Suppression of ATEM Data using Minimum Curvature Method","authors":"Wang Hao, Shengqing Xiong, Li Fei, Wanyin Wang, Feng Bin, Jifeng Zhang, Yang Min","doi":"10.1190/geo2022-0107.1","DOIUrl":"https://doi.org/10.1190/geo2022-0107.1","url":null,"abstract":"The airborne transient electromagnetic (ATEM) method is a geophysical exploration technique that offers several advantages, including rapid exploration, minimal interference from complex terrains, dense data sampling, and environmental friendliness. However, the full-frequency time-domain secondary field signal acquisition of ATEM is susceptible to various electromagnetic interferences, resulting in noise in the collected data. This noise complicates subsequent processing and interpretation. To address this issue, we have developed a one-dimensional minimum-curvature difference equation with unequal spacing based on the minimum-curvature differential equations. Through research, we propose an ATEM noise-suppression method based on multiple single-step and multiple superposition-step implicit iterative formats. In addition, we have investigated the convergence of explicit and implicit iterative schemes via Fourier spectrum analysis, and proved that the implicit iterative schemes converge. This method effectively suppresses noise in both synthetic and real ATEM data during processing. Our results show that the minimum-curvature method is particularly effective in suppressing sferic and Gaussian noises in ATEM data. As a result, our proposed method provides high-quality and reliable ATEM data for further data processing and interpretation. Moreover, the minimum-curvature method is also capable of suppressing noise in ground and ground–space transient electromagnetic data, as well as magnetotelluric data noise. Therefore, this method exhibits a wide range of potential applications.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139958560","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 acoustic wavefields from moving sources in the presence of a time-varying free surface 模拟存在时变自由表面时来自移动声源的声波场
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-19 DOI: 10.1190/geo2023-0527.1
Khalid Almuteri, Jeffrey Shragge, Paul Sava
Marine vibrators are increasingly being recognized as a viable alternative to seismic airguns for ocean-bottom acquisition due to their ability to generate more low-frequency content and their less adverse impact on marine wildlife. However, their use introduces processing challenges, such as the Doppler effect and time-dependent source-receiver offsets, which are negligible in conventional airgun acquisition. In addition, the time-varying nature of the sea surface during the multi-second acquisition time introduces further challenges for processing and inversion. To accurately account for source motion and time-varying sea surface effects in seismic data processing, we develop a reliable and robust numerical modeling tool. We use a mimetic finite-difference method in a generalized coordinate system to model the full acoustic wavefield triggered by a moving source in the presence of a time-varying sea surface. Our approach uses a time- and space-dependent coordinate transformation, which tracks the source movement and conforms to the irregular time-varying sea surface, to map an irregular physical domain in Cartesian coordinates to a regular computational domain in generalized coordinates. We formulate this coordinate transformation such that both coordinate systems conformally match below the ocean-bottom level. Numerical examples demonstrate that this approach is accurate and stable, even for an unrealistically exaggerated sea state. This computational tool is not limited to modeling, but could also be used to develop advanced processing techniques for marine vibrator data, such as imaging and inversion.
由于海洋振动器能够产生更多的低频内容,而且对海洋野生动物的不利影响较小,因此越来越多的人认为海洋振动器是海底采集地震气枪的可行替代品。然而,使用它们会带来处理方面的挑战,如多普勒效应和随时间变化的震源-接收器偏移,而这在传统气枪采集中可以忽略不计。此外,在多秒采集时间内海面的时变特性也给处理和反演带来了更多挑战。为了在地震数据处理中准确考虑震源运动和时变海面效应,我们开发了一种可靠、稳健的数值建模工具。我们使用广义坐标系中的模拟有限差分法,对存在时变海面的移动震源引发的全声波场进行建模。我们的方法使用与时间和空间相关的坐标变换,该坐标变换可跟踪声源运动并符合不规则的时变海面,从而将笛卡尔坐标中的不规则物理域映射到广义坐标中的规则计算域。我们制定的这种坐标变换使两个坐标系在海底水平面以下保持一致。数值示例表明,这种方法既精确又稳定,即使对于不切实际的夸张海况也是如此。这种计算工具不仅限于建模,还可用于开发海洋振动器数据的高级处理技术,如成像和反演。
{"title":"Modeling acoustic wavefields from moving sources in the presence of a time-varying free surface","authors":"Khalid Almuteri, Jeffrey Shragge, Paul Sava","doi":"10.1190/geo2023-0527.1","DOIUrl":"https://doi.org/10.1190/geo2023-0527.1","url":null,"abstract":"Marine vibrators are increasingly being recognized as a viable alternative to seismic airguns for ocean-bottom acquisition due to their ability to generate more low-frequency content and their less adverse impact on marine wildlife. However, their use introduces processing challenges, such as the Doppler effect and time-dependent source-receiver offsets, which are negligible in conventional airgun acquisition. In addition, the time-varying nature of the sea surface during the multi-second acquisition time introduces further challenges for processing and inversion. To accurately account for source motion and time-varying sea surface effects in seismic data processing, we develop a reliable and robust numerical modeling tool. We use a mimetic finite-difference method in a generalized coordinate system to model the full acoustic wavefield triggered by a moving source in the presence of a time-varying sea surface. Our approach uses a time- and space-dependent coordinate transformation, which tracks the source movement and conforms to the irregular time-varying sea surface, to map an irregular physical domain in Cartesian coordinates to a regular computational domain in generalized coordinates. We formulate this coordinate transformation such that both coordinate systems conformally match below the ocean-bottom level. Numerical examples demonstrate that this approach is accurate and stable, even for an unrealistically exaggerated sea state. This computational tool is not limited to modeling, but could also be used to develop advanced processing techniques for marine vibrator data, such as imaging and inversion.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140449214","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
Efficient Preconditioned Least-Squares Wave-Equation Migration 高效预处理最小二乘波方程迁移
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-19 DOI: 10.1190/geo2023-0048.1
Siamak Abolhassani, D. J. Verschuur
Since the appearance of wave-equation migration, many have tried to improve the resolution and effectiveness of this technology. Least-squares wave-equation migration is one of those attempts that tries to fill the gap between the migration assumptions and reality in an iterative manner. However, these iterations do not come cheap. A proven solution to limit the number of least-squares iterations is to correct the gradient direction within each iteration via the action of a preconditioner that approximates the Hessian inverse. However, the Hessian computation, or even the Hessian approximation computation, in large-scale seismic imaging problems involves an expensive computational bottleneck, making it unfeasible. Therefore, we propose an efficient computation of the Hessian approximation operator, in the context of one-way wave-equation migration (WEM) in the space-frequency domain. We build the Hessian approximation operator depth by depth, considerably reducing the operator size each time it is calculated. We prove the validity of our proposed method with two numerical examples. We then extend our proposal to the framework of full-wavefield migration, which is based on WEM principles but includes interbed multiples. Finally, this efficient preconditioned least-squares full-wavefield migration is successfully applied to a dataset with strong interbed multiple scattering.
自从波方程迁移技术出现以来,许多人都试图提高这项技术的分辨率和有效性。最小二乘法波方程迁移就是其中的一种尝试,它试图以迭代的方式填补迁移假设与现实之间的差距。然而,这些迭代并不便宜。限制最小二乘迭代次数的一个行之有效的解决方案是,在每次迭代过程中,通过近似海斯逆的预处理器的作用来修正梯度方向。然而,在大规模地震成像问题中,黑森计算,甚至黑森近似计算,都涉及昂贵的计算瓶颈,使其不可行。因此,我们在空间-频率域单向波方程迁移(WEM)的背景下,提出了一种高效计算黑森近似算子的方法。我们逐深建立黑森近似算子,大大减少了每次计算时算子的大小。我们用两个数值示例证明了所提方法的有效性。然后,我们将建议扩展到全波场迁移框架,全波场迁移基于 WEM 原理,但包括层间倍数。最后,我们将这种高效的预条件最小二乘全波场迁移方法成功地应用于具有强烈层间多重散射的数据集。
{"title":"Efficient Preconditioned Least-Squares Wave-Equation Migration","authors":"Siamak Abolhassani, D. J. Verschuur","doi":"10.1190/geo2023-0048.1","DOIUrl":"https://doi.org/10.1190/geo2023-0048.1","url":null,"abstract":"Since the appearance of wave-equation migration, many have tried to improve the resolution and effectiveness of this technology. Least-squares wave-equation migration is one of those attempts that tries to fill the gap between the migration assumptions and reality in an iterative manner. However, these iterations do not come cheap. A proven solution to limit the number of least-squares iterations is to correct the gradient direction within each iteration via the action of a preconditioner that approximates the Hessian inverse. However, the Hessian computation, or even the Hessian approximation computation, in large-scale seismic imaging problems involves an expensive computational bottleneck, making it unfeasible. Therefore, we propose an efficient computation of the Hessian approximation operator, in the context of one-way wave-equation migration (WEM) in the space-frequency domain. We build the Hessian approximation operator depth by depth, considerably reducing the operator size each time it is calculated. We prove the validity of our proposed method with two numerical examples. We then extend our proposal to the framework of full-wavefield migration, which is based on WEM principles but includes interbed multiples. Finally, this efficient preconditioned least-squares full-wavefield migration is successfully applied to a dataset with strong interbed multiple scattering.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140451118","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 marine controlled-source electromagnetic application using towed and seafloor-based receivers capable of mapping both seafloor and embedded massive sulfides 海洋受控源电磁应用,使用拖曳式和海底接收器,能够绘制海底和嵌入式块状硫化物地图
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-19 DOI: 10.1190/geo2023-0389.1
K. Ishizu, T. Kasaya, T. Goto, Katsuaki Koike, W. Siripunvaraporn, H. Iwamoto, Y. Kawada, Jun-Ichiro Ishibashi
Deep-sea massive sulfide deposits formed by hydrothermal fluid circulation are potential metal resources. They can exist not only as mound manifestations on the seafloor (seafloor massive sulfides) but also as embedded anomalies buried beneath the seafloor (embedded massive sulfides). The distribution of embedded massive sulfides is largely unknown, despite their expected high economic value. Recent drilling surveys have revealed a complex model suggesting embedded massive sulfides coexist beneath seafloor massive sulfides. In the coexisting case, geophysical methods are required to distinguish and map both seafloor and embedded massive sulfides for accurate resource estimation. Marine controlled-source electromagnetic (CSEM) methods are useful for mapping massive sulfides as they exhibit higher electrical conductivity compared to the surrounding host rock. However, CSEM applications capable of distinguishing and mapping both massive sulfides are lacking. We employ a towed electric dipole transmitter with two types of receivers: stationary ocean bottom electric (OBE) and short-offset towed receivers. This combination utilizes differences in sensitivity: the towed receiver data are sensitive to seafloor massive sulfides and the stationary OBE receiver data are sensitive to embedded massive sulfides. Our synthetic data example demonstrates that the combined inversion of towed and OBE data can recover resistivities and positions of both massive sulfides more accurately than the existing inversion methods using individual applications. We perform the combined inversion of measured CSEM data obtained from the middle Okinawa Trough. The inversion models demonstrate that a combined inversion can map the location and shape of embedded massive sulfides identified during drilling more accurately than the inversion of individual datasets.
由热液循环形成的深海块状硫化物矿床是潜在的金属资源。它们不仅可以以海底土丘的形式存在(海底块状硫化物),也可以以埋藏在海底下的嵌入式异常点的形式存在(嵌入式块状硫化物)。尽管嵌入式块状硫化物具有很高的经济价值,但其分布情况在很大程度上仍不为人所知。最近的钻探勘测揭示了一个复杂的模型,表明嵌入式块状硫化物与海底块状硫化物共存。在共存的情况下,需要采用地球物理方法来区分海底块状硫化物和内蕴块状硫化物并绘制地图,以进行准确的资源估算。海洋可控源电磁(CSEM)方法可用于绘制块状硫化物地图,因为与周围的主岩相比,块状硫化物具有更高的导电性。然而,目前还缺乏能够区分和绘制两种块状硫化物的 CSEM 应用。我们采用了拖曳式电偶极子发射器和两种类型的接收器:固定海底电(OBE)和短偏移拖曳接收器。这种组合利用了灵敏度上的差异:拖曳式接收器数据对海底块状硫化物敏感,而固定式 OBE 接收器数据对嵌入式块状硫化物敏感。我们的合成数据示例表明,与现有的单独应用反演方法相比,拖曳数据和 OBE 数据的组合反演能更准确地恢复两种块状硫化物的电阻率和位置。我们对冲绳海槽中部的 CSEM 测量数据进行了组合反演。反演模型表明,与单个数据集的反演相比,组合反演能更准确地绘制出钻探过程中发现的嵌入式块状硫化物的位置和形状。
{"title":"A marine controlled-source electromagnetic application using towed and seafloor-based receivers capable of mapping both seafloor and embedded massive sulfides","authors":"K. Ishizu, T. Kasaya, T. Goto, Katsuaki Koike, W. Siripunvaraporn, H. Iwamoto, Y. Kawada, Jun-Ichiro Ishibashi","doi":"10.1190/geo2023-0389.1","DOIUrl":"https://doi.org/10.1190/geo2023-0389.1","url":null,"abstract":"Deep-sea massive sulfide deposits formed by hydrothermal fluid circulation are potential metal resources. They can exist not only as mound manifestations on the seafloor (seafloor massive sulfides) but also as embedded anomalies buried beneath the seafloor (embedded massive sulfides). The distribution of embedded massive sulfides is largely unknown, despite their expected high economic value. Recent drilling surveys have revealed a complex model suggesting embedded massive sulfides coexist beneath seafloor massive sulfides. In the coexisting case, geophysical methods are required to distinguish and map both seafloor and embedded massive sulfides for accurate resource estimation. Marine controlled-source electromagnetic (CSEM) methods are useful for mapping massive sulfides as they exhibit higher electrical conductivity compared to the surrounding host rock. However, CSEM applications capable of distinguishing and mapping both massive sulfides are lacking. We employ a towed electric dipole transmitter with two types of receivers: stationary ocean bottom electric (OBE) and short-offset towed receivers. This combination utilizes differences in sensitivity: the towed receiver data are sensitive to seafloor massive sulfides and the stationary OBE receiver data are sensitive to embedded massive sulfides. Our synthetic data example demonstrates that the combined inversion of towed and OBE data can recover resistivities and positions of both massive sulfides more accurately than the existing inversion methods using individual applications. We perform the combined inversion of measured CSEM data obtained from the middle Okinawa Trough. The inversion models demonstrate that a combined inversion can map the location and shape of embedded massive sulfides identified during drilling more accurately than the inversion of individual datasets.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139958775","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
Strain field reconstruction from helical-winding fiber distributed acoustic sensing and its application in anisotropic elastic reverse time migration 螺旋缠绕光纤分布式声学传感的应变场重建及其在各向异性弹性反向时间迁移中的应用
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-02-16 DOI: 10.1190/geo2023-0354.1
Lele Zhang, Yang Zhao, Lu Liu, Fenglin Niu, Wei Wu, Chuangyang Wang, Hengyu Tang, Jingming Li, J. Zuo, Yi Yao, Yixin Wang
Optical fiber-based distributed acoustic sensing (DAS) technology has been a popular seismic acquisition tool due to its easy deployment, wide bandwidth, and dense sampling. However, the sensitivity of straight optical fiber to only single-axis strain presents challenges in fully characterizing multi-components seismic wavefields, making it difficult to use these data in elastic reverse time migration (ERTM). The helical-winding fiber receives projecting signals projected onto the fiber from all seismic strain field components and has the potential to reconstruct those strain components for ERTM imaging. Here we give detailed mathematical principles of helical fiber-based DAS with crucial parameters such as pitch angle, gauge length and rotating angle. At least six points of DAS responses are required in one or several winding periods to rebuild the strain fields within the seismic wavelength. The projecting matrix of conventional regular helical-winding fiber is singular and ill-conditioned, which results in computation challenges for the inverse of the Hessian matrix for strain component reconstruction. To tackle this problem, we develop a non-regular variant-pitch angle winding configuration for helical fiber. Our winding design is validated using the rank and condition number of the projecting matrix, which is proven as an important tool in the reconstruction of the original seismic strains. The recovered strain components from the DAS response are then used to backward propagate receiver wavefields in ERTM with an efficient P/S decoupled approach. To sum up, we develop a novel winding design of helical fiber to recover the strain fields, and then propose an efficient 3D anisotropic P/S wave-mode decomposition method for generating vector P- and S-wavefields during their propagation. Both methods are applied to build an anisotropic DAS-ERTM workflow for producing PP- and PS- images. Two synthetic examples demonstrate the effectiveness of our approach.
基于光纤的分布式声学传感(DAS)技术因其易于部署、带宽宽、采样密集等优点而成为一种流行的地震采集工具。然而,直光纤只对单轴应变敏感,这给全面描述多成分地震波场带来了挑战,使这些数据难以用于弹性反向时间迁移(ERTM)。螺旋缠绕光纤可接收投射到光纤上的所有地震应变场分量的投射信号,并有可能为 ERTM 成像重建这些应变分量。在此,我们给出了基于螺旋光纤的 DAS 的详细数学原理,其中包括俯仰角、轨距长度和旋转角度等关键参数。在一个或多个绕组周期内,至少需要六个点的 DAS 响应,才能重建地震波长内的应变场。传统常规螺旋缠绕光纤的投影矩阵是奇异和条件不良的,这给应变分量重建的赫塞斯矩阵逆计算带来了挑战。为解决这一问题,我们开发了一种非规则变距角螺旋光纤缠绕结构。我们的缠绕设计利用投影矩阵的秩和条件数进行了验证,该矩阵被证明是重建原始地震应变的重要工具。然后,从 DAS 响应中恢复的应变分量可用于在 ERTM 中通过高效的 P/S 解耦方法反向传播接收波场。总之,我们开发了一种新颖的螺旋光纤缠绕设计来恢复应变场,然后提出了一种高效的三维各向异性 P/S 波模分解方法,用于在传播过程中生成矢量 P 波场和 S 波场。这两种方法都被用于建立各向异性 DAS-ERTM 工作流程,以生成 PP 和 PS 图像。两个合成实例证明了我们方法的有效性。
{"title":"Strain field reconstruction from helical-winding fiber distributed acoustic sensing and its application in anisotropic elastic reverse time migration","authors":"Lele Zhang, Yang Zhao, Lu Liu, Fenglin Niu, Wei Wu, Chuangyang Wang, Hengyu Tang, Jingming Li, J. Zuo, Yi Yao, Yixin Wang","doi":"10.1190/geo2023-0354.1","DOIUrl":"https://doi.org/10.1190/geo2023-0354.1","url":null,"abstract":"Optical fiber-based distributed acoustic sensing (DAS) technology has been a popular seismic acquisition tool due to its easy deployment, wide bandwidth, and dense sampling. However, the sensitivity of straight optical fiber to only single-axis strain presents challenges in fully characterizing multi-components seismic wavefields, making it difficult to use these data in elastic reverse time migration (ERTM). The helical-winding fiber receives projecting signals projected onto the fiber from all seismic strain field components and has the potential to reconstruct those strain components for ERTM imaging. Here we give detailed mathematical principles of helical fiber-based DAS with crucial parameters such as pitch angle, gauge length and rotating angle. At least six points of DAS responses are required in one or several winding periods to rebuild the strain fields within the seismic wavelength. The projecting matrix of conventional regular helical-winding fiber is singular and ill-conditioned, which results in computation challenges for the inverse of the Hessian matrix for strain component reconstruction. To tackle this problem, we develop a non-regular variant-pitch angle winding configuration for helical fiber. Our winding design is validated using the rank and condition number of the projecting matrix, which is proven as an important tool in the reconstruction of the original seismic strains. The recovered strain components from the DAS response are then used to backward propagate receiver wavefields in ERTM with an efficient P/S decoupled approach. To sum up, we develop a novel winding design of helical fiber to recover the strain fields, and then propose an efficient 3D anisotropic P/S wave-mode decomposition method for generating vector P- and S-wavefields during their propagation. Both methods are applied to build an anisotropic DAS-ERTM workflow for producing PP- and PS- images. Two synthetic examples demonstrate the effectiveness of our approach.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-02-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139960986","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
An improved 27-point frequency-domain elastic-wave average-derivative method with applications to a moving point source 改进的 27 点频域弹性波平均衍射法及其在移动点源中的应用
IF 3.3 2区 地球科学 Q1 Earth and Planetary Sciences Pub Date : 2024-01-30 DOI: 10.1190/geo2023-0257.1
Hao Wang, Jingbo Chen
Exploration seismology based on moving sources like high-speed trains has attracted increasingly great attention. It is important task to perform near-surface seismic imaging and inversion using moving sources. In this context, an accurate and efficient forward-modeling engine for moving sources is crucial. We initially consider a single moving point source. To implement the numerical simulation in terms of a moving point source, two aspects need to be considered: (1) the numerical discretization strategy of the seismic-wave equation and (2) the discrete representation of the moving point source. For the 3D heterogeneous elastic-wave equation, we propose an improved frequency-domain average-derivative numerical method (ADM) by introducing the weighted average of four sets of grid points for the mass-acceleration term. Due to the numerical discretization template, the continuously moving point source is discretized as a series of fixed sources located at different grid points and excited at different times. The corresponding discrete representation of spatio-temporal variation is given by modifying the right-hand side of the resulting linear system of equations. A numerical experiment in a 3D homogeneous half-space model validates the higher computational accuracy of the improved ADM and the feasibility of the numerical simulation scheme for a moving point source. Furthermore, we test the performance of the moving-point-source numerical simulation scheme in the 3D heterogeneous overthrust model. The successful implementation of the 3D frequency-domain numerical simulation for a moving point source establishes a foundation for subsequent practical applications of moving sources in seismic imaging, to be performed in the future.
基于高速列车等移动震源的勘探地震学越来越受到人们的关注。利用移动震源进行近地表地震成像和反演是一项重要任务。在这种情况下,精确高效的动源前向建模引擎至关重要。我们首先考虑单个移动点源。要实现移动点源的数值模拟,需要考虑两个方面:(1) 地震波方程的数值离散化策略;(2) 移动点源的离散表示。针对三维异质弹性波方程,我们提出了一种改进的频域平均衍生数值方法(ADM),为质量加速项引入了四组网格点的加权平均。由于采用了数值离散模板,连续运动的点源被离散为一系列位于不同网格点、在不同时间被激发的固定源。时空变化的相应离散表示是通过修改所得到的线性方程组的右边给出的。在三维均质半空间模型中进行的数值实验验证了改进 ADM 的更高计算精度,以及针对移动点源的数值模拟方案的可行性。此外,我们还在三维异质推力模型中测试了移动点源数值模拟方案的性能。移动点源三维频域数值模拟的成功实施,为今后移动源在地震成像中的实际应用奠定了基础。
{"title":"An improved 27-point frequency-domain elastic-wave average-derivative method with applications to a moving point source","authors":"Hao Wang, Jingbo Chen","doi":"10.1190/geo2023-0257.1","DOIUrl":"https://doi.org/10.1190/geo2023-0257.1","url":null,"abstract":"Exploration seismology based on moving sources like high-speed trains has attracted increasingly great attention. It is important task to perform near-surface seismic imaging and inversion using moving sources. In this context, an accurate and efficient forward-modeling engine for moving sources is crucial. We initially consider a single moving point source. To implement the numerical simulation in terms of a moving point source, two aspects need to be considered: (1) the numerical discretization strategy of the seismic-wave equation and (2) the discrete representation of the moving point source. For the 3D heterogeneous elastic-wave equation, we propose an improved frequency-domain average-derivative numerical method (ADM) by introducing the weighted average of four sets of grid points for the mass-acceleration term. Due to the numerical discretization template, the continuously moving point source is discretized as a series of fixed sources located at different grid points and excited at different times. The corresponding discrete representation of spatio-temporal variation is given by modifying the right-hand side of the resulting linear system of equations. A numerical experiment in a 3D homogeneous half-space model validates the higher computational accuracy of the improved ADM and the feasibility of the numerical simulation scheme for a moving point source. Furthermore, we test the performance of the moving-point-source numerical simulation scheme in the 3D heterogeneous overthrust model. The successful implementation of the 3D frequency-domain numerical simulation for a moving point source establishes a foundation for subsequent practical applications of moving sources in seismic imaging, to be performed in the future.","PeriodicalId":55102,"journal":{"name":"Geophysics","volume":null,"pages":null},"PeriodicalIF":3.3,"publicationDate":"2024-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140481993","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