A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers

<p>Cross-borehole seismic tomography is a powerful tool to investigate the subsurface with a very high spatial resolution. In a set of boreholes, comprehensive three-dimensional investigations at different depths can be conducted to analyse velocity anisotropy effects due to local changes with...

Full description

Bibliographic Details
Main Authors: S. Hellmann, M. Grab, C. Patzer, A. Bauder, H. Maurer
Format: Article
Language:English
Published: Copernicus Publications 2023-07-01
Series:Solid Earth
Online Access:https://se.copernicus.org/articles/14/805/2023/se-14-805-2023.pdf
_version_ 1797770651709210624
author S. Hellmann
S. Hellmann
M. Grab
M. Grab
M. Grab
C. Patzer
A. Bauder
A. Bauder
H. Maurer
author_facet S. Hellmann
S. Hellmann
M. Grab
M. Grab
M. Grab
C. Patzer
A. Bauder
A. Bauder
H. Maurer
author_sort S. Hellmann
collection DOAJ
description <p>Cross-borehole seismic tomography is a powerful tool to investigate the subsurface with a very high spatial resolution. In a set of boreholes, comprehensive three-dimensional investigations at different depths can be conducted to analyse velocity anisotropy effects due to local changes within the medium. Especially in glaciological applications, the drilling of boreholes with hot water is cost-efficient and provides rapid access to the internal structure of the ice. In turn, movements of the subsurface such as the continuous flow of ice masses cause deformations of the boreholes and complicate a precise determination of the source and receiver positions along the borehole trajectories. Here, we present a three-dimensional inversion scheme that considers the deviations of the boreholes as additional model parameters next to the common velocity inversion parameters. Instead of introducing individual parameters for each source and receiver position, we describe the borehole trajectory with two orthogonal polynomials and only invert for the polynomial coefficients. This significantly reduces the number of additional model parameters and leads to much more stable inversion results. In addition, we also discuss whether the inversion of the borehole parameters can be separated from the velocity inversion, which would enhance the flexibility of our inversion scheme. In that case, updates of the borehole trajectories are only performed if this further reduces the overall error in the data sets. We apply this sequential inversion scheme to a synthetic data set and a field data set from a temperate Alpine glacier. With the sequential inversion, the number of artefacts in the velocity model decreases compared to a velocity inversion without borehole adjustments. In combination with a rough approximation of the borehole trajectories, for example, from additional a priori information, heterogeneities in the velocity model can be imaged similarly to an inversion with fully correct borehole coordinates. Furthermore, we discuss the advantages and limitations of our approach in the context of an inherent seismic anisotropy of the medium and extend our algorithm to consider an elliptic velocity anisotropy. With this extended version of the algorithm, we analyse the interference between a seismic anisotropy in the medium and the borehole coordinate adjustment. Our analysis indicates that the borehole inversion interferes with seismic velocity anisotropy. The inversion can compensate for such a velocity anisotropy. Based on the modelling results, we propose considering polynomials up to degree 3. For such a borehole trajectory inversion, third-order polynomials are a good compromise between a good representation of the true borehole trajectories and minimising compensation for velocity anisotropy.</p>
first_indexed 2024-03-12T21:26:15Z
format Article
id doaj.art-0d5702ebd5424801a92a2ee7bc7e497c
institution Directory Open Access Journal
issn 1869-9510
1869-9529
language English
last_indexed 2024-03-12T21:26:15Z
publishDate 2023-07-01
publisher Copernicus Publications
record_format Article
series Solid Earth
spelling doaj.art-0d5702ebd5424801a92a2ee7bc7e497c2023-07-28T07:55:24ZengCopernicus PublicationsSolid Earth1869-95101869-95292023-07-011480582110.5194/se-14-805-2023A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciersS. Hellmann0S. Hellmann1M. Grab2M. Grab3M. Grab4C. Patzer5A. Bauder6A. Bauder7H. Maurer8Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, SwitzerlandInstitute of Geophysics, ETH Zurich, Zurich, SwitzerlandLaboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, SwitzerlandInstitute of Geophysics, ETH Zurich, Zurich, SwitzerlandTerra Vermessungen AG, Othmarsingen, SwitzerlandGeological Survey of Finland (GTK), Espoo, FinlandLaboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, SwitzerlandSwiss Federal Institute for Forest, Snow and Landscape Research (WSL), Birmensdorf, SwitzerlandInstitute of Geophysics, ETH Zurich, Zurich, Switzerland<p>Cross-borehole seismic tomography is a powerful tool to investigate the subsurface with a very high spatial resolution. In a set of boreholes, comprehensive three-dimensional investigations at different depths can be conducted to analyse velocity anisotropy effects due to local changes within the medium. Especially in glaciological applications, the drilling of boreholes with hot water is cost-efficient and provides rapid access to the internal structure of the ice. In turn, movements of the subsurface such as the continuous flow of ice masses cause deformations of the boreholes and complicate a precise determination of the source and receiver positions along the borehole trajectories. Here, we present a three-dimensional inversion scheme that considers the deviations of the boreholes as additional model parameters next to the common velocity inversion parameters. Instead of introducing individual parameters for each source and receiver position, we describe the borehole trajectory with two orthogonal polynomials and only invert for the polynomial coefficients. This significantly reduces the number of additional model parameters and leads to much more stable inversion results. In addition, we also discuss whether the inversion of the borehole parameters can be separated from the velocity inversion, which would enhance the flexibility of our inversion scheme. In that case, updates of the borehole trajectories are only performed if this further reduces the overall error in the data sets. We apply this sequential inversion scheme to a synthetic data set and a field data set from a temperate Alpine glacier. With the sequential inversion, the number of artefacts in the velocity model decreases compared to a velocity inversion without borehole adjustments. In combination with a rough approximation of the borehole trajectories, for example, from additional a priori information, heterogeneities in the velocity model can be imaged similarly to an inversion with fully correct borehole coordinates. Furthermore, we discuss the advantages and limitations of our approach in the context of an inherent seismic anisotropy of the medium and extend our algorithm to consider an elliptic velocity anisotropy. With this extended version of the algorithm, we analyse the interference between a seismic anisotropy in the medium and the borehole coordinate adjustment. Our analysis indicates that the borehole inversion interferes with seismic velocity anisotropy. The inversion can compensate for such a velocity anisotropy. Based on the modelling results, we propose considering polynomials up to degree 3. For such a borehole trajectory inversion, third-order polynomials are a good compromise between a good representation of the true borehole trajectories and minimising compensation for velocity anisotropy.</p>https://se.copernicus.org/articles/14/805/2023/se-14-805-2023.pdf
spellingShingle S. Hellmann
S. Hellmann
M. Grab
M. Grab
M. Grab
C. Patzer
A. Bauder
A. Bauder
H. Maurer
A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
Solid Earth
title A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
title_full A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
title_fullStr A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
title_full_unstemmed A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
title_short A borehole trajectory inversion scheme to adjust the measurement geometry for 3D travel-time tomography on glaciers
title_sort borehole trajectory inversion scheme to adjust the measurement geometry for 3d travel time tomography on glaciers
url https://se.copernicus.org/articles/14/805/2023/se-14-805-2023.pdf
work_keys_str_mv AT shellmann aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT shellmann aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT cpatzer aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT abauder aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT abauder aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT hmaurer aboreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT shellmann boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT shellmann boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT mgrab boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT cpatzer boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT abauder boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT abauder boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers
AT hmaurer boreholetrajectoryinversionschemetoadjustthemeasurementgeometryfor3dtraveltimetomographyonglaciers