The regularizing Levenberg–Marquardt scheme for history matching of petroleum reservoirs     

Marco A. Iglesias  Clint Dawson 《Computational Geosciences》2013,17(6):1033-1053

In this paper, we study on a history matching approach that consists of finding stable approximations to the problem of minimizing the weighted least-squares functional that penalizes the misfit between the reservoir model predictions G(u) and noisy observations y ^η . In other words, we are interested in computing an approximation to the minimizer of $\frac {1}{2}\vert \vert \Gamma ^{-1/2}(y^{\eta }-G(u))\vert \vert _{Y}^{2} $ where Γ is the measurements error covariance, Y is the observation space, and X is a set of admissible parameters. This is an ill-posed nonlinear inverse problem that we address by means of the regularizing Levenberg–Marquardt scheme developed by Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010). Under certain conditions on G, the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) ensures the convergence of the scheme to stable approximations to the inverse problem. We propose an implementation of the regularizing Levenberg–Marquardt scheme that enforces prior knowledge on the geologic properties. In particular, the prior mean $\overline {u}$ is incorporated in the initial guess of the algorithm, and the prior error covariance C is enforced through the definition of the parameter space X. Our main goal is to numerically show that the proposed implementation of the regularizing Levenberg–Marquardt scheme of Hanke is a robust method capable of providing accurate estimates of the geologic properties for small noise measurements. In addition, we provide numerical evidence of the convergence and regularizing results predicted by the theory of Hanke (Inverse Probl. 13:79–95, 1997; J. Integr. Equ. Appl. 22(2):259–283, 2010) for a prototypical oil–water reservoir model. The performance for recovering the true permeability with the regularizing Levenberg–Marquardt scheme is compared to the typical approach of computing the maximum a posteriori (MAP) estimator. In particular, we compare the proposed application of the regularizing Levenberg–Marquardt (LM) scheme against the standard LM approach of Li et al. (SPE J. 8(4):328–340, 2003) and Reynolds et al. (2008) for computing the MAP. Our numerical experiments suggest that the history matching approach based on iterative regularization is robust and could potentially be used to improve further on various methodologies already proposed as effective tools for history matching in petroleum reservoirs.  相似文献   

Combined geological and gravimetric mapping and modelling for an improved understanding of observed high-resolution gravity variations: a case study for the Global Geodynamics Project (GGP) station Moxa, Germany     

Norbert Kasch  Marco Naujoks  Jonas Kley  Thomas Jahr 《International Journal of Earth Sciences》2013,102(5):1257-1270

Geodynamic observatories around the globe continuously monitor signals like gravity, tilt and strain as a function of time. However, global signals are often masked by local effects, caused by the direct surroundings of the station, including the local geological setting. This link is well established for superconducting gravimeters (SG) that observe the gravity field variations at very high resolution. An enhancement of the SG time series by the application of local correction is exemplarily shown here in a very practicable procedure for the Geodynamic Observatory Moxa, Germany. We show how the combination of geological and gravimetrical mapping and modelling around Moxa results in a significant correction of the original gravity data, as can be proven by comparison with the satellite-derived gravity field. Detailed geological mapping of the observatory surroundings, including measurements of fold axes, foliations and joints, was the basis of the present study. The fold structure in the area of interest was interpreted by geometrical 3D modelling. The complete representation of different rock types in space also provides a better understanding of the local hydro-geological situation. Using a combination of the 3D geological model with the high-resolution Bouguer map of the observatory surroundings, we developed a 3D density model. This model enables the correction of small gravity effects caused by mass variations close to the SG, in particular local hydro-geological mass changes. Thus, the procedure presented here, based on a close connection of geology and gravimetry, is a powerful tool for the reduction of local gravity effects on SG recordings. It should be applicable to SG stations worldwide, where similar hydrologically driven mass changes can be assumed.  相似文献   

Geodynamic evolution of ophiolites from Albania and Greece (Dinaric-Hellenic belt): one, two, or more oceanic basins?   总被引：1，自引：0，他引：1  

Valerio Bortolotti  Marco Chiari  Michele Marroni  Luca Pandolfi  Gianfranco Principi  Emilio Saccani 《International Journal of Earth Sciences》2013,102(3):783-811

All the geological constraints for an exhaustive reconstruction of the Triassic to Tertiary tectonic history of the southern Dinaric-Hellenic belt can be found in Albania and Greece. This article aims to schematically reconstruct this long tectonic evolution primarily based on a detailed analysis of the tectonic setting, the stratigraphy, the geochemistry, and the age of the ophiolites. In contrast to what was previously reported in the literature, we propose a new subdivision on a regional scale of the ophiolite complexes cropping out in Albania and Greece. This new subdivision includes six types of ophiolite occurrences, each corresponding to different tectonic units derived from a single obducted sheet. These units are represented by: (1) sub-ophiolite mélange, (2) Triassic ocean-floor ophiolites, (3) metamorphic soles, (4) Jurassic fore-arc ophiolites, (5) Jurassic intra-oceanic-arc ophiolites, and (6) Jurassic back-arc basin ophiolites. The overall features of these ophiolites are coherent with the existence of a single, though composite, oceanic basin located east of the Adria/Pelagonian continental margin. This oceanic basin was originated during the Middle Triassic and was subsequently (Early Jurassic) affected by an east-dipping intra-oceanic subduction. This subduction was responsible for the birth of intra-oceanic-arc and back-arc oceanic basins separated by a continental volcanic arc during the Early to Middle Jurassic. From the uppermost Middle Jurassic to the Early Cretaceous, an obduction developed, during which the ophiolites were thrust westwards firstly onto the neighboring oceanic lithosphere and then onto the Adria margin.  相似文献   

The slope aspect: A predisposing factor for landsliding?     

Marco Capitani  Adriano Ribolini  Monica Bini 《Comptes Rendus Geoscience》2013,345(11-12):427-438

The influence of slope aspect on the distribution of landslides was studied in the Milia and Roglio basins in Tuscany, Italy. For each basin, the new Tuscany region landslide inventory that was initiated in 2010 was used. The landslides were split into separate datasets based on their prevailing movement typology. To assess the results that were obtained from the different slope aspect values, maps of the lithology, slope angle, distances to streams, and distances to tectonic lineaments were included in the bivariate statistical analysis as comparison terms. For each basin, all of the geo-environmental factor maps were compared with the different landslide typologies with GIS software. Pearson's Chi² (χ²) coefficient was used to test the degree of spatial association between each predictor variable and landslide type. In addition, Cramer's V test was used to quantify the strength of the degree of association. Next, a conditional analysis was applied to all of the possible combinations that occurred between the slope aspect and other landslide-predisposing factors. Overall, the slope aspect significantly affected the distribution of superficial landslide types, but apparently not that of other landslide types.  相似文献   

Distinction between amorphous and healed planar deformation features in shocked quartz using composite color scanning electron microscope cathodoluminescence (SEM‐CL) imaging          下载免费PDF全文

Maartje F. Hamers  Gill M. Pennock  Marco Herwegh  Martyn R. Drury 《Meteoritics & planetary science》2016,51(10):1914-1931

Planar deformation features (PDFs) in quartz are one of the most reliable and most widely used forms of evidence for hypervelocity impact. PDFs can be identified in scanning electron microscope cathodoluminescence (SEM‐CL) images, but not all PDFs show the same CL behavior: there are nonluminescent and red luminescent PDFs. This study aims to explain the origin of the different CL emissions in PDFs. Focused ion beam (FIB) thin foils were prepared of specific sample locations selected in composite color SEM‐CL images and were analyzed in a transmission electron microscope (TEM). The FIB preparation technique allowed a direct, often one‐to‐one correlation between the CL images and the defect structure observed in TEM. This correlation shows that composite color SEM‐CL imaging allows distinction between amorphous PDFs on one hand and healed PDFs and basal Brazil twins on the other: nonluminescent PDFs are amorphous, while healed PDFs and basal Brazil twins are red luminescent, with a dominant emission peak at 650 nm. We suggest that the red luminescence is the result of preferential beam damage along dislocations, fluid inclusions, and twin boundaries. Furthermore, a high‐pressure phase (possibly stishovite) in PDFs can be detected in color SEM‐CL images by its blue luminescence.  相似文献   

Endogenous and nonimpact origin of the Arkenu circular structures (al‐Kufrah basin—SE Libya)     

Corrado CIGOLINI  Mario DI MARTINO  Marco LAIOLO  Diego COPPOLA  Piergiorgio ROSSETTI  Marco MORELLI 《Meteoritics & planetary science》2012,47(11):1772-1788

Abstract– The twin Arkenu circular structures (ACS), located in the al‐Kufrah basin in southeastern Libya, were previously considered as double impact craters (the “Arkenu craters”). The ACS consist of a NE (Arkenu 1) and a SW structure (Arkenu 2), with approximate diameters of about 10 km. They are characterized by two shallow depressions surrounded by concentric circular ridges and silica‐impregnated sedimentary dikes cut by local faults. Our field, petrographic, and textural observations exclude that the ACS have an impact origin. In fact, we did not observe any evidence of shock metamorphism, such as planar deformation features in the quartz grains of the collected samples, and the previously reported “shatter cones” are wind‐erosion features in sandstones (ventifacts). Conversely, the ACS should be regarded as a “paired” intrusion of porphyritic stocks of syenitic composition that inject the Nubia Formation and form a rather simple and eroded ring dike complex. Stock emplacement was followed by hydrothermal activity that involved the deposition of massive magnetite–hematite horizons (typical of iron oxide copper‐gold deposits). Their origin was nearly coeval with the development of silicified dikes in the surroundings. Plugs of tephritic‐phonolitic rocks and lamprophyres (monchiquites) inject the Nubian sandstone along conjugate fracture zones, trending NNW–SSE and NE–SW, that crosscut the structural axis of the basin.  相似文献   

Millennial-scale shifts in microtidal ecosystems during the Holocene: dynamics and drivers of change from the Po Plain coastal record (NE Italy)     

Veronica Rossi  Giulia Barbieri  Stefano Claudio Vaiani  Marco Cacciari  Luigi Bruno  Bruno Campo  Marco Marchesini  Silvia Marvelli  Alessandro Amorosi 《第四纪科学杂志》2021,36(6):961-979

Framed into a robust stratigraphic context, multivariate analyses on the Holocene palaeobiological record (pollen, benthic foraminifers, ostracods) of the Po coastal plain (NE Italy) allowed the investigation of microtidal ecosystems variability and driving parameters along a 35-km-long land–sea transect. Millennial-scale ecosystem shifts are documented by coeval changes in the meiofauna, reflecting variations in organic matter–water depth (shallow-marine environments) and degree of confinement-salinity (back-barrier settings). In-phase shifts of vegetation communities track unsteady water-table levels and river dynamics in freshwater palustrine areas. Five environmental–ecological stages followed one another crossing four tipping points that mark changes in relative sea level (RSL), climate and/or fluvial regime. At the culmination of Mediterranean RSL rise, after the 8200 event, remarkable growth of peatlands took place in the Po estuary, while low accumulation rates typified the shelf. At the transgressive–regressive turnaround (~7000 cal a bp ), the estuary turned into a delta plain with tidally influenced interdistributary embayments. River flow regime oscillations after the Climate Optimum (post-5000 cal a bp ) favoured isolation of the bays and the development of brackish wetlands surrounded by wooded peatlands. The youngest threshold (~800 cal a bp ), which led to the establishment of the modern delta, reflects a major avulsion of the Po River.  相似文献   

Towards a better understanding of the evolution of the flood risk in Mediterranean urban areas: the case of Barcelona     

Maria Cortès  Maria Carmen Llasat  Joan Gilabert  Montserrat Llasat-Botija  Marco Turco  Raül Marcos  Juan Pedro Martín Vide  Lluís Falcón 《Natural Hazards》2018,93(1):39-60

This contribution explores the evolution of the flood risk in the Metropolitan Area of Barcelona (MAB; Northeast Spain) from 1981 to 2015, and how it has been affected by changes in land use, population and precipitation. To complete this study, we analysed PRESSGAMA and INUNGAMA databases to look for all the information related to the floods and flash floods that have affected the chosen region. The “Consorcio de Compensación de Seguros”, a state insurance company for extraordinary risks, provided data on economic damage. The extreme precipitation trend was analysed by the Fabra Observatory and El Prat-Airport Observatory, and daily precipitation data were provided by the State Meteorological Agency of Spain (AEMET) and the Meteorological Service of Catalonia (SMC). Population data were obtained from the Statistical Institute of Catalonia (IDESCAT). Changes in land use were estimated from the land use maps for Catalonia corresponding to 1956, 1993, 2000, 2005 and 2009. Prevention measures like rainwater tanks and improvements to the drainage system were also been considered. The specific case of Barcelona is presented, a city recognised by United Nations International Strategy for Disaster Reduction as a model city for urban resilience to floods. The evolution of flood events in the MAB does not show any significant trend for this period. We argue that the evolution in floods can be explained, at least in part, by the lack of trend in extreme precipitation indices, and also by the improvements in flood prevention measures.  相似文献   

The onset of the Messinian salinity crisis in the deep Eastern Mediterranean basin   总被引：1，自引：0，他引：1       下载免费PDF全文

Vinicio Manzi  Rocco Gennari  Stefano Lugli  Davide Persico  Matteo Reghizzi  Marco Roveri  B. Charlotte Schreiber  Ran Calvo  Ittai Gavrieli  Zohar Gvirtzman 《地学学报》2018,30(3):189-198

Astronomical tuning of the Messinian pre‐salt succession in the Levant Basin allows for the first time the reconstruction of a detailed chronology of the Messinian salinity crisis (MSC) events in deep setting and their correlation with marginal records that supports the CIESM ( 2008 ) 3‐stage model. Our main conclusions are (1) MSC events were synchronous across marginal and deep basins, (2) MSC onset in deep basins occurred at 5.97 Ma, (3) only foraminifera‐barren, evaporite‐free shales accumulated in deep settings between 5.97 and 5.60 Ma, (4) deep evaporites (anhydrite and halite) deposition started later, at 5.60 Ma and (5) new and published ⁸⁷Sr/⁸⁶Sr data indicate that during all stages, evaporites precipitated from the same water body in all the Mediterranean sub‐basins. The wide synchrony of events and ⁸⁷Sr/⁸⁶Sr homogeneity implies inter‐sub‐basin connection during the whole MSC and is not compatible with large sea‐level fall and desiccation of the Mediterranean.  相似文献   

Modeling groundwater/surface-water interactions in an Alpine valley (the Aosta Plain,NW Italy): the effect of groundwater abstraction on surface-water resources     

Gennaro A. Stefania  Marco Rotiroti  Letizia Fumagalli  Fulvio Simonetto  Pietro Capodaglio  Chiara Zanotti  Tullia Bonomi 《Hydrogeology Journal》2018,26(1):147-162

A groundwater flow model of the Alpine valley aquifer in the Aosta Plain (NW Italy) showed that well pumping can induce river streamflow depletions as a function of well location. Analysis of the water budget showed that ～80% of the water pumped during 2 years by a selected well in the downstream area comes from the baseflow of the main river discharge. Alluvial aquifers hosted in Alpine valleys fall within a particular hydrogeological context where groundwater/surface-water relationships change from upstream to downstream as well as seasonally. A transient groundwater model using MODFLOW2005 and the Streamflow-Routing (SFR2) Package is here presented, aimed at investigating water exchanges between the main regional river (Dora Baltea River, a left-hand tributary of the Po River), its tributaries and the underlying shallow aquifer, which is affected by seasonal oscillations. The three-dimensional distribution of the hydraulic conductivity of the aquifer was obtained by means of a specific coding system within the database TANGRAM. Both head and flux targets were used to perform the model calibration using PEST. Results showed that the fluctuations of the water table play an important role in groundwater/surface-water interconnections. In upstream areas, groundwater is recharged by water leaking through the riverbed and the well abstraction component of the water budget changes as a function of the hydraulic conditions of the aquifer. In downstream areas, groundwater is drained by the river and most of the water pumped by wells comes from the base flow component of the river discharge.  相似文献