1,720,955 research outputs found
Coupled poroelastic waves and electromagnetic fields in layered media: Theory, Modeling, and Interferometric Synthesis
In this thesis, I study coupled poroelastic waves and electromagnetic fields in layered media. The focus is two-fold:1. Increase the theoretical and physical understanding of the seismo-electromagnetic phenomenon by analytically-based numerical modeling.2. Investigate the potential of seismo-electromagnetic interferometry.After presenting the governing equations that form the basis of the theoretical framework, I capture this system into a matrix-vector representation of the wave equation. I first use literature eigenvector sets, which I normalize with respect to power-flux. I then derive new, alternative power-flux normalized eigenvector sets that I prove to be numerically more stable and accurate. The eigenvector sets form the basis of the analytically-based numerical modeling code `ESSEMOD' that I developed to model seismo-electromagnetic wave/field propagation/diffusion in layered-Earth media. The alternative eigenvector set models scenarios with no seismo-electromagnetic coupling correctly, where the literature eigenvector sets fail. In addition, the alternative set properly deals with scenarios where both small amplitude signals and large amplitude signals occur in the record, whereas the literature eigenvector sets result in noise levels masking the small events. The same holds for scenarios with a small seismo-electromagnetic coupling coefficient. I design an effective global reflection scheme that properly describes the primary and multiple reflections in the models. I implement the correct boundary conditions to account for scenarios with a free-surface, and also for scenarios containing fluid/porous medium/fluid transitions. To transform all the seismo-electromagnetic source-receiver combinations in a numerically effective way back from the horizontal wavenumber-frequency domain to the space-frequency domain, I derive and implement explicit Fourier-Bessel transformations.I then validate the developed modeling code in numerous ways. First of all, I compare the results of seismo-EM layer-code modeling in a homogeneous medium with explicit homogeneous space Green's function expressions. This comparison provides a clear validation that the layer-code models the dynamic responses in homogeneous scenarios correctly. Next, I check numerical consistency by carrying out reciprocity checks. I study homogeneous space models, models containing a free-surface and models with interfaces.As a next step, I validate the modeling results of seismo-EM layer-code modeling for typical seismo-electromagnetic laboratory configurations, i.e. models containing fluid/porous medium/fluid transitions. I first compare the purely electromagnetic part of the seismo-EM layer-code with an independently developed purely electromagnetic layered-Earth code. The results match perfectly in both phase and amplitude for full transmission and pure reflection experiments, as well as for a combination of both. I then carry out a seismo-electromagnetic reciprocity test for a fluid halfspace overlying a porous medium halfspace, proving that the coupled poroelastic and electromagnetic fields are modeled consistently and yield the expected results.As a final validation step, I compare ESSEMOD with an independently developed seismo-electromagnetic layered-Earth modeling code. The results display an almost perfect match in both phase and relative amplitudes, and a constant amplitude correction factor of 4 needs to be applied to let the absolute amplitudes match.I then carry out a small feasibility test to study the potential of the seismo-electromagnetic effect for exploration purposes. I investigate different source-receiver combinations for the same model, and focus on the signal strength recorded at different distances from the target depth level. I conclude that for the source-receiver combinations studied, the electric field due to a volume injection monopole source, as well as the magnetic field due to a seismic bulk force source, yield the strongest converted signals. The receiver-distance from the target of interest plays an important role in the signal measurability. The closer the receivers to the target, the higher the signal strengths. However, when the receivers are located too close to the target, the coseismic reflected fields can mask the interface response fields that we are mainly interested in.Next, I study if nature itself can help us to overcome the very low signal-to-noise ratio of seismo-electromagnetic converted fields, by investigating the effects of thin-bed geological structures on the seismo-electromagnetic signal. To investigate the effects of bed-thinning on the seismo-electromagnetic interference patterns, I numericallysimulate seismo-electromagnetic wave propagation through horizontally layered media with different amounts and thicknesses of thin-beds. I demonstrate seismo-electromagnetic sensitivity to changes in medium parameters on a spatial scale much smaller than the seismic resolution. By simulating moving oil/water contacts duringproduction, where the oil layer is gradually being thinned, seismo-electromagnetic signals are proven very sensitive to oil/water contacts.I now explore the application of interferometric techniques to the seismo-electromagnetic system, which might eventually lead to an improved signal-to-noise ratio of the weak converted fields.I derive the theory for interferometric retrieval of 2D SH-TE seismo-electromagnetic Green's functions.Using both a circular source configuration and a line source configuration, I show that it is possible to correctly retrieve the dynamic seismo-electromagnetic 2D SH-TE response in a homogeneous medium, using seismic boundary sources only. Using seismo-EM layer-code data, I then show that it is also possible to correctly retrieve the direct shear wave-related causal coseismic field in a homogeneous medium, in both phase and amplitude. To obtain a perfect match in absolute amplitudes, I apply a single linear scaling factor. I finally carry out interferometric experiments in a model containing a single interface at 800 m depth, proving that it is possible to correctlyretrieve all 2D SH-TE causal seismic-related direct and reflected coseismic fields, as well as interface response fields, by cross-correlation interferometry, using seismic boundary sources only.These results are promising for the application of 3D seismo-electromagnetic interferometry using seismo-EM layer-code modeling, and later on, in the field.Next, I present an alternative way to effectively decompose fields into their up- and downgoing components and different field types, using recordings at multiple depth levels. I present the theory of this MDL decomposition scheme, followed by successful decomposition of synthetic elastodynamic data sets. I additionally study the implications of laterally-varying media on the horizontal wavenumber-frequency domain MDL decomposition scheme.I demonstrate successful decomposition, using an acoustic approximation and applying a combined multi-component / MDL decomposition approach, of a field data set recorded in Annerveen, in the North of the Netherlands. I address how to effectively use the MDL decomposition scheme in a unified fashion, applied to all wave phenomena including seismo-electromagnetic phenomena.I then make a step towards seismo-electromagnetic inversion, presenting an effective way to carry out a seismo-electromagnetic sensitivity analysis using resolution functions. I start by explaining the theory of resolution functions using a seismo-electromagnetic example. I define the seismo-electromagnetic resolution function for inversion for a bulk density perturbation. I demonstrate the effectiveness of this method by first carrying out a purely electromagnetic sensitivity analysis for a point perturbation in conductivity, located in an isotropic homogeneous half-space. These results are compared with literature results based on analytical homogeneous space Green's function expressions. The result using the seismo-EM layer-code is nearly identical to the literature result. The position of the scatterer is correctly resolved. At the end of this section, I present the results of the fully-coupled seismo-electromagnetic senstivity analysis for a bulk density contrast for a specific source-receiver combination, using single-frequency multi-component line data. I show that the coupled seismo-electromagnetic system is sensitive to a perturbation in bulk density and that the position of the perturbation can be correctly recovered.I finalize this thesis by discussing potential seismo-electromagnetic applications, as well as by providing a brief outlook for future research
Seismoelectric Modelling of the Flux-Normalized P-SV-TM Propagation Mode
Elastodynamic and electromagnetic processes are coupled together in saturated, porous media, by a phenomenon known as the electrokinetic effect. In horizontally layered media, the seismoelectric system, which contains the coupled elastodynamic and electromagnetic systems, can be separated into two independent modes of propagation: SH-TE and P-SV-TM. The SH-TE mode contains horizontally polarized shear waves coupled with transverse electric polarized electromagnetic waves. In the P-SV-TM mode, both fast and slow compressional waves are coupled with vertically polarized shear waves and transverse magnetic polarized electromagnetic waves. In this thesis, the P-SV-TM mode of the two-dimensional seismoelectric system was expressed in the form of both the two-way and one-way wave equations. The principle of normalizing energy flux across boundaries was applied, improving the matrix amplitude balance of the system and allowing for the implementation of one-way reciprocity theorems. We carried out full-waveform modelling of the flux-normalized P-SV-TM seismoelectric system in a 2-D fluid-saturated, horizontally-stratified, porous media. Both one-way and two-way wavefields were modelled, allowing the composition of one-way wavefields into two-way wavefields to be clearly observed. We investigated both the generation of electromagnetic fields due to the propagation of a seismic pertubation and the generation of seismic waves due to the propagation of a diffusive electromagnetic wave. Reciprocity of the wavefields was verified by applying reciprocity theorems to both one-way and two-way wave vectors. The electromagnetic field that is created when a seismic wave traverses a contrast in medium parameters is rapidly attenuated during propagation. To mitigate the decay in the amplitude of the signal with distance, we modelled a Vertical ElectroSeismic Profiling (VESP) survey, in which receivers could be placed in near proximity to the target layer. In another model, the sensitivity of the seismoelectric method to pore fluid contrasts was tested by simulating the influx of contaminants into an aquifer. It was observed that a small change in the conductivity of the aquifer led to a significant change in the strength of the electromagnetic signal that was generated at the top of the aquifer.Applied Geophysics and PetrophysicsGeotechnologyCivil Engineering and Geoscience
Seismo-Electromagnetic Thin-Bed Responses: Natural Signal Enhancements?
We study if nature can help us overcome the very low signal-to-noise ratio of seismo-electromagnetic converted fields by investigating the effects of thin-bed geological structures on the seismo-electromagnetic signal. To investigate the effects of bed thinning on the seismo-electromagnetic interference patterns, we numerically simulate seismo-electromagnetic wave propagation through horizontally layered media with different amounts and thicknesses of thin beds. We distinguish two limits of bed thickness. Below the upper limit, the package of thin beds starts acting like an “effective” medium. Below the lower limit, further thinning does not affect the seismo-electromagnetic interface response signal strength anymore. We demonstrate seismo-electromagnetic sensitivity to changes in medium parameters on a spatial scale much smaller than the seismic resolution. Increasing amounts of thin beds can cause the interface response signal strength to increase or decrease. Whether constructive or destructive interference occurs seems to be dependent on the seismo-electromagnetic coupling coefficient contrasts. When the combined result of the contrast, between upper half-space and package of thin beds and the internal thin-bed contrast, is positive, constructive interference occurs. Destructive interference occurs when the combined contrast is negative. Maximum amplitude tuning occurs for thicknesses of thin-bed packages similar to the dominant pressure and shear wavelengths. Artifacts due to model periodicity are excluded by comparing periodic media with random models. By simulating moving oil/water contacts during production, where the oil layer is gradually being thinned, seismo-electromagnetic signals are proven very sensitive to oil/water contacts. An oil layer with a thickness of <1% of the dominant shear wavelength is still recognized
Analysis of 2D homogeneous space solutions of the seismoelectric P-SV-TM mode for interferometric purposes
Seismic and electromagnetic imaging methods both provide the geophysicist with different types of medium parameters. Seismic methods are sensitive to the elastic properties of the medium, while electromagnetic methods are sensitive to the electric properties. In porous-saturated media, these two wave fields occur as a coupled system, which is known as 'seismoelectrics'. This coupling is caused by physical interactions at the grain surface boundary and is a function of several medium parameters, such as dynamic permeability. This medium parameter is valuable to the oil and gas industry, as well to the field of hydrology. By conducting a seismoelectric survey it would theoretically be possible to provide an extra control on this medium parameter. However, both practice and theory have shown that this coupling mechanism also results in a low signal-to-noise ratio. A possible solution to this problem would be to apply interferometric Green's retrieval, which is a technique based on stacking of cross-correlated data. This approach has been proved successful for the SH-TE mode in 1D. The SH-TE mode forms together with the P-SV-TM mode, the total seismoelectric system. In this thesis the first steps are taken towards the proof that this technique could also work for the P-SV-TM mode of the system. This is supported by a modelling experiment of 2D homogeneous space solutions of the seismoelectric P-SV-TM mode for different configurations. This analysis turned out that the unwanted artefacts observed in the interferometric retrieval are generated by cross-correlations between P-waves and SV-waves.IDEA League Joint Master's in Applied GeophysicsGeoscience & EngineeringCivil Engineering and Geoscience
Going Beyond Counting First Authors in Author Co-citation Analysis
The present study examines one of the fundamental aspects of author co-citation analysis (ACA) - the way co-citation
counts are defined. Co-citation counting provides the data on which all subsequent statistical analyses and mappings
are based, and we compare ACA results based on two different types of co-citation counting - the traditional type that
only counts the first one among a cited work's authors on the one hand and a non-traditional type that takes into
account the first 5 authors of a cited work on the other hand. Results indicate that the picture produced through this non-traditional author co-citation counting contains more coherent author groups and is therefore considerably clearer. However, this picture represents fewer specialties in the research field being studied than that produced through the traditional first-author co-citation counting when the same number of top-ranked authors is selected and analyzed. Reasons for these effects are discussed
Variations on the Author
“Variations on the Author” discusses two of Eduardo Coutinho’s recent films (Um Dia na Vida, from 2010, and Últimas Conversas, posthumously released in 2015) and their contribution to the general question of documentary authorship. The director’s filmography is characterized by a consistent yet self-effacing form of authorial self-inscription: Coutinho often features as an interviewer that rather than express opinions propels discourses; an interviewer that is good at listening. This mode of self-inscription characterizes him as an author who is not expressive but who is nonetheless markedly present on the screen. In Um Dia na Vida, however, Coutinho is completely absent form the image, while Últimas Conversas, on the contrary, includes a confessional prologue that moves the director from the margins to the center of his films. This article examines the ways in which these works stand out in the filmography of a director who offers new insights into the notion of cinematic authorship
Appropriate Similarity Measures for Author Cocitation Analysis
We provide a number of new insights into the methodological discussion about author cocitation analysis. We first argue that the use of the Pearson correlation for measuring the similarity between authors’ cocitation profiles is not very satisfactory. We then discuss what kind of similarity measures may be used as an alternative to the Pearson correlation. We consider three similarity measures in particular. One is the well-known cosine. The other two similarity measures have not been used before in the bibliometric literature. Finally, we show by means of an example that our findings have a high practical relevance.information science;Pearson correlation;cosine;similarity measure;author cocitation analysis
Comparison of Eigenvectors for Coupled Seismo-Electromagnetic Layered-Earth Modeling
We study the accuracy and numerical stability of three eigenvector sets for modelling the coupled poroelastic and electromagnetic layered-Earth response. We use a known eigenvector set, its flux-normalized version and a newly derived flux-normalized set. The new set is chosen such that the system is properly uncoupled when the coupling between the poroelastic and electromagnetic fields vanishes. We carry out two different numerical stability tests: the first test focuses on the internal system, eigenvector and eigenvalue consistency; the second test investigates the stability and preciseness of the flux-normalized systems by looking at identity relations. We find that the known set shows the largest deviation for both tests, whereas the new set performs best. In two additional numerical modelling experiments, these numerical inaccuracies are shown to generate numerical noise levels comparable to small signals, such as signals coming from the important interface conversion responses, especially when the coupling coefficient is small. When coupling vanishes completely, the known set does not produce proper results. The new set produces numerically stable and accurate results in all situations. We therefore strongly recommend to use this newly derived set for future layered-Earth seismo-electromagnetic modelling experiments
Electromagnetic & Seismoelectric sensitivity analysis using resolution functions
In the field of exploration geophysics various methods are applied to determine the physical properties of the subsurface of the Earth. Some of the methods most widely used are seismic and electromagnetic surveys, which are each used according to the type of information that is being sought and their ability to provide that information. The method of seismoelectrics is in that regard a promising technique because theoretically it should be sensitive to a wide range of subsurface parameters, spanning both the acoustic and the electromagnetic methods. We aim to perform a parameter sensitivity analysis for the seismoelectric problem, investigating how perturbations in different parameters affect the data and how well these perturbations can be inverted for. We will first study some elements of inverse theory with a special focus on how to construct resolution functions from the basic integral equation of scattering theory. This integral equation can be derived using both superposition and reciprocity principles. This will be shown for both the electromagnetic and the acoustic case, before moving on to the seismoelectric case; it will become clear that the latter poses additional challenges. In addition, the results for the electromagnetic case can be compared directly with the results obtained by [Slob and Mulder, 2011] who performed an electromagnetic parameter sensitivity analysis. The results of this comparison could serve as a validation of the seismoelectric forward modelling code ESSEMOD [Grobbe and Slob, 2013], which will be used to generate synthetic data for this thesis. It is concluded that the superposition principle and the reciprocity theorem provide identical expressions for a scattered field. This conclusion was the basis for deriving the complete reciprocity theorem for the seismoelectric system, which was in turn used to derive resolution functions for a perturbation in bulk density and a perturbation in the seismoelectric coupling coefficient. For the electromagnetic case, we show that the resolution functions computed, using explicit analytical Green's function solutions are identical to resolution functions computed with data produced by ESSEMOD and that this indeed serves as a validation of ESSEMOD. Using the same approach, we have successfully computed a resolution function for inversion for the coupled seismoelectric system for a perturbation in bulk density. It is recommended that this innovative result is used as the basis for the analysis of the seismoelectric sensitivity to more complex parameter-contrasts, which could serve as an assessment of the true potential of the seismoelectric method.IDEA League Joint Master's in Applied GeophysicsGeoscience and Remote SensingCivil Engineering and Geoscience
- …
