An adaptive inversion algorithm for one-dimensional magnetotelluric problems

<p>      As we know, the traditional one-dimensional (1-D) magnetotelluric (MT) regularization inversion needs the geometry model of the 1-D Earth conductivity model, i.e., the number of layers and the thickness of each layer to be given in advance and cannot be changed during the inversion. In this way, too few layers cannot approximate the 1-D conductivity model accurately, while too many layers will increase the non-uniqueness of the inversion problem and hence may result in unreasonable results. Aiming to solve this issue, an adaptive inversion algorithm has been proposed for 1-D MT problems, where the layer number and the thickness of each layer can be adjusted automatically during the inversion process. To this end, three pseudo a-posterior error estimators has been proposed to guide the adjustment of the 1D geometry model, which are based on the gradient of the data misfit term of the penalty function, the diagonal elements of the model resolution matrix, and the weighted elements of the sensitivity matrix, respectively. The inversion results of the synthetic and field data by using our proposal adaptive inversion algorithm and the traditional regularization inversion not only validate the proposed algorithm, but also show that our proposed algorithm can obtain more accurate and reasonable results than traditional one. Subsequently, the proposed algorithm will be extended for 3-D magnetotelluric inversion problems soon.</p>

Effect of Displacement Current on the Finite-Difference Modeling of Natural Source Electromagnetic Diffusion

For electromagnetic (EM) modeling based on the electric-field formulation at low frequencies, the quasi-static approximation (i.e., only the conduction current is considered and the displacement current is ignored) is commonly applied, and a small conductivity value for the air layer is chosen subjectively. Actually, in the air layer, due to the use of the small conductivity value, the quasi-static approximation is ubiquitously violated. However, the effect of the violation of the quasi-static approximation in the air on EM modeling is not well examined in the literature. In this paper, we investigate this issue by comparing the finite-difference modeling results from the calculation with the quasi-static approximation and those considering both the conduction and displacement currents. For the quasi-static approximation, the conductivity in the air is set to be different small values, and zero air conductivity is used for the modeling with both the conduction and displacement currents considered. Several simple models are designed to verify the numerical solution and study how the assigned conductivity for the air affects the modeling accuracy. One flat model and two models with topography are designed to examine the effect of the quasi-static approximation on the EM modeling results. For frequencies used in typical geophysical applications of EM diffusion, using the quasi-static approximation is able to produce accurate modeling results for models with typical earth conductivity. However, if the rough surface topography is considered, the use of the quasi-static approximation can reduce the EM modeling accuracy substantially at much lower frequencies (as low as several hundred Hz), which is probably due to the inaccurate description of EM waves in the air, and poses problems for applications based on direct EM field interpretation.

Modeling of direct current traction power system with the consideration for earth conductivity

Nowadays, one of the main tasks facing the JSC Russian Railways is the increase of volumes of freight transportation. The solution of this task is directly connected with the increase of the load on traction power system devices and their strengthening, which is primarily related to railways electrified with direct current that are characterized by high traction currents. In order to assess the possibility to strengthen a railway section and reveal «bottlenecks» in traction power supply system it is necessary to create a mathematical model of traction power system that must consider parameters of a specific section including earth conductivity that affects the values of skin resistance, potential of railway network and stray currents. The paper proposes a mathematical model of direct current traction power system that represents a combination of a catenary model and a model of railway network. The model considers parameters of traction substations, connection scheme of catenary suspensions and possible influence of rails from adjacent tracks on each other.

Solar quiet variation of the horizontal and vertical components of geomagnetic field using wavelet analysis

The solar quiet (Sq) variations of horizontal and vertical (SqH and SqZ) components of the geomagnetic field obtained from both the Northern Hemisphere and Southern Hemisphere of the International Real-Time Magnetic Observatory Network (INTERMAGNET) during solar maximum year 2001 were investigated. The results show enlargement of the SqH component of the geomagnetic field during the daytime, attributed to equatorial electrojet (EEJ) current closer to the geomagnetic equator at the electrojet stations (BNG and MBO), which are produced from large eastward flow of the current. It was observed that SqZ is positive at the southward and negative at the northward hemispheres. SqZ is amplified at HER and HBK around the daytime. Wavelet power spectrum based approach was employed to analyse the SqH, SqZ, and rate of induction (SqZ/SqH) time series in a sequence of time scaling from January to December. The higher energy of SqH and SqZ of the wavelet coefficients is noticeable at high frequency. The monthly variation rate of induction (SqZ/SqH) analyses during the Sq variations are associated with the influence of equatorwards penetration of electric fields from the field-aligned current, Earth conductivity, effect of the ocean, and ionospheric conductivity.

Quantitative influence of coast effect on geomagnetically induced currents in power grids: a case study

In recent years, several magnetic storms have disrupted the normal operation of power grids in the mid-low latitudes. Data obtained from the monitoring of geomagnetically induced currents (GIC) indicate that GIC tend to be elevated at nodes near the ocean-land interface. This paper discusses the influence of the geomagnetic coast effect on GIC in power grids based on geomagnetic data from a coastal power station on November 9, 2004. We used a three-dimensional (3D) Earth conductivity model to calculate the induced electric field using the finite element method (FEM), and compared it to a one-dimensional (1D) layered model, which could not incorporate a coastal effect. In this manner, the GIC in the Ling’ao power plant was predicted while taking the coast effect into consideration in one case and ignoring it in the other. We found that the GIC predicted by the 3D model, which took the coastal effect into consideration, showed only a 2.9% discrepancy with the recorded value, while the 1D model underestimated the GIC by 23%. Our results demonstrate that the abrupt lateral variations of Earth conductivity structures significantly influence GIC in the power grid. We can infer that high GIC may appear even at mid-low latitude areas that are subjected to the coast effect. Therefore, this effect should be taken into consideration while assessing GIC risk when power networks are located in areas with lateral shifts in Earth conductivity structures, such as the shoreline and the interfaces of different geological structures.