Forchheimer and added mass effects in the flow towards a well

In this note, we examine the flow towards a well in a confined aquifer in the presence of an interaction force defined by the sum of three terms, namely, a Darcy term (linear in the velocity), a Forchheimer term (quadratic in the velocity), and an added‐mass term (linear in the acceleration). We obtain the exact dynamic solution for the piezometric head distribution around the well and investigate the relative importance of the non‐Darcian terms. Copyright © 2013 John Wiley & Sons, Ltd.     

Declustering of Iran earthquake catalog (1983–2017) using the epidemic-type aftershock sequence (ETAS) model

The main goal of this article is to decluster Iranian plateau seismic catalog by the epidemic-type aftershock sequence (ETAS) model and compare the results with some older methods. For this purpose, Iranian plateau bounded in 24°–42°N and 43°–66°E is subdivided into three major tectonic zones: (1) North of Iran (2) Zagros (3) East of Iran. The extracted earthquake catalog had a total of 6034 earthquakes (Mw?>?4) in the time span 1983–2017. The ETAS model is an accepted stochastic approach for seismic evaluation and declustering earthquake catalogs. However, this model has not yet been used to decluster the seismic catalog of Iran. Until now, traditional methods like the Gardner and Knopoff space–time window method and the Reasenberg link-based method have been used in most studies for declustering Iran earthquake catalog. Finally, the results of declustering by the ETAS model are compared with result of Gardner and Knopoff (Bull Seismol Soc Am 64(5):1363–1367, 1974), Uhrhammer (Earthq Notes 57(1):21, 1986), Gruenthal (pers. comm.) and Reasenberg (Geophys Res 90:5479–5495, 1985) declustering methods. The overall conclusion is difficult, but the results confirm the high ability of the ETAS model for declustering Iranian earthquake catalog. Use of the ETAS model is still in its early steps in Iranian seismological researches, and more parametric studies are needed.     

Operational model for direct determination of evapotranspiration for well watered crops in Mediterranean region

Direct calculation of actual evapotranspiration ET_c based on Penman-Monteith type models gives more accurate values than indirect models, which need the determination of reference evapotranspiration and crop coefficient. However, the direct models need the measurement of weather variables above the crop, which is limiting and not easily feasible in practice. An operational version of a known ET_c direct model is described and tested. This new version is based on the determination of the weather variables collected in a standard agro-meteorological station. The original and the operational versions of the ET_c model were validated on two crops with contrasting height: soybean (0.8 m) and sweet sorghum (3 m). For soybean, ET_c calculated with the two versions gave results very similar at both hourly and daily scales. For sweet sorghum, ET_c calculated with the operational version is good at daily scale and not as good, although acceptable, at the hourly scale.     

An operational model to estimate hourly and daily crop evapotranspiration in hilly terrain: validation on wheat and oat crops

In this paper, we present an operational model to estimate the actual evapotranspiration (ET) of crops cultivated on hilly terrains. This new model has the following three characteristics: (1) ET modelling is based on a Penman?CMonteith (PM) type equation (Monteith 1965) where canopy resistance is simulated by following an approach already illustrated by Katerji and Perrier (Agronomie 3(6):513?C521, 1983); (2) the estimation of ET, by means of the PM equation, is made by using meteorological variables simulated on sloped sites as input; (3) these variables are simulated by using simple relationships linking the variables measured at a reference site on plane to the topographic characteristics of the site (slope, orientation, altitude as difference between reference, and sloped sites). This approach presents two advantages if compared with previously proposed models: Not only computation steps are greatly simplified but also error sources due to the simulation of climatic variables in sloped sites and the ET estimation are well distinguished. This model was validated at hourly and daily scales at four sites cultivated with wheat and oats offering a wide range of slope and orientation values: a reference site on plane, site 1 (9° sloping, NW orientation, 7 m from the reference site in plane), site 2 (6°, SE, 12 m) and site 3 (1°, SE, 18 m). At hourly scale, the new model performed well at all sites studied. The observed slope of the linear relationships between estimated and measured ET values ranged between 0.93 and 1.03, with coefficients of determination, r ², between 0.80 and 0.98. At daily scale, the slopes of the linear relationships between measured and estimated ET for the sites on plane and the sloped sites were practically the same (0.98?±?0.01); however, the coefficient of determination r ² observed in the site on plane was clearly greater (0.98) than that observed in the sloped sites (0.83). The presented analysis does not show any significant systematic effect of topography (slope and orientation of the plots) on the good performance of the proposed model for the ET estimation. Furthermore, we observed that coefficients of determination tend to decrease with the increase in the slope of the site, which translates into increased inaccuracy of the climatic variables simulation, in particular the net radiation, as the slope of the site increases. The proposed model allows to verify the different steps for calculating the fluxes, to identify the eventual sources of error and to make the needed corrections. For this reasons, the proposed model seems to be particularly ??operational??, i.e. a useful tool for estimating fluxes on hilly terrains.     

Fluid–rock interactions associated with regional tectonics and basin evolution

An integrated approach consisting of fracture analysis, petrography, carbon, oxygen and strontium‐isotope analyses, as well as fluid‐inclusion micro‐thermometry, led to a better understanding of the evolution of fluid–rock interactions and diagenesis of the Upper Permian to Upper Triassic carbonates of the United Arab Emirates. The deposited carbonates were first marked by extensive early dolomitization. During progressive burial, the carbonates were affected by dolomite recrystallization as well as precipitation of vug and fracture‐filling dolomite, quartz and calcite cements. After considerable burial during the Middle Cretaceous, sub‐vertical north–south oriented fractures (F1) were cemented by dolomite derived from mesosaline to hypersaline fluids. Upon the Late Cretaceous maximum burial and ophiolite obduction, sub‐vertical east–west fractures (F2) were cemented by dolomite (Dc2) and saddle dolomite (Ds) derived from hot, highly saline fluids. Then, minor quartz cement has precipitated in fractures from hydrothermal brines. Fluid‐inclusion analyses of the various diagenetic phases imply the involvement of increasingly hot (200°C) saline brines (20 to 23% NaCl eq.). Through one‐dimensional burial history numerical modelling, the maximum temperatures reached by the studied rocks are estimated to be in the range of 160 to 200°C. Tectonically‐driven flux of hot fluids and associated diagenetic products are interpreted to have initiated during the Late Cretaceous maximum burial and lasted until the Oligocene–Miocene compressional tectonics and related uplift. The circulation of such hydrothermal brines led to partial dissolution of dolomites (Dc2 and Ds) and to precipitation of hydrothermal calcite C1 in new (mainly oriented north–south; F3) and pre‐existing, reactivated fractures. The integration of the obtained data confirms that the diagenetic evolution was controlled primarily by the interplay of the burial thermal evolution of the basin and the regional tectonic history. Hence, this contribution highlights the impacts of regional tectonics and basin history on diagenetic processes, which may subsequently affect reservoir properties.     

Comparative probabilistic assessment of the hydrological performance of reconstructed and natural watersheds

The oil sands industry has committed to returning the mine sites to a productive condition. The reconstructed soil covers must have sufficient available water holding capacity (AWHC) to supply enough moisture over the growing season, to promote vegetation. In order to assess the sustainability of various soil cover alternatives, a generic, system dynamic watershed model entitled GSDW was used along with the available historical meteorological records to estimate the maximum soil moisture deficit and annual evapotranspiration fluxes. A probabilistic framework was adopted; consequently, frequency curves of the maximum annual moisture deficit values are constructed and used to assess the probability that various reconstructed and natural watersheds can provide the associated moisture demands. In general, the study showed a tendency for the reconstructed watershed to provide less moisture for evapotranspiration than natural systems. Watersheds of various soil types, layering, thicknesses and topography were studied. The gained knowledge was used to predict the possible performance of a hypothetical reclamation cover. The results indicated that the hypothetical cover performed in a similar manner to the thickest existing soil cover which confirmed a high probability of that cover to survive under the same existing climatic conditions. Moreover, this probabilistic framework was found to be useful for integrating information gained from natural watersheds (e.g. the canopy of mature natural systems and transfer the results to the reconstructed system). The results show that the canopy influenced the moisture deficit regime positively which signifies a greater possibility that reconstructed covers will adapt to vegetation type. In brief, the adopted approach enables better understanding of the response of reconstructed systems via multiple simulations of ‘what‐if’ scenarios using different soil/vegetation alternatives. Copyright © 2010 John Wiley & Sons, Ltd.     

Diagenesis of the Khuff Formation (Permian–Triassic), northern United Arab Emirates

This diagenetic study (including fieldwork, petrographic, fluid inclusion, and stable isotope investigations) deals with the outcrop of Upper Permian–Lower Triassic carbonate rocks, which are equivalent to the Khuff Formation. The studied succession, which outcrops in the Ras Al Khaimah region, northern United Arab Emirates, comprises three formations, including the Bih, the Hagil, and the Ghail formations. The study focuses on unraveling the conditions and fluid compositions encountered during diagenesis of the succession. Emphasize is also made on linking diagenesis to major stratigraphic surfaces and to highlight reservoir property evolution and heterogeneity of the studied rocks. The evolution of fluids and related diagenetic products can be summarized as follows: (1) formation of near-surface to shallow burial, fine-crystalline dolomite (dolomite matrix) through pervasive dolomitization of carbonate sediments by modified marine pore waters; (2) formation of coarse-crystalline dolomite cement by highly evolved marine pore waters (13–23 wt.% NaCl eq.) at elevated temperatures (120–208°C), and (3) calcite cementation by highly saline fluid (20–23 wt.% NaCl eq.) at high temperature (170–212°C). A final calcite cement generation has been formed by the percolation of meteoric fluids during uplift. Fracture- and vug-filling diagenetic minerals are mainly restricted to the mid-Bih breccia marker level, suggesting preferential focused fluid flow through specific stratigraphic surfaces as well as along tectonic-related structures. Reservoir properties have been evolved as result of the interplay of the original sedimentary texture and the diagenetic evolution. Porosity is higher in the Bih Formation, which is dominated by dolomitized packstones and grainstones, than in the Hagil and Ghail formations, consisting mainly of dolomitized mudstones and wackestones. Image analyses were used to quantify the visual porosity in thin sections. The highest porosity values were measured in the Bih Formation, which is characterized by significant amounts of vug- and fracture-filling cements. This feature is attributed to the increase of porosity owing to substantial dissolution of abundant intergranular and vug-filling cements. In contrast, the Hagil and Ghail formations, which consist of finer-grained rock than the Bih Formation, were less cemented, and thus, the porosity enhancement by cement dissolution was insignificant.