We present an interpretation method for the gravity anomaly of an arbitrary interface separating two homogeneous media. It consists essentially of a downward continuation of the observed anomaly and the division of the continued anomaly by a scale factor involving the density contrast between the media. The knowledge of the interface depth at isolated points is used to estimate the depth d1 of the shallowest point of the interface, the density contrast ρ between the two media, and the coefficients c1 and c2 of a first-order polynomial representing a linear trend to be removed from data. The solutions are stabilized by introducing a damping parameter in the computation of the downward-continued anomaly by the equivalent layer method. Different from other interface mapping methods using gravity data, the proposed method: (1) takes into account the presence of an undesirable linear trend in data; (2) requires just intervals for both ρ (rather than the knowledge of its true value) and coefficients c1 and c2; and (3) does not require the knowledge of the average interface depth z0. As a result of (3), the proposed method does not call for extensive knowledge of the interface depth to obtain a statistically significant estimate of z0; rather, it is able to use the knowledge of the interface depth at just a few isolated points to estimate d1, ρ, c1, and c2. Tests using synthetic data confirm that the method produces good and stable estimates as far as the established premises (smooth interface separating two homogeneous media and, at most, the presence of an unremoved linear trend in data) are not violated. If the density contrast is not uniform, the method may still be applied using Litinsky's concept of effective density. The method was applied to gravity data from Recôncavo Basin, Brazil, producing good correlations of estimated lows and terraces in the basement with corresponding known geological features.
We demonstrate the potential of gravity data to detect and to locate in-depth subtle normal faults in the basement relief of a sedimentary basin. This demonstration is accomplished by inverting the gravity data with the constraint that the estimated basement relief presents local abrupt faults and is smooth elsewhere. We inverted the gravity data from the onshore Almada Basin in northeastern Brazil, and we mapped several normal faults whose locations and plane geometries were already known from seismic imaging. The inversion method delineated well both the discontinuities with small or large slips and a sequence of step faults. Using synthetic data, we performed a systematic search of normal fault slips versus fault displacement depths to map the fault-detectable region in this space. This mapping helps to assess the ability of gravity inversion to detect normal faults. Mapping shows that normal faults with small (<0.5km), medium (about 1km), and large (about 2km) vertical slips can be detected if the maximum midpoint depths of the fault planes are smaller than 1.8, 3.8, and 6.8km, respectively.
This study presents the results from the interpretation of an extensive and recent regional two-dimensional seismic survey focused on the understanding of the salt tectonics in the western Gulf of Cadiz (GoC). Two different salt units were identified: an autochthonous salt unit of the Late Triassic or the Early Jurassic (Hettangian) and an allochthonous unit that originated from the Hettangian salt. Interpretation of the pattern of distribution of the salt in the basin allowed subdivision of the area of study into three distinct salt domains: the eastern domain characterized by the presence of a conspicuous allochthonous salt nappe (Esperança Salt), the central domain dominated by salt diapirs with mild deformation of Miocene strata and wide salt-withdrawal minibasins, and the southwestern domain where present-day tectonics induces impressive salt deformation affecting the sea floor. This complex pattern is mainly the result of the interaction of inherited basement structure, complex tectonic history, and stress regime of the basin. The intense halokinesis observed has created several salt-related trap geometries and fluid migration pathways. As the focus of worldwide exploration along passive margins is gradually shifting to deep-water regions, the western GoC has the potential to become a deep-water petroleum province in the near future.
The marine controlled-source electromagnetic (CSEM) method provides complementary information to seismic imaging in the exploration of sedimentary basins. The CSEM is mainly used for reservoir scanning and appraisal. The CSEM interpretation workflow is heavily based on inversion and forward - modeling for hypothesis testing. Until the recent past, the effectiveness of a given workflow was achieved after the drilling results, as there wasn’t any geological complex model available to serve as a benchmark. In the present paper, we describe the workflow to build up Marlim R3D, a realistic and complex geoelectric model. Marlim R3D aims to be a reference model of turbidite reservoirs of the Brazilian continental margin. Our model is based on seismic interpretation and constrained by the input of available well-log information. The workflow used is composed of seven steps: seismic and well-log dataset loading, well-tie, Vp cube construction, Vp resistivity calibration, time-depth conversion, resistivity cube construction, and quality-control check. As a result, we obtained an interpreted dataset composed by main stratigraphic horizons, pseudo-well logs, and the resistivity cubes. These elements were made freely available for research or commercial use, under the Creative Common License, at the Zenodo platform.
Synthetic data provided by geoelectric earth models are a powerful tool to evaluate a priori a controlled-source electromagnetic (CSEM) workflow effectiveness. Marlim R3D (MR3D) is an open-source complex and realistic geoelectric model for CSEM simulations of the postsalt turbiditic reservoirs at the Brazilian offshore margin. We have developed a 3D CSEM finite-difference time-domain forward study to generate the full-azimuth CSEM data set for the MR3D earth model. To that end, we fabricated a full-azimuth survey with 45 towlines striking the north–south and east–west directions over a total of 500 receivers evenly spaced at 1 km intervals along the rugged seafloor of the MR3D model. To correctly represent the thin, disconnected, and complex geometries of the studied reservoirs, we have built a finely discretized mesh of 100×100×20 m cells leading to a large mesh with a total of approximately 90 million cells. We computed the six electromagnetic field components (Ex, Ey, Ez, Hx, Hy, and Hz) at six frequencies in the range of 0.125–1.25 Hz. In our efforts to mimic noise in real CSEM data, we summed to the data a multiplicative noise with a 1% standard deviation. Both CSEM data sets (noise free and noise added), with inline and broadside geometries, are distributed for research or commercial use, under the Creative Common License, at the Zenodo platform.