A Parametric Approach for Dealing with Compositional Rounded Zeros

In this work, a parametric approach for replacing data below the detection limit, also known as rounded zeros, in compositional data sets is proposed. Compositional rounded zeros correspond to small proportions of some whole that cannot be reliably detected by the analytical instruments under given operating conditions. This kind of zeros appear frequently in the data collection process in geosciences. They must be treated in an adequate way before some multivariate analysis can be applied. Our procedure results from a modification of the Expectation-Maximization (EM) algorithm and is based on the additive log-ratio transformation. Its coherence with the nature of compositional data and with basic operations in the simplex sample space is checked. Using real data sets, we find that this approach improves other parametric and non-parametric techniques for compositional rounded zeros.     

Influence of initial density of cohesionless soil on evolution of passive earth pressure

This paper focuses on the influence of the initial void ratio on the evolution of the passive earth pressure and the formation of shear zones in a dry sand body behind a retaining wall. For the numerical simulation a rigid and very rough retaining wall undergoing a horizontal translation against the backfill is considered. The essential mechanical properties of cohesionless granular soil are described with a micro-polar hypoplastic model which takes into account stresses and couple stresses, pressure dependent limit void ratios and the mean grain size as a characteristic length. Numerical investigations are carried out with an initially medium dense and initially loose sand using a homogeneous and random distribution of the initial void ratio. The geometry of calculated shear zones is discussed and compared with a corresponding laboratory model test.     

坝基集中渗漏通道的温度场反分析与试验研究

渗漏是严重影响堤坝安全的主要形式,确定渗漏通道的位置是进行堤坝渗漏整治的前提和基础。温度场探测堤坝渗漏的研究是其中一种较便捷的方法,但是在理论上还需要进一步完善。在库水和地下水存在着温度差异的基础上,考虑堤坝发生集中渗漏后,渗漏通道中的地下水与周围土体进行热量交换,造成渗流通道附近地层温度发生改变,建立起温度场对渗漏通道的反演模型。在前人研究的基础上,考虑地温梯度的影响,可以建立有集中渗漏通道影响下的二维热传导方程,进而以现场测井的温度数据作为边界条件,对泛定方程进行求解,最终得到渗漏通道的位置。计算结果和实测结果的对比表明,通过温度场反推渗漏通道参数的方法是可行的。     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

NUMERICAL STUDY ON THE MIXING OF UNSORTED SEDIMENT PARTICLES DISCHARGED INTO A CROSS-FLOW BY MULTIPHASE PARTICLE-IN-CELL （MP-PIC） METHOD

The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.     

Gaussian beams in inhomogeneous media: A review

The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method, which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are twins, though they differ by names. The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered medium.     

Interpretation of Gravity Anomalies with the Normalized Full Gradient (NFG) Method and an Example

The Normalized Full Gradient (NFG) method which was put forward about 50 years ago has been used for downward continuation of gravity potential data, especially in the former Union of Soviet Socialist Republics. This method nullifies perturbations due to the passage of mass depth during downward continuation. The method depends on the downwards analytical continuation of normalized full gradient values of gravity data. Analytical continuation discriminates certain structural anomalies which cannot be distinguished in the observed gravity field. This method has been used in various petroleum and tectonic studies. The Trapeze method was used for the determination of Fourier coefficients during the application of this method. No other techniques for calculating these coefficients have been used. However, the Filon method was used for the determination of Fourier coefficients during the application of the NFG method in this work. This method, rather than the Trapeze method, should be preferred for indicating abnormal mass resources at the lower harmonics. In this study, the NFG method using the Filon method has been applied the first time to theoretical models of gravity profiles as example field at the Hasankale-Horasan petroleum exploration province where successful results were achieved. Hydrocarbon presence was shown on the NFG sections by the application of NFG downward continuation operations on theoretical models. Important signs of hydrocarbon structure on the NFG section for field and model data at low harmonics are obtained more effectively using this method.     

EM Monitoring of Crustal Processes Including the Use of the Network-MT Observations

There are several kinds of coupling mechanisms which can convert mechanical, chemical or thermal energies due to seismic or volcanic activities into electromagnetic energies. As a result of concentrated efforts in laboratory and theoretical research, the basic relationship between the intensity of electromagnetic sources and changes in mechanical, chemical and thermal state is becoming established. Also with the progress of the electromagnetic simulation techniques, it has been possible to evaluate in situ sensitivity. Based on this progress and also due to extensive improvement in measuring techniques, many field experiments have been performed to elucidate subsurface geophysical processes underlying the preparation stage, onset, and subsequent healing stage of earthquakes and volcanic eruptions. In volcanic studies, many studies have reported the measurement of electromagnetic signals which were successfully interpreted in terms of various driving mechanisms. Although there have been numerous reports about the existence of precursory electromagnetic signals in seismic studies, only a few of them could be successfully explained by the proposed mechanisms, whereas coseismic phenomena are often consistent with those mechanisms including the absence of detectable signals. In many cases, one or two orders of higher sensitivity were required, especially for precursory signals. Generally, electromagnetic methods are more sensitive to near-surface phenomena. It will be necessary to discriminate electromagnetic signals due to these near-surface sources, which often possess no relationship with the crustal activities. Further efforts to enhance in situ sensitivity through improvements in observation techniques and in data processing techniques are recommended. At the same time, multi-disciplinary confirmation against the validity of electromagnetic phenomena will inevitably be necessary. A Network-MT observation technique has been developed to determine large-scale deep electrical conductivity structure. In the method, a telephone line network or purpose-built long baseline cables are utilized to measure voltage differences with long electrode separations. Because of the averaging effect of the electric fields, static shift problems due to small-scale, near-surface lateral heterogeneities can be alleviated. Several field experiments revealed regional scale deep electrical conductivity structures related to slab subduction or its stagnation, which enable us to elucidate underlying physical processes caused by the slab motion. The technique can also be applied to monitor the electric potential field related to crustal activities. The annual variation of the potential field and electrical conductivity in the French Alps were interpreted to be caused by the annual variation of lake water level. The method was also used to monitor the regional scale spatio-temporal variation of the SP field and electrical conductivity before and at the onset of earthquakes and volcanic eruptions.