Abstract: Geostatistical methods provide a powerful tool to understand the complexity of data arising from Earth sciences. Since the mid 70’s, this numerical approach is widely used to understand the spatial variation of natural phenomena in various domains like Oil and Gas, Mining or Environmental Industries. Considering the huge amount of data available, standard imple- mentations of these numerical methods are not efficient enough to tackle current challenges in geosciences. Moreover, most of the software packages available for geostatisticians are de- signed for a usage on a desktop computer due to the trial and error procedure used during the interpolation. The Geological Data Management (GDM ) software package developed by the French geological survey (BRGM) is widely used to build reliable three-dimensional geological models that require a large amount of memory and computing resources. Considering the most time-consuming phase of kriging methodology, we introduce an efficient out-of-core algorithm that fully benefits from graphics cards acceleration on desktop computer. This way we are able to accelerate kriging on GPU with data 4 times bigger than a classical in-core GPU algorithm, with a limited loss of performances.

Abstract: With the advanced image, drilling and field observation technology, geological structure of reservoirs can be described in more details. A novel 3D mesh generation method for meshing reservoir models is proposed and implemented with arbitrary stratigraphical surface constraints, which ensures the detailed structure geometries and material properties of reservoirs are better described and analysed. The stratigraphic interfaces are firstly extracted and meshed, and then a tetrahedron mesh is generated in 3D with the constraints of such meshed surfaces. The proposed approach includes the following five steps: (1) extracting stratum interfaces; (2) creating a background mesh with size field on the interfaces; (3) constructing geodesic isolines from interface boundaries to the interior; (4) employing a geodesic-based approach to create surface triangles on the area between adjacent isolines and then merge them together; (5) generating tetrahedron mesh for 3D reservoirs with constraints of generated surface triangular mesh. The proposed approach has been implemented and applied to the Lawn Hill reservoir as a practical example to demonstrate its effectiveness and usefulness.

Abstract: We are developing a coupling software ppOpen-MATH/MP. ppOpen-MATH/MP is characterized by its wide applicability. This feature comes from the design that grid point correspondence and interpolation coefficients should be calculated in advance. However, calculation of these values on the unstructured grid model requires a lot of computation time in general. So, we developed new effective algorithm and program for calculating the grid point correspondence as a pre-processor of ppOpen-MATH/MP. In this article, an algorithm and performance evaluation of the program is presented in the first half, and next, an application example of ppOpen-MATH/MP, targeting atmospheric model NICAM and ocean model COCO coupling, is described.

Abstract: In order to execute various finite-difference method applications on large-scale parallel computers with a reasonable cost of computer resources, a framework using an adaptive mesh refinement (AMR) technique has been developed. AMR can realize high-resolution simulations while saving computer resources by generating and removing hierarchical grids dynamically. In the AMR framework, a dynamic domain decomposition (DDD) technique, as a dynamic load balancing method, is also implemented to correct the computational load imbalance between each process associated with parallelization. By performing a 3D AMR test simulation, it is confirmed that dynamic load balancing can be achieved and execution time can be reduced by introducing the DDD technique.