Fast hyperbolic deconvolutive Radon transform using generalized Fourier slice theorem

The hyperbolic Radon transform has a long history of applications in seismic data processing because of its ability to focus/sparsify the data in the transform domain. Recently, deconvolutive Radon transform has also been proposed with an improved time resolution which provides improved processing results. The basis functions of the (deconvolutive) Radon transform, however, are time-variant, making the classical Fourier based algorithms ineffective to carry out the required computations. A direct implementation of the associated summations in the time–space domain is also computationally expensive, thus limiting the application of the transform on large data sets. In this paper, we present a new method for fast computation of the hyperbolic (deconvolutive) Radon transform. The method is based on the recently proposed generalized Fourier slice theorem which establishes an analytic expression between the Fourier transforms associated with the data and Radon plane. This allows very fast computations of the forward and inverse transforms simply using fast Fourier transform and interpolation procedures. These canonical transforms are used within an efficient iterative method for sparse solution of (deconvolutive) Radon transform. Numerical examples from synthetic and field seismic data confirm high performance of the proposed fast algorithm for filling in the large gaps in seismic data, separating primaries from multiple reflections, and performing high-quality stretch-free stacking.     

Simulation of relationship between river discharge and sediment yield in the semi-arid river watersheds

Complex and variable nature of the river sediment yield caused many problems in estimating the long-term sediment yield and problems input into the reservoirs. Sediment Rating Curves (SRCs) are generally used to estimate the suspended sediment load of the rivers and drainage watersheds. Since the regression equations of the SRCs are obtained by logarithmic retransformation and have a little independent variable in this equation, they also overestimate or underestimate the true sediment load of the rivers. To evaluate the bias correction factors in Kalshor and Kashafroud watersheds, seven hydrometric stations of this region with suitable upstream watershed and spatial distribution were selected. Investigation of the accuracy index (ratio of estimated sediment yield to observed sediment yield) and the precision index of different bias correction factors of FAO, Quasi-Maximum Likelihood Estimator (QMLE), Smearing, and Minimum-Variance Unbiased Estimator (MVUE) with LSD test showed that FAO coefficient increases the estimated error in all of the stations. Application of MVUE in linear and mean load rating curves has not statistically meaningful effects. QMLE and smearing factors increased the estimated error in mean load rating curve, but that does not have any effect on linear rating curve estimation.     

Integration of rock physics and seismic inversion for rock typing and flow unit analysis: A case study

Rock typing and flow unit detection are more challenging in clastic reservoirs with a uniform pore system. An integrated workflow based on well logs, inverted seismic data and rock physics models is proposed and developed to address such challenges. The proposed workflow supplies a plausible reservoir model for further investigation and adds extra information. Then, this workflow has been implemented in order to define different rock types and flow units in an oilfield in the Persian Gulf, where some of these difficulties have been observed. Here, rock physics models have the leading role in our proposed workflow by providing a diagnostic framework in which we successfully differentiate three rock types with variant characteristics on the given wells. Furthermore, permeability and porosity are calculated using the available rock physics models to define several flow units. Then, we extend our investigation to the entire reservoir by means of simultaneous inversion and rock physics models. The outcomes of the study suggest that in sediments with homogeneous pore size distribution, other reservoir properties such as shale content and cementation (which have distinct effects on the elastic domain) can be used to identify rock types and flow units. These reservoir properties have more physical insights for modelling purposes and can be distinguished on seismic cube using proper rock physics models. The results illustrate that the studied reservoir mainly consists of rock type B, which is unconsolidated sands and has the characteristics of a reservoir for subsequent fluid flow unit analysis. In this regard, rock type B has been divided into six fluid units in which the first detected flow unit is considered as the cleanest unit and has the highest reservoir process speed about 4800 to 5000 mD. Here, reservoir quality decreases from flow unit 1 to flow unit 6.     

Sources of the predynastic grinding stones in the hu-semaineh region,upper Egypt,and their cultural context

Grinding stones excavated at the Predynastic site of HG in Upper Egypt were analyzed, and consist of igneous rocks (rhyolite porphyry, basalt, granite) and metamorphic rocks (marble, quartzite). These rocks were brought to the site in Predynastic times (ca. 4000–3000 B. C.) from the Dokkan Volcanics located in the Wadi Hammamat, possibly as far away as 150 km. This is further evidence for the widespread exchange network that existed in Predynastic Upper Egypt. Complex economic interaction, as evidenced by artifacts excavated at HG, was a prime force as the early state evolved in Egypt. © 1993 John Wiley & Sons, Inc.     

Assessment of Estimated Bedrock and Stream Sediment Geochemical Backgrounds in Catchment Basin Analysis

Natural Resources Research - The weighted mean and the multiple regression techniques are two methods that are employed to estimate elemental background concentration of lithologies upstream of...     

Effects of urbanization on river morphology of the Talar River,Mazandarn Province,Iran

In the present study, we investigate the effects of urbanization growth on river morphology in the downstream part of Talar River, east of Mazandaran Province, Iran. Morphological and morphometric parameters in 10 equal sub-reaches were defined along a 11.5 km reach of the Talar River after land cover maps were produced for 1955, 1968, 1994, 2005 and 2013. Land cover types changed extremely during the study period. Residential lands were found to have increased in area by about 1631%, while forest land and riparian vegetation decreased in by approximately 99.9 and 96.2%, respectively. The results of morphometric and morphological factors showed that average channel width (W) for all 11.5 km of the study river decreased by 84% during the study period, while the flow length increased by about 2.14%.     

Modeling the connectivity and intersection of hydraulically loaded cracks with in situ fractures in rock

The interaction of an advancing hydraulically loaded crack and in situ fracture network can yield highly complex patterns. We model the connectivity of cells in a finite element domain and in a fracture network by a simplicial complex data structure. The complete adjacency information between cells is determined by one level down facet and one level up cofacet neighborhood information. Combined with a disjoint set data structure, explicit algorithms are derived to efficiently track network connectivity and load transfer between independent fracture sets. We also propose an approach to regularize the application of hydraulic load to newly intersected in situ cracks to smoothen the transition of pressure on intersected cracks from ambient to hydraulic pressure and to avoid the sudden loading of the entire length of these cracks. Numerical results demonstrate the performance of crack connectivity and load transfer models and the effect of regularization model. The results show that as the angle between an incoming hydraulically loaded crack and an in situ crack increases, the effect of in situ crack shifts from slight realignment to diversion/offsetting of the loaded crack. As the angle difference approaches the normal angle, the loaded crack tends to penetrate through the in situ crack. The proposed schemes are also used for transient simulation of 2D reservoirs with multiple perforations surrounded by in situ cracks with and without a bias in the distribution of their orientation. It was shown that from 2 perforations with angles closer to in situ cracks at low loading rates to all perforations at higher loading rates can result in active hydraulic crack propagation. The h‐adaptive method of asynchronous space‐time discontinuous Galerkin method is used to exactly track complex fracture patterns in these dynamic fracture simulations.     

Distributed empirical algorithms to estimate catchment scale sediment connectivity and yield in a subtropical region

The intensity of soil loss and sediment delivery, representing hydrologic and geomorphic processes within a catchment, accelerates with rapid changes in land cover and rainfall events. An underlying component of sustainable management of water resources is an understanding of spatial and temporal variability and the adverse influences of regional parameters involved in generating sediment following widespread changes in land cover. A calibrated algorithm of soil loss coupled with a sediment delivery ratio (SDR) was applied in raster data layers to improve the capability of a combined model to estimate annual variability in sediment yields related to changes in vegetation cover identified by analyses of SPOT imagery. Four catchments in Kangaroo River State forest were assessed for annual changes in sediment yields. Two catchments were selectively logged in 2007, while the two other sites remained undisturbed. Results of SDR estimates indicated that only a small proportion of total eroded sediment from hillslopes is transported to catchment outlets. Larger SDR values were estimated in regions close to catchment outlets, and the SDR reduced sharply on hillslopes further than 200–300 m from these areas. Estimated sediment yield increased by up to 30% two years after land cover change (logging) in 2009 when more storm events were recorded, despite the moderate density of vegetation cover in 2009 having almost recovered to its initial pre‐logging (2005) condition. Rainfall had the most significant influence on streamflow and sediment delivery in all catchments, with steeply sloping areas contributing large amounts of sediment during moderate and high rainfall years in 2007 and 2009. It is concluded that the current scenario of single‐tree selection logging utilized in the study area is an acceptable and environmentally sound land management strategy for preservation of soil and water resources. Copyright © 2013 John Wiley & Sons, Ltd.     

Magmatic and post-magmatic evolution of post-collisional rare-metal bearing granite: The Neoproterozoic Homrit Akarem Granitic Intrusion,south Eastern Desert of Egypt,Arabian-Nubian Shield

The Homrit Akarem granitic intrusion (HAGI) outcrops near the western edge of the south Eastern Desert basement exposure in Egypt. It is a composite of two cogenetic intrusive bodies: an early albite granite phase shallowly emplaced at the apex of a magmatic cupola, and a later subjacent pink granite phase with a marginal zone of muscovite granite and better preservation of magmatic features. Mineral chemistry of primary biotite and garnet, together with whole-rock chemistry, identify the HAGI as a highly fractionated A-type peraluminous intrusion. The chemistry of F-dominant, Li-bearing, Fe³⁺-rich primary magmatic mica in the pink granite resembles that typically found in highly evolved Nb-Y-F pegmatites. The HAGI is the evolved product of a primary magma generated by partial melting of juvenile crust of the Arabian-Nubian Shield (ANS), emplaced along a regional strike-slip fault system that promoted its ascent. The main emplacement mechanism and evolutionary sequence of the HAGI was magmatic, although secondary minerals and textures resulting from hydrothermal fluid interactions are observed, especially at its margins. Primary columbite-(Mn) crystallized from melt and was partly replaced by secondary fluorcalciomicrolite. The high fluorine content of magmatic fluids exsolved from the intrusion is indicated by quartz-fluorite veins, greisenization, albitization, and F-bearing secondary oxide minerals. The magmatic derivation of this fluid is demonstrated by the F-dominant primary mica, a siderophyllite-polylithionite solid solution commonly known as zinnwaldite. The chemistry of zinnwaldite constrains the F/OH activity ratio and oxygen fugacity of its parental melt and thereby resolves the ambiguity between pressure and the effects of F in controlling the normative quartz content of rare-metal granites. The HAGI is less mineralized than the post-collisional rare-metal granites found further east in the south Eastern Desert, replicating a trend observed previously in the central Eastern Desert and suggesting that east-west zoning in rare metal enrichment is a persistent feature across latitudes at the western edge of the ANS.