Next Article in Journal
Classification Algorithm of 3D Pattern Film Using the Optimal Widths of a Histogram
Previous Article in Journal
Neuro-Robotic Synergy: Crafting the Secure Future of Industries in the Post Pandemic Era
Previous Article in Special Issue
A Next POI Recommendation Based on Graph Convolutional Network by Adaptive Time Patterns
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TEM Strata Inversion Imaging with IP Effect Based on Enhanced GCN by Extracting Long-Dependency Features

1
School of Information Engineering, Hubei University of Economics, Wuhan 430205, China
2
State Key Laboratory of Power Transmission Equipment and System Security and New Technology, Chongqing University, Chongqing 400044, China
3
Hubei Internet Finance Information Engineering Technology Research Center, Hubei University of Economics, Wuhan 430205, China
*
Authors to whom correspondence should be addressed.
Electronics 2023, 12(19), 4138; https://doi.org/10.3390/electronics12194138
Submission received: 25 August 2023 / Revised: 29 September 2023 / Accepted: 2 October 2023 / Published: 4 October 2023
(This article belongs to the Special Issue Mechanism and Modeling of Graph Convolutional Networks)

Abstract

:
Utilizing neural network models to inverse time-domain electromagnetic signals enables rapid acquisition of electrical structures, a non-intrusive method widely applied in geological and environmental surveys. However, traditional multi-layer perceptron (MLP) feature extraction is limited, struggling with cases involving complex electrical media with induced polarization effects, thereby limiting the inversion model’s predictive capacity. A graph-topology-based neural network model for strata electrical structure imaging with long-dependency feature extraction was proposed. We employ graph convolutional networks (GCN) for capturing non-Euclidean features like resistivity-thickness coupling and Long Short-Term Memory (LSTM) to capture long-dependency features. The LSTM compensates for GCN’s constraints in capturing distant node relationships. Using case studies with 5-strata and 9-strata resistivity models containing induced polarization effects, compared to traditional MLP networks, the proposed model utilizing time-domain features and graph-topology-based electrical structure extraction significantly improves performance. The mean absolute error in inversion misfit is reduced from 10–20% to around 2–3%.

1. Introduction

Using coils or long straight wires as signal transmitters generates time-domain-induced electromagnetic signals underground, enabling inference of the basic structure of subsurface electrical media [1,2,3,4]. This non-intrusive geophysical method, known as transient electromagnetic (TEM) detection, is extensively employed for subsurface imaging and exploration [5,6,7], particularly for low-resistivity target areas [8]. It finds broad applications in environmental monitoring [9,10], source exploration [11,12], and geological studies [13,14].
The process of inferring the electrical media structure from acquired time-domain electromagnetic induction signals is termed TEM inversion [15,16]. Strata media inversion based on single stations is commonly utilized for extensive geological surveys, often involving substantial data processing [13,17,18]. Traditional optimization methods like Newton’s gradient descent are computationally inefficient for massive TEM survey data [15,19,20,21,22,23]. The rise of electromagnetic signal inversion using traditional multi-layer perceptron (MLP) addresses this challenge, with network models progressively deepening the extraction of electromagnetic signal features and geoelectric structural characteristics per layer [24,25,26]. However, a targeted feature extraction strategy is still lacking to enhance neural network performance [27].
When subsurface electrical media exhibit induced polarization (IP) effects [28,29], charges rearrange, impacting the trend of electromagnetic induction signals within specific time windows [30,31,32,33,34]. This further exacerbates the nonlinearity of TEM inversion, posing a greater challenge to feature extraction in neural network inversion models.
For time-domain-induced electromagnetic signals, long short-term memory (LSTM) passes signals step by step via the time series, utilizing gate mechanisms to selectively retain time-domain signal features carrying crucial electrical media information [35,36]. This optimizes feature extraction from time-domain signals, enhancing the neural network model’s predictive capability in electromagnetic signal inversion. However, concerning geoelectric structural features, containing attributes like resistivity and layer thickness, traditional Convolutional Neural Network (CNN) feature extraction, which excels in image processing, struggles to accommodate the non-Euclidean nature of resistivity-layer thickness pairs [25,37]. Given the success of topological graphs in handling irregular non-Euclidean data [38] incapable of satisfying translational invariance in tasks like node classification [39,40], graph classification [41], and link prediction [42], leveraging graph neural networks to further extract electrical structural features incorporating resistivity-layer thickness pairs is a viable option.
Given these considerations, we establish a neural network inversion model for time-domain-induced electromagnetic signals based on GCN’s electrical media feature extraction. However, since each layer of GCN considers only first-order neighborhood information, it may struggle to capture complex relationships among distant nodes. To address this limitation, we use LSTM to capture long-dependency features [43,44], compensating for this shortcoming and enhancing the overall information aggregation capacity of the neural network model.
In the subsequent sections, we delve into the methodology in Section 2, reviewing the fundamental principles of TEM forward and inverse modeling. We elaborate on how the proposed neural network model employs LSTM to capture temporal dependency features from TEM-induced electromagnetic signals and uses GCN to capture spatial features of the geoelectric structure. Section 3 and Section 4 present illustrative examples using five-strata and nine-strata geoelectric models, respectively. The discussions are presented in Section 5 and conclusions are presented in Section 6.

2. TEM Inversion with IP Effect

2.1. TEM Forward and Inversion

Consider a TEM device with a central loop, with the surface located outside the observation center point and surrounded by a coil with a radius of a. The center point is considered the observation point of the station and is used to receive electromagnetic induction signals. The z-axis is positive and negative below the ground [45]. Pass a current of amplitude I through the coil and immediately turn it off. The vertical magnetic field generated by the central observation station at the time of the cut-off time of 0, Hz, can be expressed as [15,24]:
H z = I a 2 0 e u 0 z + r T E e u 0 z λ 2 u 0 J 1 λ a d λ ,
where z = 0, u0 can be approximated as λ, J1 represents a first-order Bessel function, which is discretized using the Hankel transform for numerical integration. The surface reflection coefficient rTE = (λu1)/(λ + u1). ui is obtained by the layered recursive formula [46,47]:
u i = u i u i + 1 + u i tanh ( u i h i ) u i + u i + 1 tanh ( u i h i ) u i = λ 2 k i 2 1 / 2 ,
here, N denotes the total stratum count of the geoelectric model, when i = N, uN = uN. k is the wave number, ki = (−iωμσ)1/2, ω and σ denote the angular frequency and conductivity.
The Cole-Cole model is one of the mathematical models used to describe the polarization response of materials. It is widely employed in geophysics and electromagnetic exploration. Based on the complex dielectric function, the model depicts the relationship between electrical conductivity and permittivity of the geoelectric structure at various frequencies [48]. It encompasses parameters such as the polarization strength index, relaxation time, and initial phase. The Cole-Cole model is one of the expressions of this model and can be described as [15,49]:
ρ ω = ρ 1 m 1 1 1 + i ω τ c ,
where ω denotes angular frequency, τ represents the relaxation time associated with the planning mechanism. In the equivalent circuit of the rock unit, c denotes the frequency correlation coefficient. The coefficient m controls the reactive power performance.
The time-domain-induced vertical magnetic field Hz(t) and induced electromotive force (EMF) ε can be transferred from frequency-domain induced magnetic field signals [3], and be obtained using the piecewise linear approximation method as follows [27,50]:
H z ( t ) = 2 π t 2 i = 0 N t Im [ F ω i + 1 ] ω i + 1 Im [ F ω i ] ω i ω i + 1 ω i cos ( ω i t ) cos ( ω i + 1 t ) ,
ε ( t ) = 2 μ π t 2 i = 0 N t Re [ F ω i + 1 ] Re [ F ω i ] ω i + 1 ω i cos ( ω i t ) cos ( ω i + 1 t ) ,
where Nt represents the number of time points, t denotes the time channel.

2.2. IP Effect on Inversion

When inversion is performed using the traditional gradient class method [15,51] when the initial geoelectric structure is a uniform half-space model, it is difficult to find the geoelectric structure that makes the TEM signal appear in the wrong direction (Figure 1a) during the inversion process, and the predicted geoelectric structure cannot be close to the preset structure (Figure 1b). The optimization process is stalled.

3. Our Networks Framework

3.1. Temporal Feature Capturing

The LSTM aims to tackle challenges such as gradient vanishing and the complexities associated with processing extended sequential data within recurrent neural networks. The structure of the LSTM model (Figure 2) involves an input layer, a recurrent hidden layer, and an output layer. Within the recurrent hidden layer, there exist components such as the forget gate (F), input gate (I), output gate (O), and memory cell (C). The memory cell (C) runs straight through the entire chain without any activation, preventing the gradients from vanishing when training the LSTM model using backpropagation [52].
LSTM model’s gating update process is illustrated as follows:
Step 1: The forget gate, labeled as Ft, plays a crucial role in deciding which information should be retained within the memory cell Ct−1 and what should be discarded [53]:
F t = σ W f [ h t 1 , x t ] + b f ,
Step 2: The input gate, denoted as It, is responsible for determining the new information that should be integrated into the memory cell [54,55]:
I t = σ W i [ h t 1 , x t ] + b i ,
C t = tanh W c [ h t 1 , x t ] + b c ,
Step 3: Transitioning from the prior memory cell Ct to the updated memory cell involves [44,56]:
C t = F t C t 1 + I t C t ,
Step 4: The role of the output gate, designated as Ot, is to regulate the quantity of information from the present memory cell Ct that should be directed towards the subsequent layer:
h t = O t tanh ( C t ) ,
where xt symbolizes the input to the LSTM model during time step t, while ht−1 denotes the preceding time step’s output from the recurrent hidden layer, and b denotes bias coefficients. Wf, Wi, Wo, Wc correspond to the weight matrices associated with F, I, O, and C, respectively.

3.2. Spatial Feature Capturing

We define the graph and feature matrix as follows:
Definition 1. 
We utilize a graph G = (V, E, A) to depict the topology of a geoelectric media structure, treating one stratum to be a distinct nodal point within this framework. Here, V = {v1, v2, …, vN} signifies the collection of nodal points that corresponds to the electrical layers, and N signifies the count of nodal points. E denotes the edges within the graph. The adjacency matrix A captures the connectivity among electrical strata, denoted as ARN×N. This adjacency matrix solely consists of elements 0 and 1. An element within the matrix is assigned a value of 0 if there is no linkage between two layers; conversely, it takes on a value of 1 to signify a connection.
Definition 2. 
We utilize the electrical information within the geoelectric structure as attributes for the nodes in the network, denoted as XRN×P; here, P denotes the count of attribute features for every node, corresponding to the length of temporal series. Additionally, XtRN×i refers to the time-domain electromagnetic induction signals within each time window of instance i.
Therefore, the challenge of TEM inversion for geoelectric structures m can be viewed as the acquisition of a supposing function f with the geoelectric strata topology G and the feature matrix X, as illustrated in Equation (1):
m   = f G ;   X ,
The GCN model constructs a Fourier filter, allowing it to operate on the graph’s nodes and capture spatial characteristics between them, specifically involving their first-order neighborhoods [57]. Multiple convolutional layers are then stacked to form the GCN model. Considering the adjacency matrix A, the propagation between layers in graph convolution is formulated as [58,59]:
H l + 1 = s i g m o i d D 1 / 2 A D 1 / 2 H l W l ,
where D denotes the degree matrix of A . Hl symbolizes the feature matrix for the l-th layer, while Wl pertains to the subsequent (l + 1)-th layer. The initial layer, denoted as H0, corresponds to the feature matrix X. Wl stands as the weight matrix for the l-th layer.
In this study, a 2-layer GCN model is selected for capturing spatial correlations, formulated as follows [57]:
f ( X , A ) = s i g m o i d A ReLU ( A X W 0 ) W 1 ,
A = A + I
A and D are calculated using the following formulas:
A i j = 1 , i j 0 , i = j ,
D i j = j = 1 N A i j , i = j 0 , i j ,
where A represents the connection between each node. We choose a 2-layer GCN model, and Equation (12) is rewritten as [60,61]:
H 2 = sigmoid D 1 / 2 A D 1 / 2 sigmoid D 1 / 2 A D 1 / 2 H 1 W 1 W 2 ,
We apply the GCN model to extract spatiotemporal characteristics from the time-domain-induced electromagnetic signals and the geoelectric strata structure. Illustrated in Figure 3, the model captures topological relationships between the central strata and its neighboring counterparts. It encodes both the road network topology and the geoelectric properties of the strata, thereby deriving spatial dependency relationships.

3.3. Network Inverse Frame

In TEM inversion, the induced time-domain electromagnetic response signals and the geoelectric structure exhibit complex non-Euclidean spatial patterns. The signal in each time-domain window is influenced by preceding time-domain window signals. Additionally, each geoelectric layer is interconnected with its adjacent layers. Our model comprises two primary components: a GCN and an LSTM network. We initially process the time-domain electromagnetic response signals of individual stations using the LSTM network to capture temporal features within the response signals (Figure 4). Subsequently, we employ the GCN to capture spatial correlations between different layers’ topological structures in the geoelectric structure. The electrical features within each layer are treated as signals on the graph. Ultimately, we derive the outcomes of geoelectric structure from the induced TEM responses using fully connected layers.
By feeding the time-domain electromagnetic induction signals of each observation station into the LSTM network, the two types of TEM signals containing spatial information features, namely the induced magnetic field Hz signal and the induced EMF signal ε, are processed to extract temporal features. This results in features representing early, middle, and late stages. Subsequently, these features are input into the GCN to extract geoelectric spatial features. The function is expressed as follows:
Y t = σ W H e Y H e + W H m Y H m + W H l Y H l + W v e Y v e + W v m Y v m + W v l Y v l
where Yt represents the prediction target for the t-th time interval, ⊙ denotes the element-wise multiplication, Ye, Ym, Yl represent the feature outputs for the early, middle, and late stages, respectively, extracted using the LSTM network, We, Wm, Wl are learnable parameter tensors representing the degrees of influence from the early, middle, and late stages, respectively. The subscript H and the subscript v, respectively, indicate the magnetic field signal and the induced voltage signal.

3.4. Training

We established separate training samples for each layer based on different depths and used them to train different neural networks. In constructing the training samples, we randomly sampled the resistivity information for each layer within a range of [1 × 10−1, 3.5 × 103]. For the layer thickness information in the geoelectric structure, we initially created a thickness sample library with values ranging from [10, 1000] in a logarithmic scale with a base of log2. For simulating the excitation polarization effect, we considered the Cole-Cole model and randomly sampled the parameters m, τ, and c within the range of 0–1.
We divided the dataset containing 10 million samples into training and validation sets in an 8:2 ratio. We trained the proposed neural networks for 100 epochs, with the mean squared error as the objective loss function. The training process for the neural network suitable for a five-layer geoelectric model is shown in Figure 5. To determine the values of hyperparameters such as learning rate, network depth, and regularization rate, we also conducted a sensitivity analysis on these parameters.
When the learning rate is set to 0.01, both the training loss and validation loss decrease slowly during training, in contrast to the performance when the learning rate is 0.1. This indicates that with a learning rate of 0.01, the adjustment step size for the coefficients of each neuron is too small, resulting in inefficient training of the neural network.
Additionally, we compared the training performance of the neural network with different network depths, specifically 5, 10, and 15 layers. Regardless of whether the learning rate is 0.01 or 0.1, deeper networks result in smaller fitting errors for the neural network. The performance of the proposed network model on the training and validation datasets is summarized in the Table 1.

4. Numerical Models

4.1. Model 1

To validate the precision of the GCN network model proposed in this study for transient electromagnetic inversion, we designed the multi-layer geoelectric simulation models and compared the inversion prediction capabilities of different neural network models based on TEM responses. Model 1 is set as a 5-strata geoelectric structure with a total of 9 predictive parameters.
In the setup, the TEM excitation-reception device uses a central loop and a circular loop with 20 m radius. A step turn-off current with an amplitude of 1 A is in the loop. The moment of current interruption is considered as time 0 s. The vertical component of the time-domain electromagnetic response signal is received starting from the observation station. The input signals to the neural network model are the induced magnetic field Hz and the induced EMF ε.
In the five-strata geoelectric model (Figure 6a), the resistivity of each layer is 50, 300, 100, 500, and 200 Ω·m, respectively. The thicknesses of the first four layers are 200, 250, 300, and 500 m, respectively. The bottom layer is considered an infinitely extending uniform space with infinite thickness. The third and fourth layers among the five layers are complex resistivity layers exhibiting IP effects. The Cole-Cole model represents the electrical structure of these two layers, and their frequency correlation coefficients are both set to 0.3.
From the perspective of predicted geoelectric structure (Figure 6a), all four neural network models accurately predict the ρ and thickness h of the 1-st stratum. However, significant misfits are shown in the predictive resistivity of the 2nd layer, with MLP exhibiting the largest deviation among the methods. This deviation is further magnified in the prediction of the fourth layer’s resistivity. The deviations follow a pattern of MLP having the highest, GCN and LSTM falling in the middle, and the proposed model has the smallest deviation.
Examining the accuracy of predicted polarization parameters for the third and fourth layers of this geoelectric model (Figure 6b), various methods exhibit relatively accurate restoration of the IP parameters for the third layer, while the accuracy drops for the fourth layer’s IP parameters. In the prediction of the fourth layer’s IP parameters, MLP performs the poorest, LSTM shows slight improvement, and the two GCN-based models achieve the highest accuracy in IP parameter prediction. The detailed predictions of the five-strata geoelectric model parameters by the four network models are represented in Table 2.
The accuracy of parameter predictions for the geoelectric model by different neural networks ranks from high to low as follows: our proposed model, GCN, LSTM, and MLP (Figure 7a). The inaccuracies in the geoelectric structure primarily concentrate on high-resistance layers, specifically layers with resistivities of 300 Ω·m and 500 Ω·m. In layers containing induced polarization effects, the accuracy of predicting the third-layer IP parameters generally surpasses that of predicting the fourth-layer IP parameters.
To comprehensively assess the prediction misfits of various network models, we conducted a comparative analysis of accuracy for the predicted geoelectric models’ corresponding TEM responses as well (Figure 7c,d). The induced magnetic field misfits analysis reveals that the proposed method’s prediction accuracy exhibits a substantial increase in two orders of magnitude in comparison to other methods, reducing the fitting discrepancy from 1 × 10−1 to 1 × 10−3. The advantage of the proposed method, not manifested in the magnetic field misfits, becomes apparent in the comparison of EMF misfit (Figure 7d) prior to 1 × 10−8 s. Throughout the prediction process, the comparison of mean absolute percentage errors (MAPE) and root mean square errors (RMSE) for predicted resistivity, thickness, IP parameters, induced magnetic field, and induced EMF using different network models is presented in Table 3.

4.2. Model 2

Model 2 is a nine-strata geoelectric model with a total of 17 predicted parameters. The reference comparative neural network models include MLP, LSTM, and traditional GCN. In the nine-strata geoelectric model (Figure 8a), the resistivity of each layer is, respectively, 50, 300, 100, 500, 200, 100, 300, 50, and 200 Ω·m. The thickness of the first four layers is 200, 250, 300, 500, 300, 100, 200, and 300 m. The bottom-most layer is considered to be an infinitely extending homogeneous space with infinite thickness. Layers 2, 4, 6, and 8 are complex resistivity layers exhibiting IP effects. The frequency correlation coefficients are 0.3, 0.4, 0.2, and 0.5, respectively.
From the predicted geoelectric structure (Figure 8a), it is observed that all four neural network models accurately predict the resistivity and thickness of the first, third, and fifth layers. However, for the second, fourth, seventh, and ninth layers, noticeable deviations are evident in the predicted resistivity results among different methods, with MLP exhibiting the largest deviation. For the resistivity of the fourth layer, except for the proposed method, MLP, LSTM, and GCN all show significant discrepancies. Below the fifth layer, MLP’s predicted results significantly deviate from the predefined geoelectric structure.
In the comparison of the accuracy in predicting the IP parameters (Figure 8b), it is observed that various methods achieve relatively accurate restoration of the IP parameters for the second and fourth layers. However, the accuracy decreases for predicting the IP parameters of the sixth and eighth layers. Similar to the resistivity-thickness relationship’s prediction results (Figure 8a), below the fifth layer’s structure, the predictive electric parameters of the traditional MLP model significantly deviate from the predefined geoelectric structure. Additionally, although LSTM’s prediction accuracy is higher than that of MLP, it still exhibits a noticeable misfit. The two types of GCN models exhibit the highest accuracy in predicting the IP parameters. Detailed information about the predictions of the 9-strata electric structure parameters by the four network models is shown in Table 4.
Taking into account the comprehensive comparison of the misfits of various neural networks in predicting geoelectric model parameters (Figure 9a,b), their prediction accuracy from high to low ranks is as follows: our model, GCN, LSTM, and MLP. The intervals with larger error values are primarily concentrated in high-resistivity layers, particularly for layers with resistivity values of 500 Ω·m and 300 Ω·m. For layers with induced polarization effects, MLP’s predictions show significant inaccuracies. LSTM starts exhibiting noticeable misfits from the fourth layer onwards.
We conducted a comparative analysis of the accuracy of TME responses predicted using various neural network models (Figure 9c,d). The induction magnetic field misfit shows that within the observation time window range of 1 × 10−6 to 1 × 10−3 s, the predictive accuracy of our method is significantly superior to that of the network model by two orders of magnitude, reducing the fitting discrepancy from 1 × 10−2 to 1 × 10−4 (Figure 9c). This advantage is also reflected in the misfit of the induced EMF (Figure 9d). For the late-stage induced magnetic field after 1 × 10−3 s, GCN-based network models perform slightly worse than LSTM-based models. However, this performance difference is not observed in the misfit of the induced EMF. Using MAPE and RMSE metrics, we performed a comparative analysis of the errors in predicted geoelectric model parameters and the errors in corresponding induced TEM signals predicted using different network models (Table 5).

5. Discussion

5.1. Denoising

To assess the stability and generalization ability of the proposed neural network for inverting TEM data with consideration of IP effects, we designed a case study involving a noisy 5-layer geoelectric model. Gaussian random noise at levels of 5% and 10% was added to the TEM responses of a five-layer geoelectric model, with a noise-free case serving as the reference. The resistivity values for the five layers of the geoelectric model were set as follows: ρ1 = 300 Ω·m, ρ2 = 200 Ω·m, ρ3 = 500 Ω·m, ρ4 = 50 Ω·m and ρ5 = 1000 Ω·m. The corresponding layer thicknesses were h1 = 200 m, h2 = 300 m, h3 = 300 m and h4 = 300 m. We considered a simple excitation polarization effect with variations only in the “m” coefficient, where the m values for the second and fourth layers were set to 0.3 and 0.4, respectively.
We compared the inversion performance of the proposed algorithm on noisy TEM data with noise levels of 5% and 10% (see Table 6). As the noise increased, both the MAPE and the RMSE gradually increased. To account for the noise levels, we subtracted the noise level from the MAPE, considering it as a relative MAPE error. The Re-MAPE increased continuously with the rising noise levels, indicating that as the noise grew, the neural network’s generalization ability decreased, rising from an initial 4.04% to 7.68%.
Furthermore, the RMSE reflected the level of data consistency, with larger RMSE values indicating poorer data consistency. This suggests that as noise levels increased, the stability of the neural network also decreased.
Figure 10 presents the inversion-predicted geoelectric structure and the alignment status of time-domain-induced signals generated using the five-layer model at different noise levels. For noise-free data, the proposed algorithm achieves precise reconstruction of resistivity, thickness, and induced polarization parameters for the assumed geoelectric model, resulting in a MAPE of 4.04%. For data with 5% and 10% random noise, the resistivity, thickness, and polarization parameters can still broadly characterize the underlying trends, but the precision in capturing details has diminished. Particularly, at a noise level of 10%, the depth estimation for the third layer appears shallower.

5.2. Complex Model

To further validate the network’s capability to handle complex geoelectric structures, a 3D simulation model of a seafloor massive sulfide (SMS) deposit has been established. The resistivity of the overlying layer of minerals is set to 10 Ω·m. As magmatic-hydrothermal upwelling occurs, alteration zones gradually form around the upwelling channels, resulting in electrical resistivities of 0.5, 1, and 2 Ω·m, corresponding to strong, medium, and weak alteration regions, respectively. In the intermediate region between two magmatic-hydrothermal upwelling channels, there exists an alteration zone (depicted in yellow) that does not make direct contact with the seafloor surface alongside a buried, intact sulfide sedimentary mineral deposit (depicted in black). Figure 11 displays a perspective view, top view, and sliced view of this 3D model.
The inversion results of the proposed method are compared with those of the traditional gradient inversion method, as illustrated in Figure 12. The traditional gradient inversion method employs OCCAM inversion [51], and the specific implementation and inversion termination strategy are referenced from the literature. Due to the absence of a strategy for optimizing polarization parameters, the inversion results at y = 70 m and y = 110 m in the gradient inversion show an inability to accurately recover the predefined geoelectric structure of SMS. In contrast, the inversion slices obtained using the proposed method demonstrate an effective characterization of the fundamental geoelectric structural features of magmatic-hydrothermal upwelling channels, alteration regions, and buried SMS minerals. This underscores that the TEM inversion framework based on the GCN-LSTM neural network can effectively account for the induced polarization effect in electrical media while achieving accurate geoelectric structure reconstructions.

6. Conclusions

A strata electrical imaging neural network model based on graph topology is established. This model utilizes GCN to capture non-Euclidean features, such as resistivity-thickness relationships and employs LSTM to capture long-range dependency features. The weakness of GCN in capturing distant node relationships is compensated for by LSTM. The effectiveness of the approach was validated using case studies involving five-strata and nine-strata geoelectric models. Comparative disintegration experiments included results contrasting with traditional MLP, LSTM, and GCN, further illustrating the robustness of the proposed network model. Misfit analysis indicates that, compared to the traditional MLP, the proposed method reduces prediction errors in electrical structure estimation from 10–20% down to 2–3%.

Author Contributions

Conceptualization, R.L. and Y.D.; methodology, Y.D.; software, L.G.; validation, R.L., H.T. and L.G.; investigation, R.L.; writing—original draft preparation, R.L. and Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 42104072 and 52201363.

Data Availability Statement

Data available on request due to restrictions e.g., privacy or ethical.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, S.; Huang, Q.; Zhao, L. A Deep Learning-Based Network for the Simulation of Airborne Electromagnetic Responses. Geophys. J. Int. 2022, 233, 253–263. [Google Scholar] [CrossRef]
  2. Wang, S.; Yu, Y.; Zhao, X.; Wang, Y.; Gao, Q.; Ji, Y. Characterization of Excitation Effects and Data Interpretation of Combined Time-Domain Multiwaveform Transmission Currents. IEEE Sens. J. 2023, 23, 17159–17172. [Google Scholar] [CrossRef]
  3. Li, J.; Farquharson, C.G.; Hu, X. Three Effective Inverse Laplace Transform Algorithms for Computing Time-Domain Electromagnetic responsesInverse Laplace Transform Algorithms. Geophysics 2016, 81, E113–E128. [Google Scholar] [CrossRef]
  4. Cai, H.; Liu, M.; Zhou, J.; Li, J.; Hu, X. Effective 3D-Transient Electromagnetic Inversion Using Finite-Element Method with a Parallel Direct Solver. Geophysics 2022, 87, E377–E392. [Google Scholar] [CrossRef]
  5. Cai, H.; Hu, X.; Xiong, B.; Zhdanov, M.S. Finite-Element Time-Domain Modeling of Electromagnetic Data in General Dispersive Medium Using Adaptive Padé Series. Comput. Geosci. 2017, 109, 194–205. [Google Scholar] [CrossRef]
  6. Maurya, P.K.; Christensen, F.E.; Kass, M.A.; Pedersen, J.B.; Frederiksen, R.R.; Foged, N.; Christiansen, A.V.; Auken, E. Technical Note: Efficient Imaging of Hydrological Units below Lakes and Fjords with a Floating, Transient Electromagnetic (FloaTEM) System. Hydrol. Earth Syst. Sci. 2022, 26, 2813–2827. [Google Scholar] [CrossRef]
  7. Sandersen, P.B.E.; Kallesøe, A.J.; Møller, I.; Høyer, A.-S.; Jørgensen, F.; Pedersen, J.B.; Christiansen, A.V. Utilizing the Towed Transient ElectroMagnetic Method (tTEM) for Achieving Unprecedented near-Surface Detail in Geological Mapping. Eng. Geol. 2021, 288, 106125. [Google Scholar] [CrossRef]
  8. Li, R.; Hu, X.; Xu, D.; Liu, Y.; Yu, N. Characterizing the 3D Hydrogeological Structure of a Debris Landslide Using the Transient Electromagnetic Method. J. Appl. Geophys. 2020, 175, 103991. [Google Scholar] [CrossRef]
  9. Pondthai, P.; Everett, M.E.; Micallef, A.; Weymer, B.A.; Faghih, Z.; Haroon, A.; Jegen, M. 3D Characterization of a Coastal Freshwater Aquifer in SE Malta (Mediterranean Sea) by Time-Domain Electromagnetics. Water 2020, 12, 1566. [Google Scholar] [CrossRef]
  10. Auken, E.; Foged, N.; Larsen, J.J.; Lassen, K.V.T.; Maurya, P.K.; Dath, S.M.; Eiskjær, T.T. tTEM—A Towed Transient Electromagnetic System for Detailed 3D Imaging of the Top 70 m of the Subsurface. Geophysics 2019, 84, E13–E22. [Google Scholar] [CrossRef]
  11. Cheng, M.; Yang, D.; Luo, Q. Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion. Minerals 2023, 13, 34. [Google Scholar] [CrossRef]
  12. Grombacher, D.; Maurya, P.K.; Lind, J.C.; Lane, J.; Auken, E. Rapid Mapping of Hydrological Systems in Tanzania Using a Towed Transient Electromagnetic System. Groundwater 2022, 60, 35–46. [Google Scholar] [CrossRef]
  13. Lin, J.; Chen, J.; Liu, F.; Zhang, Y. The Helicopter Time-Domain Electromagnetic Technology Advances in China. Surv. Geophys. 2021, 42, 585–624. [Google Scholar] [CrossRef]
  14. Maurya, P.K.; Christiansen, A.V.; Pedersen, J.; Auken, E. High Resolution 3D Subsurface Mapping Using a Towed Transient Electromagnetic System - tTEM: Case Studies. Near Surf. Geophys. 2020, 18, 249–259. [Google Scholar] [CrossRef]
  15. Li, R.; Hu, X.; Li, J. Pseudo-3D Constrained Inversion of Transient Electromagnetic Data for a Polarizable SMS Hydrothermal System in the Deep Sea. Stud. Geophys. Et Geod. 2018, 62, 512–533. [Google Scholar] [CrossRef]
  16. Yang, H.; Cai, H.; Liu, M.; Xiong, Y.; Long, Z.; Li, J.; Hu, X. Three-dimensional Inversion of Semi-airborne Transient Electromagnetic Data Based on Finite Element Method. Near Surf. Geophys. 2022, 20, 661–678. [Google Scholar] [CrossRef]
  17. Wu, Q.; Ma, Y.; Li, D.; Wang, Y.; Ji, Y. Denoising Algorithm of Ground-Airborne Time-Domain Electromagnetic Method Based on Variational Bayesian-Based Adaptive Kalman Filter (VBAKF). J. Appl. Geophys. 2022, 202, 104674. [Google Scholar] [CrossRef]
  18. Barsukov, P.O.; Fainberg, E.B. Mapping Bedrock Topography and Moraine Deposits by Transient Electromagnetic Sounding: Oslo Graben, Norway. Near Surf. Geophys. 2020, 18, 123–133. [Google Scholar] [CrossRef]
  19. Bai, P.; Vignoli, G.; Hansen, T.M. 1D Stochastic Inversion of Airborne Time-Domain Electromagnetic Data with Realistic Prior and Accounting for the Forward Modeling Error. Remote Sens. 2021, 13, 3881. [Google Scholar] [CrossRef]
  20. Kong, R.; Hu, X.; Cai, H. Three-dimensional joint inversion of gravity and magnetic data using Gramian constraints and Gauss-Newton method. Chin. J. Geophys. 2023, 66, 3493–3513. [Google Scholar]
  21. Wang, X.; Cai, H.; Liu, L.; Revil, A.; Hu, X. Three-Dimensional Inversion of Long-Offset Transient Electromagnetic Method over Topography. Minerals 2023, 13, 908. [Google Scholar] [CrossRef]
  22. Belliveau, P.; Haber, E. Achieving Depth Resolution with Gradient Array Survey Data through Transient Electromagnetic Inversion. In Proceedings of the SEG International Exposition and Annual Meeting, Dallas, TX, USA, 16–21 October 2016; pp. 857–861. [Google Scholar] [CrossRef]
  23. Murzina, E.V.; Pospeev, A.V.; Seminskiy, I.K.; Buddo, I.V.; Nemtseva, D.B.; Emelianov, V.S.; Agafonov, Y.A. Features of High-Density Transient Electromagnetic Sounding Data Inversion under Oil and Gas Prospecting in the Nepa-Botuoba Anteclise Area. Jour 2022, 45, 137–151. [Google Scholar] [CrossRef]
  24. Li, R.; Zhuang, Q.; Yu, N.; Li, R.; Zhang, H. Improved Hybrid Particle Swarm Optimizer with Sine-Cosine Acceleration Coefficients for Transient Electromagnetic Inversion. Curr. Bioinform. 2021, 16, 60–76. [Google Scholar] [CrossRef]
  25. Li, R.; Yu, N.; Wang, X.; Liu, Y.; Cai, Z.; Wang, E. Model-Based Synthetic Geoelectric Sampling for Magnetotelluric Inversion with Deep Neural Networks. IEEE Trans. Geosci. Remote Sens. 2020, 60, 1–14. [Google Scholar] [CrossRef]
  26. Asif, M.R.; Bording, T.S.; Maurya, P.K.; Zhang, B.; Fiandaca, G.; Grombacher, D.J.; Christiansen, A.V.; Auken, E.; Larsen, J.J. A Neural Network-Based Hybrid Framework for Least-Squares Inversion of Transient Electromagnetic Data. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–10. [Google Scholar] [CrossRef]
  27. Li, R.; Wu, X.; Tian, H.; Yu, N.; Wang, C. Hybrid Memetic Pretrained Factor Analysis-Based Deep Belief Networks for Transient Electromagnetic Inversion. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–20. [Google Scholar] [CrossRef]
  28. Kozhevnikov, N.O.; Antonov, E.Y. Aftereffects in the Transient Electromagnetic Method: Inductively Induced Polarization. Russ. Geol. Geophys. 2021, 62, 1440–1448. [Google Scholar] [CrossRef]
  29. Grombacher, D.; Auken, E.; Foged, N.; Bording, T.; Foley, N.; Doran, P.T.; Mikucki, J.; Dugan, H.A.; Garza-Giron, R.; Myers, K.; et al. Induced Polarization Effects in Airborne Transient Electromagnetic Data Collected in the McMurdo Dry Valleys, Antarctica. Geophys. J. Int. 2021, 226, 1574–1583. [Google Scholar] [CrossRef]
  30. Wu, C.; Zou, C.; Wu, T.; Shen, L.; Zhou, J.; Tao, C. Experimental Study on the Detection of Metal Sulfide under Seafloor Environment Using Time Domain Induced Polarization. Mar. Geophys. Res. 2021, 42, 17. [Google Scholar] [CrossRef]
  31. Liu, H.; Zhao, X.; Yu, Y.; Qiu, S.; Ji, Y. 3-D Full-Waveform Modeling and Analysis of Induced Polarization and Magnetic Viscosity Effect in Time-Domain Electromagnetic Method. IEEE Trans. Geosci. Remote Sens. 2023, 61, 1–14. [Google Scholar] [CrossRef]
  32. He, Z.; Cai, H.; Li, S.; Xian, J.; Hu, X. Extracting IP Parameters of Rock Samples Using Machine Learning. Geophys. J. Int. 2023, 235, 862–878. [Google Scholar] [CrossRef]
  33. Ageev, V.V. On the Nature of Negative Values of Induced Polarization in the Transient Electromagnetic Process When Working with an Axial Electrical Installation on the Sea Shelf. Izv. Phys. Solid Earth 2023, 59, 328–336. [Google Scholar] [CrossRef]
  34. Aigner, L.; Katona, T.; Michel, H.; Ahmed, A.; Hermans, T.; Flores Orozco, A. Quantification of Clay Content Using the Transient Electromagnetic and Spectral Induced Polarization Method. In Proceedings of the EGU General Assembly 2021, Online, 19–30 April 2021. [Google Scholar]
  35. Fedeli, A.; Schenone, V.; Pastorino, M.; Randazzo, A. An LSTM Based Strategy for Data and Model Calibration in Subsurface Electromagnetic Imaging. In Proceedings of the 2022 IEEE Conference on Antenna Measurements and Applications (CAMA), Guangzhou, China, 14–17 December 2022; pp. 1–3. [Google Scholar]
  36. Wu, F.; Fan, M.; Liu, W.; Liang, B.; Zhou, Y. An Efficient Time-Domain Electromagnetic Algorithm Based on LSTM Neural Network. IEEE Antennas Wirel. Propag. Lett. 2021, 20, 1322–1326. [Google Scholar] [CrossRef]
  37. Guo, L.; Song, G.; Wu, H. Complex-Valued Pix2pix—Deep Neural Network for Nonlinear Electromagnetic Inverse Scattering. Electronics 2021, 10, 752. [Google Scholar] [CrossRef]
  38. Hu, K.; Wu, J.; Li, Y.; Lu, M.; Weng, L.; Xia, M. FedGCN: Federated Learning-Based Graph Convolutional Networks for Non-Euclidean Spatial Data. Mathematics 2022, 10, 1000. [Google Scholar] [CrossRef]
  39. Wu, Q.; Zhao, W.; Li, Z.; Wipf, D.P.; Yan, J. NodeFormer: A Scalable Graph Structure Learning Transformer for Node Classification. Adv. Neural Inf. Process. Syst. 2022, 35, 27387–27401. [Google Scholar]
  40. Chen, D.; Lin, Y.; Zhao, G.; Ren, X.; Li, P.; Zhou, J.; Sun, X. Topology-Imbalance Learning for Semi-Supervised Node Classification. In Proceedings of the Advances in Neural Information Processing Systems, Online, 6–14 December 2021; Curran Associates, Inc.: Red Hook, NY, USA, 2021; Volume 34, pp. 29885–29897. [Google Scholar]
  41. Chen, Y.; Coskunuzer, B.; Gel, Y. Topological Relational Learning on Graphs. In Proceedings of the Advances in Neural Information Processing Systems, Online, 6–14 December 2021; Curran Associates, Inc.: Red Hook, NY, USA, 2021; Volume 34, pp. 27029–27042. [Google Scholar]
  42. Chen, J.; He, H.; Wu, F.; Wang, J. Topology-Aware Correlations Between Relations for Inductive Link Prediction in Knowledge Graphs. In Proceedings of the AAAI Conference on Artificial Intelligence, Vancouver, BC, Canada, 2–9 February 2021; Volume 35, pp. 6271–6278. [Google Scholar] [CrossRef]
  43. Behera, R.K.; Jena, M.; Rath, S.K.; Misra, S. Co-LSTM: Convolutional LSTM Model for Sentiment Analysis in Social Big Data. Inf. Process. Manag. 2021, 58, 102435. [Google Scholar] [CrossRef]
  44. Bi, J.; Zhang, X.; Yuan, H.; Zhang, J.; Zhou, M. A Hybrid Prediction Method for Realistic Network Traffic with Temporal Convolutional Network and LSTM. IEEE Trans. Autom. Sci. Eng. 2022, 19, 1869–1879. [Google Scholar] [CrossRef]
  45. Li, J.; Wang, X.; Hu, X.; Yuan, W.; Wang, P. Land-Based TEM Data Processing: From Turn-off Ramp to Full Waveform. Geophysics, 2022; Preprints. [Google Scholar]
  46. Li, R.; Gao, L.; Yu, N.; Li, J.; Liu, Y.; Wang, E.; Feng, X. Memetic Strategy of Particle Swarm Optimization for One-Dimensional Magnetotelluric Inversions. Mathematics 2021, 9, 519. [Google Scholar] [CrossRef]
  47. Nabighian, M.; projecteditor; Corbett, J. Electromagnetic Methods in Applied Geophysics; Society of Exploration Geophysicists: Houston, TX, USA, 1988. [Google Scholar]
  48. Wei, W.; Cai, J.; Hu, X.; Han, Q. An Electrical Conductivity Model for Fractal Porous Media: Fractal and Electrical Conductivity. Geophys. Res. Lett. 2015, 42, 4833–4840. [Google Scholar] [CrossRef]
  49. Madsen, L.M.; Fiandaca, G.; Auken, E.; Christiansen, A.V. Time-Domain Induced Polarization—An Analysis of Cole–Cole Parameter Resolution and Correlation Using Markov Chain Monte Carlo Inversion. Geophys. J. Int. 2017, 211, 1341–1353. [Google Scholar] [CrossRef]
  50. Li, R.; Zhang, H.; Zhuang, Q.; Li, R.; Chen, Y. BP Neural Network and Improved Differential Evolution for Transient Electromagnetic Inversion. Comput. Geosci. 2020, 137, 104434. [Google Scholar] [CrossRef]
  51. Li, R.; Yu, N.; Wang, E.; Wang, X.; Mao, L.; Sun, Z.; Wu, X. Airborne Transient Electromagnetic Simulation: Detecting Geoelectric Structures for HVdc Monopole Operation. IEEE Geosci. Remote Sens. Mag. 2021, 10, 274–288. [Google Scholar] [CrossRef]
  52. Shi, X.; Chen, Z.; Wang, H.; Yeung, D.; Wong, W.; Woo, W. Convolutional LSTM Network: A Machine Learning Approach for Precipitation Nowcasting. In Advances in Neural Information Processing Systems; The MIT Press: Cambridge, MA, USA, 13 June 2015. [Google Scholar]
  53. Greff, K.; Srivastava, R.K.; Koutnik, J.; Steunebrink, B.R.; Schmidhuber, J. LSTM: A Search Space Odyssey. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 2222–2232. [Google Scholar] [CrossRef]
  54. Sherstinsky, A. Fundamentals of Recurrent Neural Network (RNN) and Long Short-Term Memory (LSTM) Network. Phys. D Nonlinear Phenom. 2020, 404, 132306. [Google Scholar] [CrossRef]
  55. Gers, F.A.; Schmidhuber, J.; Cummins, F. Learning to Forget: Continual Prediction with LSTM. Neural Comput. 2000, 12, 2451–2471. [Google Scholar] [CrossRef]
  56. Deepa, N.; Prabakeran, S.; Devi, T. (NMRNN-LSTM)—Novel Modified RNN with Long and Short-Term Memory Unit in Healthcare and Big Data Applications. In Proceedings of the 2022 International Conference on Advancements in Smart, Secure and Intelligent Computing (ASSIC), Bhubaneswar, India, 19–20 November 2022; pp. 1–6. [Google Scholar] [CrossRef]
  57. Zhao, L.; Song, Y.; Zhang, C.; Liu, Y.; Wang, P.; Lin, T.; Deng, M.; Li, H. T-GCN: A Temporal Graph Convolutional Network for Traffic Prediction. IEEE Trans. Intell. Transp. Syst. 2020, 21, 3848–3858. [Google Scholar] [CrossRef]
  58. Hu, J.; Li, X.; Coleman, K.; Schroeder, A.; Ma, N.; Irwin, D.J.; Lee, E.B.; Shinohara, R.T.; Li, M. SpaGCN: Integrating Gene Expression, Spatial Location and Histology to Identify Spatial Domains and Spatially Variable Genes by Graph Convolutional Network. Nat. Methods 2021, 18, 1342–1351. [Google Scholar] [CrossRef] [PubMed]
  59. He, X.; Deng, K.; Wang, X.; Li, Y.; Zhang, Y.; Wang, M. LightGCN: Simplifying and Powering Graph Convolution Network for Recommendation. In Proceedings of the 43rd International ACM SIGIR Conference on Research and Development in Information Retrieval, Virtual, 25–30 July 2020; pp. 639–648. [Google Scholar] [CrossRef]
  60. Li, G.; Knoop, V.L.; Van Lint, H. Multistep Traffic Forecasting by Dynamic Graph Convolution: Interpretations of Real-Time Spatial Correlations. Transp. Res. Part C: Emerg. Technol. 2021, 128, 103185. [Google Scholar] [CrossRef]
  61. Cui, Z.; Henrickson, K.; Ke, R.; Wang, Y. Traffic Graph Convolutional Recurrent Neural Network: A Deep Learning Framework for Network-Scale Traffic Learning and Forecasting. IEEE Trans. Intell. Transport. Syst. 2020, 21, 4883–4894. [Google Scholar] [CrossRef]
Figure 1. TEM inversion results considering IP-induced polarization effect. (a) represents the predicted TEM signal corresponding to the inversion iteration step, and (b) represents the predicted geoelectric structure corresponding to the inversion iteration step.
Figure 1. TEM inversion results considering IP-induced polarization effect. (a) represents the predicted TEM signal corresponding to the inversion iteration step, and (b) represents the predicted geoelectric structure corresponding to the inversion iteration step.
Electronics 12 04138 g001
Figure 2. LSTM network architecture diagram.
Figure 2. LSTM network architecture diagram.
Electronics 12 04138 g002
Figure 3. Assuming the central node (purple) represents the current geoelectric layer. The green nodes depict the connections to the current layer. We acquire spatial features by obtaining the topological relationships between the current geoelectric layer and its surrounding layers.
Figure 3. Assuming the central node (purple) represents the current geoelectric layer. The green nodes depict the connections to the current layer. We acquire spatial features by obtaining the topological relationships between the current geoelectric layer and its surrounding layers.
Electronics 12 04138 g003
Figure 4. GCN-LSTM Model Framework.
Figure 4. GCN-LSTM Model Framework.
Electronics 12 04138 g004
Figure 5. Training processes of the proposed neural networks correspond to different learning rates, network depths, and regularization rates. (a) represents the training loss decay; (b) represents the loss decay of valid dataset.
Figure 5. Training processes of the proposed neural networks correspond to different learning rates, network depths, and regularization rates. (a) represents the training loss decay; (b) represents the loss decay of valid dataset.
Electronics 12 04138 g005
Figure 6. The predictions of different neural network models for the five-strata geoelectric structure based on TEM responses. (a) represents the depth-resistivity relationship chart for the assumed and predicted five-strata geoelectric structures. (b) represents the depth-IP (Induced Polarization) parameter relationship chart for the assumed and predicted five-strata geoelectric structures.
Figure 6. The predictions of different neural network models for the five-strata geoelectric structure based on TEM responses. (a) represents the depth-resistivity relationship chart for the assumed and predicted five-strata geoelectric structures. (b) represents the depth-IP (Induced Polarization) parameter relationship chart for the assumed and predicted five-strata geoelectric structures.
Electronics 12 04138 g006
Figure 7. Misfit analysis of predicting the five-strata geoelectric structure using different neural network models based on TEM responses. (a) Depth and resistivity misfits of the predicted five-strata geoelectric structure, (b) Depth and IP parameter misfits of the predicted five-layer geoelectric structure, (c) Induced magnetic field misfits in the predicted TEM responses, (d) Induced EMF misfits in the predicted TEM responses.
Figure 7. Misfit analysis of predicting the five-strata geoelectric structure using different neural network models based on TEM responses. (a) Depth and resistivity misfits of the predicted five-strata geoelectric structure, (b) Depth and IP parameter misfits of the predicted five-layer geoelectric structure, (c) Induced magnetic field misfits in the predicted TEM responses, (d) Induced EMF misfits in the predicted TEM responses.
Electronics 12 04138 g007
Figure 8. The predictions of different neural network models for the nine-strata geoelectric structure based on TEM responses. (a) represents the depth-resistivity relationship chart for the assumed and predicted nine-strata geoelectric structures. (b) represents the depth-IP (Induced Polarization) parameter relationship chart for the assumed and predicted nine-strata geoelectric structures.
Figure 8. The predictions of different neural network models for the nine-strata geoelectric structure based on TEM responses. (a) represents the depth-resistivity relationship chart for the assumed and predicted nine-strata geoelectric structures. (b) represents the depth-IP (Induced Polarization) parameter relationship chart for the assumed and predicted nine-strata geoelectric structures.
Electronics 12 04138 g008
Figure 9. Misfit analysis of predicting the nine-strata geoelectric structure using different neural network models based on TEM responses. (a) Depth and resistivity misfits of the predicted nine-strata geoelectric structure, (b) Depth and IP parameter misfits of the predicted nine-layer geoelectric structure, (c) Induced magnetic field misfits in the predicted TEM responses, (d) Induced EMF misfits in the predicted TEM responses.
Figure 9. Misfit analysis of predicting the nine-strata geoelectric structure using different neural network models based on TEM responses. (a) Depth and resistivity misfits of the predicted nine-strata geoelectric structure, (b) Depth and IP parameter misfits of the predicted nine-layer geoelectric structure, (c) Induced magnetic field misfits in the predicted TEM responses, (d) Induced EMF misfits in the predicted TEM responses.
Electronics 12 04138 g009
Figure 10. Inversion outcomes of the five-strata geoelectric model with an intermediate polarization layer. (a) represent the theoretical geoelectric model and inversed geoelectric model; (b) represent the predictive polarization parameters and the theoretical parameters; (c) represent the HZ data fitting curve.
Figure 10. Inversion outcomes of the five-strata geoelectric model with an intermediate polarization layer. (a) represent the theoretical geoelectric model and inversed geoelectric model; (b) represent the predictive polarization parameters and the theoretical parameters; (c) represent the HZ data fitting curve.
Electronics 12 04138 g010
Figure 11. A 3D simulation model of a seafloor massive sulfide (SMS) deposit is depicted here. With the upwelling of magmatic-hydrothermal fluids (in red), gradual alteration of the surrounding rock layers occurs, resulting in the formation of alteration zones. These alteration zones, progressing from the inner to outer regions, exhibit varying degrees of alteration—strong, moderate, and weak (in green). In the intermediate space between two magmatic-hydrothermal upwelling channels, there is an alteration zone (in yellow) that does not make direct contact with the seafloor surface, along with a buried intact sulfide sedimentary mineral deposit (in black).
Figure 11. A 3D simulation model of a seafloor massive sulfide (SMS) deposit is depicted here. With the upwelling of magmatic-hydrothermal fluids (in red), gradual alteration of the surrounding rock layers occurs, resulting in the formation of alteration zones. These alteration zones, progressing from the inner to outer regions, exhibit varying degrees of alteration—strong, moderate, and weak (in green). In the intermediate space between two magmatic-hydrothermal upwelling channels, there is an alteration zone (in yellow) that does not make direct contact with the seafloor surface, along with a buried intact sulfide sedimentary mineral deposit (in black).
Electronics 12 04138 g011aElectronics 12 04138 g011b
Figure 12. Comparison between the gradient inversion and the SMS model of the proposed GCN-LSTM networks inversion.
Figure 12. Comparison between the gradient inversion and the SMS model of the proposed GCN-LSTM networks inversion.
Electronics 12 04138 g012
Table 1. Performance of the dataset under different neural network architectures.
Table 1. Performance of the dataset under different neural network architectures.
The Stack DepthLoss
TrainValid
GCN + LSTM151050.080.22
GCN9630.121.86
LSTM64212.943.17
MLP151050.162.04
Table 2. Comparative details of four network models’ predictions for the five-strata geoelectric model parameters.
Table 2. Comparative details of four network models’ predictions for the five-strata geoelectric model parameters.
rho1rho2rho3rho4rho5
Synthetic50300100500200
GCN + LSTM49.72294.0199.01516.67198.07
GCN46.19299.70104.13474.36211.40
LSTM46.39305.93108.02543.94188.85
MLP48.96341.54107.87426.10210.80
h1h2h3h4-
Synthetic200250300500
GCN + LSTM197.80246.51289.21488.01
GCN182.96244.69270.20472.07
LSTM199.31243.80301.43476.49
MLP212.36260.88304.71457.72
c1c2c3c4c5
Synthetic--0.30.3-
GCN + LSTM--0.2930.294-
GCN--0.2790.286-
LSTM--0.2800.272-
MLP--0.2880.335-
Note: rhoi, hi, ci represent the resistivity, thickness, and Cole-Cole model frequency correction coefficient of the i-th layer, respectively.
Table 3. Misfit details of four network models’ predictions for the five-strata geoelectric model parameters.
Table 3. Misfit details of four network models’ predictions for the five-strata geoelectric model parameters.
HzVzrhoThickness (h)c
MAPERMSEMAPERMSEMAPERMSEMAPERMSEMAPERMSE
GCN + LSTM2.8010.313.6012.912.092.482.612.937.841.45
GCN4.846.8014.230.675.866.272.603.2319.103.51
LSTM7.2328.15.83312.944.905.845.546.4425.324.29
MLP19.8672.921.7260.448.179.016.558.1520.543.90
Table 4. Comparative details of four network models’ predictions for the nine-strata geoelectric model parameters.
Table 4. Comparative details of four network models’ predictions for the nine-strata geoelectric model parameters.
rho1rho2rho3rho4rho5rho6rho7rho8rho9
Synthetic5030010050020010030050200
GCN + LSTM49.51294.54101.09520.51208.18100.92294.9851.77198.85
GCN54.36319.12104.57467.58194.4193.78270.0748.16207.98
LSTM49.32311.6695.14450.98201.2995.59326.7754.06195.71
MLP51.42341.5990.57453.96190.50112.99290.1646.60179.12
h1h2h3h4h5h6h7h8-
Synthetic200250300500300100200300-
GCN + LSTM208.09238.33300.97510.82290.3898.37193.75294.66-
GCN205.01252.15296.34478.74300.10105.23210.50304.56-
LSTM180.99258.57320.23547.15273.4299.01203.30311.20-
MLP193.83240.60266.80490.26263.24103.44170.66306.59-
c1c2c3c4c5c6c7c8c9
Synthetic-0.3-0.4-0.2-0.5-
GCN + LSTM-0.30-0.40-0.19-0.48-
GCN-0.31-0.40-0.18-0.47-
LSTM-0.31-0.39-0.18-0.55-
MLP-0.28-0.41-0.20-0.46-
Note: rhoi, hi, ci represent the resistivity, thickness, and Cole-Cole model frequency correction coefficient of the i-th layer, respectively.
Table 5. Misfit details of four network models’ predictions for the nine-strata geoelectric model parameters.
Table 5. Misfit details of four network models’ predictions for the nine-strata geoelectric model parameters.
HzVzrhoThickness (h)c
MAPERMSEMAPERMSEMAPERMSEMAPERMSEMAPERMSE
GCN + LSTM1.131.661.783.611.561.862.122.330.891.41
GCN9.4011.0410.4014.504.535.176.537.192.353.81
LSTM11.9112.9917.0422.516.316.762.002.673.205.14
MLP10.9311.7924.6942.508.7910.055.135.723.125.52
Table 6. Inversion accuracy with noise data.
Table 6. Inversion accuracy with noise data.
HKH-Type ModelRMSE (%)MAE (%)Re-MAE
no noise14.074.044.04
5% noise20.319.914.91
10% noise31.5917.687.68
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, R.; Di, Y.; Tian, H.; Gan, L. TEM Strata Inversion Imaging with IP Effect Based on Enhanced GCN by Extracting Long-Dependency Features. Electronics 2023, 12, 4138. https://doi.org/10.3390/electronics12194138

AMA Style

Li R, Di Y, Tian H, Gan L. TEM Strata Inversion Imaging with IP Effect Based on Enhanced GCN by Extracting Long-Dependency Features. Electronics. 2023; 12(19):4138. https://doi.org/10.3390/electronics12194138

Chicago/Turabian Style

Li, Ruiheng, Yi Di, Hao Tian, and Lu Gan. 2023. "TEM Strata Inversion Imaging with IP Effect Based on Enhanced GCN by Extracting Long-Dependency Features" Electronics 12, no. 19: 4138. https://doi.org/10.3390/electronics12194138

APA Style

Li, R., Di, Y., Tian, H., & Gan, L. (2023). TEM Strata Inversion Imaging with IP Effect Based on Enhanced GCN by Extracting Long-Dependency Features. Electronics, 12(19), 4138. https://doi.org/10.3390/electronics12194138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop