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...
Main Authors: | , , , , |
---|---|
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 |