Ernesto Bonomi - Profile on Academia.edu (original) (raw)
Papers by Ernesto Bonomi
From the 3D phase-shift extrapolation of individual phase-space components, we derive for prestac... more From the 3D phase-shift extrapolation of individual phase-space components, we derive for prestacked data a vector relation between kh, the horizontal offset wavenumber, and kz, the vertical one. While in 2D this relation depends only on the tangent of the scattering angle θ and not on the structural deep, in 3D, in general, this is no longer true. The resulting vector formula takes into account the orientation of the scattering plane containing the slowness vectors ps and pr, one describing the down-going wave and the other the up-going one.
Simulation numérique et mécanique statistique : extension et étude de quelques problèmes d'ingénierie
Le contenu de ce travail, rassemble en une serie d’articles –publies, ou en cours de publication-... more Le contenu de ce travail, rassemble en une serie d’articles –publies, ou en cours de publication- touche aussi bien a la physique, aux telecommunications qu’a l’optimisation combinatoire. Les idees utilisees s’articulent autour de deux disciplines : la simulation numerique et la mecanique statistique des systemes, aussi bien physiques que d’ingenierie. Nous entendons par systeme, un ensemble d’elements en interaction entre eux et avec l’environnement qui les entoure. Le point de vue adopte est celui du modele om le probleme est ramene a un niveau tel, qu’une analyse mathematique puisse etre envisagee, ou tout au moins formulee. Il s’est agi pour nous d’analyser l’ordre statistique qui prend place dans un espace des configurations dans lequel les interactions font evoluer les differents elements. Le cadre de depart a ete microscopique et l’utilisation d’une description complete, imposee jusque dans les plus petits details du systeme, nous a amenes tout naturellement a definir des pro...
Computers in Physics, 1998
The Datacrossing DSS: A Data-GRID Based Decision Support System for Groundwater Management
2006 Second IEEE International Conference on e-Science and Grid Computing (e-Science'06), 2006
... the technologies we adopted to design distributed GIS applications, providing original soluti... more ... the technologies we adopted to design distributed GIS applications, providing original solutions for an ... have played a key role in the development of an effective, authoritative, interoperable ... groundwater model to support decision-making in the field of sustainable water resources ...
This document presents the general problem of electromagnetic inversion in the frequency domain. ... more This document presents the general problem of electromagnetic inversion in the frequency domain. The 1-dimensional analysis of a stratified medium is performed at different frequencies; the data measured on the surface (the admittance or the apparent electric ...
PSPI for High Performance Echo-Reconstructive Imaging
Ernesto Bonomi (ernesto@ crs4. it) obtained his Ph. D. in Theoretical Physics at the University o... more Ernesto Bonomi (ernesto@ crs4. it) obtained his Ph. D. in Theoretical Physics at the University of Paris XI, Orsay and is currently the head of geophysics at CRS4. His earlier professional activities included work in statistical mechanics and computational physics; ...
Methods for Optimization in Presence of Relative Minima: A Conjugate-Direction Method
Seismic imaging of medical ultrasound data: Towards in vivo applications
Europhysics Letters
This study presents the application to medical ultrasonography of a high value-added imaging meth... more This study presents the application to medical ultrasonography of a high value-added imaging method developed in the context of geophysical exploration and based on the undulatory description of the physical process. Currently, the most serious limitation of conventional medical ultrasound systems is the implicit assumption of constant-velocity media stemming from the widespread use of geometrical ray-based imaging algorithms and leading to phase aberration phenomena in the case typical of the human body of two or more tissues with different velocities. As a result, imaging of targets underlying tissue layers can be severely degraded. To address this problem, the proposed algorithm is designed expressly for variable-velocity media, providing high-resolution images by implementing the solution of the one-way wave equation. From a macro-velocity model, the illuminated structures are imaged by correlating, in the frequency domain and for each source element, the simulated emitted wavef...
The SmartGeo Portal: A retrospective
Although time-domain depth migration techniques have been successfully ported to run on modern ha... more Although time-domain depth migration techniques have been successfully ported to run on modern hardware accelerators, their ultimate obstacle is the I/O overhead present during the imaging step. Frequency-domain depth migration algorithms overcome this limitation and can exploit the full potential of new computing technologies. In particular, our implementation of Phase Shift Plus Interpolation (PSPI) method is characterized by fast running time, good quality results under lowsignal-to-noise ratio conditions and excellent results for steep dips. We provide a novel computational dataflow scheme to perform acceleration of PSPI on a generic dataflow engine. We present speedup results obtained on the state-of-the-art dataflow technology for synthetic VTI datasets. Our measurements indicate that a dataflow approach can achieve high speedups despite larger and larger computational domains, increased complexity of the anisotropic approach and the I/O overhead during angle-gathers calculation.
Physical Review Applied, 2020
Conventional medical ultrasound systems normally implement ray-based imaging algorithms, such as ... more Conventional medical ultrasound systems normally implement ray-based imaging algorithms, such as delay-and-sum beamforming, whose severest limitation derives from the implicit assumption of constantvelocity media. As a result, in the case of two or more tissues with different velocities, the image of the underlying targets appears strongly degraded both in placement and in resolution. The proposed ultrasound-imaging strategy, a value-added application of concepts developed in the context of seismic prospecting, avoids this restriction by relying on the undulatory description of the physical process and not on the geometric one. Echoes are sensed in the synthetic aperture configuration by the transducer, whose elements sequentially emit a nearly spherical wave front that covers the whole region of interest. Given a macrovelocity model, the recorded echoes become the initial condition for the downward propagator of the ultrasound wavefield. The processing is in three steps: time-reverse propagation of the sensed echoes, forward simulation of the wavefield emitted by the source element, and partial imaging by computing the zero-lag temporal correlation of the two wavefields to detect the scattering structures. The final reconstruction is obtained by stacking all the partial results on a single image. Laboratory tests, performed on experimental data acquired on a physical phantom, with and without an aberrant layer, prove the effectiveness of the proposed method even in the case of vertical and lateral velocity variations, with images of impressive spatial resolution and highly accurate target placement.
Amplitude-Preserving Monte Carlo 3D Prestack Migration
64th EAGE Conference & Exhibition, 2002
Cloud computing is establishing worldwide as a new high performance computing paradigm that offer... more Cloud computing is establishing worldwide as a new high performance computing paradigm that offers formidable possibilities to industry and science. The presented cloud-computing portal, part of the Grida3 project, provides an innovative approach to seismic data processing by combining open-source state-of-the-art processing software and cloud-computing technology, making possible the effective use of distributed computation and data management with administratively distant resources. We substituted the user-side demanding hardware and software requirements by remote access to high-performance grid-computing facilities. As a result, data processing can be done quasi in real-time being ubiquitously controlled via Internet by a user-friendly web-browser interface. Besides the obvious advantages over locally installed seismic-processing packages, the presented cloud-computing solution creates completely new possibilities for scientific education, collaboration, and presentation of repr...
3D Zero-Offset CRS Stack for Narrow-Azimuth Data - Formulation and Examples
EAGE/SEG Summer Research Workshop - Processing and imaging of seismic data, 2003
The common-reflection-surface (CRS) stack is a macro-model independent seismic reflection imaging... more The common-reflection-surface (CRS) stack is a macro-model independent seismic reflection imaging method that has proved to accurately detect and characterise meaningful events in the pre-stack data. As any 3D data-driven method, it requires an azimuthal distribution of traces as uniform as possible, in order to properly infer necessary information from the data. Its application to 3D datasets with narrow-azimuth acquisitions requires some simplified assumptions of the underlying theory. We present an approximated formulation and a related processing method that still represents an improvement with respect to conventional techniques currently employed. The advantages of the method when performing stacking velocity analysis are demonstrated on a synthetic example
Earth Science Informatics, 2009
Shared Resources Manager for Environmental Data Analysis and Applications) is a multidisciplinary... more Shared Resources Manager for Environmental Data Analysis and Applications) is a multidisciplinary project designed to deliver an integrated system to forge solutions to some environmental challenges such as the constant increase of polluted sites, the sustainability of natural resources usage and the forecast of extreme meteorological events. The GRIDA3 portal is mainly based on Web 2.0 technologies and EnginFrame framework. The portal, now at an advanced stage of development, provides end-users with intuitive Web-interfaces and tools that simplify job submission to the underneath computing resources. The framework manages the user authentication and authorization, then controls the action and job execution into the grid computing environment, collects the results and transforms them into an useful format on the client side. The GRIDA3 Portal framework will provide a problemsolving platform allowing, through appropriate access policies, the integration and the sharing of skills, resources and tools located at multiple sites across federated domains.
Journal of Geophysics and Engineering, 2017
The present article outlines an innovative derivation of the well-known CRS traveltime formula. T... more The present article outlines an innovative derivation of the well-known CRS traveltime formula. This is made possible by approximating the isochrones tangent to a small reflecting element embedded in a 3D homogeneous auxiliary medium. The arising paraxial traveltime formula is then parametrised by six attributes, each one characterising geometrically the wavefront of the reference normal-incident ray at the emergence point x 0 of the datum plane. For a layered medium, the assignment of the six attributes for each azimuthal direction around x 0 describes locally the emerging wavefront and establishes the mapping between reflecting elements of the two media, the auxiliary and the real, which respond, however, with different traveltimes. This 3D homeomorphism provides the time correction that gives rise to the CRS traveltime formula with eight attributes. For time imaging applications, the traveltime profile must be numerically shaped by improving iteratively the value of the eight attributes, so as to intercept, without the need of a velocity model, the largest number of coherent data in the volume of seismic traces gathered in the midpoint-offset domain.
Pspi: Streamlining 3D Echo-Reconstructive Imaging
Subsurface depth imaging� being as it is the outcome of repeated steps of �3D seismic data migrat... more Subsurface depth imaging� being as it is the outcome of repeated steps of �3D seismic data migration� requires Gbytes of data which must be reduced� transformed� visualized and interpreted to obtain meaningful information�. Severe performance requirements oblige high performance computing hardware and techniques� enormous effort has also historically gone into simplifying the migration model so as to reduce the cost of the operation while retaining the essential features of the signal propagation� One such model leads to the phase-shift-plus-interpolation �PSPI� algorithm� a highquality method for seismic migration� By optimizing the algorithm to reduce calculation requirements and exploiting its decoupling in the frequency domain for concurrency in parallel implementations� we have created in PSPI a cost-effective method�.
3D spectral reverse time migration with no‐wraparound absorbing conditions
SEG Technical Program Expanded Abstracts 1998, 1998
Acoustic earth imaging by seismic migration
Lecture Notes in Physics
Geophysical methods of prospecting for oil and gas resources are based on observations of elastic... more Geophysical methods of prospecting for oil and gas resources are based on observations of elastic wave fields, treated to a first approximation as acoustic waves. Emitted pressure impulses penetrate the crust and then are backscattered to the surface where they are ...
Stochastic lattice gases and parallel flow
Transport Theory and Statistical Physics, 1994
From the 3D phase-shift extrapolation of individual phase-space components, we derive for prestac... more From the 3D phase-shift extrapolation of individual phase-space components, we derive for prestacked data a vector relation between kh, the horizontal offset wavenumber, and kz, the vertical one. While in 2D this relation depends only on the tangent of the scattering angle θ and not on the structural deep, in 3D, in general, this is no longer true. The resulting vector formula takes into account the orientation of the scattering plane containing the slowness vectors ps and pr, one describing the down-going wave and the other the up-going one.
Simulation numérique et mécanique statistique : extension et étude de quelques problèmes d'ingénierie
Le contenu de ce travail, rassemble en une serie d’articles –publies, ou en cours de publication-... more Le contenu de ce travail, rassemble en une serie d’articles –publies, ou en cours de publication- touche aussi bien a la physique, aux telecommunications qu’a l’optimisation combinatoire. Les idees utilisees s’articulent autour de deux disciplines : la simulation numerique et la mecanique statistique des systemes, aussi bien physiques que d’ingenierie. Nous entendons par systeme, un ensemble d’elements en interaction entre eux et avec l’environnement qui les entoure. Le point de vue adopte est celui du modele om le probleme est ramene a un niveau tel, qu’une analyse mathematique puisse etre envisagee, ou tout au moins formulee. Il s’est agi pour nous d’analyser l’ordre statistique qui prend place dans un espace des configurations dans lequel les interactions font evoluer les differents elements. Le cadre de depart a ete microscopique et l’utilisation d’une description complete, imposee jusque dans les plus petits details du systeme, nous a amenes tout naturellement a definir des pro...
Computers in Physics, 1998
The Datacrossing DSS: A Data-GRID Based Decision Support System for Groundwater Management
2006 Second IEEE International Conference on e-Science and Grid Computing (e-Science'06), 2006
... the technologies we adopted to design distributed GIS applications, providing original soluti... more ... the technologies we adopted to design distributed GIS applications, providing original solutions for an ... have played a key role in the development of an effective, authoritative, interoperable ... groundwater model to support decision-making in the field of sustainable water resources ...
This document presents the general problem of electromagnetic inversion in the frequency domain. ... more This document presents the general problem of electromagnetic inversion in the frequency domain. The 1-dimensional analysis of a stratified medium is performed at different frequencies; the data measured on the surface (the admittance or the apparent electric ...
PSPI for High Performance Echo-Reconstructive Imaging
Ernesto Bonomi (ernesto@ crs4. it) obtained his Ph. D. in Theoretical Physics at the University o... more Ernesto Bonomi (ernesto@ crs4. it) obtained his Ph. D. in Theoretical Physics at the University of Paris XI, Orsay and is currently the head of geophysics at CRS4. His earlier professional activities included work in statistical mechanics and computational physics; ...
Methods for Optimization in Presence of Relative Minima: A Conjugate-Direction Method
Seismic imaging of medical ultrasound data: Towards in vivo applications
Europhysics Letters
This study presents the application to medical ultrasonography of a high value-added imaging meth... more This study presents the application to medical ultrasonography of a high value-added imaging method developed in the context of geophysical exploration and based on the undulatory description of the physical process. Currently, the most serious limitation of conventional medical ultrasound systems is the implicit assumption of constant-velocity media stemming from the widespread use of geometrical ray-based imaging algorithms and leading to phase aberration phenomena in the case typical of the human body of two or more tissues with different velocities. As a result, imaging of targets underlying tissue layers can be severely degraded. To address this problem, the proposed algorithm is designed expressly for variable-velocity media, providing high-resolution images by implementing the solution of the one-way wave equation. From a macro-velocity model, the illuminated structures are imaged by correlating, in the frequency domain and for each source element, the simulated emitted wavef...
The SmartGeo Portal: A retrospective
Although time-domain depth migration techniques have been successfully ported to run on modern ha... more Although time-domain depth migration techniques have been successfully ported to run on modern hardware accelerators, their ultimate obstacle is the I/O overhead present during the imaging step. Frequency-domain depth migration algorithms overcome this limitation and can exploit the full potential of new computing technologies. In particular, our implementation of Phase Shift Plus Interpolation (PSPI) method is characterized by fast running time, good quality results under lowsignal-to-noise ratio conditions and excellent results for steep dips. We provide a novel computational dataflow scheme to perform acceleration of PSPI on a generic dataflow engine. We present speedup results obtained on the state-of-the-art dataflow technology for synthetic VTI datasets. Our measurements indicate that a dataflow approach can achieve high speedups despite larger and larger computational domains, increased complexity of the anisotropic approach and the I/O overhead during angle-gathers calculation.
Physical Review Applied, 2020
Conventional medical ultrasound systems normally implement ray-based imaging algorithms, such as ... more Conventional medical ultrasound systems normally implement ray-based imaging algorithms, such as delay-and-sum beamforming, whose severest limitation derives from the implicit assumption of constantvelocity media. As a result, in the case of two or more tissues with different velocities, the image of the underlying targets appears strongly degraded both in placement and in resolution. The proposed ultrasound-imaging strategy, a value-added application of concepts developed in the context of seismic prospecting, avoids this restriction by relying on the undulatory description of the physical process and not on the geometric one. Echoes are sensed in the synthetic aperture configuration by the transducer, whose elements sequentially emit a nearly spherical wave front that covers the whole region of interest. Given a macrovelocity model, the recorded echoes become the initial condition for the downward propagator of the ultrasound wavefield. The processing is in three steps: time-reverse propagation of the sensed echoes, forward simulation of the wavefield emitted by the source element, and partial imaging by computing the zero-lag temporal correlation of the two wavefields to detect the scattering structures. The final reconstruction is obtained by stacking all the partial results on a single image. Laboratory tests, performed on experimental data acquired on a physical phantom, with and without an aberrant layer, prove the effectiveness of the proposed method even in the case of vertical and lateral velocity variations, with images of impressive spatial resolution and highly accurate target placement.
Amplitude-Preserving Monte Carlo 3D Prestack Migration
64th EAGE Conference & Exhibition, 2002
Cloud computing is establishing worldwide as a new high performance computing paradigm that offer... more Cloud computing is establishing worldwide as a new high performance computing paradigm that offers formidable possibilities to industry and science. The presented cloud-computing portal, part of the Grida3 project, provides an innovative approach to seismic data processing by combining open-source state-of-the-art processing software and cloud-computing technology, making possible the effective use of distributed computation and data management with administratively distant resources. We substituted the user-side demanding hardware and software requirements by remote access to high-performance grid-computing facilities. As a result, data processing can be done quasi in real-time being ubiquitously controlled via Internet by a user-friendly web-browser interface. Besides the obvious advantages over locally installed seismic-processing packages, the presented cloud-computing solution creates completely new possibilities for scientific education, collaboration, and presentation of repr...
3D Zero-Offset CRS Stack for Narrow-Azimuth Data - Formulation and Examples
EAGE/SEG Summer Research Workshop - Processing and imaging of seismic data, 2003
The common-reflection-surface (CRS) stack is a macro-model independent seismic reflection imaging... more The common-reflection-surface (CRS) stack is a macro-model independent seismic reflection imaging method that has proved to accurately detect and characterise meaningful events in the pre-stack data. As any 3D data-driven method, it requires an azimuthal distribution of traces as uniform as possible, in order to properly infer necessary information from the data. Its application to 3D datasets with narrow-azimuth acquisitions requires some simplified assumptions of the underlying theory. We present an approximated formulation and a related processing method that still represents an improvement with respect to conventional techniques currently employed. The advantages of the method when performing stacking velocity analysis are demonstrated on a synthetic example
Earth Science Informatics, 2009
Shared Resources Manager for Environmental Data Analysis and Applications) is a multidisciplinary... more Shared Resources Manager for Environmental Data Analysis and Applications) is a multidisciplinary project designed to deliver an integrated system to forge solutions to some environmental challenges such as the constant increase of polluted sites, the sustainability of natural resources usage and the forecast of extreme meteorological events. The GRIDA3 portal is mainly based on Web 2.0 technologies and EnginFrame framework. The portal, now at an advanced stage of development, provides end-users with intuitive Web-interfaces and tools that simplify job submission to the underneath computing resources. The framework manages the user authentication and authorization, then controls the action and job execution into the grid computing environment, collects the results and transforms them into an useful format on the client side. The GRIDA3 Portal framework will provide a problemsolving platform allowing, through appropriate access policies, the integration and the sharing of skills, resources and tools located at multiple sites across federated domains.
Journal of Geophysics and Engineering, 2017
The present article outlines an innovative derivation of the well-known CRS traveltime formula. T... more The present article outlines an innovative derivation of the well-known CRS traveltime formula. This is made possible by approximating the isochrones tangent to a small reflecting element embedded in a 3D homogeneous auxiliary medium. The arising paraxial traveltime formula is then parametrised by six attributes, each one characterising geometrically the wavefront of the reference normal-incident ray at the emergence point x 0 of the datum plane. For a layered medium, the assignment of the six attributes for each azimuthal direction around x 0 describes locally the emerging wavefront and establishes the mapping between reflecting elements of the two media, the auxiliary and the real, which respond, however, with different traveltimes. This 3D homeomorphism provides the time correction that gives rise to the CRS traveltime formula with eight attributes. For time imaging applications, the traveltime profile must be numerically shaped by improving iteratively the value of the eight attributes, so as to intercept, without the need of a velocity model, the largest number of coherent data in the volume of seismic traces gathered in the midpoint-offset domain.
Pspi: Streamlining 3D Echo-Reconstructive Imaging
Subsurface depth imaging� being as it is the outcome of repeated steps of �3D seismic data migrat... more Subsurface depth imaging� being as it is the outcome of repeated steps of �3D seismic data migration� requires Gbytes of data which must be reduced� transformed� visualized and interpreted to obtain meaningful information�. Severe performance requirements oblige high performance computing hardware and techniques� enormous effort has also historically gone into simplifying the migration model so as to reduce the cost of the operation while retaining the essential features of the signal propagation� One such model leads to the phase-shift-plus-interpolation �PSPI� algorithm� a highquality method for seismic migration� By optimizing the algorithm to reduce calculation requirements and exploiting its decoupling in the frequency domain for concurrency in parallel implementations� we have created in PSPI a cost-effective method�.
3D spectral reverse time migration with no‐wraparound absorbing conditions
SEG Technical Program Expanded Abstracts 1998, 1998
Acoustic earth imaging by seismic migration
Lecture Notes in Physics
Geophysical methods of prospecting for oil and gas resources are based on observations of elastic... more Geophysical methods of prospecting for oil and gas resources are based on observations of elastic wave fields, treated to a first approximation as acoustic waves. Emitted pressure impulses penetrate the crust and then are backscattered to the surface where they are ...
Stochastic lattice gases and parallel flow
Transport Theory and Statistical Physics, 1994