{"title":"各向同性弹性迭代逆时偏移的无矩阵多参数下降和共轭梯度法","authors":"W. Mulder","doi":"10.3997/2214-4609.202112725","DOIUrl":null,"url":null,"abstract":"Multi-parameter inversion of linear systems appears in many problems. The focus here is on isotropic elastic iterative reverse-time migration for three position-dependent subsurface model parameters, which amounts to data fitting of processed seismic data with synthetics from the Born approximation of the elastic wave equation. In that case, the matrix of the linear system is the hessian. As it is impractical to form, a matrix-free formulation is needed, which is readily derived for the gradient descent method. For single-parameter inversion, the conjugate-gradient (CG) method is generally more efficient than simple descent. However, the multiple-parameter CG method has a significantly higher cost than the descent method. Here, first a matrix-free data-domain reformulation is derived. Then, its performance is compared to the simple descent method to see of its faster convergence justifies the higher cost. A comparison on a marine 2-D toy problem with a salt body and sea-bottom receivers shows that the multiple-parameter descent method wins in terms of efficiency if the number of iterations is limited and that the single-parameter CG method is even faster.","PeriodicalId":143998,"journal":{"name":"82nd EAGE Annual Conference & Exhibition","volume":"206 1","pages":"0"},"PeriodicalIF":0.0000,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"A matrix-free reformulation of the multi-parameter descent and conjugate-gradient method for isotropic elastic iterative reverse-time migration\",\"authors\":\"W. Mulder\",\"doi\":\"10.3997/2214-4609.202112725\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Multi-parameter inversion of linear systems appears in many problems. The focus here is on isotropic elastic iterative reverse-time migration for three position-dependent subsurface model parameters, which amounts to data fitting of processed seismic data with synthetics from the Born approximation of the elastic wave equation. In that case, the matrix of the linear system is the hessian. As it is impractical to form, a matrix-free formulation is needed, which is readily derived for the gradient descent method. For single-parameter inversion, the conjugate-gradient (CG) method is generally more efficient than simple descent. However, the multiple-parameter CG method has a significantly higher cost than the descent method. Here, first a matrix-free data-domain reformulation is derived. Then, its performance is compared to the simple descent method to see of its faster convergence justifies the higher cost. A comparison on a marine 2-D toy problem with a salt body and sea-bottom receivers shows that the multiple-parameter descent method wins in terms of efficiency if the number of iterations is limited and that the single-parameter CG method is even faster.\",\"PeriodicalId\":143998,\"journal\":{\"name\":\"82nd EAGE Annual Conference & Exhibition\",\"volume\":\"206 1\",\"pages\":\"0\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"1900-01-01\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"82nd EAGE Annual Conference & Exhibition\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/10.3997/2214-4609.202112725\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"82nd EAGE Annual Conference & Exhibition","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.3997/2214-4609.202112725","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
A matrix-free reformulation of the multi-parameter descent and conjugate-gradient method for isotropic elastic iterative reverse-time migration
Multi-parameter inversion of linear systems appears in many problems. The focus here is on isotropic elastic iterative reverse-time migration for three position-dependent subsurface model parameters, which amounts to data fitting of processed seismic data with synthetics from the Born approximation of the elastic wave equation. In that case, the matrix of the linear system is the hessian. As it is impractical to form, a matrix-free formulation is needed, which is readily derived for the gradient descent method. For single-parameter inversion, the conjugate-gradient (CG) method is generally more efficient than simple descent. However, the multiple-parameter CG method has a significantly higher cost than the descent method. Here, first a matrix-free data-domain reformulation is derived. Then, its performance is compared to the simple descent method to see of its faster convergence justifies the higher cost. A comparison on a marine 2-D toy problem with a salt body and sea-bottom receivers shows that the multiple-parameter descent method wins in terms of efficiency if the number of iterations is limited and that the single-parameter CG method is even faster.