Inverse analysis of stochastic moment equations for transient flow in randomly heterogeneous media     

Monica Riva  Alberto Guadagnini  Shlomo P. Neuman  Emanuela Bianchi Janetti  Bwalya Malama 《Advances in water resources》2009

We present a nonlinear stochastic inverse algorithm that allows conditioning estimates of transient hydraulic heads, fluxes and their associated uncertainty on information about hydraulic conductivity (K) and hydraulic head (h  ) data collected in a randomly heterogeneous confined aquifer. Our algorithm is based on Laplace-transformed recursive finite-element approximations of exact nonlocal first and second conditional stochastic moment equations of transient flow. It makes it possible to estimate jointly spatial variations in natural log-conductivity (Y=lnK)$(Y=\ln K)$, the parameters of its underlying variogram, and the variance–covariance of these estimates. Log-conductivity is parameterized geostatistically based on measured values at discrete locations and unknown values at discrete “pilot points”. Whereas prior values of Y at pilot point are obtained by generalized kriging, posterior estimates at pilot points are obtained through a maximum likelihood fit of computed and measured transient heads. These posterior estimates are then projected onto the computational grid by kriging. Optionally, the maximum likelihood function may include a regularization term reflecting prior information about Y. The relative weight assigned to this term is evaluated separately from other model parameters to avoid bias and instability. We illustrate and explore our algorithm by means of a synthetic example involving a pumping well. We find that whereas Y and h can be reproduced quite well with parameters estimated on the basis of zero-order mean flow equations, all model quality criteria identify the second-order results as being superior to zero-order results. Identifying the weight of the regularization term and variogram parameters can be done with much lesser ambiguity based on second- than on zero-order results. A second-order model is required to compute predictive error variances of hydraulic head (and flux) a posteriori. Conditioning the inversion jointly on conductivity and hydraulic head data results in lesser predictive uncertainty than conditioning on conductivity or head data alone.  相似文献   

Joint Inversion of Body-Wave Arrival Times and Surface-Wave Dispersion for Three-Dimensional Seismic Structure Around SAFOD   总被引：2，自引：0，他引：2  

Haijiang Zhang  Monica Maceira  Philippe Roux  Clifford Thurber 《Pure and Applied Geophysics》2014,171(11):3013-3022

We incorporate body-wave arrival time and surface-wave dispersion data into a joint inversion for three-dimensional P-wave and S-wave velocity structure of the crust surrounding the site of the San Andreas Fault Observatory at Depth. The contributions of the two data types to the inversion are controlled by the relative weighting of the respective equations. We find that the trade-off between fitting the two data types, controlled by the weighting, defines a clear optimal solution. Varying the weighting away from the optimal point leads to sharp increases in misfit for one data type with only modest reduction in misfit for the other data type. All the acceptable solutions yield structures with similar primary features, but the smaller-scale features change substantially. When there is a lower relative weight on the surface-wave data, it appears that the solution over-fits the body-wave data, leading to a relatively rough V _s model, whereas for the optimal weighting, we obtain a relatively smooth model that is able to fit both the body-wave and surface-wave observations adequately.  相似文献   

Inorganic synthesis of Fe–Ca–Mg carbonates at low temperature     

Christopher S. Romanek  Concepcin Jimnez-Lpez  Alejandro Rodriguez Navarro  Monica Snchez-Romn  Nita Sahai  Max Coleman 《Geochimica et cosmochimica acta》2009,73(18):5361-5376

A set of free-drift experiments was undertaken to synthesize carbonates of mixed cation content (Fe, Ca, Mg) from solution at 25 and 70 °C to better understand the relationship between the mineralogy and composition of these phases and the solutions from which they precipitate. Metastable solid solutions formed at 25 °C which are not predicted from the extrapolation of higher temperature equilibrium assemblages; instead, solids formed that were intermediary in chemical composition to known magnesite–siderite and dolomite solid solutions. A calcite–siderite solid solution precipitated at 25 °C, with the percentage of CaCO₃ in the solid being proportional to the aqueous Ca/Fe ratio of the solution, while Mg was excluded from the crystal structure except at relatively high aqueous Mg/Ca and Mg/Fe ratios and a low Ca content. Alternatively, at 70 °C Mg was the predominant cation of the solid solutions. These results are compatible with the hypothesis that the relative dehydration energies of Fe, Ca and Mg play an important role in the formation of mixed cation carbonates in nature.  相似文献   

地磁数据分析的一些思考   总被引：3，自引：1，他引：2       下载免费PDF全文

Laurentiu ASIMOPOLOS  Agata Monica PESTINA  Natalia Silvia ASIMOPOLOS 《地球物理学报》2010,53(3):765-772

通过若干实例,本文讨论了地磁数据分析中的一些方法.这些方法包括时间序列的多项式拟合,相关性分析,频谱分析和小波分析等.在Matlab和AutoSignal软件下实现了以上方法,并基于这些方法研究了地磁场形态和地磁现象的频谱特征.我们分析了Surlari观测台（余纬45.3°,经度26.3°）的地磁时间序列,并与Ottawa（余纬44.6°,〖JP2〗经度284.5°）,Canberra（余纬125.3°,经度149.3°）, Kakioka（余纬53.8°,经度140.2°）,Vernadsky〖JP〗（余纬155.3°,经度295.7°）等台站的结果进行了对比分析,讨论了Fourier变换时间窗和小波基函数选取的影响,有限长时间序列引起的边缘效应,以及小波尺度和Fourier频率的关系等问题.我们采用谱分析和小波分析处理了2007年两个台站的地磁数据,研究了地磁场形态和地磁现象的频谱特征,并讨论了Fourier变换和小波变换的优缺点.  相似文献   

Effects of uncertainty of lithofacies, conductivity and porosity distributions on stochastic interpretations of a field scale tracer test   总被引：2，自引：2，他引：0  

Monica Riva  Laura Guadagnini  Alberto Guadagnini 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(7):955-970

We investigate the importance of selecting two different methodologies for the determination of hydraulic conductivity from available grain-size distributions on the stochastic modeling of the depth-averaged breakthrough curve observed during a forced-gradient tracer test experiment. The latter was performed in the Lauswiesen alluvial aquifer, located near the city of Tübingen, Germany, by injecting NaBr into a well at a distance of about 50 m from a pumping well. We also examine the joint effect of the choice of the transport model adopted to describe solute transport at the site and the way the spatial distribution of porosity is assessed. In the absence of direct measurements of porosity, we consider: (a) the model used by Riva et al. (J Contam Hydrol 88:92–118, 2006; J Contam Hydrol 101:1–13, 2008), which relates the natural logarithms of effective porosity and conductivity through an empirical, experimentally-based, linear relationship derived for a nearby experimental site; and (b) a model based on a commonly used relationship linking the total porosity to the coefficient of uniformity of grain size distributions. Transport is described in terms of a purely advective process and/or by including mass exchange processes between mobile and immobile regions. Modeling of flow and transport is performed within a Monte Carlo framework, upon conceptualizing the aquifer as a random composite medium. Our results indicate that the model adopted to describe the correlation between conductivity and porosity and the way grain-sieve information are incorporated to depict the heterogeneous distribution of hydraulic conductivity can have relevant effects in the interpretation of the data at the site. All the conceptual models employed to describe the structural heterogeneity of the system and transport features can reasonably reproduce the global characteristics of the experimental depth-averaged breakthrough curve. Specific details, such as the peak concentration and the time of first arrival, can be better reproduced by a double porosity transport model when a correlation between conductivity and porosity based on grain size information at the site is considered. The best prediction of the late-time behavior of the measured breakthrough curves, in terms of the observed heavy tailing, is offered by directly linking porosity distribution to the spatial variability of particle size information.  相似文献   

Evolution and depositional environments of the Eberswalde fan delta, Mars   总被引：2，自引：0，他引：2  

Monica Pondrelli  Angelo Pio Rossi 《Icarus》2008,197(2):429-451

The Eberswalde crater and its contributing basins have been analyzed in detail in order to reconstruct the geological evolution of the water-related landforms with particular focus on the Eberswalde delta-like feature. Based on a complex strata organization characterized by a topset-foreset-bottomset geometry, typical of delta progradation on Earth, we interpret the Eberswalde feature to be a fan delta associated with a lacustrine system. Depositional sub-environments have been recognized and mapped and the sedimentary processes discussed. A sequence stratigraphy approach has been used to evaluate the system, which we interpret to result from three depositional sequences. These sequences suggest relative water level fluctuations and a longer trend over time towards decreasing water content inside the basin.  相似文献   

BACI design reveals the decline of the seagrass Posidonia oceanica induced by anchoring     

Montefalcone M  Chiantore M  Lanzone A  Morri C  Albertelli G  Nike Bianchi C 《Marine pollution bulletin》2008,56(9):1637-1645

  相似文献   

Scenarios of Earthquake-Generated Tsunamis for the Italian Coast of the Adriatic Sea     

Mara Monica Tiberti  Stefano Lorito  Roberto Basili  Vanja Kastelic  Alessio Piatanesi  Gianluca Valensise 《Pure and Applied Geophysics》2008,165(11-12):2117-2142

We calculated the expected impact on the Italian coast of the Adriatic Sea of a large set of tsunamis resulting from potential earthquakes generated by major fault zones. Our approach merges updated knowledge on the regional tectonics and scenario-like calculations of expected tsunami impact. We selected six elongated potential source zones. For each of them we determined a Maximum Credible Earthquake and the associated Typical Fault, described by its size, geometry and kinematics. We then let the Typical Fault float along strike of its parent source zone and simulated all tsunamis it could generate. Simulations are based on the solution of the nonlinear shallow water equations through a finite-difference technique. For each run we calculated the wavefields at specified simulation times and the maximum water height field (above mean sea level), then generated travel-time maps and maximum water height profiles along the target coastline. Maxima were also classified in a three-level code of expected tsunami threat. We found that the southern portion of Apulia facing Albania and the Gargano promontory are especially prone to the tsunami threat. We also found that some bathymetric features are crucial in determining the focalization-defocalization of tsunami energy. We suggest that our results be taken into account in the design of early-warning strategies.  相似文献   

Estimation of the behavior factor of existing RC-MRF buildings     

Marco Vona  Monica Mastroberti 《地震工程与工程振动(英文版)》2018,17(1):191-204

In recent years, several research groups have studied a new generation of analysis methods for seismic response assessment of existing buildings. Nevertheless, many important developments are still needed in order to define more reliable and effective assessment procedures. Moreover, regarding existing buildings, it should be highlighted that due to the low knowledge level, the linear elastic analysis is the only analysis method allowed. The same codes (such as NTC2008, EC8) consider the linear dynamic analysis with behavior factor as the reference method for the evaluation of seismic demand. This type of analysis is based on a linear-elastic structural model subject to a design spectrum, obtained by reducing the elastic spectrum through a behavior factor. The behavior factor (reduction factor or q factor in some codes) is used to reduce the elastic spectrum ordinate or the forces obtained from a linear analysis in order to take into account the non-linear structural capacities. The behavior factors should be defined based on several parameters that influence the seismic nonlinear capacity, such as mechanical materials characteristics, structural system, irregularity and design procedures. In practical applications, there is still an evident lack of detailed rules and accurate behavior factor values adequate for existing buildings. In this work, some investigations of the seismic capacity of the main existing RC-MRF building types have been carried out. In order to make a correct evaluation of the seismic force demand, actual behavior factor values coherent with force based seismic safety assessment procedure have been proposed and compared with the values reported in the Italian seismic code, NTC08.  相似文献   

Illuminating Subsurface Microbial Water Quality Patterns Using Adenosine Triphosphate and Dynamic Time Warping Approaches     

Alex H. S. Chik  Monica B. Emelko  Alfred P. Blaschke  Jack F. Schijven 《Ground Water Monitoring & Remediation》2020,40(3):55-66

Aquifer microbial water quality evaluations are often performed by collecting groundwater samples from monitoring wells. While samples collected from continuously pumped sources are seldom disputed as representative of the aquifer, natural biofilm present in the vicinity of well screens may introduce unwanted microbial artefacts in monitoring wells that are only periodically sampled. The need for well water purging to obtain samples void of these artefacts has been widely recognized. However, purging methods are not standardized; many approaches presume that physico-chemical water quality stability achieved through the removal of 3 to 5 well volumes is indicative of the stability of target analytes. Using a data set collected from a shallow unconfined aquifer in Southern Ontario, Canada, the need for using dedicated approaches that account for the time-dependent nature of microbial water quality changes was demonstrated. Specifically, the utility of adenosine triphosphate (ATP) as a rapid, field-ready biochemical indicator of microbial water quality stability was investigated. This work shows that ATP concentrations reflect time-limited (bio)colloid transport processes that are consistent with other microbial water quality parameters monitored, but different from commonly measured physical and chemical water quality indicators of well purging adequacy. ATP concentrations occasionally fluctuated even after 3 or 4 h of purging, indicating that microbial artefacts attributable to biofilms in the vicinity of the well screen can still persist. The recurrence of characteristic ATP patterns in each well was systematically examined through the novel application of dynamic time warping (DTW), a nonparametric time series analysis approach. These patterns are believed to be linked with seasonal hydrogeological conditions, which warrant consideration in the design and interpretation of subsurface microbial water quality investigations.  相似文献