scholarly journals Reframing the magnetotelluric phase tensor for monitoring applications: improved accuracy and precision in strike determinations.

2021 ◽  
Author(s):  
Ana Gabriela Bravo-Osuna ◽  
Enrique Gómez-Treviño ◽  
Olaf Josafat Cortés-Arroyo ◽  
Néstor Fernando Delgadillo-Jáuregui ◽  
Rocío Fabiola Arellano-Castro

Abstract The magnetotelluric method is increasingly being used to monitor electrical resistivity changes in the subsurface. One of the preferred parameters derived from the surface impedance is the strike direction, which is very sensitive to changes in the direction of the subsurface electrical current flow. The preferred method for estimating the strike changes is that provided by the phase tensor because it is immune to galvanic distortions. However, it is also a fact that the associated analytic formula is unstable for noisy data, something that limits its applicability for monitoring purposes, because in general this involves comparison of two or more very similar data sets. One of the issues is that the noise complicates the distribution of estimates between the four quadrants. This can be handled by sending all values to the same quadrant by adding or subtracting the appropriate amount. This is justified by showing that the analytic formula is also a least squares solution. This is equivalent to define penalty functions for the matrix of eigenvalues and then select the minima numerically. Contrary to the analytic formula this numerical approach can be generalized to compute strikes using windows of any number of periods, thus providing tradeoffs between variance and resolution. The performance of the proposed approach is illustrated by its application to synthetic data and to real data from a monitoring array in the Cerro Prieto geothermal field, México.

2020 ◽  
Author(s):  
Ana Gabriela Bravo-Osuna ◽  
Enrique Gómez-Treviño ◽  
Olaf Josafat Cortés-Arroyo ◽  
Néstor Fernando Delgadillo-Jáuregui ◽  
Rocío Fabiola Arellano-Castro

Abstract The magnetotelluric method is increasingly being used to monitor electrical resistivity changes in the subsurface. One of the preferred parameters derived from the surface impedance is the strike direction, which is very sensitive to changes in the direction of the subsurface electrical current flow. The preferred method for estimating the strike changes is that provided by the phase tensor because it is immune to galvanic distortions. However, it is also a fact that the associated analytic formula is unstable for noisy data, something that limits its applicability for monitoring purposes, because in general this involves comparison of two or more very similar data sets. One of the issues is that the noise complicates the distribution of estimates between the four quadrants. This can be handled by sending all values to the same quadrant by adding or subtracting the appropriate amount. This is justified by showing that the analytic formula is also a least squares solution. This is equivalent to define penalty functions for the matrix of eigenvalues and then select the minima numerically. Contrary to the analytic formula this numerical approach can be generalized to compute strikes using windows of any number of periods, thus providing tradeoffs between variance and resolution. The performance of the proposed approach is illustrated by its application to synthetic data and to real data from a monitoring array in the Cerro Prieto geothermal field, México.


2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Ana G. Bravo-Osuna ◽  
Enrique Gómez-Treviño ◽  
Olaf J. Cortés-Arroyo ◽  
Nestor F. Delgadillo-Jauregui ◽  
Rocío F. Arellano-Castro

AbstractThe magnetotelluric method is increasingly being used to monitor electrical resistivity changes in the subsurface. One of the preferred parameters derived from the surface impedance is the strike direction, which is very sensitive to changes in the direction of the subsurface electrical current flow. The preferred method for estimating the strike changes is that provided by the phase tensor because it is immune to galvanic distortions. However, it is also a fact that the associated analytic formula is unstable for noisy data, something that limits its applicability for monitoring purposes, because in general this involves comparison of two or more very similar datasets. One of the issues is that the noise complicates the distribution of estimates between the four quadrants. This can be handled by sending all values to the same quadrant by adding or subtracting the appropriate amount. This is justified by showing that the analytic formula is also a least squares solution. This is equivalent to define penalty functions for the matrix of eigenvalues and then select the minima numerically. Contrary to the analytic formula, this numerical approach can be generalized to compute strikes using windows of any number of periods, thus providing tradeoffs between variance and resolution. The performance of the proposed approach is illustrated by its application to synthetic data and to real data from a monitoring array in the Cerro Prieto geothermal field, México.


2020 ◽  
Author(s):  
Ana Gabriela Bravo-Osuna ◽  
Enrique Gómez-Treviño ◽  
Olaf Josafat Cortés-Arroyo ◽  
Néstor Fernando Delgadillo-Jáuregui ◽  
Rocío Fabiola Arellano-Castro

Abstract The magnetotelluric method is increasingly being used to monitor electrical resistivity changes in the subsurface. One of the preferred parameters derived from the surface impedance is the strike direction, which is very sensitive to changes in the direction of the subsurface electrical current flow. The preferred method for estimating the strike changes is that provided by the phase tensor because it is immune to galvanic distortions. However, it is also a fact that the associated analytic formula is unstable for noisy data, something that limits its applicability for monitoring purposes, because in general this involves comparison of two or more very similar data sets. On the other hand, the classical Swift’s approach for strike is very stable for noisy data but it is severely affected by galvanic distortions. In this paper we impose the criterion of Swift’s approach to the phase tensor. Rather than developing an analytical formula we optimize numerically the same criterion. This stabilizes the estimation of strike by relaxing an exact condition to an optimal condition in the presence of noise. This has the added benefit that it can be applied to windows of several periods, thus providing tradeoffs between variance and resolution. The performance of the proposed approach is illustrated by its application to synthetic data and to real data from a monitoring array in the Cerro Prieto geothermal field, México.


2020 ◽  
Author(s):  
Ana Gabriela Bravo-Osuna ◽  
Enrique Gómez-Treviño ◽  
Olaf Josafat Cortés-Arroyo ◽  
Néstor Fernando Delgadillo-Jáuregui ◽  
Rocío Fabiola Arellano-Castro

Abstract The magnetotelluric method is increasingly being used to monitor electrical resistivity changes in the subsurface. One of the preferred parameters derived from the surface impedance is the strike direction, which is very sensitive to changes in the direction of the subsurface electrical current flow. The preferred method for estimating the strike changes is that provided by the phase tensor because it is immune to galvanic distortions. However, it is also a fact that the associated analytic formula is unstable for noisy data, something that limits its applicability for monitoring purposes, because in general this involves comparison of two or more very similar data sets. On the other hand, the classical Swift’s approach for strike is very stable for noisy data but it is severely affected by galvanic distortions. In this paper we impose the criterion of Swift’s approach to the phase tensor. Rather than developing an analytical formula we optimize numerically the same criterion. This stabilizes the estimation of strike by relaxing an exact condition to an optimal condition in the presence of noise. This has the added benefit that it can be applied to windows of several periods, thus providing tradeoffs between variance and resolution. The performance of the proposed approach is illustrated by its application to synthetic data and to real data from a monitoring array in the Cerro Prieto geothermal field, México.


2020 ◽  
Vol 223 (3) ◽  
pp. 1565-1583
Author(s):  
Hoël Seillé ◽  
Gerhard Visser

SUMMARY Bayesian inversion of magnetotelluric (MT) data is a powerful but computationally expensive approach to estimate the subsurface electrical conductivity distribution and associated uncertainty. Approximating the Earth subsurface with 1-D physics considerably speeds-up calculation of the forward problem, making the Bayesian approach tractable, but can lead to biased results when the assumption is violated. We propose a methodology to quantitatively compensate for the bias caused by the 1-D Earth assumption within a 1-D trans-dimensional Markov chain Monte Carlo sampler. Our approach determines site-specific likelihood functions which are calculated using a dimensionality discrepancy error model derived by a machine learning algorithm trained on a set of synthetic 3-D conductivity training images. This is achieved by exploiting known geometrical dimensional properties of the MT phase tensor. A complex synthetic model which mimics a sedimentary basin environment is used to illustrate the ability of our workflow to reliably estimate uncertainty in the inversion results, even in presence of strong 2-D and 3-D effects. Using this dimensionality discrepancy error model we demonstrate that on this synthetic data set the use of our workflow performs better in 80 per cent of the cases compared to the existing practice of using constant errors. Finally, our workflow is benchmarked against real data acquired in Queensland, Australia, and shows its ability to detect the depth to basement accurately.


1995 ◽  
Vol 38 (3-4) ◽  
Author(s):  
F. Batini ◽  
M. Caputo ◽  
R. Console

In this paper we model the geometry of a seismic source as a dislocation occurring on an elemental flat fault in an arbitrary direction with respect to the fault plane. This implies the use of a fourth parameter in addition to the three usual ones describing a simple double couple mechanism. We applied the radiation pattern obtained from the theory to a computer code written for the inversion of the observation data (amplitudes and polarities of the first onsets recorded by a network of stations). It allows the determination of the fault mechanism gener- alized in the above mentioned way. The computer code was verified on synthetic data and then applied to real data recorded by the seismic network operated by the Ente Nazionale per l'Energia Elettrica (ENEL), monitoring the geothermal field of Larderello. The experimental data show that for some events the source mechanism exhibits a significant dipolar component. However, due to the high standard deviation of the amplitude data, F-test applied to the results of the analysis shows that only for two events the confidence level for the general- ized model exceeds 90%.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. W15-W30 ◽  
Author(s):  
Gary F. Margrave ◽  
Michael P. Lamoureux ◽  
David C. Henley

We have extended the method of stationary spiking deconvolution of seismic data to the context of nonstationary signals in which the nonstationarity is due to attenuation processes. As in the stationary case, we have assumed a statistically white reflectivity and a minimum-phase source and attenuation process. This extension is based on a nonstationary convolutional model, which we have developed and related to the stationary convolutional model. To facilitate our method, we have devised a simple numerical approach to calculate the discrete Gabor transform, or complex-valued time-frequency decomposition, of any signal. Although the Fourier transform renders stationary convolution into exact, multiplicative factors, the Gabor transform, or windowed Fourier transform, induces only an approximate factorization of the nonstationary convolutional model. This factorization serves as a guide to develop a smoothing process that, when applied to the Gabor transform of the nonstationary seismic trace, estimates the magnitude of the time-frequency attenuation function and the source wavelet. By assuming that both are minimum-phase processes, their phases can be determined. Gabor deconvolution is accomplished by spectral division in the time-frequency domain. The complex-valued Gabor transform of the seismic trace is divided by the complex-valued estimates of attenuation and source wavelet to estimate the Gabor transform of the reflectivity. An inverse Gabor transform recovers the time-domain reflectivity. The technique has applications to synthetic data and real data.


Author(s):  
P.L. Nikolaev

This article deals with method of binary classification of images with small text on them Classification is based on the fact that the text can have 2 directions – it can be positioned horizontally and read from left to right or it can be turned 180 degrees so the image must be rotated to read the sign. This type of text can be found on the covers of a variety of books, so in case of recognizing the covers, it is necessary first to determine the direction of the text before we will directly recognize it. The article suggests the development of a deep neural network for determination of the text position in the context of book covers recognizing. The results of training and testing of a convolutional neural network on synthetic data as well as the examples of the network functioning on the real data are presented.


2021 ◽  
Vol 22 (1) ◽  
Author(s):  
João Lobo ◽  
Rui Henriques ◽  
Sara C. Madeira

Abstract Background Three-way data started to gain popularity due to their increasing capacity to describe inherently multivariate and temporal events, such as biological responses, social interactions along time, urban dynamics, or complex geophysical phenomena. Triclustering, subspace clustering of three-way data, enables the discovery of patterns corresponding to data subspaces (triclusters) with values correlated across the three dimensions (observations $$\times$$ × features $$\times$$ × contexts). With increasing number of algorithms being proposed, effectively comparing them with state-of-the-art algorithms is paramount. These comparisons are usually performed using real data, without a known ground-truth, thus limiting the assessments. In this context, we propose a synthetic data generator, G-Tric, allowing the creation of synthetic datasets with configurable properties and the possibility to plant triclusters. The generator is prepared to create datasets resembling real 3-way data from biomedical and social data domains, with the additional advantage of further providing the ground truth (triclustering solution) as output. Results G-Tric can replicate real-world datasets and create new ones that match researchers needs across several properties, including data type (numeric or symbolic), dimensions, and background distribution. Users can tune the patterns and structure that characterize the planted triclusters (subspaces) and how they interact (overlapping). Data quality can also be controlled, by defining the amount of missing, noise or errors. Furthermore, a benchmark of datasets resembling real data is made available, together with the corresponding triclustering solutions (planted triclusters) and generating parameters. Conclusions Triclustering evaluation using G-Tric provides the possibility to combine both intrinsic and extrinsic metrics to compare solutions that produce more reliable analyses. A set of predefined datasets, mimicking widely used three-way data and exploring crucial properties was generated and made available, highlighting G-Tric’s potential to advance triclustering state-of-the-art by easing the process of evaluating the quality of new triclustering approaches.


2021 ◽  
Vol 40 (3) ◽  
pp. 1-12
Author(s):  
Hao Zhang ◽  
Yuxiao Zhou ◽  
Yifei Tian ◽  
Jun-Hai Yong ◽  
Feng Xu

Reconstructing hand-object interactions is a challenging task due to strong occlusions and complex motions. This article proposes a real-time system that uses a single depth stream to simultaneously reconstruct hand poses, object shape, and rigid/non-rigid motions. To achieve this, we first train a joint learning network to segment the hand and object in a depth image, and to predict the 3D keypoints of the hand. With most layers shared by the two tasks, computation cost is saved for the real-time performance. A hybrid dataset is constructed here to train the network with real data (to learn real-world distributions) and synthetic data (to cover variations of objects, motions, and viewpoints). Next, the depth of the two targets and the keypoints are used in a uniform optimization to reconstruct the interacting motions. Benefitting from a novel tangential contact constraint, the system not only solves the remaining ambiguities but also keeps the real-time performance. Experiments show that our system handles different hand and object shapes, various interactive motions, and moving cameras.


Sign in / Sign up

Export Citation Format

Share Document