Multi‐dimensional analyses of the SEAM controlled source electromagnetic data—the story of a blind test of interpretation workflows

Using a subset of the SEG Advanced Modeling Program Phase I controlled‐source electromagnetic data, we apply our standard controlled‐source electromagnetic interpretation workflows to delineate a simulated hydrocarbon reservoir. Experience learned from characterizing such a complicated model offers us an opportunity to refine our workflows to achieve better interpretation quality. The exercise proceeded in a blind test style, where the interpreting geophysicists did not know the true resistivity model until the end of the project. Rather, the interpreters were provided a traditional controlled‐source electromagnetic data package, including electric field measurements, interpreted seismic horizons, and well log data. Based on petrophysical analysis, a background resistivity model was established first. Then, the interpreters started with feasibility studies to establish the recoverability of the prospect and carefully stepped through 1D, 2.5D, and 3D inversions with seismic and well log data integrated at each stage. A high‐resistivity zone is identified with 1D analysis and further characterized with 2.5D inversions. Its lateral distribution is confirmed with a 3D anisotropic inversion. The importance of integrating all available geophysical and petrophysical data to derive more accurate interpretation is demonstrated.     

1D joint multi‐offset inversion of time‐domain marine controlled source electromagnetic data

The accurate estimation of sub‐seafloor resistivity features from marine controlled source electromagnetic data using inverse modelling is hindered due to the limitations of the inversion routines. The most commonly used one‐dimensional inversion techniques for resolving subsurface resistivity structures are gradient‐based methods, namely Occam and Marquardt. The first approach relies on the smoothness of the model and is recommended when there are no sharp resistivity boundaries. The Marquardt routine is relevant for many electromagnetic applications with sharp resistivity contrasts but subject to the appropriate choice of a starting model. In this paper, we explore the ability of different 1D inversion schemes to derive sub‐seafloor resistivity structures from time domain marine controlled source electromagnetic data measured along an 8‐km‐long profile in the German North Sea. Seismic reflection data reveal a dipping shallow amplitude anomaly that was the target of the controleld source electromagnetic survey. We tested four inversion schemes to find suitable starting models for the final Marquardt inversion. In this respect, as a first scenario, Occam inversion results are considered a starting model for the subsequent Marquardt inversion (Occam–Marquardt). As a second scenario, we employ a global method called Differential Evolution Adaptive Metropolis and sequentially incorporate it with Marquardt inversion. The third approach corresponds to Marquardt inversion introducing lateral constraints. Finally, we include the lateral constraints in Differential Evolution Adaptive Metropolis optimization, and the results are sequentially utilized by Marquardt inversion. Occam–Marquardt may provide accurate estimation of the subsurface features, but it is dependent on the appropriate conversion of different multi‐layered Occam model to an acceptable starting model for Marquardt inversion, which is not straightforward. Employing parameter spaces, the Differential Evolution Adaptive Metropolis approach can be pertinent to determine Marquardt a priori information; nevertheless, the uncertainties in Differential Evolution Adaptive Metropolis optimization will introduce some inaccuracies in Marquardt inversion results. Laterally constrained Marquardt may be promising to resolve sub‐seafloor features, but it is not stable if there are significant lateral changes of the sub‐seafloor structure due to the dependence of the method to the starting model. Including the lateral constraints in Differential Evolution Adaptive Metropolis approach allows for faster convergence of the routine with consistent results, furnishing more accurate estimation of a priori models for the subsequent Marquardt inversion.     

Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data

Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.     

Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan

Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.     

Repeatability analysis of land time‐lapse seismic data: CO2CRC Otway pilot project case study

Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO₂ sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO₂ sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO₂ sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO₂ sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.     

Burying receivers for improved time‐lapse seismic repeatability: CO2CRC Otway field experiment

4D seismic is widely used to remotely monitor fluid movement in subsurface reservoirs. This technique is especially effective offshore where high survey repeatability can be achieved. It comes as no surprise that the first 4D seismic that successfully monitored the CO₂ sequestration process was recorded offshore in the Sleipner field, North Sea. In the case of land projects, poor repeatability of the land seismic data due to low S/N ratio often obscures the time‐lapse seismic signal. Hence for a successful on shore monitoring program improving seismic repeatability is essential. Stage 2 of the CO2CRC Otway project involves an injection of a small amount (around 15,000 tonnes) of CO₂/CH₄ gas mixture into a saline aquifer at a depth of approximately 1.5 km. Previous studies at this site showed that seismic repeatability is relatively low due to variations in weather conditions, near surface geology and farming activities. In order to improve time‐lapse seismic monitoring capabilities, a permanent receiver array can be utilised to improve signal to noise ratio and hence repeatability. A small‐scale trial of such an array was conducted at the Otway site in June 2012. A set of 25 geophones was installed in 3 m deep boreholes in parallel to the same number of surface geophones. In addition, four geophones were placed into boreholes of 1–12 m depth. In order to assess the gain in the signal‐to‐noise ratio and repeatability, both active and passive seismic surveys were carried out. The surveys were conducted in relatively poor weather conditions, with rain, strong wind and thunderstorms. With such an amplified background noise level, we found that the noise level for buried geophones is on average 20 dB lower compared to the surface geophones. The levels of repeatability for borehole geophones estimated around direct wave, reflected wave and ground roll are twice as high as for the surface geophones. Both borehole and surface geophones produce the best repeatability in the 30–90 Hz frequency range. The influence of burying depth on S/N ratio and repeatability shows that significant improvement in repeatability can be reached at a depth of 3 m. The level of repeatability remains relatively constant between 3 and 12 m depths.     

Seismic time‐lapse imaging using interferometric least‐squares migration: Case study

Seismic time‐lapse surveys are susceptible to repeatability errors due to varying environmental conditions. To mitigate this problem, we propose the use of interferometric least‐squares migration to estimate the migration images for the baseline and monitor surveys. Here, a known reflector is used as the reference reflector for interferometric least‐squares migration, and the data are approximately redatumed to this reference reflector before imaging. This virtual redatuming mitigates the repeatability errors in the time‐lapse migration image. Results with synthetic and field data show that interferometric least‐squares migration can sometimes reduce or eliminate artifacts caused by non‐repeatability in time‐lapse surveys and provide a high‐resolution estimate of the time‐lapse change in the reservoir.     

3D joint inversion of magnetotelluric and airborne tipper data: a case study from the Morrison porphyry Cu–Au–Mo deposit,British Columbia,Canada

Z‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.     

Time‐lapse electrical resistivity imaging of solute transport in a karst conduit

The use of electrical resistivity surveys to locate karst conduits has shown mixed success. However, time‐lapse electrical resistivity imaging combined with salt injection improves conduit detection and can yield valuable insight into solute transport behaviour. We present a proof of concept above a known karst conduit in the Kentucky Horse Park (Lexington, Kentucky). A salt tracer solution was injected into a karst window over a 45‐min interval, and repeat resistivity surveys were collected every 20 min along a 125‐m transect near a monitoring well approximately 750 m downgradient from the injection site. In situ fluid conductivity measurements in the well peaked at approximately 25% of the initial value about 3 h after salt injection. Time‐lapse electrical resistivity inversions show two broad zones at the approximate conduit depth where resistivity decreased and then recovered in general agreement with in situ measurements. Combined salt injection and electrical resistivity imaging are a promising tool for locating karst conduits. The method is also useful for gaining insight into conduit geometry and could be expanded to include multiple electrical resistivity transects. Copyright © 2015 John Wiley & Sons, Ltd.     

Navigating marine electromagnetic transmitters using dipole field geometry

The marine controlled source electromagnetic (CSEM) technique has been adopted by the hydrocarbon industry to characterize the resistivity of targets identified from seismic data prior to drilling. Over the years, marine controlled source electromagnetic has matured to the point that four‐dimensional or time lapse surveys and monitoring could be applied to hydrocarbon reservoirs in production, or to monitor the sequestration of carbon dioxide. Marine controlled source electromagnetic surveys have also been used to target shallow resistors such as gas hydrates. These novel uses of the technique require very well constrained transmitter and receiver geometry in order to make meaningful and accurate geologic interpretations of the data. Current navigation in marine controlled source electromagnetic surveys utilize a long base line, or a short base line, acoustic navigation system to locate the transmitter and seafloor receivers. If these systems fail, then rudimentary navigation is possible by assuming the transmitter follows in the ship's track. However, these navigational assumptions are insufficient to capture the detailed orientation and position of the transmitter required for both shallow targets and repeat surveys. In circumstances when acoustic navigation systems fail we propose the use of an inversion algorithm that solves for transmitter geometry. This algorithm utilizes the transmitter's electromagnetic dipole radiation pattern as recorded by stationary, close range (<1000 m), receivers in order to model the geometry of the transmitter. We test the code with a synthetic model and validate it with data from a well navigated controlled source electromagnetic survey over the Scarborough gas field in Australia.     

Three‐dimensional inversion of towed streamer electromagnetic data

A towed streamer electromagnetic system capable of simultaneous seismic and electromagnetic data acquisition has recently been developed and tested in the North Sea. We introduce a 3D inversion methodology for towed streamer electromagnetic data that includes a moving sensitivity domain. Our implementation is based on the 3D integral equation method for computing responses and Fréchet derivatives and uses the re‐weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present two model studies relevant to hydrocarbon exploration in the North Sea. First, we demonstrate the ability of a towed electromagnetic system to detect and characterize the Harding field, a medium‐sized North Sea hydrocarbon target. We compare our 3D inversion of towed streamer electromagnetic data with 3D inversion of conventional marine controlled‐source electromagnetic data and observe few differences between the recovered models. Second, we demonstrate the ability of a towed streamer electromagnetic system to detect and characterize the Peon discovery, which is representative of an infrastructure‐led shallow gas play in the North Sea. We also present an actual case study for the 3D inversion of towed streamer electromagnetic data from the Troll field in the North Sea and demonstrate our ability to image all the Troll West Oil and Gas Provinces and the Troll East Gas Province. We conclude that 3D inversion of data from the current generation of towed streamer electromagnetic systems can adequately recover hydrocarbon‐bearing formations to depths of approximately 2 km. We note that by obviating the need for ocean‐bottom receivers, the towed streamer electromagnetic system enables electromagnetic data to be acquired over very large areas in frontier and mature basins for higher acquisition rates and relatively lower cost than conventional marine controlled‐source electromagnetic methods.     

Time‐lapse seismic imaging using regularized full‐waveform inversion with a prior model: which strategy?

Full‐waveform inversion is an appealing technique for time‐lapse imaging, especially when prior model information is included into the inversion workflow. Once the baseline reconstruction is achieved, several strategies can be used to assess the physical parameter changes, such as parallel difference (two separate inversions of baseline and monitor data sets), sequential difference (inversion of the monitor data set starting from the recovered baseline model) and double‐difference (inversion of the difference data starting from the recovered baseline model) strategies. Using synthetic Marmousi data sets, we investigate which strategy should be adopted to obtain more robust and more accurate time‐lapse velocity changes in noise‐free and noisy environments. This synthetic application demonstrates that the double‐difference strategy provides the more robust time‐lapse result. In addition, we propose a target‐oriented time‐lapse imaging using regularized full‐waveform inversion including a prior model and model weighting, if the prior information exists on the location of expected variations. This scheme applies strong prior model constraints outside of the expected areas of time‐lapse changes and relatively less prior constraints in the time‐lapse target zones. In application of this process to the Marmousi model data set, the local resolution analysis performed with spike tests shows that the target‐oriented inversion prevents the occurrence of artefacts outside the target areas, which could contaminate and compromise the reconstruction of the effective time‐lapse changes, especially when using the sequential difference strategy. In a strongly noisy case, the target‐oriented prior model weighting ensures the same behaviour for both time‐lapse strategies, the double‐difference and the sequential difference strategies and leads to a more robust reconstruction of the weak time‐lapse changes. The double‐difference strategy can deliver more accurate time‐lapse variation since it can focus to invert the difference data. However, the double‐difference strategy requires a preprocessing step on data sets such as time‐lapse binning to have a similar source/receiver location between two surveys, while the sequential difference needs less this requirement. If we have prior information about the area of changes, the target‐oriented sequential difference strategy can be an alternative and can provide the same robust result as the double‐difference strategy.     

Imaging a 3D geological structure from HEM, airborne magnetic and ground ERT data in Kalat-e-Reshm area, Iran

A set of geophysical data collected in an area in Iran are analyzed to check the validity of a geological map that was prepared in connection to a mineral prospecting project and also to image the spatial electrical resistivity distribution. The data set includes helicopter electromagnetic (HEM), airborne magnetic and ground electrical resistivity measurement. Occam approach was used to invert the HEM data to model the resistivity using a layered earth model with fixed thicknesses. The algorithm is based on a nonlinear inverse problem in a least-squares sense.The algorithm was tested on a part of an HEM dataset acquired with a DIGHEM helicopter EM system at Kalat-e-Reshm, Semnan in Iran. The area contains a resistive porphyry andesite that is covered by Eocene sedimentary units. The results are shown as resistivity sections and maps confirming the existence of an arc like resistive structure in the survey area. The resistive andesite seems to be thicker than it is indicated in the geological maps. The results are compared with the reduced to the pole (RTP) airborne magnetic anomaly field data as well as with two ground resistivity profiles. We found reasonable correlations between the HEM 1D resistivity models and 2D models from electrical resistivity tomography (ERT) inversions. A 3D visualization of the 1D models along all flight lines provided a useful tool for the study of spatial variations of the resistivity structure in the investigation area.     

A comparison of receiver technologies in borehole MMR and EM surveys

A series of downhole magnetometric resistivity (DHMMR) and downhole electromagnetic (DHEM) surveys were conducted near Broken Hill, New South Wales, Australia, and at Zinkgruvan, Sweden, to determine how probe and receiver equipment choices affect the amount of noise visible in borehole MMR and EM data. Noise analyses performed on the data, using the standard deviation to gauge the relative noise levels between different probes and receiver systems, indicate that high noise levels in MMR data result primarily from the use of a three‐component EM probe, which has a reduced effective area and hence a higher noise floor compared with a single‐component EM probe. High noise, attributable to cultural sources such as nearby power lines, in either MMR or EM data can be reduced through the use of full‐waveform, multipurpose receiver systems. These systems allow for the use of tapered stacking, which is a more effective method of eliminating coherent noise associated with power‐line transients than the boxcar‐stacking method used by traditional receiver systems.     

Sharp spatially constrained inversion with applications to transient electromagnetic data

Time‐domain electromagnetic data are conveniently inverted by using smoothly varying 1D models with fixed vertical discretization. The vertical smoothness of the obtained models stems from the application of Occam‐type regularization constraints, which are meant to address the ill‐posedness of the problem. An important side effect of such regularization, however, is that horizontal layer boundaries can no longer be accurately reproduced as the model is required to be smooth. This issue can be overcome by inverting for fewer layers with variable thicknesses; nevertheless, to decide on a particular and constant number of layers for the parameterization of a large survey inversion can be equally problematic. Here, we present a focusing regularization technique to obtain the best of both methodologies. The new focusing approach allows for accurate reconstruction of resistivity distributions using a fixed vertical discretization while preserving the capability to reproduce horizontal boundaries. The formulation is flexible and can be coupled with traditional lateral/spatial smoothness constraints in order to resolve interfaces in stratified soils with no additional hypothesis about the number of layers. The method relies on minimizing the number of layers of non‐vanishing resistivity gradient, instead of minimizing the norm of the model variation itself. This approach ensures that the results are consistent with the measured data while favouring, at the same time, the retrieval of horizontal abrupt changes. In addition, the focusing regularization can also be applied in the horizontal direction in order to promote the reconstruction of lateral boundaries such as faults. We present the theoretical framework of our regularization methodology and illustrate its capabilities by means of both synthetic and field data sets. We further demonstrate how the concept has been integrated in our existing spatially constrained inversion formalism and show its application to large‐scale time‐domain electromagnetic data inversions.     

Field test of sub‐basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data

The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section.     

A radio magnetotelluric study to evaluate the extents of a limestone quarry in Estonia

Electromagnetic signals from distant radio transmitters in the frequency range 15–250 kHz were measured to model an electrical resistivity structure beneath 7 profiles in the vicinity of the Karinu limestone quarry in Estonia with the aim to map the extent of the economically exploitable limestone. The resistivity models from a 2D inversion of determinant resistivity and phase values using an Occam type of regularization contained reasonably accurate information about the geometry, namely depth to the top and the bottom of the target high‐resistivity limestone. The resistivity models correlated well with existing geological evidences as well as information from closely located boreholes. However, the sharp lithological boundaries seen in the boreholes were not resolved exactly in the resistivity models. This is probably because of the smoothing regularization used in the inversion process. Combined use of borehole data together with resistivity models resulted in two major geological interpretations; a) towards the western part of the existing limestone quarry there is a NNW to NS striking fault, covered by post‐glacial sediments, b) a potential cost‐effective exploitable area containing high quality highly resistive limestone is located south of the existing quarry. This case study shows the applicability of the reasonably fast radio magnetotelluric (RMT) method for the exploration of near‐surface resources.     

Seismic time‐lapse effects of solution salt mining – a feasibility study

This article addresses the question whether time‐lapse seismic reflection techniques can be used to follow and quantify the effects of solution salt mining. Specifically, the production of magnesium salts as mined in the north of the Netherlands is considered. The use of seismic time‐lapse techniques to follow such a production has not previously been investigated. For hydrocarbon production and CO₂ storage, time‐lapse seismics are used to look at reservoir changes mainly caused by pressure and saturation changes in large reservoirs, while for solution mining salt is produced from caverns with a limited lateral extent, with much smaller production volumes and a fluid (brine) replacing a solid (magnesium salt). In our approach we start from the present situation of the mine and then study three different production scenarios, representing salt production both in vertical and lateral directions of the mine. The present situation and future scenarios have been transformed into subsurface models that were input to an elastic finite‐difference scheme to create synthetic seismic data. These data have been analysed and processed up to migrated seismic images, such that time‐lapse analyses of intermediate and final results could be done. From the analyses, it is found that both vertical and lateral production is visible well above the detection threshold in difference data, both at pre‐imaging and post‐imaging stages. In quantitative terms, an additional production of the mine of 6 m causes time‐shifts in the order of 2 ms (pre‐imaging) and 4 ms (post‐imaging) and amplitude changes of above 20% in the imaged sections. A laterally oriented production causes even larger amplitude changes at the edge of the cavern due to replacement of solid magnesium salt with brine introducing a large seismic contrast. Overall, our pre‐imaging and post‐imaging time‐lapse analysis indicates that the effects of solution salt mining can be observed and quantified on seismic data. The effects seem large enough to be observable in real seismic data containing noise.     

Improved normalization of time‐lapse seismic data using normalized root mean square repeatability data to improve automatic production and seismic history matching in the Nelson field

Updating of reservoir models by history matching of 4D seismic data along with production data gives us a better understanding of changes to the reservoir, reduces risk in forecasting and leads to better management decisions. This process of seismic history matching requires an accurate representation of predicted and observed data so that they can be compared quantitatively when using automated inversion. Observed seismic data is often obtained as a relative measure of the reservoir state or its change, however. The data, usually attribute maps, need to be calibrated to be compared to predictions. In this paper we describe an alternative approach where we normalize the data by scaling to the model data in regions where predictions are good. To remove measurements of high uncertainty and make normalization more effective, we use a measure of repeatability of the monitor surveys to filter the observed time‐lapse data. We apply this approach to the Nelson field. We normalize the 4D signature based on deriving a least squares regression equation between the observed and synthetic data which consist of attributes representing measured acoustic impedances and predictions from the model. Two regression equations are derived as part of the analysis. For one, the whole 4D signature map of the reservoir is used while in the second, 4D seismic data is used from the vicinity of wells with a good production match. The repeatability of time‐lapse seismic data is assessed using the normalized root mean square of measurements outside of the reservoir. Where normalized root mean square is high, observations and predictions are ignored. Net: gross and permeability are modified to improve the match. The best results are obtained by using the normalized root mean square filtered maps of the 4D signature which better constrain normalization. The misfit of the first six years of history data is reduced by 55 per cent while the forecast of the following three years is reduced by 29 per cent. The well based normalization uses fewer data when repeatability is used as a filter and the result is poorer. The value of seismic data is demonstrated from production matching only where the history and forecast misfit reductions are 45% and 20% respectively while the seismic misfit increases by 5%. In the best case using seismic data, it dropped by 6%. We conclude that normalization with repeatability based filtering is a useful approach in the absence of full calibration and improves the reliability of seismic data.     

Electromagnetic mapping of electrical conductivity beneath the Columbia basalts1

Sedimentary rocks beneath the Columbia River Basalt Group are recognized as having potential for oil and gas production, but the overlying layered basalts effectively mask seismic reflections from the underlying sediments. Four electromagnetic (EM) methods have been applied on profiles crossing Boylston Ridge, a typical east–west trending anticline of the Yakima Fold Belt, in an attempt to map the resistivity interface between the basalts and the sediments and to map variations in structure and resistivity within the sediments. The EM surveys detected strong variations in resistivity within the basalts, and in particular the continuous magnetotelluric array profiling (EMAP) revealed resistivity lows beneath the surface anticlines. These low resistivity zones probably coincide with fracturing in the core of the anticlines and they appear to correlate well with similar zones of low seismic velocity observed on a nearby seismic profile. The controlled-source EM surveys (in-loop transient, long-offset transient, and variable-offset frequency-domain) were designed in anticipation of relatively uniform high resistivity basalts, and were found to have been seriously distorted by the intrabasalt conductors discovered in the field. In particular, the resistivity sections derived from 1D inversions were found to be inconsistent and misleading. The EMAP survey provided the most information about the subsurface resistivity distribution, and was certainly the most cost-effective. However, both controlled-source and EMAP surveys call for accurate 2D or 3D inversion to accommodate the geological objectives of this project.