Validation of an electroseismic and seismoelectric modeling code, for layered earth models, by the explicit homogeneous space solutions

We have developed an analytically based, energy fluxnormalized numerical modeling code (ESSEMOD), capable of modeling the wave propagation of all existing ElectroSeismic and SeismoElectric source-receiver combinations in horizontally layered configurations. We compare the results of several of these modeled source-receiver combinations in a homogeneous medium with explicitly derived homogeneous space Green’s function solutions, in order to be able to validate the results of ESSEMOD both in arrival times and amplitudes. Especially the amplitudes are important due to the fact that the main reason seismoelectric phenomena are not yet used in industry, are the weak amplitudes of these phenomena. Here we show that ESSEMOD correctly models the wave propagation of components of the electric field generated by different components of bulk forces, as well as the particle velocity fields generated by a bulk force source and an electric current source. We are capable of validating both amplitudes and arrival times of the results of ESSEMOD for all electroseismic and seismoelectric source-receiver combinations in homogeneous media. Herewith, we reduce uncertainty in our modeling results (also for heterogeneous scenarios) and can get better insights in which parameters affect the amplitudes most. In addition, we show that ESSEMOD is capable of modeling reciprocal source-receivers combinations correctly, implicitly indicating correct modeling of both geometrical configurations (source located above or below the receiver level). ESSEMOD can now be used for comparison with and validation of existing seismoelectric layered earth numerical modeling codes. Afterwards, ESSEMOD can be used for validation of existing seismoelectric finite-element and finite-difference codes.