Yinfeng Wang, Jilong Lu, Ying Shi, Ning Wang, Liguo Han
{"title":"利用分数粘弹性波方程高精度模拟瑞利波","authors":"Yinfeng Wang, Jilong Lu, Ying Shi, Ning Wang, Liguo Han","doi":"10.3390/fractalfract7120880","DOIUrl":null,"url":null,"abstract":"The propagation of Rayleigh waves is usually accompanied by dispersion, which becomes more complex with inherent attenuation. The accurate simulation of Rayleigh waves in attenuation media is crucial for understanding wave mechanisms, layer thickness identification, and parameter inversion. Although the vacuum formalism or stress image method (SIM) combined with the generalized standard linear solid (GSLS) is widely used to implement the numerical simulation of Rayleigh waves in attenuation media, this type of method still has its limitations. First, the GSLS model cannot split the velocity dispersion and amplitude attenuation term, thus limiting its application in the Q-compensated reverse time migration/full waveform inversion. In addition, GSLS-model-based wave equation is usually numerically solved using staggered-grid finite-difference (SGFD) method, which may result in the numerical dispersion due to the harsh stability condition and poses complexity and computational burden. To overcome these issues, we propose a high-accuracy Rayleigh-waves simulation scheme that involves the integration of the fractional viscoelastic wave equation and vacuum formalism. The proposed scheme not only decouples the amplitude attenuation and velocity dispersion but also significantly suppresses the numerical dispersion of Rayleigh waves under the same grid sizes. We first use a homogeneous elastic model to demonstrate the accuracy in comparison with the analytical solutions, and the correctness for a viscoelastic half-space model is verified by comparing the phase velocities with the dispersive images generated by the phase shift transformation. We then simulate several two-dimensional synthetic models to analyze the effectiveness and applicability of the proposed method. The results show that the proposed method uses twice as many spatial step sizes and takes 0.6 times that of the GSLS method (solved by the SGFD method) when achieved at 95% accuracy.","PeriodicalId":12435,"journal":{"name":"Fractal and Fractional","volume":"26 7","pages":""},"PeriodicalIF":3.6000,"publicationDate":"2023-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"High-Accuracy Simulation of Rayleigh Waves Using Fractional Viscoelastic Wave Equation\",\"authors\":\"Yinfeng Wang, Jilong Lu, Ying Shi, Ning Wang, Liguo Han\",\"doi\":\"10.3390/fractalfract7120880\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"The propagation of Rayleigh waves is usually accompanied by dispersion, which becomes more complex with inherent attenuation. The accurate simulation of Rayleigh waves in attenuation media is crucial for understanding wave mechanisms, layer thickness identification, and parameter inversion. Although the vacuum formalism or stress image method (SIM) combined with the generalized standard linear solid (GSLS) is widely used to implement the numerical simulation of Rayleigh waves in attenuation media, this type of method still has its limitations. First, the GSLS model cannot split the velocity dispersion and amplitude attenuation term, thus limiting its application in the Q-compensated reverse time migration/full waveform inversion. In addition, GSLS-model-based wave equation is usually numerically solved using staggered-grid finite-difference (SGFD) method, which may result in the numerical dispersion due to the harsh stability condition and poses complexity and computational burden. To overcome these issues, we propose a high-accuracy Rayleigh-waves simulation scheme that involves the integration of the fractional viscoelastic wave equation and vacuum formalism. The proposed scheme not only decouples the amplitude attenuation and velocity dispersion but also significantly suppresses the numerical dispersion of Rayleigh waves under the same grid sizes. We first use a homogeneous elastic model to demonstrate the accuracy in comparison with the analytical solutions, and the correctness for a viscoelastic half-space model is verified by comparing the phase velocities with the dispersive images generated by the phase shift transformation. We then simulate several two-dimensional synthetic models to analyze the effectiveness and applicability of the proposed method. The results show that the proposed method uses twice as many spatial step sizes and takes 0.6 times that of the GSLS method (solved by the SGFD method) when achieved at 95% accuracy.\",\"PeriodicalId\":12435,\"journal\":{\"name\":\"Fractal and Fractional\",\"volume\":\"26 7\",\"pages\":\"\"},\"PeriodicalIF\":3.6000,\"publicationDate\":\"2023-12-12\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Fractal and Fractional\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://doi.org/10.3390/fractalfract7120880\",\"RegionNum\":2,\"RegionCategory\":\"数学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"MATHEMATICS, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Fractal and Fractional","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.3390/fractalfract7120880","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
High-Accuracy Simulation of Rayleigh Waves Using Fractional Viscoelastic Wave Equation
The propagation of Rayleigh waves is usually accompanied by dispersion, which becomes more complex with inherent attenuation. The accurate simulation of Rayleigh waves in attenuation media is crucial for understanding wave mechanisms, layer thickness identification, and parameter inversion. Although the vacuum formalism or stress image method (SIM) combined with the generalized standard linear solid (GSLS) is widely used to implement the numerical simulation of Rayleigh waves in attenuation media, this type of method still has its limitations. First, the GSLS model cannot split the velocity dispersion and amplitude attenuation term, thus limiting its application in the Q-compensated reverse time migration/full waveform inversion. In addition, GSLS-model-based wave equation is usually numerically solved using staggered-grid finite-difference (SGFD) method, which may result in the numerical dispersion due to the harsh stability condition and poses complexity and computational burden. To overcome these issues, we propose a high-accuracy Rayleigh-waves simulation scheme that involves the integration of the fractional viscoelastic wave equation and vacuum formalism. The proposed scheme not only decouples the amplitude attenuation and velocity dispersion but also significantly suppresses the numerical dispersion of Rayleigh waves under the same grid sizes. We first use a homogeneous elastic model to demonstrate the accuracy in comparison with the analytical solutions, and the correctness for a viscoelastic half-space model is verified by comparing the phase velocities with the dispersive images generated by the phase shift transformation. We then simulate several two-dimensional synthetic models to analyze the effectiveness and applicability of the proposed method. The results show that the proposed method uses twice as many spatial step sizes and takes 0.6 times that of the GSLS method (solved by the SGFD method) when achieved at 95% accuracy.
期刊介绍:
Fractal and Fractional is an international, scientific, peer-reviewed, open access journal that focuses on the study of fractals and fractional calculus, as well as their applications across various fields of science and engineering. It is published monthly online by MDPI and offers a cutting-edge platform for research papers, reviews, and short notes in this specialized area. The journal, identified by ISSN 2504-3110, encourages scientists to submit their experimental and theoretical findings in great detail, with no limits on the length of manuscripts to ensure reproducibility. A key objective is to facilitate the publication of detailed research, including experimental procedures and calculations. "Fractal and Fractional" also stands out for its unique offerings: it warmly welcomes manuscripts related to research proposals and innovative ideas, and allows for the deposition of electronic files containing detailed calculations and experimental protocols as supplementary material.