PHoToNs–A parallel heterogeneous and threads oriented code for cosmological N-body simulation

We introduce a new code for cosmological simulations, PHo To Ns, which incorporates features for performing massive cosmological simulations on heterogeneous high performance computer(HPC) systems and threads oriented programming. PHo To Ns adopts a hybrid scheme to compute gravitational force, with the conventional Particle-Mesh(PM) algorithm to compute the long-range force,the Tree algorithm to compute the short range force and the direct summation Particle-Particle(PP) algorithm to compute gravity from very close particles. A self-similar space filling a Peano-Hilbert curve is used to decompose the computing domain. Threads programming is advantageously used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. PHo To Ns scales well and efficiency of the PP kernel achieves68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement.     

Observed daily large-scale rainfall patterns during BOBMEX-1999

A daily rainfall dataset and the corresponding rainfall maps have been produced by objective analysis of rainfall data. The satellite estimate of rainfall and the raingauge values are merged to form the final analysis. Associated with epochs of monsoon these rainfall maps are able to show the rainfall activities over India and the Bay of Bengal region during the BOBMEX period. The intra-seasonal variations of rainfall during BOBMEX are also seen using these data. This dataset over the oceanic region compares well with other available popular datasets like GPCP and CMAP. Over land this dataset brings out the features of monsoon in more detail due to the availability of more local raingauge stations.     

Damping characteristics of friction damped braced frame and its effectiveness in the mega-sub controlled structure system

Based on energy dissipation and structural control principle, a new structural configuration, called the mega- sub controlled structure (MSCS) with friction damped braces (FDBs), is first presented. Meanwhile, to calculate the damping coefficient in the slipping state a new analytical method is proposed. The damping characteristics of one-storey friction damped braced frame (FDBF) are investigated, and the influence of the structural parameters on the energy dissipation and the practical engineering design are discussed. The nonlinear dynamic equations and the analytical model of the MSCS with FDBs are established. Three building structures with different structural configurations, which were designed with reference to the conventional mega-sub structures such as used in Tokyo City Hall, are comparatively investigated. The results illustrate that the structure presented in the paper has excellent dynamic properties and satisfactory control effectiveness.     

Wave Propagation Modeling and Amplitude-Variation-with-Offset Response in a Fractured Coalbed

Coalbed methane can be detected employing the amplitudevariation-with-offset technique. However, there are two issues in applying this technique to a coalbed: strong azimuthal anisotropy resulting from high-density fractures, and the seismic response being composed of many or several individual reflections within the coalbed. To overcome these difficulties, we present an exact solution for reflections in extensive dilatancy anisotropy media. First, we build a three-layer model and simulate the wave propagation in this model. Then we derive an exact P- and converted S-wave reflection coefficient equation based on boundary conditions. Finally, substituting given model parameters into the exact equation, we obtain the variation in the reflection coefficient with incidence angle. The results show that the fracture factors, wavelet frequency and thickness of the coalbed have different effects on the reflection coefficient. Furthermore, we create a synthetic seismogram by forward calculation, and the result fits well with results of the exact equation.     

Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.     

Hybridization of the probability perturbation method with gradient information

Geostatistically based history-matching methods make it possible to devise history-matching strategies that will honor geologic knowledge about the reservoir. However, the performance of these methods is known to be impeded by slow convergence rates resulting from the stochastic nature of the algorithm. It is the purpose of this paper to introduce a method that integrates qualitative gradient information into the probability perturbation method to improve convergence. The potential of the proposed method is demonstrated on a synthetic history-matching example. The results indicate that inclusion of qualitative gradient information improves the performance of the probability perturbation method.     

A level-set corrector to an adaptive multiscale permeability prediction

We consider the inverse problem of permeability estimation for two-phase porous-media flow. The novel approach is based on regularization by zonation, where the geometry and size of the regions are chosen adaptively during the optimization procedure. To achieve this, we have utilized level-set functions to represent the permeability. The available data are sparsely distributed in space; hence, it is reasonable to confine the estimation to coarse-scale structures. The level-set approach is able to alter the boundaries between regions of different permeability without strict restrictions on their shape; however, when the data are sparse, a reasonable initial guess for the permeability is required. For this task, we use adaptive multiscale permeability estimation, which has the potential of identifying main permeability variations. These are described by a piecewise constant function, where the constant values are attained on rectangular zones. In the current work, we develop a level-set corrector strategy, assuming adaptive multiscale permeability estimation as a predictor.     

Types of sedimentary environment in alluvial sediments distinguished on the basis of its chemical constitution: the example of the lower course of the Obra river (Western Poland)

Research concerning chemical constitution of alluvial sediments was done in the lower course of the Obra river (Western Poland). The fragment of vertical profile, which consisted of various alluvial sediments (fine sands, peats, and sandy silts) was chosen for detailed analysis. The main research problem was to determine if lithology and chemical constitution of alluvial deposits are interconnected in a distinct way within studied section of the Obra river valley, and, if changes of chemical elements’ concentration could be used to illustrate depositional processes, which take place in riverbed and floodplain. Concentrations of Fe, Mn, Cu, Zn, Ca, Mg, and K were determined in collected sediment samples. Next, geochemical groups of alluvial sediments were distinguished using cluster analysis. Investigated changes of chemical elements’ concentration show that there is a distinct border between the environment of organic sediments, which marks the place of former functioning of the Obra riverbed, and sandy silts, which were deposited within floodplain during floods. Besides, floodplain sediments, riverbed sediments, and reductive environment within peat deposits were singled out on the basis of the analysis.