Contribution of electrical tomography methods in geotechnical investigations at Mavropigi lignite open pit mine,Northern Greece

In this paper, the application of 2D and 3D electrical resistivity methods in geotechnical investigations is explored through a case study in Northern Greece. These two methods were employed at a lignite surface mining operation where fracture zones and discontinuities have been recently observed close to the pit boundaries. The main aim of the geophysical survey was to estimate the inclination of the contact between the Neogene and Schist/Carbonate formations near the southern limits of the pit, as well as to estimate the thickness of the carbonate rocks on top of the Schist formations to evaluate the stability of the southern slopes. Synthetic data were initially generated to help plan an efficient electrical tomography survey, in a region with complex geology and irregular terrain. Three configurations (Wenner–Schlumberger and dipole–dipole or pole–dipole) proved essential in such conditions and helped improving the resolution of the resistivity section. The sections were then calibrated by boreholes. Finally, the geophysical survey provided invaluable data regarding the geometry of the bedrock and possible faults, which was essential for the slope stability calculations.     

ISORROPIA: A New Thermodynamic Equilibrium Model for Multiphase Multicomponent Inorganic Aerosols

A computationally efficient and rigorous thermodynamic model that predicts the physical state and composition of inorganic atmospheric aerosol is presented. One of the main features of the model is the implementation of mutual deliquescence of multicomponent salt particles, which lowers the deliquescence point of the aerosol phase.The model is used to examine the behavior of four types of tropospheric aerosol (marine, urban, remote continental and non-urban continental), and the results are compared with the predictions of two other models currently in use. The results of all three models were generally in good agreement. Differences were found primarily in the mutual deliquescence humidity regions, where the new model predicted the existence of water, and the other two did not. Differences in the behavior (speciation and water absorbing properties) between the aerosol types are pointed out. The new model also needed considerably less CPU time, and always shows stability and robust convergence.     

Numerical Simulation of the 2011 Tohoku Tsunami Based on a New Transient FEM Co-seismic Source: Comparison to Far- and Near-Field Observations

In this work, we simulate the 2011 M9 Tohoku-Oki tsunami using new coseismic tsunami sources based on inverting onshore and offshore geodetic data, using 3D Finite Element Models (FEM). Such FEMs simulate elastic dislocations along the plate boundary interface separating the stiff subducting Pacific Plate from the relatively weak forearc and volcanic arc of the overriding Eurasian plate. Due in part to the simulated weak forearc materials, such sources produce significant shallow slip (several tens of meters) along the updip portion of the rupture near the trench. To assess the accuracy of the new approach, we compare observations and numerical simulations of the tsunami's far- and near-field coastal impact for: (i) one of the standard seismic inversion sources (UCSB; Shao et al. 2011); and (ii) the new FEM sources. Specifically, results of numerical simulations for both sources, performed using the fully nonlinear and dispersive Boussinesq wave model FUNWAVE-TVD, are compared to DART buoy, GPS tide gauge, and inundation/runup measurements. We use a series of nested model grids with varying resolution (down to 250 m nearshore) and size, and assess effects on model results of the latter and of model physics (such as when including dispersion or not). We also assess the effects of triggering the tsunami sources in the propagation model: (i) either at once as a hot start, or with the spatiotemporal sequence derived from seismic inversion; and (ii) as a specified surface elevation or as a more realistic time and space-varying bottom boundary condition (in the latter case, we compute the initial tsunami generation up to 300 s using the non-hydrostatic model NHWAVE). Although additional refinements are expected in the near future, results based on the current FEM sources better explain long wave near-field observations at DART and GPS buoys near Japan, and measured tsunami inundation, while they simulate observations at distant DART buoys as well or better than the UCSB source. None of the sources, however, are able to explain the largest runup and inundation measured between 39.5° and 40.25°N, which could be due to insufficient model resolution in this region (Sanriku/Ria) of complex bathymetry/topography, and/or to additional tsunami generation mechanisms not represented in the coseismic sources (e.g., splay faults, submarine mass failure). This will be the object of future work.     

Dimensionality reduction for production optimization using polynomial approximations

The objective of this paper is to introduce a novel paradigm to reduce the computational effort in waterflooding global optimization problems while realizing smooth well control trajectories amenable for practical deployments in the field. In order to overcome the problems of slow convergence and non-smooth impractical control strategies, often associated with gradient-free optimization (GFO) methods, we introduce a generalized approach which represent the controls by smooth polynomial approximations either by a polynomial function or by a piecewise polynomial interpolation, which we denote as function control method (FCM) and interpolation control method (ICM), respectively. Using these approaches, we aim to optimize the coefficients of the selected functions or the interpolation points in order to represent the well-control trajectories along a time horizon. Our results demonstrate significant computational savings, due to a substantial reduction in the number of control parameters, as we seek the optimal polynomial coefficients or the interpolation points to describe the control trajectories as opposed to directly searching for the optimal control values (bottom hole pressure) at each time interval. We demonstrate the efficiency of the method on two and three-dimensional models, where we found the optimal variables using a parallel dynamic-neighborhood particle swarm optimization (PSO). We compared our FCM-PSO and ICM-PSO to the traditional formulation solved by both gradient-free and gradient-based methods. In all comparisons, both FCM and ICM show very good to superior performances.     

Risk provision using field spectroscopy to identify spectral regions for the detection of defects in flexible pavements

Natural and physical hazards accelerate the deterioration of asphalted surfaces. Climatic factors are unavoidable and can affect the properties of asphalt mixtures, making them weaker and less durable. Thus, continuous monitoring of bituminous surfaces is something that can reduce the risks of public health. Remote sensing techniques have become an effective, noninvasive method for early detection of damaged asphalt pavements. This paper outlines a range of different remote sensing methodologies that can be used to monitor asphalt road pavements. This is complemented by the use of field spectroscopy for the examination of asphalt pavements of varying age and conditions. The results of the study found spectral differences regarding asphalt defects, such as physical cracking, patched cracking and polishing. These spectral changes were examined through “in-band” simulation analysis of the Landsat 7 ETM+ sensor, using appropriate relative spectral response filters, concluding that the ratio band 5/band 1 can be used to distinguish asphalt pavements of different date of construction and condition.