Impact of tailings dam failure on spatial features of copper contamination (Mazraeh mine area,Iran)

Mining activities pose a potential risk of metal contamination around mining sites. On May 6, 2010, a tailings dam failure of the Mazraeh copper mine near Ahar in East Azerbaijan province, Iran, released vast amounts of mine wastes. To better understand the magnitude of copper contamination in the waste-affected soils, it is important to assess the spatial distribution of soil copper content at unsampled points. A total of 30 soil samples and their surficial sediments together with the 6 uncontaminated control samples (0–10 and 10–30 cm) were collected along the stream flow that joined Ahar-Chai River. Some of soil properties as well as total copper concentration were determined in all samples. The mean value of the latter in the surface contaminated soils was found to be approximately two times more than controls. Furthermore, the mean concentration of copper in the surface loaded material was 10 times more than the soils. High copper concentrations were observed in surficial sediments of the soils near the broken tailings dam. The Inverse Distance Weighting (IDW) method was employed in data analysis. The spherical and Gaussian semivariogram models were properly fitted to the data of copper contents in soils and surficial sediments.     

Study on the improving method for gas production prediction in tight clastic reservoir

Relative permeability and resistivity index are important parameters for the exploration and development in a tight sandstone gas field. In the splitting method which uses permeability (K), reservoir thickness (H), and relative permeability (K), briefly referred to as the KHK splitting method, the accuracy of the relative permeability is crucial. According to the relationship between resistivity index and relative permeability of the Mesozoic Lower Safa Formation at Obaiyed Field in the Western Desert of Egypt, we improved the split method and made it more in line with the real situation by adopting Pairoys’ model which is more suitable to our study area. In this paper, we use a radial basis function (RBF) to establish the relationship between logging data and the gas production split point to point in tight sandstone gas reservoirs. To compare with the result by support vector regression (SVR), our method is better as indicated by mean absolute error values. In order to solve the problem that the relative permeability is difficult to obtain in the well logging evaluation, we also provide a convenient method and application example.     

Preservation of organic matter in soils of a climo-biosequence in the Main Range of Peninsular Malaysia

Limited information is available about factors of soil organic carbon(SOC) preservation in soils along a climo-biosequence. The objective of this study was to evaluate the role of soil texture and mineralogy on preservation of SOC in the topsoil and subsoil along a climo-biosequence in the Main Range of Peninsular Malaysia. Soil samples from the A and B-horizons of four representative soil profiles were subjected to particle-size fractionation and mineralogical analyses including X-ray diffraction and selective dissolution. The proportion of SOC in the 250-2000 μm fraction(SOC associated with coarse sand) decreased while the proportion of SOC in the 53 μm fraction(SOC associated with clay and silt)increased with depth. This reflected the importance of the fine mineral fractions of the soil matrix for SOC storage in the subsoil. Close relationships between the content of SOC in the 53 μm fraction and the content of poorly crystalline Fe oxides [oxalate-extractable Fe(Fe_o) – pyrophosphate-extractable Fe(Fe_p)] and poorly crystalline inorganic forms of Al [oxalateextractable Al(Al_o) – pyrophosphate-extractable Al(Al_p)] in the B-horizon indicated the importance of poorly crystalline Fe oxides and poorly crystalline aluminosilicates for the preservation of SOC in the Bhorizon. The increasing trend of Fe_o-Fe_p and Al_o-Al_p over elevation suggest that the importance of poorly crystalline Fe oxides and poorly crystalline aluminosilicates for the preservation of SOC in the Bhorizon increased with increasing elevation. This study demonstrates that regardless of differences in climate and vegetation along the studied climobiosequence, preservation of SOC in the subsoil depends on clay mineralogy.     

A new simplified approach for assessing nonlinear seismic response of plan-asymmetric structures considering soil-structure interaction

An innovative approximate method is presented to consider the plan asymmetry, nonlinear structural behaviour and soil-structure interaction (SSI) effects simultaneously. The proposed method so-called Flexible base 2DMPA (F2MPA) is an extension of 2 degrees of freedom modal pushover analysis (2DMPA) approach to consider foundation flexibility in seismic response analysis of plan asymmetric structures which itself were developed based on Uncoupled Modal Response History Analysis method for inelastic fixed-base asymmetric structures. In F2MPA for each mode shape using 2DMPA procedure, the elastic and inelastic properties of 2DOF modal systems corresponding to the fixed-base structure are initially derived. Then in each time step, displacements and inelastic restoring forces of the superstructure are computed from modal equations of the flexibly-supported structure. In each time step, the nonlinear secant stiffness matrix corresponding to the n-th MDOF modal equations of soil-structure system is updated using the corresponding modal 2DOF system of fixed-base structure. To update the transformed modal stiffness matrix of the SSI system, this matrix is partitioned and it is assumed that the non-linear variation of the superstructure can be estimated from the variation of modal stiffness matrix of the fixed-base structure. Accuracy of the proposed method was verified on an 8-story asymmetric-plan building under different seismic excitations. The results obtained from F2MPA method were compared with those obtained by nonlinear response history analysis of the asymmetric soil-structure system as a reference response. It was shown that the proposed approach could predict the results of the nonlinear time history analysis with a good accuracy. The main advantage of F2MPA is that this method is much less time-consuming and useful for the practical aims such as massive analysis of a nonlinear structure under different records with multiple intensity levels.     

Seismic reliability assessment of a steel moment-resisting frame with two different ductility levels using a cloud analysis approach

A cloud method for generating percentile engineering demand parameter versus intensity measure(EDP-IM) curves of a structure subjected to a set of synthetic ground motions is presented. To this end, an ensemble of synthetic ground motions based on available real ones is generated. This is done by using attenuation relationships, duration and suitable Gutenberg-Richter relations attributed to the considered seismic hazard at a given site by estimating a suitable distribution of magnitude and site to source distance. The study aims to clarify the significance of the duration and frequency content on the seismic performance of structures, which were not considered in developing percentile incremental dynamic analysis(IDA) curves. The collapse probabilities of two steel moment-resisting frames with different ductility levels generated by IDA and the proposed cloud method are compared at different intensity levels. When compared with conventional IDA, the suggested cloud analysis(SCA) methodology with the same run number of dynamic analyses was able to develop response hazard curves that were more consistent with site-specific seismic hazards. Eliminating the need to find many real records by generating synthetic records consistent with site-specific seismic hazards from a few available recorded ground motions is another advantage of using this scheme over the IDA method..     

Prediction of Remediation of a Heterogeneous Aquifer: A Case Study

Contaminant plumes whose characteristic length is smaller than the horizontal integral scale of the hydraulic conductivity, K, are abundant in shallow, phreatic aquifers. In such cases, the aquifer can be regarded as layered, with K being only a function of the vertical coordinate. The heterogeneity of K has a critical role upon the efficiency of remediation of such sites, for example, by Pump and Treat schemes. The expected efficiency is a random variable, with uncertainty. Quantifying this uncertainty can be of great importance to decision making. In this study, we focus on a case study in the coastal aquifer of Israel and compare two different approaches for constructing realizations of K: continuous and indicator. We observe a significant difference between the constructed realizations, which results in a considerable difference in the predicted remediation efficiency and its uncertainty. Furthermore, we study the effect of conditioning the realizations by a rather limited number of K data points. We find that the conditioning results in a major reduction of the uncertainty. In addition, we compare the results of the transport model to a simplified semi‐analytical solution that is based on assuming radial flow. We find a good agreement with the three‐dimensional numerical model. This result illustrates that the simplified solution can be used for prediction of the remediation efficiency when the flow at the plume vicinity can be regarded as radial.     

S-system theory applied to array-based GNSS ionospheric sensing

The GPS carrier-phase and code data have proven to be valuable sources of measuring the Earth’s ionospheric total electron content (TEC). With the development of new GNSSs with multi frequency data, many more ionosphere-sensing combinations of different precision can be formed as input of ionospheric modelling. We present the general way of interpreting such combinations through an application of S-system theory and address how their precision propagates into that of the unbiased TEC solution. Presenting the data relevant to TEC determination, we propose the usage of an array of GNSS antennas to improve the TEC precision and to expedite the rather long observational time-span required for high-precision TEC determination.     

Effects of rotational components of earthquake on seismic response of arch concrete dams

The present study deals with dynamic analysis of arch concrete dams, taking rotational components of earthquakes into account. A modified methodology was used to evaluate the rotational components of the earthquake. The translational components of the earthquake have been used in to obtain the rotational components of the earthquake, based on the intersecting isotropic elastic wave propagation. Two rotational components of Taft, Tabas and San-Fernando earthquakes are evaluated based on the translational components of the earthquakes and considering frequency dependencies of incident angle and wave velocity. Finally, dynamic analyses of Morrow Point Dam are presented to evaluate the effects of combined translational and rotational components on the seismic response of the dam. Various conditions of reservoirs, including full and empty state, are considered in the analyses. Fluid–structure interaction was completely taken into account. It was realized that incorporating rotational components increased the maximum compressive and tensile stresses in both empty and full reservoir analyses. Distribution of maximum tensile stresses is very sensitive to the rotational components of the earthquake. Also, it can be concluded that the segregated effect of the rocking component on the response of concrete dams is more effective than the sole effect of the torsional component.     

Sparse Bayesian linearized amplitude-versus-angle inversion

The technique of seismic amplitude-versus-angle inversion has been widely used to estimate lithology and fluid properties in seismic exploration. The amplitude-versus-angle inversion problem is intrinsically ill-posed and generally stabilized by the use of L2-norm regularization methods but with drawback of smoothing important boundaries between adjacent layers. In this study, we propose a sparse Bayesian linearized solution for amplitude-versus-angle inversion problem to preserve the sharp geological interfaces. In this regard, a priori constraint term with two regularization functions is presented: the sparse constraint regularization and the low-frequency model information. In addition, to obtain high-resolution reflectivity estimation, the model parameters decorrelation technique combined with dipole decomposition method is employed. We validate the applicability of the presented method by both synthetic and real seismic data from the Gulf of Mexico. The accuracy improvement of the presented method is also confirmed by comparing the results with the commonly used Bayesian linearized amplitude-versus-angle inversion.     

A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.