HPC and edge elements for geophysical electromagnetic problems: an overview

In Finite Element Methods for solving electromagnetic problems, the use of Nédélec Elements has become very popular. In fact, Nédélec Elements are often said to be a cure to many difficulties that are encountered, particularly eliminating spurious solutions, and are claimed to yield accurate results. In this paper, we present our first steps in developing a Nédélec Elements code for simulation of geophysical electromagnetic problems and first ideas on how implement the key issues of Edge Elements in an efficient way on HPC platforms

HPC and edge elements for geophysical electromagnetic problems: an overview

In Finite Element Methods for solving electromagnetic problems, the use of Nédélec Elements has become very popular. In fact, Nédélec Elements are often said to be a cure to many difficulties that are encountered, particularly eliminating spurious solutions, and are claimed to yield accurate results. In this paper, we present our first steps in developing a Nédélec Elements code for simulation of geophysical electromagnetic problems and first ideas on how implement the key issues of Edge Elements in an efficient way on HPC platforms

Elastic full-waveform inversion of vertical seismic profile data acquired with distributed acoustic sensors

Distributed acoustic sensing (DAS) is a rapidly developing technology particularly useful for the acquisition of vertical seismic profile (VSP) surveys. DAS data are increasingly used for seismic imaging, but not for estimating rock properties. We have developed a workflow for estimating elastic properties of the subsurface using full-waveform inversion (FWI) of DAS VSP data. Whereas conventional borehole geophones usually measure three components of particle velocity, DAS measures a single quantity, which is an approximation of the strain or strain rate along the fiber. Standard FWI algorithms are developed for particle velocity data, and hence their application to DAS data requires conversion of these data to particle velocity along the fiber. This conversion can be accomplished by a specially designed filter. Field measurements show that the conversion result is close to vertical particle velocity as measured by geophones. Elastic time-domain FWI of a synthetic multioffset VSP data set for a vertical well shows that the inversion of the vertical component alone is sufficient to recover elastic properties of the subsurface. Application of the proposed workflow to a multioffset DAS data set acquired at the CO2CRC Otway Project site in Victoria, Australia, reveals salient subhorizontal layering consistent with the known geology of the site. The inverted VP model at the well location matches the upscaled VP log with a correlation coefficient of 0.85

Interpretation of deep directional resistivity measurements acquired in high-angle and horizontal wells using 3-D inversion

The interpretation of resistivity measurements acquired in high-angle and horizontal wells is a critical technical problem in formation evaluation. We develop an efficient parallel 3-D inversion method to estimate the spatial distribution of electrical resistivity in the neighbourhood of a well from deep directional electromagnetic induction measurements. The methodology places no restriction on the spatial distribution of the electrical resistivity around arbitrary well trajectories. The fast forward modelling of triaxial induction measurements performed with multiple transmitter-receiver configurations employs a parallel direct solver. The inversion uses a pre-conditioned gradient-based method whose accuracy is improved using the Wolfe conditions to estimate optimal step lengths at each iteration. The large transmitter-receiver offsets, used in the latest generation of commercial directional resistivity tools, improve the depth of investigation to over 30 m from the wellbore. Several challenging synthetic examples confirmthe feasibility of the full 3-D inversion-based interpretations for these distances, hence enabling the integration of resistivity measurements with seismic amplitude data to improve the forecast of the petrophysical and fluid properties. Employing parallel direct solvers for the triaxial induction problems allows for large reductions in computational effort, thereby opening the possibility to invert multiposition 3-D data in practical CPU times

Edge-elements for geophysical electromagnetic problems: a new implementation challenge

Peer ReviewedPostprint (author's final draft

Genome-wide sequence analyses of ethnic populations across Russia

The Russian Federation is the largest and one of the most ethnically diverse countries in the world, however no centralized reference database of genetic variation exists to date. Such data are crucial for medical genetics and essential for studying population history. The Genome Russia Project aims at filling this gap by performing whole genome sequencing and analysis of peoples of the Russian Federation. Here we report the characterization of genome-wide variation of 264 healthy adults, including 60 newly sequenced samples. People of Russia carry known and novel genetic variants of adaptive, clinical and functional consequence that in many cases show allele frequency divergence from neighboring populations. Population genetics analyses revealed six phylogeographic partitions among indigenous ethnicities corresponding to their geographic locales. This study presents a characterization of population-specific genomic variation in Russia with results important for medical genetics and for understanding the dynamic population history of the world's largest country