News and Announcements
ARTICLE
Trigonal meshes in diffraction tomography with optimum regularization: an application forcarbon sequestration monitoring
E.T.F. SANTOS1 ,
J.M. HARRIS2 ,
A. BASSREI3 ,
J.C. COSTA4
Show Less
1
Federal Center of Technological Education of Bahia (CEFET-BA), Salvador, BA, Brazil. eduardot@cefetba.br,
2
Department of Geophysics, Stanford University, Stanford, CA, U.S.A.,
3
Institute of Physics & Research Center in Geophysics and Geology, Federal University of Bahia, Caixa Postal 1001, 40001-970 Salvador, BA, Brazil.,
4
Faculty of Geophysics, Federal University of Pará, Belém, PA, Brazil.,
JSE 2009, 18(2), 135–156;
Submitted: 9 June 2025 | Revised: 9 June 2025 | Accepted: 9 June 2025 | Published: 9 June 2025
© 2025 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract
Santos, E.T.F., Harris, J.M., Bassrei, A. and Costa, J.C., 2009. Trigonal meshes in diffraction tomography with optimum regularization: an application for carbon sequestration monitoring. Journal of Seismic Exploration, 18: 135-156. Diffraction tomography is an inversion technique that provides the reconstruction of a subsurface velocity field from scattered acoustic field data. High-resolution imaging conventionally requires estimation of a large number of parameters. A trigonal mesh is applied in the study described in this paper, in order to strongly reduce the number of parameters. Thus, instead of a velocity-estimate for each cell in a regular grid, the velocity is estimated only at triangle vertices, which act as control points for the interpolation of velocity field within each triangle. Regularization is required to avoid sharp artifacts due to trigonal elements. A synthetic model is adopted to test the feasibility of the proposed method for reservoir monitoring.
Keywords
inverse problems
diffraction tomography
regularization
trigonal meshes
CO2 sequestration
References
- Ajo-Franklin, J.B., Urban, J.A. and Harris, J.M., 2006. Using resolution-constrained adaptive
- meshes for traveltime tomography. J. Seismic Explor., 14: 371-392.
- Bassrei, A. and Rodi, W.L., 1993. Regularization and inversion of linear geophysical data. 3rd
- Internat. Congr. Brazil. Geophys. Soc., Rio de Janeiro, Brazil, I: 111-116.
- Belge, M., Kilmer, M.E. and Miller, E.L., 2002. Efficient determination of multiple regularization
- parameters in a generalized L-curve framework. Inverse Problems, 18: 1161-1183.
- Bottema, O., 1982. On the area of a triangle in barycentric coordinates. Crux Mathematic., 8:
- 228-231.
- 156 SANTOS, HARRIS, BASSREI & COSTA
- Devaney, A.J., 1984. Geophysical diffraction tomography.Inst. Electr. Electron. Engin. Transact.
- Geosci. Remote Sens., 22: 3-13.
- Hansen, P.C., 1992. Analysis of discrete ill-posed problems by means of the I-curve. Soc. Industr.
- Appl. Mathemat. Rev., 34: 561-580.
- Hansen, P.C., 1998. Rank-Deficient and Discrete Ill-Posed Problems. Soc. Industr. Appl.
- Mathemat, Philadelphia.
- Hansen, P.C. and O’Leary, D.P., 1993. The use of the L-curve in the regularization of discrete
- ill-posed problems. Soc. Industr. Appl. Mathemat. J. Scientif. Computat., 14: 1487-1503.
- Harris, J.M., 1987. Diffraction tomography with discrete arrays of sources and receivers. Inst.
- Electr. Electron. Engin. Transact. Geosci. Remote Sens., 25: 448-455.
- Kilmer, M.E., O’Leary, D.P., 2001. Choosing Regularization Parameters in Iterative methods for
- ill-posed problems. Soc. Industr. Appl. Mathemat. J. Matrix Analys. Applic., 22:
- 1204-1221.
- Lanczos, C., 1961. Linear Differential Operators. Van Nostrand, London.
- Levenberg, K., 1944. A method for the solution of certain non-linear problems in least-squares.
- Quart. Appl. Mathemat., 2: 164-168.
- Lo, T.-W. and Inderwiesen, P.L., 1994. Fundamentals of Seismic Tomography. Geophysical
- Monograph Series, SEG, Tulsa, OK.
- Marquardt, D.W., 1963. An algorithm for least-squares of nonlinear parameters. J. Soc. Industr.
- Appl. Mathemat., 11: 431-441.
- Penrose, R., 1955. A generalized inverse for matrices. Proc. Cambridge Philosoph. Soc., 51:
- 406-413.
- Persson, P.-O. and Strang, G., 2004. A simple mesh generator in MATLAB. Soc. Industr. Appl.
- Mathemat. Rev., 46: 329-345.
- Reiter, T.D. and Rodi, W., 1996. Nonlinear waveform tomography applied to crosshole seismic
- data. Geophysics, 61: 902-913.
- Rocha Filho, A.A., Harris, J.M. and Bassrei, A., 1996. A simple matrix formulation diffraction
- tomography algorithm. 39th Congr. Brasil. Geol., Salvador, BA, Brasil, 2: 312-315.
- Santos, E.T.F., Bassrei, A. and Costa, J., 2006. Evaluation of L-curve and @-curve approaches for
- the selection of regularization parameter in anisotropic traveltime tomography. J. Seismic
- Explor., 15: 245-272.
- Santos, E.T.F. and Bassrei, A., 2007. L- and @-curve approaches for the selection of regularization
- parameter in geophysical diffraction tomography. Comput. Geosci., 33: 618-629.
- Thompson, D.R., Rodi, W. and Tokséz, M.N., 1994. Nonlinear seismic diffraction tomography
- using minimum structure constraints. J. Acoust. Soc. Am., 95: 324-330.
- Titterington, D.M., 1985. General structure of regularization procedures in image reconstruction.
- Astron. Astrophys., 144: 381-387.
- Twomey, S., 1963. On the numerical solution of Fredholm integral equations of the first kind by
- the inversion of the linear system produced by quadrature. J. Assoc. Comput. Mach., 10:
- 97-101.
- Wu, R.-S. and Toks6z, M.N., 1987. Diffraction tomography and multisource holography applied
- to seismic imaging. Geophysics, 52: 11-25.
