{"title":"基于高斯过程回归的扩散-粘性波动方程参数反演","authors":"Zhaowei Bai, Haixia Zhao, Shaoru Wang","doi":"10.1093/jge/gxad085","DOIUrl":null,"url":null,"abstract":"Abstract The diffusive-viscous wave (DVW) equation is used to characterize the relationship between frequency-dependent seismic responses and saturated fluids by incorporating the frictional dissipation and viscous damping to the scalar wave equation. Simultaneous inversion of three model parameters in DVW equation is essential for seismic interpretations. Traditional inversion methods require continuous forward modeling updates, resulting in low computational efficiency. Moreover, the traditional methods have limitations in simultaneously inverting multi-parameters of wave equations such as DVW equation, usually fixing one parameter to invert the other two parameters. Gaussian process regression (GPR) is a kernel-based non-parametric probabilistic model that introduces prior variables through Gaussian processes (GP). We present a method for the inversion of the three parameters (velocity, diffusive and viscous attenuation coefficients) of the DVW equation based on GPR. The procedure consists of initially implementing the central finite difference approximation to discretize the DVW equation in the time domain. Subsequently, a Gaussian prior is provided on two snapshots of the DVW equation to obtain the corresponding kernel functions. Furthermore, the hyperparameters in kernel functions and the three model parameters are simultaneously trained by minimizing the negative logarithmic marginal likelihood with few training samples while incorporating the underlying physics in terms of encoding the DVW equation into the kernel functions. It is worth noting that it is the first time to implement three-parameter simultaneous inversion based on DVW equation. The numerical examples in homogeneous, layered and heterogeneous media demonstrate the effectiveness of this method.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":"47 1","pages":"0"},"PeriodicalIF":1.6000,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Parameter inversion of the diffusive-viscous wave equation based on Gaussian process regression\",\"authors\":\"Zhaowei Bai, Haixia Zhao, Shaoru Wang\",\"doi\":\"10.1093/jge/gxad085\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Abstract The diffusive-viscous wave (DVW) equation is used to characterize the relationship between frequency-dependent seismic responses and saturated fluids by incorporating the frictional dissipation and viscous damping to the scalar wave equation. Simultaneous inversion of three model parameters in DVW equation is essential for seismic interpretations. Traditional inversion methods require continuous forward modeling updates, resulting in low computational efficiency. Moreover, the traditional methods have limitations in simultaneously inverting multi-parameters of wave equations such as DVW equation, usually fixing one parameter to invert the other two parameters. Gaussian process regression (GPR) is a kernel-based non-parametric probabilistic model that introduces prior variables through Gaussian processes (GP). We present a method for the inversion of the three parameters (velocity, diffusive and viscous attenuation coefficients) of the DVW equation based on GPR. The procedure consists of initially implementing the central finite difference approximation to discretize the DVW equation in the time domain. Subsequently, a Gaussian prior is provided on two snapshots of the DVW equation to obtain the corresponding kernel functions. Furthermore, the hyperparameters in kernel functions and the three model parameters are simultaneously trained by minimizing the negative logarithmic marginal likelihood with few training samples while incorporating the underlying physics in terms of encoding the DVW equation into the kernel functions. It is worth noting that it is the first time to implement three-parameter simultaneous inversion based on DVW equation. The numerical examples in homogeneous, layered and heterogeneous media demonstrate the effectiveness of this method.\",\"PeriodicalId\":54820,\"journal\":{\"name\":\"Journal of Geophysics and Engineering\",\"volume\":\"47 1\",\"pages\":\"0\"},\"PeriodicalIF\":1.6000,\"publicationDate\":\"2023-10-19\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Geophysics and Engineering\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/10.1093/jge/gxad085\",\"RegionNum\":3,\"RegionCategory\":\"地球科学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q3\",\"JCRName\":\"GEOCHEMISTRY & GEOPHYSICS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Geophysics and Engineering","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1093/jge/gxad085","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
Parameter inversion of the diffusive-viscous wave equation based on Gaussian process regression
Abstract The diffusive-viscous wave (DVW) equation is used to characterize the relationship between frequency-dependent seismic responses and saturated fluids by incorporating the frictional dissipation and viscous damping to the scalar wave equation. Simultaneous inversion of three model parameters in DVW equation is essential for seismic interpretations. Traditional inversion methods require continuous forward modeling updates, resulting in low computational efficiency. Moreover, the traditional methods have limitations in simultaneously inverting multi-parameters of wave equations such as DVW equation, usually fixing one parameter to invert the other two parameters. Gaussian process regression (GPR) is a kernel-based non-parametric probabilistic model that introduces prior variables through Gaussian processes (GP). We present a method for the inversion of the three parameters (velocity, diffusive and viscous attenuation coefficients) of the DVW equation based on GPR. The procedure consists of initially implementing the central finite difference approximation to discretize the DVW equation in the time domain. Subsequently, a Gaussian prior is provided on two snapshots of the DVW equation to obtain the corresponding kernel functions. Furthermore, the hyperparameters in kernel functions and the three model parameters are simultaneously trained by minimizing the negative logarithmic marginal likelihood with few training samples while incorporating the underlying physics in terms of encoding the DVW equation into the kernel functions. It is worth noting that it is the first time to implement three-parameter simultaneous inversion based on DVW equation. The numerical examples in homogeneous, layered and heterogeneous media demonstrate the effectiveness of this method.
期刊介绍:
Journal of Geophysics and Engineering aims to promote research and developments in geophysics and related areas of engineering. It has a predominantly applied science and engineering focus, but solicits and accepts high-quality contributions in all earth-physics disciplines, including geodynamics, natural and controlled-source seismology, oil, gas and mineral exploration, petrophysics and reservoir geophysics. The journal covers those aspects of engineering that are closely related to geophysics, or on the targets and problems that geophysics addresses. Typically, this is engineering focused on the subsurface, particularly petroleum engineering, rock mechanics, geophysical software engineering, drilling technology, remote sensing, instrumentation and sensor design.