Joint inversion of surface and borehole 4D gravity data for continuous characterization of fluid contact movement

Author(s):  
Richard A. Krahenbuhl ◽  
Yaoguo Li
Keyword(s):  
2020 ◽  
Author(s):  
W. Soyer ◽  
R. Mackie ◽  
S. Hallinan ◽  
F. Miorelli ◽  
A. Pavesi
Keyword(s):  

2021 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino

<p>We present results from a joint inversion study of new seismic and gravity data to constrain a 2D high-resolution image of one of the most prominent geophysical anomalies of the European Alps: the Ivrea geophysical body (IGB). Our work exploits both new data and multidisciplinary a priori constraints, to better resolve the shallow crustal structure in the Ivrea-Verbano zone (IVZ), where the IGB is known to reach anomalously shallow depths and partially outcrop at the surface.</p><p>A variety of previous studies, ranging from gravity surveys to vintage refraction seismics and recent local earthquake tomographies (Solarino et al. 2018, Diehl et al. 2009), provide comprehensive but spatially sparse information on the IGB structure, which we aim at investigating at higher resolution, along a linear profile crossing the IVZ. To this purpose, we deployed 10 broadband seismic stations (MOBNET pool, IG CAS Prague), 5 km spaced along a linear West-East profile, along Val Sesia and crossing Lago Maggiore. This network operated for 27 months and allowed us to produce a new database of ca. 1000 seismic high-quality receiver functions (RFs). In addition, we collected new gravity data in the IVZ, with a data coverage of 1 gravity point every 1-2 km along the seismic profile. The newly collected data was used to set up an inversion scheme, in which RFs and gravity anomalies are jointly used to constrain the shape and the physical property contrasts across the IGB interface.</p><p>We model the IGB as a single interface between far-field constraints, whose geometry is defined by the coordinates of four nodes which may vary in space, and  density and V<sub>S</sub> shear-wave velocity contrasts associated with the interface itself, varying independently. A Markov chain Monte Carlo (MCMC) sampling method with Metropolis-Hastings selection rule was implemented to efficiently explore the model space, directing the search towards better fitting areas.</p><p>For each model, we perform ray-tracing and RFs migration using the actual velocity structure both for migration and computation of synthetic RFs, to be compared with the observations via cross-correlation of the migration images. Similarly, forward gravity modelling for a 2D density distribution is implemented and the synthetic gravity anomaly is compared with the observations along the profile. The joint inversion performance is the product of these two misfits.</p><p>The inversion results show that the IGB reaches the shallowest depths in the western part of the profile, preferentially locating the IGB interface between 3 and 7 km depth over a horizontal distance of ca. 20 km (between Boccioleto and Civiasco, longitudes 8.1 and 8.3). Within this segment, the shallowest point reaches up to 1 km below sea level. The found density and velocity contrasts are in agreement with rock physics properties of various units observed in the field and characterized in earlier studies.</p>


2020 ◽  
Vol 8 (4) ◽  
pp. SS47-SS62
Author(s):  
Thibaut Astic ◽  
Dominique Fournier ◽  
Douglas W. Oldenburg

We have carried out petrophysically and geologically guided inversions (PGIs) to jointly invert airborne and ground-based gravity data and airborne magnetic data to recover a quasi-geology model of the DO-27 kimberlite pipe in the Tli Kwi Cho (also referred to as TKC) cluster. DO-27 is composed of three main kimberlite rock types in contact with each other and embedded in a granitic host rock covered by a thin layer of glacial till. The pyroclastic kimberlite (PK), which is diamondiferous, and the volcanoclastic kimberlite (VK) have anomalously low density, due to their high porosity, and weak magnetic susceptibility. They are indistinguishable from each other based upon their potential-field responses. The hypabyssal kimberlite (HK), which is not diamondiferous, has been identified as highly magnetic and remanent. Quantitative petrophysical signatures for each rock unit are obtained from sample measurements, such as the increasing density of the PK/VK unit with depth and the remanent magnetization of the HK unit, and are represented as a Gaussian mixture model (GMM). This GMM guides the PGI toward generating a 3D quasi-geology model with physical properties that satisfies the geophysical data sets and the petrophysical signatures. Density and magnetization models recovered individually yield volumes that have physical property combinations that do not conform to any known petrophysical characteristics of the rocks in the area. A multiphysics PGI addresses this problem by using the GMM as a coupling term, but it puts a volume of the PK/VK unit at a location that is incompatible with geologic information from drillholes. To conform to that geologic knowledge, a fourth unit is introduced, PK-minor, which is petrophysically and geographically distinct from the main PK/VK unit. This inversion produces a quasi-geology model that presents good structural locations of the diamondiferous PK unit and can be used to provide a resource estimate or decide the locations of future drillholes.


2015 ◽  
Author(s):  
Dengguo Zhou ◽  
Weizhong Wang ◽  
Jie Zhang ◽  
Daniel R.H. O'Connell
Keyword(s):  

2020 ◽  
Vol 222 (2) ◽  
pp. 845-860
Author(s):  
Yann Ziegler ◽  
Sébastien B Lambert ◽  
Ibnu Nurul Huda ◽  
Christian Bizouard ◽  
Séverine Rosat

SUMMARY The rotational motions of the internal Earth layers induce resonances in the Earth nutations and tidal gravimetric response to external luni-solar gravitational forcings. The characterization of these resonances is a mean of investigating the deep Earth properties since their amplitudes and frequencies depend on a few fundamental geophysical parameters. In this work, we focus on the determination of the free core nutation and free inner core nutation periods and quality factors from the Bayesian inversion of VLBI and gravimetric data. We make a joint inversion of data from both techniques and show that, even if the results are only slightly different from the inversion of VLBI data alone, such approach may be valuable in the future if the accuracy of gravimetric data increases. We also briefly discuss the polar motion resonance, which is related to the Chandler Wobble as seen from the diurnal frequency band. Our overall estimates of the FCN period and quality factor, TFCN = (−430.2, −429.8) solar days and QFCN = (15 700, 16 700), respectively, are in good agreement with other studies, albeit slightly different for unclear reasons. Despite some concerns about the detection and characterization of the FICN, it seems that we could also successfully estimate its period, TFICN = (+600, +1300) solar days, and give a loose estimate of the upper bound on its quality factor.


2016 ◽  
Vol 64 (5) ◽  
pp. 1655-1672 ◽  
Author(s):  
Katarzyna Miernik ◽  
Adrian Bogacz ◽  
Adam Kozubal ◽  
Tomasz Danek ◽  
Marek Wojdyła

Sign in / Sign up

Export Citation Format

Share Document