BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Kees Wapenaar

Kees Wapenaar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.geo-phphysics.app-pheess.IVeess.SPphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

On the relation between the propagator matrix and the Marchenko focusing function

The propagator matrix "propagates" a full wave field from one depth level to another, accounting for all propagation angles and evanescent waves. The Marchenko focusing function forms the nucleus of data-driven Marchenko redatuming and imaging schemes, accounting for internal multiples. These seemingly different concepts appear to be closely related to each other. With this insight, the strong aspects of the propagator matrix (such as the handling of evanescent waves) can be transferred to the focusing function. Vice-versa, the propagator matrix inherits from the focusing function that it can be retrieved from the reflection response, which reduces its sensitivity to the subsurface model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Geothermal Well Doublet for Research and Heat Supply of the TU Delft Campus

A geothermal well doublet, designed for two primary aims of research and commercial heat supply, is planned to be constructed on the campus of Delft University of Technology. The plans include a comprehensive research programme, including installation of a wide range of instruments alongside a comprehensive logging and coring programme and an extensive surface monitoring network. The wells will be cored, with samples from all representative geological units down to the reservoir. An extensive suite of well-logs is planned to provide detailed information on the properties of the various units. Fibre-optic cables will be installed in both wells all the way down to the reservoir section, which is anticipated to be at 2200m, in the Lower Cretaceous Delft Sandstone. The Delft sandstone is well-known as a reservoir rock for natural gas in the West Netherlands Basin. The wells will be operated by a commercial entity, but the infrastructure is designed and explicitly installed as a research infrastructure. As such it will become part of the European EPOS (European Plate Observing System, https://www.epos-ip.org), such that accessibility and data availability will not be limited to TU Delft researchers. The university has made a decision in principle to start drilling no later than 2020.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A modified Marchenko method to retrieve the wave field inside a layered metamaterial from reflection measurements at the surface

With the Marchenko method it is possible to retrieve the wave field inside a medium from its reflection response at the surface. To date, this method has predominantly been applied to naturally occurring materials. In this paper we extend the Marchenko method for applications in metamaterials with, in the low-frequency limit, effective negative constitutive parameters. We illustrate the method with a numerical example, which shows that the method works well for vertically propagating plane waves and for dipping plane waves with small horizontal slownesses. The proposed method has potential applications for example in non-destructive testing of layered materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new role for adaptive filters in Marchenko equation-based methods for the attenuation of internal multiples

We have seen many developments in Marchenko equation-based methods for internal multiple attenuation in the past years. Starting from a wave-equation based method that required a smooth velocity model, there are now Marchenko equation-based methods that do not require any model information or user-input. In principle, these methods accurately predict internal multiples. Therefore, the role of the adaptive filter has changed for these methods. Rather than needing an aggressive adaptive filter to compensate for inaccurate internal multiple predictions, only a conservative adaptive filter is needed to compensate for minor amplitude and/or phase errors in the internal multiple predictions caused by imperfect acquisition and preprocessing of the input data. We demonstate that a conservative adaptive filter can be used to improve the attenuation of internal multiples when applying a Marchenko multiple elimination (MME) method to a 2D line of streamer data. In addition, we suggest that an adaptive filter can be used as a feedback mechanism to improve the preprocessing of the input data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A proposal for Marchenko-based target-oriented full waveform inversion

The Hessian matrix plays an important role in correct interpretation of the multiple scattered wave fields inside the FWI frame work. Due to the high computational costs, the computation of the Hessian matrix is not feasible. Consequently, FWI produces overburden related artifacts inside the target zone model, due to the lack of the exact Hessian matrix. We have shown here that Marchenko-based target-oriented Full Waveform Inversion can compensate the need of Hessian matrix inversion by reducing the non-linearity due to overburden effects. This is achieved by exploiting Marchenko-based target replacement to remove the overburden response and its interactions with the target zone from residuals and inserting the response of the updated target zone into the response of the entire medium. We have also shown that this method is more robust with respect to prior information than the standard gradient FWI. Similarly to standard Marchenko imaging, the proposed method only requires knowledge of the direct arrival time from a focusing point to the surface and the reflection response of the medium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An acoustic imaging method for layered non-reciprocal media

Given the increasing interest for non-reciprocal materials, we propose a novel acoustic imaging method for layered non-reciprocal media. The method we propose is a modification of the Marchenko imaging method, which handles multiple scattering between the layer interfaces in a data-driven way. We start by reviewing the basic equations for wave propagation in a non-reciprocal medium. Next, we discuss Green's functions, focusing functions, and their mutual relations, for a non-reciprocal horizontally layered medium. These relations form the basis for deriving the modified Marchenko method, which retrieves the wave field inside the non-reciprocal medium from reflection measurements at the boundary of the medium. With a numerical example we show that the proposed method is capable of imaging the layer interfaces at their correct positions, without artefacts caused by multiple scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Correcting for imperfectly sampled data in the iterative Marchenko scheme

The Marchenko method retrieves the responses to virtual sources in the subsurface, accounting for all orders of multiples. The method is based on two integral representations for focusing and Green's functions. In discretized form these integrals are represented by finite summations over the acquisition geometry. Consequently, the method requires ideal geometries of regularly sampled and co-located sources and receivers. However, a recent study showed that this restriction can, in theory, be relaxed by deconvolving the irregularly-sampled results with certain point spread functions (PSFs).The results are then reconstructed as if they were acquired using a perfect geometry. Here, the iterative Marchenko scheme is adapted in order to include these PSFs; thus, showing how imperfect sampling can be accounted for in practical situations. Next, the new methodology is tested on a 2D numerical example. The results show clear improvement between the proposed scheme and the standard iterative scheme. By removing the requirement for perfect geometries the Marchenko method can be more widely applied to field data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-driven internal multiple elimination applications using imperfectly sampled reflection data

We consider reflection data that have been subsampled by 70% and use Point-Spread-Functions to reconstruct the original data. The subsampled, original and reconstructed reflection data are used to image the medium of interest with the Marchenko method. The image obtained using the subsampled data shows artifacts caused by internal multiples, which are eliminated when the original and reconstructed data are used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Elastodynamic single-sided homogeneous Green's function representation: Theory and numerical examples

The homogeneous Green's function is the difference between an impulse response and its time-reversal. According to existing representation theorems, the homogeneous Green's function associated with source-receiver pairs inside a medium can be computed from measurements at a boundary enclosing the medium. However, in many applications such as seismic imaging, time-lapse monitoring, medical imaging, non-destructive testing, etc., media are only accessible from one side. A recent development of wave theory has provided a representation of the homogeneous Green's function in an elastic medium in terms of wavefield recordings at a single (open) boundary. Despite its single-sidedness, the elastodynamic homogeneous Green's function representation accounts for all orders of scattering inside the medium. We present the theory of the elastodynamic single-sided homogeneous Green's function representation and illustrate it with numerical examples for 2D laterally-invariant media. For propagating waves, the resulting homogeneous Green's functions match the exact ones within numerical precision, demonstrating the accuracy of the theory. In addition, we analyse the accuracy of the single-sided representation of the homogeneous Green's function for evanescent wave tunnelling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Employing internal multiples in time-lapse seismic monitoring, using the Marchenko method

Time-lapse seismic monitoring aims at resolving changes in a producing reservoir from changes in the reflection response. When the changes in the reservoir are very small, the changes in the seismic response can become too small to be reliably detected. In theory, multiple reflections can be used to improve the detectability of traveltime changes: a wave that propagates several times down and up through a reservoir layer will undergo a larger time shift due to reservoir changes than a primary reflection. Since we are interested in monitoring very local changes (usually in a thin reservoir layer), it would be advantageous if we could identify the reservoir-related internal multiples in the complex reflection response of the entire subsurface. We introduce a Marchenko-based method to isolate these multiples from the complete reflection response and illustrate the potential of this method with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Marchenko-based target replacement, accounting for all orders of multiple reflections

In seismic monitoring, one is usually interested in the response of a changing target zone, embedded in a static inhomogeneous medium. We introduce an efficient method which predicts reflection responses at the earth's surface for different target-zone scenarios, from a single reflection response at the surface and a model of the changing target zone. The proposed process consists of two main steps. In the first step, the response of the original target zone is removed from the reflection response, using the Marchenko method. In the second step, the modelled response of a new target zone is inserted between the overburden and underburden responses. The method fully accounts for all orders of multiple scattering and, in the elastodynamic case, for wave conversion. For monitoring purposes, only the second step needs to be repeated for each target-zone model. Since the target zone covers only a small part of the entire medium, the proposed method is much more efficient than repeated modelling of the entire reflection response.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reciprocity and representation theorems for flux- and field-normalised decomposed wave fields

We consider wave propagation problems in which there is a preferred direction of propagation. To account for propagation in preferred directions, the wave equation is decomposed into a set of coupled equations for waves that propagate in opposite directions along the preferred axis. This decomposition is not unique. We discuss flux-normalised and field-normalised decomposition in a systematic way, analyse the symmetry properties of the decomposition operators and use these symmetry properties to derive reciprocity theorems for the decomposed wave fields, for both types of normalisation. Based on the field-normalised reciprocity theorems, we derive representation theorems for decomposed wave fields. In particular we derive double- and single-sided Kirchhoff-Helmholtz integrals for forward and backward propagation of decomposed wave fields. The single-sided Kirchhoff-Helmholtz integrals for backward propagation of field-normalised decomposed wave fields find applications in reflection imaging, accounting for multiple scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Marchenko method for evanescent waves

With the Marchenko method, Green's functions in the subsurface can be retrieved from seismic reflection data at the surface. State-of-the-art Marchenko methods work well for propagating waves but break down for evanescent waves. This paper discusses a first step towards extending the Marchenko method for evanescent waves and analyses its possibilities and limitations. In theory both the downward and upward decaying components can be retrieved. The retrieval of the upward decaying component appears to be very sensitive to model errors, but the downward decaying component, including multiple reflections, can be retrieved in a reasonably stable and accurate way. The reported research opens the way to develop new Marchenko methods that can handle refracted waves in wide-angle reflection data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unified matrix-vector wave equation, reciprocity and representations

The matrix-vector wave equation is a compact first-order differential equation. It was originally used for the analysis of elastodynamic plane waves in laterally invariant media. It has been extended by various authors for laterally varying media. Other authors derived a similar formalism for other wave phenomena. This paper starts with a unified formulation of the matrix-vector wave equation for 3D inhomogeneous, dissipative media. The wave vector, source vector and operator matrix are specified in the appendices for acoustic, quantum mechanical, electromagnetic, elastodynamic, poroelastodynamic, piezoelectric and seismoelectric waves. It is shown that the operator matrix obeys unified symmetry relations for all these wave phenomena. Next, unified matrix-vector reciprocity theorems of the convolution and correlation type are derived, utilising the symmetry properties of the operator matrix. These theorems formulate mathematical relations between two wave states in the same spatial domain. A unified wave field representation is obtained by replacing one of the states in the convolution-type reciprocity theorem by a Green's state. By replacing both states in the correlation-type reciprocity theorem by Green's states, a unified representation of the homogeneous Green's matrix is obtained. Applications of the unified reciprocity theorems and representations for forward and inverse wave problems are briefly indicated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unified wave field retrieval and imaging method for inhomogeneous non-reciprocal media

Acoustic imaging methods often ignore multiple scattering. This leads to false images in cases where multiple scattering is strong.Marchenko imaging has recently been introduced as a data-driven way to deal with internal multiple scattering. Given the increasing interest in non-reciprocal materials, both for acoustic and electromagnetic applications, we propose to modify the Marchenko method for imaging such materials. We formulate a unified wave equation for non-reciprocal materials, exploiting the similarity between acoustic and electromagnetic wave phenomena. This unified wave equation forms the basis for deriving reciprocity theorems that interrelate wave fields in a non-reciprocal medium and its complementary version. Next, we reformulate these theorems for downgoing and upgoing wave fields. From these decomposed reciprocity theorems we derive representations of the Green's function inside the non-reciprocal medium, in terms of the reflection response at the surface and focusing functions inside the medium and its complementary version. These representations form the basis for deriving a modified version of the Marchenko method to retrieve the wave field inside a non-reciprocal medium and to form an image, free from artefacts related to multiple scattering. We illustrate the proposed method at the hand of the numerically modeled reflection response of a horizontally layered medium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Virtual acoustics in inhomogeneous media with single-sided access

A virtual acoustic source inside a medium can be created by emitting a time-reversed point-source response from the enclosing boundary into the medium. However, in many practical situations the medium can be accessed from one side only. In those cases the time-reversal approach is not exact. Here, we demonstrate the experimental design and use of complex focusing functions to create virtual acoustic sources and virtual receivers inside an inhomogeneous medium with single-sided access. The retrieved virtual acoustic responses between those sources and receivers mimic the complex propagation and multiple scattering paths of waves that would be ignited by physical sources and recorded by physical receivers inside the medium. The possibility to predict complex virtual acoustic responses between any two points inside an inhomogeneous medium, without needing a detailed model of the medium, has large potential for holographic imaging and monitoring of objects with single-sided access, ranging from photoacoustic medical imaging to the monitoring of induced-earthquake waves all the way from the source to the earth's surface.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Virtual sources and receivers in the real Earth -- considerations for practical applications

To enhance monitoring of the subsurface, virtual sources and receivers inside the subsurface can be created from seismic reflection data at the surface of the Earth using the Marchenko method. The response between these virtual sources and receivers can be obtained through the use of homogeneous Green's function retrieval. A homogeneous Green's function is a superposition of a Green's function and its time-reversal. The main aim of this paper is to obtain accurate homogeneous Green's functions from field data. Classical homogeneous Green's function retrieval requires an unrealistic enclosing recording surface, however, by using a recently proposed single-sided retrieval scheme, this requirement can be avoided. We first demonstrate the principles of using the single-sided representation on synthetic data and show that different source signatures can be taken into account. Because the Marchenko method is sensitive to recording limitations of the reflection data, we study five cases of recording limitations with synthetic data and demonstrate their effects on the final result. Finally, the method is demonstrated on a pre-processed field dataset which fulfills the requirements for applying the single-sided Green's function retrieval scheme. The scheme has the potential to be used in future applications, such as source localization.

0 citations0 reviewsNo code linkedOpen access