{"title":"利用快速扫描法计算二维一般倾斜横向各向同性介质中的行进时间","authors":"Yongming Lu , Wei Zhang , Jianfeng Zhang","doi":"10.1016/j.cageo.2024.105724","DOIUrl":null,"url":null,"abstract":"<div><p>Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.</p></div>","PeriodicalId":55221,"journal":{"name":"Computers & Geosciences","volume":"193 ","pages":"Article 105724"},"PeriodicalIF":4.2000,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Calculating traveltimes in 2D general tilted transversely isotropic media using fast sweeping method\",\"authors\":\"Yongming Lu , Wei Zhang , Jianfeng Zhang\",\"doi\":\"10.1016/j.cageo.2024.105724\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><p>Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.</p></div>\",\"PeriodicalId\":55221,\"journal\":{\"name\":\"Computers & Geosciences\",\"volume\":\"193 \",\"pages\":\"Article 105724\"},\"PeriodicalIF\":4.2000,\"publicationDate\":\"2024-09-13\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Computers & Geosciences\",\"FirstCategoryId\":\"89\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0098300424002073\",\"RegionNum\":2,\"RegionCategory\":\"地球科学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers & Geosciences","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0098300424002073","RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
Calculating traveltimes in 2D general tilted transversely isotropic media using fast sweeping method
Traveltime calculations play an important role in the field of exploration seismology, such as traveltime tomography and seismic imaging and so on. Seismic anisotropy poses a challenge for traveltime calculation, because anisotropic eikonal solvers are more complex than the isotropic counter part. To solve the eikonal equations in 2D tilted transversely isotropic (TTI) media, we have developed a fast algorithm combine with fast sweeping method to compute the first arrival traveltimes of quasi-P (qP)-, quasi-SV (qSV)-, and quasi-SH(qSH)-waves. For the qP- and qSV-waves, we analyzed the quartic coupled slowness surface equation derived from the Christoffel equation. Then, we constructed a local solver to relate traveltime and slowness. We found that in the local solver, one component of the slowness vector is known and the corresponding slowness equation is monotonic. This provides a strong basis for the fast iterative algorithm we proposed, where we use the Newton method to solve the qP- and qSV-wave slowness equation to determine the related traveltimes. For the qSH wave, the slowness equation is quadratic and simple to solve. Numerical experiments demonstrate that the proposed method can obtain accurate traveltimes for simple and complicated 2D TTI models.
期刊介绍:
Computers & Geosciences publishes high impact, original research at the interface between Computer Sciences and Geosciences. Publications should apply modern computer science paradigms, whether computational or informatics-based, to address problems in the geosciences.