Invariant TE and TM impedances in the marine magnetotelluric method

2020 ◽  
Vol 221 (1) ◽  
pp. 163-177 ◽  
Author(s):  
A M Montiel-Álvarez ◽  
J M Romo ◽  
S Constable ◽  
E Gómez-Treviño

SUMMARY The magnetotelluric (MT) impedance tensor has a nil diagonal when one of the axes of the coordinate system coincides with the strike of a 2-D structure. In general, real data are full tensors either because of 3-D effects or measurements not aligned to the geological strike. The usual practice to adapt the field tensor to the 2-D assumption is to rotate to a new system of coordinates. In most cases, there is no single angle of rotation that warranties that the diagonal elements become zeros as in the ideal 2-D case. Even maximizing the off-diagonal elements does not necessarily produce a nil diagonal. Consequently, the 2-D inversions proceed by neglecting whatever there is left in the diagonals. In this work, we explore an alternative that places no constraints on direction but assures a nil diagonal. We use two rotational invariants that compact the four elements of the tensor into only two and reduce in 2-D to the TE and TM impedances. These are obtained readily by solving a quadratic equation. We explore four different scenarios: (1) using the invariants, (2) rotating the tensor perpendicular to the profile, (3) rotating to the average maximum orientation for each station and (4) maximizing the off-diagonal elements of the tensor for each site, frequency to frequency. These approaches were applied to 3-D synthetic and field data. The field data correspond to two marine magnetotelluric surveys in the Gulf of California. In one of them, there is no information on the instrument orientation because the compasses failed. In this case, the rotational invariants come handy to overcome the problem. In the other survey, there was orientation information and the 2-D inversions illustrate the better performance of the invariants relative to the traditional approaches.

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. E37-E57
Author(s):  
Shunguo Wang ◽  
Steven Constable ◽  
Valeria Reyes-Ortega ◽  
Hormoz Jahandari ◽  
Colin Farquharson ◽  
...  

Two-dimensional marine magnetotelluric (MT) observations are useful for offshore geologic studies, such as natural resource exploration, fault mapping, fluid estimation at subduction zones, and the delineation of the lithosphere-asthenosphere boundary beneath the seafloor. Earth structures are often assumed to be two dimensional, which allows MT data to be decomposed into a transverse electric (TE) mode and a transverse magnetic (TM) mode. The 2D assumption can effectively reduce acquisition and computational costs. However, offline 3D effects and other problems such as lack/failure of the compass on instruments are often encountered, making it difficult to decompose data into the TE and TM modes. In these cases, 2D inversion may be misleading or may not provide an acceptable misfit to the marine MT observations. Thus, we have developed a 2D determinant inversion to the marine MT method to mitigate these difficulties, implemented in the MARE2DEM code, and we tested its utility using synthetic examples and a field example. In the synthetic examples, the determinant inversion demonstrates an ability to overcome 3D effects caused by 3D anomalies and bathymetry. With confidence from the synthetic tests, we interpreted real data acquired in the Gulf of California, Mexico, where not only is the bathymetry 3D in nature, but the external compasses failed to record the orientation. The field data can not only be fit to a reasonable misfit with a determinant inversion, but the resolved conductive zones also have a good correlation with known faults. A comparison between the resistivity model from the field data and a seismic reflection section shows that a previously interpreted fault, the Wagner Fault, should be shifted 5 km toward the southwest and made slightly steeper. Thus, the implementation of the determinant inversion may provide a new approach for using problematic 2D data.


2021 ◽  
Vol 11 (11) ◽  
pp. 5025
Author(s):  
David González-Peña ◽  
Ignacio García-Ruiz ◽  
Montserrat Díez-Mediavilla ◽  
Mª. Isabel Dieste-Velasco ◽  
Cristina Alonso-Tristán

Prediction of energy production is crucial for the design and installation of PV plants. In this study, five free and commercial software tools to predict photovoltaic energy production are evaluated: RETScreen, Solar Advisor Model (SAM), PVGIS, PVSyst, and PV*SOL. The evaluation involves a comparison of monthly and annually predicted data on energy supplied to the national grid with real field data collected from three real PV plants. All the systems, located in Castile and Leon (Spain), have three different tilting systems: fixed mounting, horizontal-axis tracking, and dual-axis tracking. The last 12 years of operating data, from 2008 to 2020, are used in the evaluation. Although the commercial software tools were easier to use and their installations could be described in detail, their results were not appreciably superior. In annual global terms, the results hid poor estimations throughout the year, where overestimations were compensated by underestimated results. This fact was reflected in the monthly results: the software yielded overestimates during the colder months, while the models showed better estimates during the warmer months. In most studies, the deviation was below 10% when the annual results were analyzed. The accuracy of the software was also reduced when the complexity of the dual-axis solar tracking systems replaced the fixed installation.


Geophysics ◽  
1989 ◽  
Vol 54 (4) ◽  
pp. 497-507 ◽  
Author(s):  
Jorge W. D. Leão ◽  
João B. C. Silva

We present a new approach to perform any linear transformation of gridded potential field data using the equivalent‐layer principle. It is particularly efficient for processing areas with a large amount of data. An N × N data window is inverted using an M × M equivalent layer, with M greater than N so that the equivalent sources extend beyond the data window. Only the transformed field at the center of the data window is computed by premultiplying the equivalent source matrix by the row of the Green’s matrix (associated with the desired transformation) corresponding to the center of the data window. Since the inversion and the multiplication by the Green’s matrix are independent of the data, they are performed beforehand and just once for given values of N, M, and the depth of the equivalent layer. As a result, a grid operator for the desired transformation is obtained which is applied to the data by a procedure similar to discrete convolution. The application of this procedure in reducing synthetic anomalies to the pole and computing magnetization intensity maps shows that grid operators with N = 7 and M = 15 are sufficient to process large areas containing several interfering sources. The use of a damping factor allows the computation of meaningful maps even for unstable transformations in the presence of noise. Also, an equivalent layer larger than the data window takes into account part of the interfering sources so that a smaller damping factor is employed as compared with other damped inversion methods. Transformations of real data from Xingú River Basin and Amazon Basin, Brazil, demonstrate the contribution of this procedure for improvement of a preliminary geologic interpretation with minimum a priori information.


2013 ◽  
Vol 2013 ◽  
pp. 1-11
Author(s):  
Changhong Hu ◽  
Shufen Liu ◽  
Ramana Reddy ◽  
Sumitra Reddy ◽  
Mingyang Liu

A model named KaM_CRK is proposed, which can supply the clustered and ranked knowledge to the users on different contexts. By comparing the attributes of contexts and JANs, our findings indicate that our model can accumulate the JANs, whose attributes are similar with the user’s contexts, together. By applying the KaM_CLU algorithm and Centre rank strategy into the KaM_CRK model, the model boosts a significant promotion on the accuracy of provision of user's knowledge. By analyzing the users' behaviors, the dynamic coefficient BehaviorFis first presented in KaM_CLU. Compared to traditional approaches of K_means and DBSCAN, the KaM_CLU algorithm does not need to initialize the number of clusters. Additionally, its synthetic results are more accurate, reasonable, and fit than other approaches for users. It is known from our evaluation through real data that our strategy performs better on time efficiency and user's satisfaction, which will save by 30% and promote by 5%, respectively.


Geophysics ◽  
1975 ◽  
Vol 40 (5) ◽  
pp. 790-797 ◽  
Author(s):  
Larry R. Grillot

A method of analyzing the tensor nature of magnetotelluric fields from short independent sections of filtered MT data is presented in which: (1) horizontal field components are filtered at chosen periods; (2) short sections which show a well defined signal over 3–4 cycles are selected from the records; (3) the polarizations of these short record sections are plotted; and (4) Fourier coefficients calculated for section pairs showing differing magnetic field polarizations are then used to calculate the elements of the impedance tensor. Results of this analysis performed on MT field data show that the use of MT signal pairs yields uniform principal resistivities for calculations using several different section pairs of the same period. These results also illustrate internal consistency of the resistivity estimates over the complete period range of the data (30–7500 sec). A brief comparison of the MT pairs calculations and results derived from a least‐squares approach at a selected period shows good agreement between the two methods. The overall results suggest that the use of independent signal pairs in the determination of the tensor impedance could complement more traditional methods in some cases.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. G15-G23
Author(s):  
Andrea Vitale ◽  
Domenico Di Massa ◽  
Maurizio Fedi ◽  
Giovanni Florio

We have developed a method to interpret potential fields, which obtains 1D models by inverting vertical soundings of potential field data. The vertical soundings are built through upward continuation of potential field data, measured on either a profile or a surface. The method assumes a forward problem consisting of a volume partitioned in layers, each of them homogeneous and horizontally finite, but with the density changing versus depth. The continuation errors, increasing with the altitude, are automatically handled by determining the coefficients of a third-order polynomial function of the altitude. Due to the finite size of the source volume, we need a priori information about the total horizontal extent of the volume, which is estimated by boundary analysis and optimized by a Markov chain process. For each sounding, a 1D inverse problem is independently solved by a nonnegative least-squares algorithm. Merging of the several inverted models finally yields approximate 2D or 3D models that are, however, shown to generate a good fit to the measured data. The method is applied to synthetic models, producing good results for either perfect or continued data. Even for real data, i.e., the gravity data of a sedimentary basin in Nevada, the results are interesting, and they are consistent with previous interpretation, based on 3D gravity inversion constrained by two gamma-gamma density logs.


2019 ◽  
Author(s):  
Mike W.-L. Cheung

Conventional meta-analytic procedures assume that effect sizes are independent. When effect sizes are non-independent, conclusions based on these conventional models can be misleading or even wrong. Traditional approaches, such as averaging the effect sizes and selecting one effect size per study, are usually used to remove the dependence of the effect sizes. These ad-hoc approaches, however, may lead to missed opportunities to utilize all available data to address the relevant research questions. Both multivariate meta-analysis and three-level meta-analysis have been proposed to handle non-independent effect sizes. This paper gives a brief introduction to these new techniques for applied researchers. The first objective is to highlight the benefits of using these methods to address non-independent effect sizes. The second objective is to illustrate how to apply these techniques with real data in R and Mplus. Researchers may modify the sample R and Mplus code to fit their data.


2008 ◽  
Vol 23 (1-2) ◽  
pp. 45
Author(s):  
G. Verdugo-Díaz ◽  
M. O. Albáñez-Lucero ◽  
R. Cervantes Duarte

Se analizó la variabilidad de la productividad primaria en Bahía de La Paz, durante noviembre de 1997 y 2000, enero de 1998 y febrero de 2001. Se identificó Agua del Golfo de California en mayor proporción; Agua Superficial Ecuatorial y Agua Subsuperficial Subtropical. En noviembre de 1997 se registró el valor promedio máximo de temperatura (~27.62 °C) y en febrero el valor mínimo (~19.16 °C). El índice de Simpson (Ö) indicó que durante noviembre de 1997 (~286.6 J m-3) la columna de agua presentó mayor estratificación, mientras que en noviembre de 1998 (~60.4 J m-3) y febrero de 2001 (~94.5 J m-3) se encontró mezclada. En noviembre de 1997 se registraron bajos valores de nitratos (~1.52 µM), nitritos (~0.04 µM) y clorofila a (~0.44 mg Cla m-3). En febrero de 2001 se registraron valores mayores de nitratos (~7.46 µM), nitritos (~0.76 µM) y clorofila a (~1.07 mg Cla m-3). La productividad primaria superficial e integrada presentaron sus promedios máximos en noviembre de 2000 (~5.09 mg C m-3 h-1 y 75.54. mg m-2 h-1) y mínimos en noviembre de 1997 (~2.93 mg C m-3 h-1 y ~55.29 mg m-2 h-1), respectivamente. De acuerdo a investigaciones anteriores, los valores reportados de productividad son característicos de la temporada fría de la bahía con excepción de noviembre de 1997, que presentó influencia de El Niño. Primary productivity stimation during autumm-winter in Bahía de La Paz, B.C.S., Mexico We analyzed the variability of primary productivity in Bahía de La Paz during November 1997 and 2000, January 1998 and February 2001. Water was identified in the Gulf of California in higher proportion; Surface Water and Water Subsuperficial Equatorial Subtropical. In November 1997 we saw the average maximum temperature (~ 27.62 ° C) in February and the lowest value (~ 19.16 ° C). The index of Simpson (Ö) indicated that during November 1997 (~ 286.6 J-3 m) column of water with greater stratification, while in November 1998 (~ 60.4 m J-3) and February 2001 (~ 94.5 J m-3) was mixed. In November 1997 there were low values of nitrate (~ 1.52 µM), nitrites (~ 0.04 µM) and chlorophyll a (Cla ~ 0.44 mg m-3). In February 2001, there were larger values of nitrate (~ 7.46 µM), nitrites (~ 0.76 µM) and chlorophyll a (Cla ~ 1.07 mg m-3). The primary productivity and integrated surface presented their highest averages in November 2000 (C ~ 5.09 mg m-3 h-1 and 75.54. mg m-2 h-1) and minimum in November 1997 (C ~ 2.93 mg m-3 h-1 and ~ 55.29 mg m-2 h-1), respectively. According to previous research, reported productivity values are characteristic of the cold season of the bay with exception of November 1997, which showed influence of El Niño.


2020 ◽  
Author(s):  
David Martínez van Dorth ◽  
Luca D'Auria ◽  
Juanjo Ledo Fernández ◽  
Perla Piña-Varas ◽  
Federico Di Paolo ◽  
...  

<p>The magnetotelluric method (MT) is a geophysical technique that provides high resolution information of the electrical resistivity of the subsurface geological structures by measuring the natural variations of the electromagnetic field recorded on the surface. Among the numerous applications, it can be used to map the presence of fluid reservoirs and localize significant structural contrasts that could be related to the presence of a geothermal or volcanic system. However, the interference of the anthropogenic noise during the MT measurements could affect significantly the correct interpretation of the collected data.</p><p>For this reason, in order to evaluate the effect of data contamination by anthropogenic sources, we analyzed the data registered by a continuous recording magnetotelluric station located inside the caldera of Las Cañadas (Tenerife, Spain). The instrumentation consisted of an ADU-08e, equipped with EPF-06 electrodes and MFS-06 magnetic coils. Two electric (Ex, Ey) and three magnetic (Hx, Hy, Hz) components have been recorded. This geophysical station was installed by the Instituto Volcanológico de Canarias (INVOLCAN), with purposes of volcano monitoring, on June 2019 and since then it has been recording data daily in the frequency range of 0.001 – 1000s.</p><p>On September 29 (2019) a significant electric blackout took place in the entire island of Tenerife in which, during approximately 6 hours the electricity supply was completely shut down. This situation represented a clear opportunity to obtain raw data almost free of anthropogenic contamination and it could help to quantify the effect of the anthropogenic noise in the MT measurements performed in a densely urbanized area as Tenerife. The first results show the clear change at 13:11:39 local time (GMT) in which both the electrical and magnetic components evidenced a pronounced change in their temporal pattern. Moreover, the comparison of the impedance tensor components between the previous hours and during the blackout reveals a noticeable difference for periods higher than 1 s.</p>


2014 ◽  
Vol 496-500 ◽  
pp. 2852-2856
Author(s):  
Jia Min Zhou ◽  
Jin Rui Guo ◽  
Ren Er Yang ◽  
Wei Hua Li

This paper studies the integer solution of the Markowitz mean-variance model. To avoid solving the quadratic equation with constraints directly, it uses geometric linear approximate method, and gives a practical and effective calculation method. Then it conducts the corresponding calculations regarding the real data, and reaches an optimal solution while the time of calculations is largely reduced, compared to the direct way and the algebraic approach.


Sign in / Sign up

Export Citation Format

Share Document