scholarly journals Multiphysics data modeling and imaging for exploration in the southern Rocky Mountains

2018 ◽  
Vol 6 (3) ◽  
pp. SG59-SG78 ◽  
Author(s):  
Maria Soledad Velasco ◽  
David Alumbaugh ◽  
Emmanuel Schnetzler

We carried out a multidata geophysics study in southern Colorado to explore for [Formula: see text] reservoirs in an area where seismic imaging is very limited due to the mountainous terrain, the presence of high-velocity volcanic rocks, and difficulty in obtaining land access permits. We have developed a modeling/interpretation methodology using ground magnetotelluric data as well as airborne magnetic and electromagnetic data combined with public domain gravity data and existing well and seismic data. We used the integration of these data sets to produce a series of 2D and 3D geophysical models that reveal basin architecture previously poorly defined through the analysis of limited seismic and well data alone. We found that this type of analysis aids in decreasing uncertainty in the interpreted geologic cross sections and a better understanding of the structural complexities of the region. Through the application of machine learning methods, we are also able to integrate several data sets into a mathematical framework resulting in a predictive model of spatial [Formula: see text] distribution. The integration of the interpretations from all data sets, predictive analytics results, and knowledge of [Formula: see text] production, allows us to delineate areas of interest for further exploration.

2021 ◽  
Author(s):  
Mohammad Shehata ◽  
Hideki Mizunaga

<p>Long-period magnetotelluric and gravity data were acquired to investigate the US cordillera's crustal structure. The magnetotelluric data are being acquired across the continental USA on a quasi-regular grid of ∼70 km spacing as an electromagnetic component of the National Science Foundation EarthScope/USArray Program. International Gravimetreique Bureau compiled gravity Data at high spatial resolution. Due to the difference in data coverage density, the geostatistical joint integration was utilized to map the subsurface structures with adequate resolution. First, a three-dimensional inversion of each data set was applied separately.</p><p>The inversion results of both data sets show a similarity of structure for data structuralizing. The individual result of both data sets is resampled at the same locations using the kriging method by considering each inversion model to estimate the coefficient. Then, the Layer Density Correction (LDC) process's enhanced density distribution was applied to MT data's spatial expansion process. Simple Kriging with varying Local Means (SKLM) was applied to the residual analysis and integration. For this purpose, the varying local means of the resistivity were estimated using the corrected gravity data by the Non-Linear Indicator Transform (NLIT), taking into account the spatial correlation. After that, the spatial expansion analysis of MT data obtained sparsely was attempted using the estimated local mean values and SKLM method at the sections where the MT survey was carried out and for the entire area where density distributions exist. This research presents the integration results and the stand-alone inversion results of three-dimensional gravity and magnetotelluric data.</p>


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. B1-B12 ◽  
Author(s):  
Juliane Hübert ◽  
Benjamin M. Lee ◽  
Lijuan Liu ◽  
Martyn J. Unsworth ◽  
Jeremy P. Richards ◽  
...  

We have evaluated results from a study combining airborne z-axis tipper electromagnetic (ZTEM) and ground-based magnetotelluric (MT) data to image an epithermal system in British Columbia. The spatially coincident use of these two methods allowed for a direct comparison of both data sets in the overlapping frequency band and showed that both measurements were consistent. Inversion of just the ZTEM data suffered from the lack of electric field amplitude information, which could be provided by the MT data. Three-dimensional inversion modeling of the two individual data sets was performed. Models of electrical resistivity derived from both data sets were consistent and could be correlated with the geological and structural setting of the mineralization. Gold is associated with disseminated pyrite and marcasite in quartz-sericite-altered felsic volcanic rocks and intrusions, especially near the contact with mafic volcanic rocks and a late diorite intrusion. The quartz-sericite alteration yields a conductivity anomaly, relative to the more resistive mafic country rocks. Although ZTEM and MT do not possess the resolution of the geologic model derived from borehole data, our model agrees well with a regional assessment of the deposit.


2015 ◽  
Vol 3 (4) ◽  
pp. SAC99-SAC106 ◽  
Author(s):  
Irina Filina ◽  
Nicholas Delebo ◽  
Gopal Mohapatra ◽  
Clayton Coble ◽  
Gary Harris ◽  
...  

A 3D gravity model was developed in the western Gulf of Mexico in the East Breaks and Alaminos Canyon protraction areas. This model integrated 3D seismic, gravity, and well data; it was constructed in support of a proprietary seismic reprocessing project and was updated iteratively with seismic. The gravity model was built from seismic horizons of the bathymetry, salt layers, and the acoustic basement; however, the latter was only possible to map in seismic data during the latest iterations. In addition, a deep layer representing the Moho boundary was derived from gravity and constrained by public-domain refraction data. A 3D density distribution was derived from the seismic velocity volume using a modified Gardner equation. The modification comprised imposing a depth dependency on the Gardner coefficient, which is constant in the classic Gardner equation. The modified coefficient was derived from well data in the study area and public-domain velocity-density data sets. The forward-calculated gravity response of the composed density model was then compared with the observed gravity field, and the mismatch was analyzed jointly by a seismic interpreter and a gravity modeler. Adjustments were then made to the gravity model to ensure that the resultant salt model was geologically reasonable and supported by gravity, seismic, and well data sets. The output of the gravity modeling was subsequently applied to the next phase of seismic processing. Overall, this integration resulted in a more robust salt model, which has led to significant improvements in subsalt seismic imaging. The analysis of the regional trend in the observed gravity field suggested that a stretched continental crust underlay our seismic reprocessing area, with an oceanic-continental transition zone located to the southeast of our reprocessing region.


Geophysics ◽  
2013 ◽  
Vol 78 (5) ◽  
pp. B317-B328 ◽  
Author(s):  
Anjana K. Shah ◽  
Paul A. Bedrosian ◽  
Eric D. Anderson ◽  
Karen D. Kelley ◽  
James Lang

We combined aeromagnetic, induced polarization, magnetotelluric, and gravity surveys as well as drillhole geologic, alteration, magnetic susceptibility, and density data for exploration and characterization of the Cu-Au-Mo Pebble porphyry deposit. This undeveloped deposit is almost completely concealed by postmineralization sedimentary and volcanic rocks, presenting an exploration challenge. Individual geophysical methods primarily assist regional characterization. Positive chargeability and conductivity anomalies are observed over a broad region surrounding the deposit, likely representing sulfide minerals that accumulated during multiple stages of hydrothermal alteration. The mineralized area occupies only a small part of the chargeability anomaly because sulfide precipitation was not unique to the deposit, and mafic rocks also exhibit strong chargeability. Conductivity anomalies similarly reflect widespread sulfides as well as water-saturated glacial sediments. Mineralogical and magnetic susceptibility data indicate magnetite destruction primarily within the Cu-Au-Mo mineralized area. The magnetic field does not show a corresponding anomaly low but the analytic signal does in areas where the deposit is not covered by postmineralization igneous rocks. The analytic signal shows similar lows over sedimentary rocks outside of the mineralized area, however, and cannot uniquely distinguish the deposit. We find that the intersection of positive chargeability anomalies with analytic signal lows, indicating elevated sulfide concentrations but low magnetite at shallow depths, roughly delineates the deposit where it is covered only by glacial sediments. Neither chargeability highs nor analytic signal lows are present where the deposit is covered by several hundred meters of sedimentary and volcanic rocks, but a 3D resistivity model derived from magnetotelluric data shows a corresponding zone of higher conductivity. Gravity data highlight geologic features within the deposit, including shallow diorite sills that locally contain higher-grade mineralization. The results thus show ways in which an integrated survey approach might be used to distinguish zones of potentially economic mineralization.


2021 ◽  
Vol 95 (2) ◽  
Author(s):  
Mirjam Bilker-Koivula ◽  
Jaakko Mäkinen ◽  
Hannu Ruotsalainen ◽  
Jyri Näränen ◽  
Timo Saari

AbstractPostglacial rebound in Fennoscandia causes striking trends in gravity measurements of the area. We present time series of absolute gravity data collected between 1976 and 2019 on 12 stations in Finland with different types of instruments. First, we determine the trends at each station and analyse the effect of the instrument types. We estimate, for example, an offset of 6.8 μgal for the JILAg-5 instrument with respect to the FG5-type instruments. Applying the offsets in the trend analysis strengthens the trends being in good agreement with the NKG2016LU_gdot model of gravity change. Trends of seven stations were found robust and were used to analyse the stabilization of the trends in time and to determine the relationship between gravity change rates and land uplift rates as measured with global navigation satellite systems (GNSS) as well as from the NKG2016LU_abs land uplift model. Trends calculated from combined and offset-corrected measurements of JILAg-5- and FG5-type instruments stabilized in 15 to 20 years and at some stations even faster. The trends of FG5-type instrument data alone stabilized generally within 10 years. The ratio between gravity change rates and vertical rates from different data sets yields values between − 0.206 ± 0.017 and − 0.227 ± 0.024 µGal/mm and axis intercept values between 0.248 ± 0.089 and 0.335 ± 0.136 µGal/yr. These values are larger than previous estimates for Fennoscandia.


2020 ◽  
Vol 94 (9) ◽  
Author(s):  
Lars E. Sjöberg

Abstract As the KTH method for geoid determination by combining Stokes integration of gravity data in a spherical cap around the computation point and a series of spherical harmonics suffers from a bias due to truncation of the data sets, this method is based on minimizing the global mean square error (MSE) of the estimator. However, if the harmonic series is increased to a sufficiently high degree, the truncation error can be considered as negligible, and the optimization based on the local variance of the geoid estimator makes fair sense. Such unbiased types of estimators, derived in this article, have the advantage to the MSE solutions not to rely on the imperfectly known gravity signal degree variances, but only the local error covariance matrices of the observables come to play. Obviously, the geoid solution defined by the local least variance is generally superior to the solution based on the global MSE. It is also shown, at least theoretically, that the unbiased geoid solutions based on the KTH method and remove–compute–restore technique with modification of Stokes formula are the same.


2020 ◽  
Vol 221 (3) ◽  
pp. 1542-1554 ◽  
Author(s):  
B C Root

SUMMARY Current seismic tomography models show a complex environment underneath the crust, corroborated by high-precision satellite gravity observations. Both data sets are used to independently explore the density structure of the upper mantle. However, combining these two data sets proves to be challenging. The gravity-data has an inherent insensitivity in the radial direction and seismic tomography has a heterogeneous data acquisition, resulting in smoothed tomography models with de-correlation between different models for the mid-to-small wavelength features. Therefore, this study aims to assess and quantify the effect of regularization on a seismic tomography model by exploiting the high lateral sensitivity of gravity data. Seismic tomography models, SL2013sv, SAVANI, SMEAN2 and S40RTS are compared to a gravity-based density model of the upper mantle. In order to obtain similar density solutions compared to the seismic-derived models, the gravity-based model needs to be smoothed with a Gaussian filter. Different smoothening characteristics are observed for the variety of seismic tomography models, relating to the regularization approach in the inversions. Various S40RTS models with similar seismic data but different regularization settings show that the smoothening effect is stronger with increasing regularization. The type of regularization has a dominant effect on the final tomography solution. To reduce the effect of regularization on the tomography models, an enhancement procedure is proposed. This enhancement should be performed within the spectral domain of the actual resolution of the seismic tomography model. The enhanced seismic tomography models show improved spatial correlation with each other and with the gravity-based model. The variation of the density anomalies have similar peak-to-peak magnitudes and clear correlation to geological structures. The resolvement of the spectral misalignment between tomographic models and gravity-based solutions is the first step in the improvement of multidata inversion studies of the upper mantle and benefit from the advantages in both data sets.


2015 ◽  
Vol 24 (07) ◽  
pp. 1550050 ◽  
Author(s):  
E. Matsinos ◽  
G. Rasche

In a previous paper, we reported the results of a partial-wave analysis (PWA) of the pion–nucleon (πN) differential cross-sections (DCSs) of the CHAOS Collaboration and came to the conclusion that the angular distribution of their π+p data sets is incompatible with the rest of the modern (meson factory) database. The present work, re-addressing this issue, has been instigated by a number of recent improvements in our analysis, namely regarding the inclusion of the theoretical uncertainties when investigating the reproduction of experimental data sets on the basis of a given "theoretical" solution, modifications in the parametrization of the form factors of the proton and of the pion entering the electromagnetic part of the πN amplitude, and the inclusion of the effects of the variation of the σ-meson mass when fitting the ETH model of the πN interaction to the experimental data. The new analysis of the CHAOS DCSs confirms our earlier conclusions and casts doubt on the value for the πN Σ term, which Stahov, Clement and Wagner have extracted from these data.


2018 ◽  
Vol 7 (1) ◽  
pp. 1-20 ◽  
Author(s):  
Aleksandra Osika ◽  
Małgorzata Wistuba ◽  
Ireneusz Malik

Abstract The aim of the study is to reconstruct the development of landslide relief in the Kamienne Mountains (Central Sudetes, SW Poland) based on a DEM from LiDAR data. Analyses of relief and geological maps in ArcGIS 10.5 and of slope cross-sections in Surfer 14 allowed to distinguish different types of landslide relief, developed in latites and trachybasalts lying above claystones and mudstones. The types vary from small, poorly visible landslides to vast landslides with complex relief. They were interpreted as consecutive stages of geomorphic evolution of hillslope-valley topography of the study area. Two main schemes have been established which explain the development of landslide slopes in the Kamienne Mts: (1) upslope, from the base of the slope towards the mountain ridge and (2) downslope, beginning on the top of the mountain ridge. The direction of landslide development depends on the thickness of volcanic rocks in relation to underlying sedimentary rocks. When the latter appear only in the lowest part of the slope, landslides develop upslope. If sedimentary rocks dominate on the slope and volcanic rocks form only its uppermost part, landslides develop downslope. The results show that landsliding leads to significant modifications of relief of the study area, including complete degradation of mountain ridges.


2009 ◽  
Vol 2 (1) ◽  
pp. 87-98 ◽  
Author(s):  
C. Lerot ◽  
M. Van Roozendael ◽  
J. van Geffen ◽  
J. van Gent ◽  
C. Fayt ◽  
...  

Abstract. Total O3 columns have been retrieved from six years of SCIAMACHY nadir UV radiance measurements using SDOAS, an adaptation of the GDOAS algorithm previously developed at BIRA-IASB for the GOME instrument. GDOAS and SDOAS have been implemented by the German Aerospace Center (DLR) in the version 4 of the GOME Data Processor (GDP) and in version 3 of the SCIAMACHY Ground Processor (SGP), respectively. The processors are being run at the DLR processing centre on behalf of the European Space Agency (ESA). We first focus on the description of the SDOAS algorithm with particular attention to the impact of uncertainties on the reference O3 absorption cross-sections. Second, the resulting SCIAMACHY total ozone data set is globally evaluated through large-scale comparisons with results from GOME and OMI as well as with ground-based correlative measurements. The various total ozone data sets are found to agree within 2% on average. However, a negative trend of 0.2–0.4%/year has been identified in the SCIAMACHY O3 columns; this probably originates from instrumental degradation effects that have not yet been fully characterized.


Sign in / Sign up

Export Citation Format

Share Document