Next Article in Journal
LGTCN: A Spatial–Temporal Traffic Flow Prediction Model Based on Local–Global Feature Fusion Temporal Convolutional Network
Previous Article in Journal
Three-Dimensional Assessment of the Effects of Kinesio Taping on Mandibular Condylar Volume with Skeletal Class II Malocclusions
Previous Article in Special Issue
Testing the Efficacy of Indirect Methods on Characterization of Sedimentary Basins by Correlation of Direct Data and Geophysical Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CSEM Optimization Using the Correspondence Principle

by
Adriany Valente
1,†,
Deivid Nascimento
2,*,† and
Jessé Costa
1,3,†
1
Faculty of Geophysics, Pará Federal University, Av. Augusto Correa, 01, Belém 66070-110, PA, Brazil
2
PETROBRAS, Petrobras Research Center—CENPES, Rio de Janeiro 21941-915, RJ, Brazil
3
INCT-GP/CNPq—National Institute of Science and Technology for Petroleum Geophysics—CNPq, Brasília 70070-010, DF, Brazil
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2024, 14(19), 8846; https://doi.org/10.3390/app14198846
Submission received: 12 July 2024 / Revised: 29 August 2024 / Accepted: 3 September 2024 / Published: 1 October 2024
(This article belongs to the Special Issue Advances in Geophysical Exploration)

Abstract

:
Traditionally, 3D modeling of marine controlled-source electromagnetic (CSEM) data (in the frequency domain) involves high-memory demand, requiring solving a large linear system for each frequency. To address this problem, we propose to solve Maxwell’s equations in a fictitious dielectric medium with time-domain finite-difference methods, with the support of the correspondence principle. As an advantage of this approach, we highlight the possibility of its implementation for execution with GPU accelerators, in addition to multi-frequency data modeling with a single simulation. Furthermore, we explore using the correspondence principle to the inversion of CSEM data by calculating the gradient of the least-squares objective function employing the adjoint-state method to establish the relationship between adjoint fields in a conductive medium and their counterparts in the fictitious dielectric medium, similar to the approach used in forward modeling. We validate this method through 2D inversions of three synthetic CSEM datasets, computed for a simple model consisting of two resistors in a conductive medium, a model adapted from a CSEM modeling and inversion package, and the last one based on a reference model of turbidite reservoirs on the Brazilian continental margin. We also evaluate the differences between the results of inversions using the steepest descent method and our proposed momentum method, comparing them with the limited-memory BFGS (Broyden–Fletcher–Goldfarb–Shanno) algorithm (L-BFGS-B). In all experiments, we use smoothing by model reparameterization as a strategy for regularizing and stabilizing the iterations throughout the inversions. The results indicate that, although it requires more iterations, our modified momentum method produces the best models, which are consistent with results from the L-BFGS-B algorithm and require less storage per iteration.

1. Introduction

The controlled source electromagnetic method (CSEM) is a consolidated method in the off-shore industry [1,2]. It is used not only to reduce the exploratory risk [3], helping to decide well placement, but also in applications aimed at reservoir monitoring [4,5], given the significant role of electrical resistivity in characterizing rock porosity and fluid saturation.
Based on measurements of the electric and magnetic fields made by receivers on the seabed (while a ship towed the submerged electromagnetic source), CSEM allows for obtaining a model of electrical conductivity in the subsurface through an inversion process, which iteratively minimizes the misfit between observed and modeled data.
Several possible approaches exist to perform forward modeling in the context of the inversion of CSEM data [6]. Some of these approaches are based on solving Maxwell’s equations for conductive media directly in the frequency domain and require solving a linear system for each frequency. Furthermore, the Gauss–Newton approximation for the model update is common during inverse problem iterations [7,8], which usually require solving large linear systems, especially in 3D problems. Another critical aspect of CSEM inversion is the regularization strategy applied to stabilize the iterative process [9,10,11,12]. We propose modeling CSEM data through the numerical time-domain finite-difference solution of Maxwell’s equations for a fictitious dielectric medium, following the prescription of [13]. The correspondence principle allows the evaluation of the fields for the corresponding conductive medium in the frequency domain [14]. In the context of inversion, we present the formulation for gradient calculation, obtained using the adjoint-state method [15,16], based on adjoint modeling following the same formulation as the forward modeling. This is the approach presented by [17] but with the correction to the equations that define the mapping between the fictitious and diffusive domains. To stabilize the inverse problem, we use model reparameterization smoothing [18,19], which does not require trial-and-error selection of a regularization parameter.
We also evaluate two optimization algorithms: the steepest descent (using a simple line search strategy) and a momentum method [20], parameterized according to the Adam algorithm [21], known as a low memory demand method. We also compared the results with those obtained using L-BFGS-B [22,23], which is traditionally used in minimization problems. Our numerical experiments show that both algorithms successfully fit the observations and correctly localize the inversion targets, but our momentum approach has low storage demand.

2. Methodology

2.1. Csem Forward Modeling

Maxwell’s equations (in the frequency domain) for a conductive medium at low frequencies are as follows:
× H ( x , ω ) + σ ( x ) E ( x , ω ) = J ( x , ω ) ,
× E ( x , ω ) i ω μ H ( x , ω ) = 0 ,
where E ( x , ω ) and H ( x , ω ) are the electric and magnetic field vectors, respectively, σ ( x ) represents the conductivity tensor, μ is the magnetic permeability, and J ( x , t ) is the electric current density vector, which represents the injection source distribution.
In the same context, Maxwell’s equations in the time domain are given as follows:
× H ( x , t ) + σ ( x ) E ( x , t ) = J ( x , t ) ,
× E ( x , t ) + μ t H ( x , t ) = 0 ,
since we can obtain Equations (1) and (2) from Equations (3) and (4) using the Fourier transform:
F ( ω ) = + d t f ( t ) e i ω t ,
f ( t ) = 1 2 π + d ω F ( ω ) e i ω t .
In a non-conductive dielectric medium, Maxwell’s equations in the time domain are:
× H ( x , t ) + ϵ ( x ) E ( x , t ) t = J ( x , t ) ,
× E ( x , t ) + μ H ( x , t ) t = 0 ,
and in the frequency domain they are the following:
× H ( x , ω ) i ω ϵ ( x ) E ( x , ω ) = J ( x , ω ) ,
× E ( x , ω ) i ω μ H ( x , ω ) = 0 ,
where ϵ ( x ) is the dielectric permittivity tensor.
Defining the fictitious dielectric permittivity,
σ x 2 ω 0 ϵ x ,
we can rewrite Equations (1) and (2):
× H ( x , ω ) i ω ϵ ( x ) E ( x , ω ) = J ( x , ω ) ,
× E ( x , ω ) i ω μ H ( x , ω ) = 0 ,
where
ω ( 1 + i ) ω ω 0 ,
E ( x , ω ) E ( x , ω ) ,
H ( x , ω ) i ω 2 ω 0 H ( x , ω ) ,
J ( x , ω ) i ω 2 ω 0 J ( x , ω ) .
Comparing with Equations (9) and (10), we obtain the equations in the fictitious time domain:
× H ( x , t ) + ϵ ( x ) E ( x , t ) t = J ( x , t ) ,
× E ( x , t ) + μ H ( x , t ) t = 0 .
The mapping defined in Equations (14)–(17) warrants this construction, which establishes a direct relationship between the electromagnetic fields in a conductive medium and their corresponding in a fictitious dielectric model, as established by the correspondence principle [14]. We solve Equations (18) and (19) using finite difference, following the prescriptions of [13], and obtain the components E i of the electric field in the conductive medium using Equation (15):
E i ( x , ω ) = 0 T d t E i ( x , t ) e i ω t = 0 T d t E i ( x , t ) e ω ω 0 t e i ω ω 0 t ,
whose numerical evaluation is stable in consequence of the exponential decay imposed by the term e ω ω 0 t , regardless of the duration of the signal recorded during the simulation in the fictitious time domain.
The injection of the source component J ( x , t ) uses the band-limited spatial Delta function approximation:
J ( x s , t ) = δ ( x x s ) J ( t ) ,
with a signature spectrum given as follows:
J ( ω ) = 2 ω 0 i ω 0 T d t J ( t ) e ω ω 0 t e i ω ω 0 t .
In practice, the choice of the source signature for modeling in the fictitious domain is flexible since the data obtained after performing the CSEM survey corresponds to the field normalized by the source at each frequency.

2.2. Gradient Calculation for CSEM Inversion

The estimate of subsurface resistivity is an ill-posed non-linear problem and inversion of CSEM data usually applies interactive algorithms, where update of the resistivity model is a step along the gradient of the objective function relative to the model parameters. The adjoint state method allows for the efficient computation of the gradient [16].
For a conductive model, σ , the objective function for inversion of CSEM data is the least squares mismatch between the observed and modeled data:
χ ( E , σ ) = s ω r W E ( x r , x s ) E o b s ( x r , ω ; x s ) E ( x r , ω ; x s | σ ) 2 ,
where index s indicates the summation over all sources, index r indicates the summation over the receiver position, and index ω indicates the summation over all recorded frequencies. The weighting factors W E penalize the residuals between the modeled and observed field components [24]. The observed electric field at the receiver position, denoted by x r , due to a source located at x s , for a given frequency, ω , is represented by the vector E o b s ( x r , x s , ω ) . The field vector E ( x r , x s , ω ) represents the corresponding modeled field for a given conductivity model, σ .
We derived the gradient of χ ( E , σ ) relative to the conductivity model through the adjoint state method, which is an extension of the Lagrange multipliers method. The starting point is the extended Lagrangian, containing the sum of the least-squares data misfit plus the state equations that enforce the modeled fields and medium properties should honor the Maxwell equations, i.e.,
L ( E , H , σ , λ E , λ H ) = s ω r W E ( x r , ω ; x s ) E o b s ( x r , ω ; x s ) E ( x r , ω ; x s | σ ) 2 + s ω Ω d Ω λ E * ( u ) × H ( u ) + σ ( x ) E + J ( u ) + s ω Ω d Ω λ H * ( u ) × E ( u ) i ω μ H ( u )
where ( z ) indicates the real part of its complex argument and u ( x , ω ; x s ) . λ E ( x , ω ; x s ) , and λ H ( x , ω ; x s ) are the adjoint-state fields.
The stationarity of the first variation of the Lagrangian functional relative to the electromagnetic fields produces the system of equations for the adjoint-state fields [15,16], λ E and λ H ,
× λ H * ( u ) + σ ( x ) λ E * ( u ) = r δ ( x x r ) W E ( u ) E o b s ( u ) E ( u ) ,
× λ E * ( u ) i ω μ λ H * ( u ) = 0 ,
and the gradient of the objective function relative to the conductivity tensor,
χ σ x = s ω λ E * ( x , ω ; x s ) E ( x , ω ; x s ) .
Therefore, to evaluate the gradient of the objective function, we need to solve the direct electromagnetic fields equation’s and the adjoint-state fields equation’s.
To solve the system Equations (25) and (26), we apply the correspondence principle proposed by [13].
Our derivation demonstrates that the corresponding adjoint-field equations in the fictitious dielectric medium in the frequency domain are as follows:
× λ H ( x , ω ) + ( i ω ) ϵ ( x ) λ E ( x , ω ) = r δ ( x x r ) J ( x , ω ) ,
× λ E ( x , ω ) + ( i ω ) μ λ H ( x , ω ) = 0 .
The relationship between the adjoint-state equations in a conductive medium and the corresponding adjoint-state equations in a fictitious dielectric medium are:
ω ( 1 + i ) ω ω 0 ,
ϵ ( x ) σ ( x ) / ( 2 ω 0 ) ,
λ E ( x , ω ) λ E * ( x , ω ) ,
λ H ( x , ω ) i ω 2 ω 0 λ H * ( x , ω ) ,
J ( x , ω ) i ω 2 ω 0 Δ E * ( x , ω ) .
The Equations from (30) to (34) map the time-domain adjoint-state fields to the corresponding frequency-domain adjoint-state fields in the conductive medium. The system of Equations (28) and (29) in the time domain is:
× λ H ( x , t ) + ϵ ( x ) t λ E ( x , t ) = r δ ( x x r ) J ( x , t ) ,
× λ E ( x , t ) + μ t λ H ( x , t ) = 0 .
To solve these equations through finite-difference, we use the algorithm proposed by [13]. The current density vector at each receiver position is as follows (34):
J ( x r , ω ) = 0 T d t J ( x r , t ) e ω ω 0 t e i ω ω 0 t .
Finally, the adjoint-field λ E * ( x , ω ) is evaluated on the fly at each finite-difference time update, applying the correspondence relation in Equation (32):
λ E * ( x , ω ) = 0 T λ E ( x , t ) e ω ω 0 t e i ω ω 0 t d t .

2.3. Adjoint Source

The source injected at a receiver during adjoint modeling in the fictitious time domain must satisfy the following equation:
Δ E * ( x r , ω ) = 2 ω 0 i ω 0 T d t J ( x r , t ) e ω ω 0 t e i ω ω 0 t .
Obtaining the adjoint source (in fictitious time domain) from Δ E * ( x r , ω ) is an underdetermined problem, which allows a high degree of freedom in the choice of J ( x r , t ) [13]. For example, one can assume an estimated pulse, J ^ ( x r , t ) , formed from the sum of pulses corresponding to the second derivative of a Gaussian:
J ^ ( x r , t ) = m = 1 N ω A m t 2 Γ ( t τ m ) ,
Γ ( t τ m ) = β m π e β m ( t τ m ) 2 ,
where N ω is the number of frequencies, and β m = π 3 / τ m 2 .
With this choice, the adjoint source depends on the 2 N ω parameters, A m and τ m , determined by minimizing the functional defined by the difference between the spectra of the residuals and the sources (calculated from Equation (40)) in the least-squares sense.
Following another strategy, we chose to estimate the adjoint source by solving the linear system corresponding to the discrete form of Equation (39), with equations given as follows:
Δ E * ( x r , ω q ) = n = 1 N t A q n J ( x r , t n ) ,
A q n = Δ t 2 ω 0 i ω q e ω q ω 0 t n e i ω q ω 0 t n ,
where q is the index associated with the frequencies, i.e., q { 1 , , N ω } , and N t is the number of time samples. To stabilize the system solution, we added the minimum norm regularization with coefficient λ [25], and we also used preconditioning via reparameterization:
J = P m ,
where the preconditioning operator P has the form of a Ricker pulse with peak frequency f p , expressed as follows:
P n k = 1 2 π 2 f p 2 t n t k 2 e π 2 f p 2 t n t k 2 .
Finally, we obtain the adjoint source (in the fictitious time domain) by solving the linear system:
AP λ I m = Δ E * 0 .
We obtain the solution through the conjugate gradient algorithm [19].

2.4. Gradient Preconditioning

The gradient calculated according to Equation (27) is subjected to a preconditioning process [26], used to compensate for its attenuation with depth. We implement two different approaches. The first approach, preconditioning 1, involves calculating the ratio of the gradient to the squared modulus of the electric field in each model sample, requiring calibration to avoid amplifying the gradient at the edges. In practice, at each frequency, the division is performed by the sum of the normalized field with a constant:
g p r e c o n d = g E max ( E ) + δ ,
where E = E x 2 + E y 2 + E z 2 , and δ can be chosen by evaluating the square modulus of the field in the region of interest.
The second strategy (referred to as preconditioning 2) involves essentially applying a depth gain factor by multiplying the gradient by an exponential [27]:
g p r e c o n d = g × e k z / δ ,
where k is a dimensionless parameter (close to 1.0 ), z is the depth related to the seabed, and δ is the skin-depth [28] at the frequency of the gradient calculation.

2.5. Model Regularization

We used the model reparameterization method proposed by [18] to promote the model’s smoothness. In this case, we assume the following:
σ = S 2 S 1 p ,
where S 1 and S 2 are exponential filters [29] that depend on the size of the previously defined smoothing window and operate in the directions where the model parameters vary.
We can obtain the gradient relative to the parameter vector, p ( x ) , using the chain rule:
χ p = χ σ σ p = S 2 S 1 χ σ .

2.6. Steepest Descent and Momentum Methods

CSEM inversion aims to obtain a subsurface conductivity model representative of the investigated region and honors the measurements of electromagnetic fields collected by receivers on the seabed. It is an iterative process that seeks to minimize the objective function defined in Equation (23) through successive steps in the parameter space, using the direction opposite to the gradient of the objective function calculated from the model obtained in each iteration. This update form is the main feature of steepest descent methods [20], which leads to the following expression:
σ k + 1 = σ k + α d k ,
where α is the step size, determined by a simple backtracking line search, and the update direction, d k , is given by the preconditioned gradient, g k , of χ ( σ ) relative to the model conductivity:
d k = g k g k .
We use the most straightforward line search strategy, starting with a step size equal to 1.0 in Equation (51). We calculate the value of the objective function using the updated model, and if there is a decrease, we proceed to the next iteration. Otherwise, we discard the update and compute a new model with a step size that is half of the previous step. This process iterates until the objective function decreases or we reach the maximum number of iterations allowed, in which case the inversion stops. In our experiments, this maximum number is 20.
Another inversion strategy used involves updating the model with a momentum vector, s , defined from the gradient of the objective function and the momentum calculated in the previous iteration as follows:
σ k + 1 = σ k α s ^ k ,
s ^ = 1 β 2 k 1 β 1 k s k ; s k = β 1 s k 1 + 1 β 1 d k ,
where we adopt the parameters β 1 = 0.9 e β 2 = 0.999 , according to the usual prescription of the Adam algorithm [21,30], formulated initially with a second momentum (which we do not use here). Indeed, our momentum update is very similar to the formulation of the Adam algorithm, where we essentially remove the term involving the second momentum from the equation. Furthermore, although it is not common practice, we combine momentum with the line search strategy already adopted in the context of steepest descent in an attempt to obtain more accurate models.

3. Numerical Experiments

3.1. Gradient Calculation

We validate our implementation to calculate the gradient of the objective function using the data from a 3D isotropic model measuring 18 km in the x direction, 10 km in the y direction, and 3 km deep (in the z direction), inspired by [31]. This model simulates a 2D distribution of conductivities consisting of a resistive body with a resistivity of 25 ohm.m, in a homogeneous background medium with 1 ohm.m, and a water layer of 800 m, with its central section shown in Figure 1. We also use a uniform grid with cell dimensions of 100 m (vertical) × 200 m × 200 m (in this work, we use a 3D modeling code, and to ensure compatibility with the 2D inversion codes, we adopt the strategy of replicating the models obtained during the optimization process along the third dimension).
With the model illustrated in Figure 1, we simulate, using reciprocity, the survey along its central section to use the horizontal component of the electric field (at frequencies of 0.25 Hz, 0.50 Hz, and 0.75 Hz) measured by the receivers on the seafloor with a regular spacing of 1 km, from coordinates 6 km to 12 km. The sources, located 100 m above the receivers, are distributed at intervals of 200 m, between 1.1 km and 16.9 km. Figure 2 illustrates the data associated with the central receiver at the three frequencies (since we use reciprocity, which allows the interchange between sources and receivers, the data illustrated in Figure 2 correspond to the field measured at the source positions for a source attributed to the position of the central receiver).
The estimation of the adjoint source involves solving the system of Equation (46) for the residual (right-hand side of the equation) at each offset, such that at the end of the process, the conjugate gradient routine provides a pulse in the (fictitious) time domain (vector m ) as the solution. An important aspect of the approach we adopted is related to the need to fit the residuals at each frequency separately, in contrast to the more intuitive strategy of simultaneously fitting the residuals at the frequencies where the gradient is to be calculated, as we can see in [17]. The adoption of this approach is essential to obtain the results that will be presented below.
Next, we analyze the accuracy of our adjoint source estimation, considering the residuals between the observed data and the data calculated with a model consisting only of the 1 ohm.m background present in the model shown in Figure 1. Figure 3 illustrates the estimated pulses (at each frequency) for source 60 and the data from receiver number 4 and their amplitude spectra, along which we indicate the residuals (in the frequency domain) adjusted during obtaining adjoint sources.
We also illustrate the gradient calculation relative to the vertical component of the conductivity tensor based on data for receiver 1 at a frequency of 0.50 Hz. Figure 4 shows the absolute value (in log scale) of the E z component of the electric field normalized by the source and the field λ E z , computed using forward and adjoint modeling, respectively, as well as their product, calculated according to Equation (27). To improve the visualization, we apply muting in a depth range corresponding to five samples below the seafloor.
The “complete” gradient is obtained by summing the gradients calculated from the data of each receiver. Figure 5 shows the results for frequencies of 0.25, 0.50, and 0.75 Hz, which have similar characteristics in shape but different depth ranges. Figure 6 illustrates the effect of the preconditioning strategies discussed in Section 2.4, where we show the result of their application to the gradient obtained with the 0.25 Hz data. As we can see, precondition 1 (performed with δ = 10 4 ) positions the anomaly deeper while precondition 2 (performed with k = 1.0 ) increases the depth range of the anomaly.

3.2. Inversion Results

In this section, we will discuss the inversion results of the data obtained under three different scenarios. The first is similar to the model used in Section 3.1, consisting of two resistors located at varying depths in a homogeneous background medium. The second is an adaptation of one of the demonstration models from the MARE2DEM modeling and inversion package (available at http://mare2dem.ucsd.edu acessed on 11 December 2023) [32], and the third corresponds to an adaptation of the synthetic model of the Marlim field (on the Brazilian continental margin) [33], constructed from seismic sections and information from the main horizons in the area. The optimizers used will be the steepest descent and the momentum method described in Section 2.6, as well as L-BFGS-B [22] (computed using a limited-memory matrix using the previous five gradients).

3.2.1. Model with Two Resistors

The model with two resistors has the same background as in Section 3.1, with the same resistivity and water layer thickness. The difference is that there are two resistors (not just one) located at depths of 1.8 km and 2.1 km (Figure 7a). The data correspond to the simulated horizontal electric field measured between 3 and 15 km, with sources every 200 m, in the interval of 1.1 km to 16.9 km. The data from each receiver were restricted to offsets between 1 and 6 km. Figure 7b illustrates the gradient at the frequency of 0.25 Hz, relative to the vertical conductivity, computed in a homogeneous medium with a resistivity of 1 ohm.m. As we can see, there is consistency in the regions where the gradient is prominent, considering the position of the resistors in the model used to generate the observed data. The higher value in the left portion is a consequence of the shallower depth of the resistor located in that region.
The inversions used data at a frequency of 0.25 Hz on a mesh with a grid interval of 100 m (along the vertical) and 200 m (along the horizontal) and a homogeneous initial model with a resistivity of 1 ohm.m. Below, we present the best results obtained with each tested method. In inversions with the steepest descent and momentum, we used preconditioning 1 ( δ = 10 4 ), and in inversion with L-BFGS-B, we used preconditioning 2 ( k = 1.0 ). Figure 8 shows the curves of the objective function’s evolution during the inversions, highlighting the greater efficiency of the steepest descent method, which takes just over 20 iterations to reach a value comparable to that obtained by the momentum inversion after 60 iterations. The L-BFGS-B method is similar to the steepest descent, although it takes more iterations to reach the minimum of the objective function value.
The models obtained from the inversions are similar (Figure 9). They are indicated by the two resistive bodies at different depths, consistent with the resistors’ positions in the model used to produce the observed data. In the cases of the steepest descent and the momentum method, we observe that the deeper resistor is slightly above the desired position, unlike the model obtained with the L-BFGS-B method, in which this resistor is better positioned.

3.2.2. MARE Model

The model we will call MARE consists of an adaptation of the isotropic conductivity model available in one of the versions of the MARE2DEM package for modeling and inversion (the MARE2DEM package [32] is a CSEM 2.5D data modeling and inversion software suite based on finite elements. The code was initially developed as part of the SEMC consortium (Seafloor Electromagnetic Methods Consortium) at Scripps Institution of Oceanography, University of California, San Diego, USA). It simulates the presence of two bodies: a smaller one with a resistivity of 50 ohm.m and a larger one, to the left of the first, of 4 ohm.m. These resistors are embedded in a medium formed by conductive layers and an extended structure with 100 ohm.m, located in the lower portion of the model. The adaptation involved replacing the bathymetry present in the original model with a horizontal seabed at a depth of 1.2 km, as we can see in Figure 10.
The acquisition used in the simulation of observed data consisted of 25 receivers positioned on the seabed, regularly spaced between 1 and 25 km, and sources located 50 m above the receivers, from coordinate 0 to coordinate 26 km, with sampling every 200 m. We computed the synthetic data using the MARE2DEM modeling code [34]. The field used corresponds to the inline electric field. The receiver data were limited to offsets greater than 1 km and those with field amplitude greater than 10 15 V/m. The inversions used 0.25 Hz data.
To perform the inversions, we limited the vertical extension of the model to a depth of 6 km on a grid with a spacing of 100 m (vertical) by 200 m (horizontal), as seen in the initial model (shown in Figure 11a), constructed by smoothing the background formed by the conductive layers present in the reference model. Figure 11b illustrates the gradient relative to vertical conductivity at a frequency of 0.25 Hz, calculated with the model presented in Figure 11a. There is a strong correlation between the regions where the gradient stands out and the locations where the resistors are present in the reference model that produced the observed data.
Next, we present the best results obtained with each methodology presented. In the case of the steepest descent with line search (Figure 12a) and the inversion with momentum (Figure 12b), we used the gradient with precondition 2. In contrast, in the inversion with L-BFGS-B (Figure 12c), the best result was obtained using precondition 1. Figure 13 shows the convergence curves for the objective function. As we can see, the steepest descent and L-BFGS-B present the best convergence behavior compared to the inversion with the momentum algorithm. The estimated models are very similar, especially concerning the smaller and more resistive body on the right, which is well-defined by the three inversions (with resistivity values slightly closer to the reference model in the results obtained with momentum and L-BFGS-B). Regarding the larger resistor on the left, it is possible to notice that there is an overestimation of resistivity in the three results, with a slightly pronounced effect in the result of the inversion with the steepest descent and L-BFGS-B.

3.2.3. Marlim Model

The model we will refer to as Marlim is an adaptation of the MR3D resistivity model constructed by [33] from 3D seismic and well-log data from the Marlim field (Brazilian continental margin). The features of interest in using electromagnetic (EM) methods include a heterogeneous reservoir embedded in a conductive background and a salt body with a resistivity of 1000 ohm.m. It is a good representation of the type of scenario where the use of the CSEM method might be required, which was also described by [35] when discussing the forward modeling performed from the MR3D model.
The modifications made relative to the original MR3D model consist of removing the effects of bathymetry and using only the vertical resistivity of the original model (in this case, we considered the model to be isotropic). Figure 14 illustrates the central section of the model (Figure 14a) from which we simulated the observed data through modeling that used a grid of 20 m (vertical) by 100 m (horizontal) and recorded the inline electric field measured by 20 receivers regularly arranged along the line that coincides with the section in Figure 14a (in the interval from 11 km to 30 km). The sources are distributed between 0 to 41 km every 200 m. We discarded data with offsets smaller than 1 km and amplitudes smaller than 10 15 V/m, and the inversions used single-frequency data at 0.25 Hz. The initial model is homogeneous, with a resistivity of 2.8 ohm.m (which corresponds to the average background resistivity in the region where the reservoir is in the Marlim model), the same used in calculating the gradient relative to vertical conductivity (Figure 14b) on a mesh with a grid spacing of 100 m (along the vertical) and 200 m (along the horizontal). Note the asymmetry in the gradient shape associated with the location of the resistor in the reference model.
Using the data and the initial model described, the inversions performed produced the results that will be illustrated below. Using the steepest descent method with gradient precondition 1, we obtained the model shown in Figure 15a. In the inversion with momentum, where we used precondition 2, we obtained the model in Figure 15b, and the inversion with L-BFGS-B (performed using the gradient with precondition 1), resulting in the model shown in Figure 15c. All results indicate the presence of a resistor in the region where the reservoir is located in the Marlim model, each of them presenting its own characteristics. The models obtained with the steepest descent and L-BFGS-B are similar, as they are characterized by an anomaly with higher resistivity in the case of the L-BFGS-B inversion. The result obtained with momentum is better distributed, covering a larger reservoir extent. The curves showing the evolution of the objective function during the inversions are shown in Figure 16, indicating a slightly more effective reduction rate in the case of L-BFGS-B, especially in the first iterations. In contrast, the others show similar behavior, with the momentum method reaching a slightly lower residual value.

4. Discussion

The results illustrated in the previous section consolidate the use of the adjoint-state method and the correspondence principle to obtain the gradient (relative to the model conductivities) of the objective function defined in Equation (23). The fictitious domain modeling strategy, already established in the context of solving the forward problem, also proved efficient in obtaining the adjoint fields. For the calculation of the gradient consistent with the perturbation observed in the data due to the presence of a resistor in the model, we opted for the approach based on obtaining the adjoint sources for each frequency separately, unlike other authors who address the topic, as seen in the work of [36].
We validated the quality of the gradient calculated with the proposed methodology, preconditioning, and regularization using model parameterization through the inversions of simulated data from three different models: the model consisting of two resistors in a conductive background, the MARE model, an adaptation of the example model from the MARE2DEM package, and the Marlim model, adapted from the MR3D model, constructed by [33]. We used three optimization methods: steepest descent, our proposed modified momentum algorithm, and L-BFGS-B.
Regardless of the dataset we used and the inversion methodology employed, we observed interesting aspects in the inversion results, which produced models with well-positioned resistors at depth and without artifacts near the receiver positions, commonly seen in CSEM inversion results. We attribute these achievements to the use of gradient preconditioning and smoothing by reparameterization.
For the model with two resistors, all three methods successfully obtained results that included the bodies at different depths, with the positioning of the deeper resistor being slightly better in the case of L-BFGS-B (Figure 9). From the perspective of the number of iterations needed for the objective function to reach its minimum value, L-BFGS-B was slightly slower than the steepest descent, with the momentum inversion requiring the most iterations (Figure 8).
The inversions for the MARE model were obtained using synthetic data produced with the MARE2DEM package, and the analysis of the objective function evolution curves during the inversions indicated slower convergence for the momentum optimization, which reached its lowest level after iteration 30, in contrast to the performance of steepest descent and L-BFGS-B, which did so at around iterations 20 and 15, respectively (Figure 13). The results obtained were similar but presented a subtle difference in the resistivity of the larger resistor, which was slightly overestimated in some cases (Figure 12). This feature was somewhat more pronounced in the model obtained with the steepest descent and less noticeable in the model obtained with momentum.
In the inversions for the Marlim model, we observed similar behavior in the convergence curve of the objective function during the first iterations with the steepest descent and momentum; the L-BFGS-B showed a sharper reduction followed by stagnation that lasted for about ten iterations, after which the objective function began to decrease again (Figure 16). The models obtained with the three optimizers differed (Figure 15). While the steepest descent result indicated the presence of a modest and localized resistor, the result obtained with momentum was more consistent with the overall geometry of the reservoir. The result produced with L-BFGS-B was somewhat similar to that obtained with the steepest descent but slightly better as it identified a larger and more resistive resistor.
The analysis of the results produced with the different optimization methodologies indicates that the momentum method tends to fit the observed data slightly better compared to the steepest descent method (as observed in the experiments conducted with the data from the two-resistor model and the MARE model). In the inversion for the Marlim model, this approach not only fit the data better but also recovered a model that was more consistent with the reference model.

5. Conclusions

In this work, we implemented CSEM inversion without solving large linear systems, which can be extremely useful in the context of 3D problems, where this task might become prohibitive. Our implementation used the correspondence principle to derive the forward modeling and gradient calculation. To improve the conditioning of the inversion, we used gradient preconditioning and model reparameterization smoothing. As a result, we obtained very satisfactory models from the inversion of data produced from three different scenarios, using the steepest descent and the momentum algorithm and comparing them to the results from L-BFGS-B. In the two-resistor and the MARE models, the superior performance of the steepest descent and L-BFGS-B was evident, as they required fewer iterations to achieve results similar to those obtained by the momentum method (which needed more iterations). In the case of the Marlim model, the performances were comparable, with better models achieved by the momentum method and L-BFGS-B. Our proposal for using momentum is based on its update from two parameters, similar to what is done in the formulation of the Adam algorithm. However, unlike this one, we used only one momentum term. As we observed, the results from the algorithm we proposed are comparable to the models produced with L-BFGS-B, with the advantage of requiring less memory consumption and fitting the observed data slightly better.

Author Contributions

Methodology, A.V., D.N. and J.C.; validation, A.V., D.N. and J.C.; numerical experiments, A.V., D.N. and J.C.; writing—review and editing, A.V., D.N. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-Brasil (CAPES)-grant number 88882.445206/2019-01 (A.V.), by the Conselho Nacional do Desenvolvimento Científico e Tecnológico-CNPq-grant number 313011/2022-2 and the Instituto Nacional de Ciência e Tecnologia-Geofísica do Petróleo-INCT-GP (J.C.) and by PETROBRAS (D.N.).

Data Availability Statement

The data used in this research can be requested directly from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Constable, S. Ten years of marine CSEM for hydrocarbon exploration. Geophysics 2010, 75, 75A67–75A81. [Google Scholar] [CrossRef]
  2. Menezes, P.T.L.; Ferreira, S.M.; Correa, J.L.; Menor, E.N. Twenty Years of CSEM Exploration in the Brazilian Continental Margin. Minerals 2023, 13, 870. [Google Scholar] [CrossRef]
  3. Buonora, M.; Correa, J.; Martins, L.; Menezes, P.; Pinho, E.; Crepaldi, J.; Ribas, M.; Ferreira, S.; Freitas, R. mCSEM data interpretation for hydrocarbon exploration: A fast interpretation workflow for drilling decision. Interpretation 2014, 2, SH1–SH11. [Google Scholar] [CrossRef]
  4. Menezes, P.T.L.; Correa, J.L.; Alvim, L.M.; Viana, A.R.; Sansonowski, R.C. Time-Lapse CSEM Monitoring: Correlating the Anomalous Transverse Resistance with SoPhiH Maps. Energies 2021, 14, 7159. [Google Scholar] [CrossRef]
  5. Ettayebi, M.; Wang, S.; Landrø, M. Time-Lapse 3D CSEM for Reservoir Monitoring Based on Rock Physics Simulation of the Wisting Oil Field Offshore Norway. Sensors 2023, 23, 7197. [Google Scholar] [CrossRef]
  6. Alumbaugh, D.L.; Newman, G.A.; Prevost, L.; Shadid, J.N. Three-dimensional wideband electromagnetic modeling on massively parallel computers. Radio Sci. 1996, 31, 1–23. [Google Scholar] [CrossRef]
  7. da Piedade, A.A.; Régis, C.; Nunes, C.M.B.; da Silva, H.F. Computational cost comparison between nodal and vector finite elements in the modeling of controlled source electromagnetic data using a direct solver. Comput. Geosci. 2021, 156, 104901. [Google Scholar] [CrossRef]
  8. Crepaldi, J.L.S.; Buonora, M.P.P.; Figueiredo, I. Fast marine CSEM inversion in the CMP domain using analytical derivatives. Geophysics 2011, 76, F303–F313. [Google Scholar] [CrossRef]
  9. Hansen, K.; Panzner, M.; Shantsev, D.; Mittet, R. TTI inversion of marine CSEM data. In SEG Technical Program Expanded Abstracts 2016; Society of Exploration Geophysicists: Houston, TX, USA, 2016; pp. 1014–1018. [Google Scholar]
  10. Meju, M.A.; Mackie, R.L.; Miorelli, F.; Saleh, A.S.; Miller, R.V. Structurally tailored 3D anisotropic controlled-source electromagnetic resistivity inversion with cross-gradient criterion and simultaneous model calibration. Geophysics 2019, 84, E387–E402. [Google Scholar] [CrossRef]
  11. Hoversten, G.M.; Mackie, R.L.; Hua, Y. Reexamination of controlled-source electromagnetic inversion at the Lona prospect, Orphan Basin, Canada. Geophysics 2021, 86, E157–E170. [Google Scholar] [CrossRef]
  12. Cai, H.; Long, Z.; Lin, W.; Li, J.; Lin, P.; Hu, X. 3D multinary inversion of controlled-source electromagnetic data based on the finite-element method with unstructured mesh. Geophysics 2021, 86, E77–E92. [Google Scholar] [CrossRef]
  13. Mittet, R. High-order finite-difference simulations of marine CSEM surveys using a correspondence principle for wave and diffusion fields FDTD simulation of marine CSEM surveys. Geophysics 2010, 75, F33–F50. [Google Scholar] [CrossRef]
  14. de Hoop, A.T. A general correspondence principle for time-domain electromagnetic wave and diffusion fields. Geophys. J. Int. 1996, 127, 757–761. [Google Scholar] [CrossRef]
  15. Plessix, R.E. A review of the adjoint-state method for computing the gradient of a functional with geophysical applications. Geophys. J. Int. 2006, 167, 495–503. [Google Scholar] [CrossRef]
  16. Chavent, G. Nonlinear Least Squares for Inverse Problems—Theoretical Foundations and Step-by-Step Guide for Applications; Number 1 in Scientific Computation; Springer: Dordrecht, The Netherlands, 2010; p. 360. [Google Scholar] [CrossRef]
  17. Valente, A.; Nascimento, D.; Costa, J.C. CSEM inversion using the correspondence principle: The adjoint-state approach. In Proceedings of the 17th International Congress of the Geophysical Society, Online, 8–11 November 2021; Sociedade Brasileira de Geofísica: Rio de Janeiro, Brazil, 2021. [Google Scholar]
  18. Harlan, W.S. Regularization by model reparameterization. Citeseer 1995.
  19. Claerbout, J.F.; Fomel, S. Image Estimation by Example: Geophysical Soundings Image Construction—Multidimensional Autoregression; Stanford University: Stanford, CA, USA, 2003. [Google Scholar]
  20. Kochenderfer, M.J.; Wheeler, T.A. Algorithms for Optimization; Massachusetts Institute of Technology: Cambridge, MA, USA, 2019. [Google Scholar]
  21. Kingma, D.P.; Ba, J. Adam: A method for stochastic optimization. arXiv 2014, arXiv:1412.6980. [Google Scholar]
  22. Zhu, C.; Byrd, R.H.; Lu, P.; Nocedal, J. Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization. ACM Trans. Math. Softw. (Toms) 1997, 23, 550–560. [Google Scholar] [CrossRef]
  23. Nguyen, A.K.; Nordskag, J.I.; Wiik, T.; Bjørke, A.K.; Boman, L.; Pedersen, O.M.; Ribaudo, J.; Mittet, R. Comparing large-scale 3D Gauss–Newton and BFGS CSEM inversions. In SEG Technical Program Expanded Abstracts 2016; Society of Exploration Geophysicists: Houston, TX, USA, 2016; pp. 872–877. [Google Scholar]
  24. Plessix, R.E.; Mulder, W. Resistivity imaging with controlled-source electromagnetic data: Depth and data weighting. Inverse Probl. 2008, 24, 034012. [Google Scholar] [CrossRef]
  25. Menke, W. Geophysical Data Analysis: Discrete Inverse Theory; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  26. Shin, C.; Jang, S.; Min, D.J. Improved amplitude preservation for prestack depth migration by inverse scattering theory. Geophys. Prospect. 2001, 49, 592–606. [Google Scholar] [CrossRef]
  27. Plessix, R.; van der Sman, P. 3D CSEM modeling and inversion in complex geologic settings. In SEG Technical Program Expanded Abstracts 2007; Society of Exploration Geophysicists: Houston, TX, USA, 2007; pp. 589–593. [Google Scholar] [CrossRef]
  28. Chave, A.D.; Constable, S.C.; Edwards, R.N. Electrical Exploration Methods for the Seafloor. In Electomagnetic Methods in Applied Geophysics, Vol. 2, Application; Nabighian, M.N., Ed.; SEG: Houston, TX, USA, 1987; Volume 2, Investigations in Geophysics, Chapter 12; pp. 931–966. [Google Scholar]
  29. Hale, D. My Favorite Ten-Line Computer Program. 2012. Available online: https://inside.mines.edu/~dhale/notebook.html (accessed on 30 May 2021).
  30. Liu, L.; Yang, B.; Zhang, Y.; Xu, Y.; Peng, Z.; Wang, F. Solving Electromagnetic Inverse Problem Using Adaptive Gradient Descent Algorithm. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5902415. [Google Scholar] [CrossRef]
  31. Støren, T.; Zach, J.; Maaø, F. Gradient calculations for 3D inversion of CSEM data using a fast finite-difference time-domain modelling code. In Proceedings of the 70th EAGE Conference and Exhibition incorporating SPE EUROPEC 2008, Rome, Italy, 9–12 June 2008. [Google Scholar]
  32. Key, K. MARE2DEM: A 2-D inversion code for controlled-source electromagnetic and magnetotelluric data. Geophys. J. Int. 2016, 207, 571–588. [Google Scholar] [CrossRef]
  33. Carvalho, B.R.; Menezes, P.T.L. Marlim R3D: A realistic model for CSEM simulations-phase I: Model building. Braz. J. Geol. 2017, 47, 633–644. [Google Scholar] [CrossRef]
  34. Key, K.; Ovall, J. A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modelling. Geophys. J. Int. 2011, 186, 137–154. [Google Scholar] [CrossRef]
  35. Correa, J.L.; Menezes, P.T.L. MR3D: A realistic model for CSEM simulations—Phase II: The CSEM dataset. In SEG Technical Program Expanded Abstracts 2018; SEG: Houston, TX, USA, 2018; pp. 959–963. [Google Scholar] [CrossRef]
  36. Yang, P. 3D fictitious wave domain CSEM inversion by adjoint source estimation. Comput. Geosci. 2023, 180, 105441. [Google Scholar] [CrossRef]
Figure 1. Resistivity model used to obtain the synthetic data to validate the gradient calculation. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative.
Figure 1. Resistivity model used to obtain the synthetic data to validate the gradient calculation. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative.
Applsci 14 08846 g001
Figure 2. (a) Amplitude and (b) phase curves of the in-line electric field at frequencies 0.25 Hz, 0.50 Hz, and 0.75 Hz associated with the central receiver in the data from the model shown in Figure 1.
Figure 2. (a) Amplitude and (b) phase curves of the in-line electric field at frequencies 0.25 Hz, 0.50 Hz, and 0.75 Hz associated with the central receiver in the data from the model shown in Figure 1.
Applsci 14 08846 g002
Figure 3. (a) Estimated adjoint source pulses derived from residuals fitted separately for each frequency; (b) the yellow circles indicate the amplitude spectra of the pulses, along with the residuals at each frequency.
Figure 3. (a) Estimated adjoint source pulses derived from residuals fitted separately for each frequency; (b) the yellow circles indicate the amplitude spectra of the pulses, along with the residuals at each frequency.
Applsci 14 08846 g003
Figure 4. (a) Absolute values (log scale) of the E z component of the electric field, normalized by the source, computed by forward modeling with the source at the position of receiver 1 (at frequency of 0.50 Hz); (b) absolute values (log scale) of the λ E z field obtained by adjoint modeling with the sources estimated from data from receiver 1 (at frequency of 0.50 Hz); (c) gradient calculated from the product of the fields from forward modeling and the adjoint modeling.
Figure 4. (a) Absolute values (log scale) of the E z component of the electric field, normalized by the source, computed by forward modeling with the source at the position of receiver 1 (at frequency of 0.50 Hz); (b) absolute values (log scale) of the λ E z field obtained by adjoint modeling with the sources estimated from data from receiver 1 (at frequency of 0.50 Hz); (c) gradient calculated from the product of the fields from forward modeling and the adjoint modeling.
Applsci 14 08846 g004
Figure 5. Sum of the gradients associated with each receptor, obtained from the product between forward modeling and adjoint modeling fields. (a) Gradient at frequency of 0.25 Hz; (b) gradient at frequency of 0.50 Hz; (c) gradient at the frequency of 0.75 Hz.
Figure 5. Sum of the gradients associated with each receptor, obtained from the product between forward modeling and adjoint modeling fields. (a) Gradient at frequency of 0.25 Hz; (b) gradient at frequency of 0.50 Hz; (c) gradient at the frequency of 0.75 Hz.
Applsci 14 08846 g005
Figure 6. (a) Gradient relative to the vertical conductivity (at frequency of 0.25 Hz), preconditioned with precondition 1 (obtained with δ = 10 4 ); (b) gradient relative to vertical conductivity (at frequency of 0.25 Hz), preconditioned with precondition 2 (obtained with k = 1.0 ).
Figure 6. (a) Gradient relative to the vertical conductivity (at frequency of 0.25 Hz), preconditioned with precondition 1 (obtained with δ = 10 4 ); (b) gradient relative to vertical conductivity (at frequency of 0.25 Hz), preconditioned with precondition 2 (obtained with k = 1.0 ).
Applsci 14 08846 g006
Figure 7. (a) Resistivity model formed by two 25 ohm.m bodies at depths of 1.8 and 2.1 km in a homogeneous medium with 1 ohm.m resistivity. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative. (b) Gradient relative to the vertical conductivity obtained from 0.25 Hz data in a homogeneous medium with a resistivity of 1 ohm.m.
Figure 7. (a) Resistivity model formed by two 25 ohm.m bodies at depths of 1.8 and 2.1 km in a homogeneous medium with 1 ohm.m resistivity. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative. (b) Gradient relative to the vertical conductivity obtained from 0.25 Hz data in a homogeneous medium with a resistivity of 1 ohm.m.
Applsci 14 08846 g007
Figure 8. Curves showing the evolution of the objective function during the inversions for the model with two resistors, using the steepest descent method (black curve), momentum (blue curve), and the L-BFGS-B (red curve).
Figure 8. Curves showing the evolution of the objective function during the inversions for the model with two resistors, using the steepest descent method (black curve), momentum (blue curve), and the L-BFGS-B (red curve).
Applsci 14 08846 g008
Figure 9. Results of the inversions with (a) the steepest descent, (b) the momentum method, and (c) the L-BFGS-B method uses the observed data simulated from the model with two resistors. Black boxes represent the edge of the resistors.
Figure 9. Results of the inversions with (a) the steepest descent, (b) the momentum method, and (c) the L-BFGS-B method uses the observed data simulated from the model with two resistors. Black boxes represent the edge of the resistors.
Applsci 14 08846 g009
Figure 10. The MARE model used to simulate observed data in the inversions. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative.
Figure 10. The MARE model used to simulate observed data in the inversions. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative.
Applsci 14 08846 g010
Figure 11. (a) Initial model used in the inversions of data modeled from the MARE model. (b) Gradient relative to vertical conductivity at a frequency of 0.25 Hz, obtained from the initial model used in the inversions.
Figure 11. (a) Initial model used in the inversions of data modeled from the MARE model. (b) Gradient relative to vertical conductivity at a frequency of 0.25 Hz, obtained from the initial model used in the inversions.
Applsci 14 08846 g011
Figure 12. Results of the inversions with (a) the steepest descent, (b) the momentum method, and (c) the L-BFGS-B method, performed using observed data simulated from the MARE model. Black outlines indicate the boundaries of the resistors.
Figure 12. Results of the inversions with (a) the steepest descent, (b) the momentum method, and (c) the L-BFGS-B method, performed using observed data simulated from the MARE model. Black outlines indicate the boundaries of the resistors.
Applsci 14 08846 g012
Figure 13. Curves showing the evolution of the objective function during the inversions for the MARE model, using the steepest descent method (black curve), momentum (blue curve), and the L-BFGS-B (red curve).
Figure 13. Curves showing the evolution of the objective function during the inversions for the MARE model, using the steepest descent method (black curve), momentum (blue curve), and the L-BFGS-B (red curve).
Applsci 14 08846 g013
Figure 14. (a) Central section of the Marlim model used to simulate the observed data utilized in the inversions. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative. (b) Gradient (at a frequency of 0.25 Hz) relative to vertical conductivity, calculated from the homogeneous medium with a resistivity of 2.8 ohm.m.
Figure 14. (a) Central section of the Marlim model used to simulate the observed data utilized in the inversions. Triangles represent the receivers, and stars represent the sources used in the simulation. Not all sources are represented, and their height relative to the seafloor is illustrative. (b) Gradient (at a frequency of 0.25 Hz) relative to vertical conductivity, calculated from the homogeneous medium with a resistivity of 2.8 ohm.m.
Applsci 14 08846 g014
Figure 15. Resistivity models obtained from the inversions using (a) steepest descent, (b) momentum, and (c) L-BFGS-B. The black outline represents the boundary of the resistor present in the Marlim model.
Figure 15. Resistivity models obtained from the inversions using (a) steepest descent, (b) momentum, and (c) L-BFGS-B. The black outline represents the boundary of the resistor present in the Marlim model.
Applsci 14 08846 g015
Figure 16. Curves showing the evolution of the objective function during the inversions for the Marlim model, using the steepest descent method (black curve), momentum (blue curve), and L-BFGS-B (red curve).
Figure 16. Curves showing the evolution of the objective function during the inversions for the Marlim model, using the steepest descent method (black curve), momentum (blue curve), and L-BFGS-B (red curve).
Applsci 14 08846 g016
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

Valente, A.; Nascimento, D.; Costa, J. CSEM Optimization Using the Correspondence Principle. Appl. Sci. 2024, 14, 8846. https://doi.org/10.3390/app14198846

AMA Style

Valente A, Nascimento D, Costa J. CSEM Optimization Using the Correspondence Principle. Applied Sciences. 2024; 14(19):8846. https://doi.org/10.3390/app14198846

Chicago/Turabian Style

Valente, Adriany, Deivid Nascimento, and Jessé Costa. 2024. "CSEM Optimization Using the Correspondence Principle" Applied Sciences 14, no. 19: 8846. https://doi.org/10.3390/app14198846

APA Style

Valente, A., Nascimento, D., & Costa, J. (2024). CSEM Optimization Using the Correspondence Principle. Applied Sciences, 14(19), 8846. https://doi.org/10.3390/app14198846

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