Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

Improved regional groundwater flow modeling using drainage features: a case study of the central northern karst aquifer system of Puerto Rico (USA)

In northern Puerto Rico (USA), subsurface conduit networks with unknown characteristics, and surface features such as springs, rivers, lagoons and wetlands, drain the coastal karst aquifers. In this study, drain lines connecting sinkholes and springs are used to improve the developed regional model by simulating the drainage effects of conduit networks. Implemented in an equivalent porous media (EPM) approach, the model with drains is able to roughly reproduce the spring discharge hydrographs in response to rainfall. Hydraulic conductivities are found to be scale dependent and significantly increase with higher test radius, indicating scale dependency of the EPM approach. Similar to other karst regions in the world, hydraulic gradients are steeper where the transmissivity is lower approaching the coastline. This study enhances current understanding of the complex flow patterns in karst aquifers and suggests that using a drainage feature improves modeling results where available data on conduit characteristics are minimal.     

Detrital zircon provenance analysis in the Zagros Orogen,SW Iran: implications for the amalgamation history of the Neo-Tethys

International Journal of Earth Sciences - The Zagros Orogen developed as a result of Arabia–Eurasia collision. New in situ detrital zircon U–Pb and Hf isotopic analyses from a Cenozoic...     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Brittleness Effect on Rock Fatigue Damage Evolution

The damage evolution mechanism of rocks is one of the most important aspects in studying of rock fatigue behavior. Fatigue damage evolution of three rock types (onyx marble, sandstone and soft limestone) with different brittleness were considered in the present study. Intensive experimental tests were conducted on the chosen rock samples and acoustic emission (AE) sensors were used in some of them to monitor the fracturing process. Experimental tests indicated that brittleness strongly influences damage evolution of rocks in the course of static and dynamic loading. AE monitoring revealed that micro-crack density induced by the applied loads during different stages of the failure processes increases as rock brittleness increases. Also, results of fatigue tests on the three rock types indicated that the rock with the most induced micro-cracks during loading cycles has the least fatigue life. Furthermore, the condition of failure surfaces of the studied rocks samples, subjected to dynamic and static loading, were evaluated and it was concluded that the roughness of failure surfaces is influenced by loading types and rock brittleness. Dynamic failure surfaces were rougher than static ones and low brittle rock demonstrate a smoother failure surface compared to high brittle rock.     

Synoptic analysis of dust storms in the Middle East

Dust storm in the Middle East and south-west Asia is a natural hazard and the Tigris-Euphrates alluvial plain has been recognized as the main dust source in this area. In this study, more than 60 dust storms that occurred during the period 2003–2011 are investigated on the basis of MODIS satellite images, and 12 of the dust storms are selected for synoptic analysis using the NCEP-NCAR Reanalysis Data. The potential dust sources in the Middle East and south-west Asian region (20°E to 80°E, 5°N to 50°N) are analyzed and used in the synoptic analysis. Dust storms in the region can be grouped into two main categories, i.e., the Shamal dust storms and the frontal dust storms. Synoptic systems, associated with the two categories, are distinguished and the frequency of the patterns is identified. For 68% of the Shamal dust storms, a high pressure system is situated between 0°E to 30°E and 27°N to 45°N, and a low pressure system between 50°E to 70°E and 23°N to 43°N. For 86% of the frontal dust storms, a high is located between 51°E to 67°E and 18°N to 33°N and a low between 28°E to 48°E and 32°N to 43°N. Three main patterns for Shamal dust storms are identified, which represent about 60% of the Shamal dust storms. This analysis confirms that the Shamal is related to the anticyclones located over northern Africa to Eastern Europe and the monsoon trough over Iraq, southern Iran, Pakistan and the Indian Subcontinent. The analysis also shows that the main dust sink for the frontal dust storms in Tigris and Euphrates alluvial plain extends from center of Iraq to west and center of Iran and, in most severe cases, to northern Iran and the southern coast of the Caspian Sea.