Iterated extended Kalman filter method for time-lapse seismic full-waveform inversion     

Kjersti Solberg Eikrem  Geir Nævdal  Morten Jakobsen 《Geophysical Prospecting》2019,67(2):379-394

Time-lapse seismic data is useful for identifying fluid movement and pressure and saturation changes in a petroleum reservoir and for monitoring of CO₂ injection. The focus of this paper is estimation of time-lapse changes with uncertainty quantification using full-waveform inversion. The purpose of also estimating the uncertainty in the inverted parameters is to be able to use the inverted seismic data quantitatively for updating reservoir models with ensemble-based methods. We perform Bayesian inversion of seismic waveform data in the frequency domain by combining an iterated extended Kalman filter with an explicit representation of the sensitivity matrix in terms of Green functions (acoustic approximation). Using this method, we test different strategies for inversion of the time-lapse seismic data with uncertainty. We compare the results from a sequential strategy (making a prior from the monitor survey using the inverted baseline survey) with a double difference strategy (inverting the difference between the monitor and baseline data). We apply the methods to a subset of the Marmousi2 P-velocity model. Both strategies performed well and relatively good estimates of the monitor velocities and the time-lapse differences were obtained. For the estimated time-lapse differences, the double difference strategy gave the lowest errors.  相似文献   

Noble gas record and cosmic‐ray exposure history of the new CO3 chondrite Moss––Comparison with Lancé and other CO chondrite falls     

Rainer BARTOSCHEWITZ  Ulrich OTT  Luitgard FRANKE  Siegfried HERRMANN  Yukio YAMAMOTO  Keisuke NAGAO  Morten BILET  Thomas GRAU 《Meteoritics & planetary science》2010,45(8):1380-1391

Abstract– The Moss meteorite is the first CO chondrite fall after a time period of 70 yr and the least terrestrially contaminated member of its group. Its cosmic‐ray exposure (CRE) age (T₃ ～ 13.5 Ma; T₂₁ ～ 14.6 Ma) is distinct among CO chondrites and, within witnessed falls is the shortest after Lancé, which we have reanalyzed. Gas retention ages are approximately 3.95 × 10⁹ yr (U/Th‐He) and approximately 4.43 × 10⁹ yr (K/Ar), respectively. Trapped Ar, Kr, and Xe are present in Moss in abundances typical for CO chondrites, with “planetary” elemental and isotopic compositions. Presence of HL‐xenon from presolar diamonds is observed in the stepwise release analysis of Lancé. It may also be present in Moss, but it is difficult to ascertain in single‐step bulk analyses. It follows from our new data combined with a survey of the literature that the abundance of trapped gases in CO chondrites is not a good indicator of their petrological subtype.  相似文献   

Weathering rates from top to bottom in a carbonate environment     

Bernard Bourdon  Sarah Bureau  Morten B. Andersen  Eric Pili  Amélie Hubert 《Chemical Geology》2009,258(3-4):275-287

This study investigates U-series, Sr isotopes, major and trace elements in a chalk aquifer system located in Eastern France. Soil and rock samples were collected along depth profiles down to 45 m in four localities as an attempt to investigate the weathering processes in the soil, the unsaturated zone and the saturated zone of the aquifer. Interstitial water was extracted from soils and rocks by a centrifugation technique. U-series offer a powerful tool to calculate weathering rates because the relative mobility of the U- and Th-isotopes can be precisely measured and it does not require the determination of a reference state as in other approaches. As expected, the data show very large mobile element depletion in the soil with large ²³⁰Th excess relative to ²³⁸U, while the rocks show more limited but not insignificant mobile element depletion. The U-series data have been used to constrain weathering rates based on a 1-D reactive transport model. Weathering rates in the near surface are about 10–100 times faster than at depth. However, when integrated over the depth of the cores, including the unsaturated and the saturated zones, this underground weathering represents more than 30% of the total weathering flux, assuming congruent dissolution of carbonates. The (²³⁴U/²³⁸U) ratios in interstitial water are consistent with solid samples showing ²³⁴U depletion near the surface and an excess ²³⁴U at depth. A leaching experiment performed on chalk shows that the excess ²³⁴U in natural waters percolating through carbonate rocks results both from preferential ²³⁴U leaching and direct recoil in the interstitial water. A new approach was used to derive the recoil ejection factor based on BET measurements and the fractal dimension of chalk surface. Consideration of preferential leaching and recoil allows a more accurate modeling of weathering rates.  相似文献   

Experimental evidence for U-U fractionation during granite weathering with implications for U/U in natural waters     

Morten B. Andersen  Yigal Erel 《Geochimica et cosmochimica acta》2009,73(14):4124-136

The daughter to parent (²³⁴U/²³⁸U) activity ratio in natural waters is often out of secular radioactive equilibrium. The major reason for this disequilibrium is related to the energetic α-decay of ²³⁸U and differential release of ²³⁴U relative to ²³⁸U. This disequilibrium originates from (1) preferential release of more loosely bound ²³⁴U from damaged mineral lattice sites or; (2) direct recoil of ²³⁴Th into surrounding media from near mineral surface boundaries, however, it is unclear which of the two mechanisms is most important in nature. To better quantify the effects of preferential release of ²³⁴U, two continuous laboratory granite leaching experiments conducted over 1100 h were performed. The leachates were characterized by declining U concentrations with time and (²³⁴U/²³⁸U) initially greater than unity (up to 1.15), which changed to below unity during leaching (∼0.95). The early elevated (²³⁴U/²³⁸U) suggests that additional ²³⁴U is released into solution by preferential release of ²³⁴U from mineral phases. However, the excess ²³⁴U constitutes a finite pool of easy leachable ²³⁴U and the (²³⁴U/²³⁸U) values become lower than unity when this pool is used up. A model based on first-order kinetics, dissolution rates and preferential release of ²³⁴U from damaged lattice sites was developed and is able to quantitatively predict the observed pattern of (²³⁴U/²³⁸U) values and U concentrations for the two granite leaching experiments. Extending the modeling to longer time scales more comparable to natural systems shows that the production of waters with high (²³⁴U/²³⁸U) ratios can be achieved in two distinct regimes (1) slow weathering where the rate of directly recoiled ²³⁴U near mineral surfaces into waters is high; (2) fast weathering where the role of incipient chemical weathering and preferential release of loosely bound ²³⁴U are important. The model is able to explain apparent opposite correlations between physical erosion rates and (²³⁴U/²³⁸U) in waters and it provides a new framework that will be useful for examining weathering regimes, their timescales and their coupling with physical erosion.  相似文献   

Megacrysts and mafic–ultramafic xenolith-bearing ignimbrites from Sirwa Volcano, Morocco: phase petrology and thermobarometry     

M. Bondi  L. Morten  P. Nimis  P. L. Rossi  C. A. Tranne 《Mineralogy and Petrology》2002,75(3-4):203-221

Summary A suite of clinopyroxene and amphibole megacrysts and mafic–ultramafic xenoliths are present in ignimbritic rocks of trachybasaltic–andesitic composition from the Sirwa volcanic district, Morocco. The stumpy prismatic and sometimes euhedral clinopyroxene megacrysts are Ti–Al-rich diopsides with mg values in the range 0.82–0.87 and Ca/(Ca + Mg) ratios in the range 0.53–0.54. The prismatic, elongated amphibole megacrysts are calcic kaersutites–kaersutites with a narrow mg range (0.66–0.68). The xenoliths are represented by gabbroic and pyroxenitic types. In the gabbroic xenoliths two distinct textural types can be distinguished: medium-sized granular and banded. The granular type is characterized by the mineral assemblage Pl + Amph + Spl + Ilm + Ap. The banded type is distinct for the absence of Ilm and the presence of Cpx and Opx and shows alternating bands enriched in Pl and Amph, respectively. The megacrysts and, probably, the xenoliths are considered not cognate with the present host rocks since the calculated liquids in equilibrium with clinopyroxene and amphibole megacrysts over a wide range of physical conditions have different trace and rare earth element contents. The observed phase relations and thermobarometric calculations indicate that the megacrysts and xenoliths crystallized from their parent melts at P ≥ 10 kbar and T ≤ 1160 °C, i.e., in the upper mantle or near the crust-mantle boundary. A deep ( ≥ 30 km) magmatic chamber, where the megacrysts and xenoliths originated, and a shallow volcanic chamber, energetically activated up to explosive conditions by injection of deep-originated melts, is suggested to explain the occurrence of high-pressure megacrysts and xenoliths in the Sirwa volcanic explosive products. Received October 8, 2000; revised version accepted September 9, 2001  相似文献   

Late Holocene palaeoceanography in Van Mijenfjorden, Svalbard   总被引：2，自引：0，他引：2  

Morten Hald  Torbjörn Dahlgren  Tor-Eirik Olsen  Erland Lebesbye 《Polar research》2001,20(1):23-35

Based on detailed stratigraphical analysis of sediment cores spanning the last ca. 4000 calendar years, we reconstruct the palaeoceanograhic changes in the fiord Van Mijenfjorden, western Svalbard. Benthic foraminiferal δ¹⁸O indicate a gradual reduction in bottom water salinities between 2200 BC and 500 BC. This reduction was probably mainly a function of reduced inflow of oceanic water to the fiord, due to isostatic shallowing of the outer fiord sill. Stable salinity conditions prevailed between 500 BC and. 1300 AD. After the onset of a major glacial surge of the tidewater Paulabreen (Paula Glacier) system (PGS) around 1300 AD, there was a foraminiferal faunal change towards glacier proximal conditions, associated with a slight bottom water salinity depletion. During a series of glacial surges occuring from 1300 AD up the present salinity in the fiord has further decreased, corresponding to a δ¹⁸O depletion of 0.5 %o. This salinity decrease corresponds to the period when the PGS lost an equivalent of 30 – 40 % of its present ice volume, mainly through calving in the fiord.  相似文献   

An ecosystem modeling approach to predicting cod recruitment     

Einar Svendsen  Morten Skogen  Paul Budgell  Geir Huse  Jan Erik Stiansen  Bjrn dlandsvik  Frode Vikeb  Lars Asplin  Svein Sundby 《Deep Sea Research Part II: Topical Studies in Oceanography》2007,54(23-26):2810

The Norwegian Ecological Model (NORWECOM) biophysical model system implemented with the ROMS ocean circulation model has been run to simulate conditions over the last 25 years for the North Atlantic. Modeled time series of water volume fluxes, primary production, and drift of cod larvae through their modeled ambient temperature fields have been analyzed in conjunction with VPA estimated time series of 3-year-old cod recruits in the Barents Sea. Individual time series account for less than 50% of the recruitment variability; however, a combination of simulated flow of Atlantic water into the Barents Sea and local primary production accounts for 70% of the variability with a 3-year lead. The associated regression predicts increased recruitment between 2007 and 2008 from about 450–700 million individuals with a standard error of nearly 150 million.  相似文献   

Kinematics of extreme waves in deep water   总被引：2，自引：0，他引：2  

John Grue  Didier Clamond  Morten Huseby  Atle Jensen 《Applied Ocean Research》2003,25(6):355-366

The velocity profiles under crest of a total of 62 different steep wave events in deep water are measured in laboratory using particle image velocimetry. The waves take place in the leading unsteady part of a wave train, focusing wave fields and random wave series. Complementary fully nonlinear theoretical/numerical wave computations are performed. The experimental velocities have been put on a nondimensional form in the following way: from the wave record (at a fixed point) the (local) trough-to-trough period, T_TT and the maximal elevation above mean water level, η_m of an individual steep wave event are identified. The local wavenumber, k and an estimate of the wave slope, ε are evaluated from ω²/(gk)=1+ε², where ω=2π/T_TT and g denotes the acceleration of gravity. A reference fluid velocity, is then defined. Deep water waves with a fluid velocity up to 75% of the estimated wave speed are measured. The corresponding kη_m is 0.62. A strong collapse of the nondimensional experimental velocity profiles is found. This is also true with the fully nonlinear computations of transient waves. There is excellent agreement between the present measurements and previously published Laser Doppler Anemometry data. A surprising result, obtained by comparison, is that the nondimensional experimental velocities fit with the exponential profile, i.e. e^ky, y the vertical coordinate, with y=0 in the mean water level.  相似文献   

Petrogenesis of Cogenetic Silica-Oversaturated and -Undersaturated Syenites by Periodic Recharge in a Crustally Contaminated Magma Chamber: the Kangerlussuaq Intrusion, East Greenland     

Riishuus  Morten S.; Peate  David W.; Tegner  Christian; Wilson  J. Richard; Brooks  C. Kent 《Journal of Petrology》2008,49(3):493-522

The Palaeogene Kangerlussuaq Intrusion (50 Ma) of East Greenlanddisplays concentric zonation from quartz-rich nordmarkite (quartzsyenite) at the margin, through pulaskite, to foyaite (nephelinesyenite) in the centre; modal layering and igneous laminationare locally developed but there are no internal intrusive contacts.This is an apparent violation of the phase relations in Petrogeny'sResidua System. We propose that this intrusion is layered, gradingfrom quartz syenite at the bottom to nepheline syenite at thetop. Mineral and whole-rock major and trace element data andSr–Nd–Hf–Pb isotope data are presented thatprovide constraints on the petrogenesis of the intrusion. Radiogenicisotope data indicate a continuously decreasing crustal componentfrom the quartz nordmarkites (⁸⁷Sr/⁸⁶Sr = 0·7061; _Ndi= 2·3; _Hfi = 5·2; ²⁰⁶Pb/²⁰⁴Pb_meas = 16·98)to the foyaites (⁸⁷Sr/⁸⁶Sr = 0·7043–0·7044;_Ndi = 3·8–4·9; _Hfi = 10·7–11·1;²⁰⁶Pb/²⁰⁴Pb_meas = 17·78–17·88); the foyaitesare dominated by a mantle isotopic signature. The average Mg-numberof amphibole cores becomes increasingly primitive, varying from26·4 in the nordmarkites to 57·4 in the pulaskites.Modal layering, feldspar lamination and the presence of hugebasaltic xenoliths derived from the chamber roof, now restingon the transient chamber floor, demonstrate bottom-upwards crystallization.The intrusion cannot, therefore, have formed in a system closedto magmatic recharge. The lack of gneissic xenoliths in thenordmarkites suggests that most contamination took place deeperin the crust. In the proposed model, the nordmarkitic magmaformed during crustal assimilation in the roof zone of a large,silica-undersaturated alkali basaltic/basanitic, stratifiedmagma chamber, prior to emplacement in the uppermost crust.The more primitive syenites, terminating with foyaite at thetop of the intrusion, formed as a consequence of repeated rechargeof the Kangerlussuaq Intrusion magma chamber by tapping lesscontaminated, more primitive phonolitic melt from deeper partsof the underlying chamber during progressive armouring of theplumbing system. KEY WORDS: Kangerlussuaq; East Greenland; syenite; crustal contamination; magma mixing  相似文献   

Evaluation of probabilistic flow predictions in sewer systems using grey box models and a skill score criterion   总被引：1，自引：1，他引：0  

Fannar ?rn Thordarson  Anders Breinholt  Jan Kloppenborg M?ller  Peter Steen Mikkelsen  Morten Grum  Henrik Madsen 《Stochastic Environmental Research and Risk Assessment (SERRA)》2012,26(8):1151-1162

  相似文献