Modelling discontinuous well log signal to identify lithological boundaries via wavelet analysis: An example from KTB borehole data

Identification of sharp and discontinuous lithological boundaries from well log signal stemming from heterogeneous subsurface structures assumes a special significance in geo-exploration studies. Well log data acquired from various geological settings generally display nonstationary/nonlinear characteristics with varying wavelengths and frequencies. Modelling of such complex well-log signals using the conventional signal processing techniques either fails to catch-up abrupt boundaries or at the best, do not provide precise information on insidious lithological discontinuities. In this paper, we have proposed a new wavelet transform-based algorithm to model the abrupt discontinuous changes from well log data by taking care of nonstationary characteristics of the signal. Prior to applying the algorithm on the geophysical well data, we analyzed the distribution of wavelet coefficients using synthetic signal generated by the first order nonstationary auto-regressive model and then applied the method on actual well log dataset obtained from the KTB bore hole, Germany. Besides identifying the formation of layered boundaries, the underlying method also maps some additional formation boundaries, which were hitherto undetected at the KTB site. The results match well with known geological lithostratigraphy and will be useful for constraining the future model of KTB bore hole data.     

Stable clustering of offshore downhole data using a combined k-means and Gaussian mixture modelling approach

Marine Geophysical Research -     

Influence of magnetic effect,anisotropy, irregularity,initial stress and heterogeneity on propagation of SH-wave in an irregular pre-stressed magnetoelastic monoclinic sandwiched layer

Monoclinic materials viz. quartz, lithium niobate and lithium tantalate are among the most abundant materials, finding numerous applications throughout the technological world. Moreover, the presence of irregularity, initial stresses, anisotropy and heterogeneity in a material medium is obvious. These facts motivate the study of magnetoelastic SH-wave propagation in an irregular monoclinic sandwiched layer between a heterogeneous isotropic layer and an isotropic half space, all under initial stress. The heterogeneity in the uppermost layer is caused due to exponential variation in rigidity, density and initial stress in terms of space variable pointing vertically downward. The dispersion relation has been obtained using first-order perturbation technique. The substantial effect of wave number, anisotropy, irregularity, width ratio of the layers, horizontal compressive/tensile initial stresses, heterogeneity and monoclinic-magnetoelastic coupling parameter associated with sandwiched layer on phase velocity of SH-wave has been studied and depicted by means of graph. Comparative study made for the case when pre-stressed irregular sandwiched layer is monoclinic-magnetoelastic to the case when it is isotropic magnetoelastic layer is one of the major highlights of the current study.     

Temporal variation of b value associated with M ��4 earthquakes in the reservoir-triggered seismic environment of the Koyna�CWarna region, Western India

It is generally found that the b values associated with reservoir-triggered seismicity (RTS) are higher than the regional b values in the frequency magnitude relation of earthquakes. In the present study, temporal and spatial variation of b value is investigated using a catalog of 3,000 earthquakes from August 2005 through December 2010 for the Koyna?CWarna region in Western India, which is a classical site of RTS globally. It is an isolated (30?×?20?km²) zone of seismicity where earthquakes of up to M ??5 are found to occur during phases of loading and unloading of the Koyna and Warna reservoirs situated 25?km apart. For the Warna region, it is found that low b values of 0.6?C0.9 are associated with earthquakes of M ??4 during the loading phase. The percentage correlation of the occurrence of an M????4 earthquake with a low b value outside the 1?? or 2?? level is as high as 78?%. A drastic drop in the b value of about 50?% being reported for an RTS site may be an important precursory parameter for short-term earthquake forecast in the future.     

Mapping of rock types using a joint approach by combining the multivariate statistics,self-organizing map and Bayesian neural networks: an example from IODP 323 site

Modeling and classification of the subsurface lithology is very important to understand the evolution of the earth system. However, precise classification and mapping of lithology using a single framework are difficult due to the complexity and the nonlinearity of the problem driven by limited core sample information. Here, we implement a joint approach by combining the unsupervised and the supervised methods in a single framework for better classification and mapping of rock types. In the unsupervised method, we use the principal component analysis (PCA), K-means cluster analysis (K-means), dendrogram analysis, Fuzzy C-means (FCM) cluster analysis and self-organizing map (SOM). In the supervised method, we use the Bayesian neural networks (BNN) optimized by the Hybrid Monte Carlo (HMC) (BNN-HMC) and the scaled conjugate gradient (SCG) (BNN-SCG) techniques. We use P-wave velocity, density, neutron porosity, resistivity and gamma ray logs of the well U1343E of the Integrated Ocean Drilling Program (IODP) Expedition 323 in the Bering Sea slope region. While the SOM algorithm allows us to visualize the clustering results in spatial domain, the combined classification schemes (supervised and unsupervised) uncover the different patterns of lithology such of as clayey-silt, diatom-silt and silty-clay from an un-cored section of the drilled hole. In addition, the BNN approach is capable of estimating uncertainty in the predictive modeling of three types of rocks over the entire lithology section at site U1343. Alternate succession of clayey-silt, diatom-silt and silty-clay may be representative of crustal inhomogeneity in general and thus could be a basis for detail study related to the productivity of methane gas in the oceans worldwide. Moreover, at the 530 m depth down below seafloor (DSF), the transition from Pliocene to Pleistocene could be linked to lithological alternation between the clayey-silt and the diatom-silt. The present results could provide the basis for the detailed study to get deeper insight into the Bering Sea’ sediment deposition and sequence.     

Earthquake chemical precursors in groundwater: a review

We review changes in groundwater chemistry as precursory signs for earthquakes. In particular, we discuss pH, total dissolved solids (TDS), electrical conductivity, and dissolved gases in relation to their significance for earthquake prediction or forecasting. These parameters are widely believed to vary in response to seismic and pre-seismic activity. However, the same parameters also vary in response to non-seismic processes. The inability to reliably distinguish between changes caused by seismic or pre-seismic activities from changes caused by non-seismic activities has impeded progress in earthquake science. Short-term earthquake prediction is unlikely to be achieved, however, by pH, TDS, electrical conductivity, and dissolved gas measurements alone. On the other hand, the production of free hydroxyl radicals (?OH), subsequent reactions such as formation of H₂O₂ and oxidation of As(III) to As(V) in groundwater, have distinctive precursory characteristics. This study deviates from the prevailing mechanical mantra. It addresses earthquake-related non-seismic mechanisms, but focused on the stress-induced electrification of rocks, the generation of positive hole charge carriers and their long-distance propagation through the rock column, plus on electrochemical processes at the rock-water interface.     

Investigation of Hemispherical Variations of Soft X-Ray Solar Flares during Solar Cycles 21 to 24

Solar System Research - In the current investigation we have studied the distribution as well as the asymmetry of solar X-ray flares during the period 1976–2017 which corresponds to solar...     

Indicators of Great Basin rangeland health

Early-warning indicators of rangeland health can be used to estimate the functional integrity of a site and may allow sustainable management of desert rangelands. The utility of several vegetation canopy-based indicators of rangeland health at 32 Great Basin rangeland locations was investigated. The indicators were originally developed in rangelands of the Chihuahuan Desert. Soil resources are lost through wind and water-driven erosion mainly from areas unprotected by plant canopies (i.e. bare soil). Study sites in Idaho had the smallest bare patches, followed by sites in Oregon. The more arid Great Basin Sagebrush Zone sites in Utah had the largest bare patches. Several vegetational indicators including percent cover by vegetation, percent cover by life-form, percent cover by sagebrush, and percent cover by resilient species were negatively related to mean bare patch size and are potential indicators of Great Basin rangeland condition. Plant community composition and the range of bare patch sizes were different at sites in the three locations in Idaho, Oregon and Utah. Therefore, expected indicator values are location specific and should not be extrapolated to other locations. The condition of study sites were often ranked differently by different indicators. Therefore, the condition of rangeland sites should be evaluated using several indicators.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

Effect of Corrugation and Reinforcement on the Dispersion of SH-wave Propagation in Corrugated Poroelastic Layer Lying over a Fibre-reinforced Half-space

The presence of porosity and reinforcement in a medium is an important factor affecting seismic wave propagation and plays vital role in many geophysical prospects. Also, the presence of salt and ore deposits, mountains, basins, mountain roots, etc. is responsible for the existence of corrugated boundary surfaces of constituent layers. Such facts brought motivation for the present paper which deals with the propagation of SH-wave in a heterogeneous fluid-saturated poroelastic layer with corrugated boundaries lying over an initially stressed fibre-reinforced elastic halfspace. Closed form of dispersion relation has been obtained and is found in well agreement to classical Love wave equation for isotropic case. The effect of corrugation, wave number, undulation, position parameter, horizontal compressive/tensile initial stress and heterogeneity on phase velocity has been analysed through numerical computation and graphical illustration. Moreover, comparative study exploring the effect of presence and absence of reinforcement in half-space on dispersion curve is the major highlight of the current study.