Quantifying fracture geometry with X-ray tomography: Technique of Iterative Local Thresholding (TILT) for 3D image segmentation

This paper presents a new method—the Technique of Iterative Local Thresholding (TILT)—for processing 3D X-ray computed tomography (xCT) images for visualization and quantification of rock fractures. The TILT method includes the following advancements. First, custom masks are generated by a fracture-dilation procedure, which significantly amplifies the fracture signal on the intensity histogram used for local thresholding. Second, TILT is particularly well suited for fracture characterization in granular rocks because the multi-scale Hessian fracture (MHF) filter has been incorporated to distinguish fractures from pores in the rock matrix. Third, TILT wraps the thresholding and fracture isolation steps in an optimized iterative routine for binary segmentation, minimizing human intervention and enabling automated processing of large 3D datasets. As an illustrative example, we applied TILT to 3D xCT images of reacted and unreacted fractured limestone cores. Other segmentation methods were also applied to provide insights regarding variability in image processing. The results show that TILT significantly enhanced separability of grayscale intensities, outperformed the other methods in automation, and was successful in isolating fractures from the porous rock matrix. Because the other methods are more likely to misclassify fracture edges as void and/or have limited capacity in distinguishing fractures from pores, those methods estimated larger fracture volumes (up to 80 %), surface areas (up to 60 %), and roughness (up to a factor of 2). These differences in fracture geometry would lead to significant disparities in hydraulic permeability predictions, as determined by 2D flow simulations.     

A multiobjective steepest descent method with applications to optimal well control

Multiobjective optimization deals with mathematical optimization problems where two or more objective functions (cost functions) are to be optimized (maximized or minimized) simultaneously. In most cases of interest, the objective functions are in conflict, i.e., there does not exist a decision (design) vector (vector of optimization variables) at which every objective function takes on its optimal value. The solution of a multiobjective problem is commonly defined as a Pareto front, and any decision vector which maps to a point on the Pareto front is said to be Pareto optimal. We present an original derivation of an analytical expression for the steepest descent direction for multiobjective optimization for the case of two objectives. This leads to an algorithm which can be applied to obtain Pareto optimal points or, equivalently, points on the Pareto front when the problem is the minimization of two conflicting objectives. The method is in effect a generalization of the steepest descent algorithm for minimizing a single objective function. The steepest-descent multiobjective optimization algorithm is applied to obtain optimal well controls for two example problems where the two conflicting objectives are the maximization of the life-cycle (long-term) net-present-value (NPV) and the maximization of the short-term NPV. The results strongly suggest the multiobjective steepest-descent (MOSD) algorithm is more efficient than competing multiobjective optimization algorithms.     

The Blasting Test and Blasting Vibration Monitoring of Vertical Crater Retreat Mining Method in the Luohe Iron Mine

Because the land requisition and demolishing became difficult more and more, the mining scheme of Luohe iron mine was changed from caving method to filling method. In order to ensure the safety of the residence and the underground tunnel cavern within the mobile belt of the underground mining, the Luohe iron mine did the blasting test on the vertical crater retreat mining method and blasting vibration monitoring. The blasting experiments use common emulsified oil explosives and non-electric initiation system. The way of caved ore adopts the cutting groove and bench side. The NUBOX-6016 intelligent vibration monitor was chosen in the blasting vibration monitoring. Twice experiments on the blasting vibration monitoring were done on the surface or in the underground refuge cavern. The first test select the three monitoring points on the ground and the second select two monitoring points on the ground and a monitoring point in the underground refuge cavern. The blasting vibration monitoring data were conducted by the regression analysis in the Sodev’s empirical formula. The vibration attenuation formula about the underground blasting vibration transmitting in three directions is derived. The blasting test on the vertical crater retreat mining method and the blasting vibration were analyzed. It is estimated if the vibration damage possibly the surface buildings and related facilities of mine.     

Complexity measurement of regional groundwater resources system using improved Lempel-Ziv complexity algorithm

In this paper, the self-adaptive artificial fish swarm algorithm (SAAFSA) is used to optimize the coarse graining of segment numbers, which are used in the Lempel-Ziv complexity algorithm. This approach improves the Lempel-Ziv complexity (LZC) algorithm of equal probability coarse graining. As a case study, the complexities of monthly series of groundwater depth were analyzed at seven farms in the Hongxinglong Administration. GIS technology was used to create spatial distributions of monthly groundwater depths. A projection pursuit model based on the SAAFSA was established and used for complexity attribution analysis at selected farms with different degrees of complexity. The three selected farms, Hongqiling, 852, and Youyi, each represent a certain degree of complexity. Further analysis shows that precipitation, evaporation, temperature, and human activities are the primary factors that cause complexity variations in local groundwater depth. The results reveal the evolution of the complexity characteristics of local groundwater depth and provide scientific evidence for the need to effectively allocate regional water resources. Additionally, the proposed method can be applied in complexity analyses of other hydrologic features, as well as in research regarding nonlinear time series in economic, engineering, medical, and signal analyses.     

Consistent seasonal snow cover depth and duration variability over the Western Himalayas (WH)

Precipitation in solid form, i.e., snow, during winter season over the Western Himalayas (WH) leads to the build-up of seasonal snow cover. Seasonal snow cover build-up (snow cover depth and duration) largely depends on atmospheric variables such as temperature, precipitation, radiation, wind, etc. Integrated (combined) influence of atmospheric variables on seasonal snow cover gets reflected in terms of spatial and temporal variability in seasonal snow cover build-up pattern. Hence spatial and temporal variability of seasonal snow cover build-up can serve as a good indicator of climate change in high altitude mountainous regions like the WH. Consistent seasonal snow cover depth and duration, delay days and early melt days of consistent seasonal snow cover at 11 stations spread across different mountain ranges over the WH were analyzed. Mean, maximum and percentiles (25th, 50th, 75th, 90th and 95th) of consistent seasonal snow cover depth and duration show decline over the WH in the recent past 2–3 decades. Consistent seasonal snow cover is found to melt early and snow cover build-up pattern is found to show changes over the WH. Decline in consistent seasonal snow cover depth, duration and changing snow cover build-up pattern over the WH in recent decades indicate that WH has undergone considerable climate change and winter weather patterns are changing in the WH.     

Performance of seismic arrays in the presence of weathering layer variations

Near-surface layer variations can degrade the desired response of seismic arrays which are typically used to attenuate horizontally traveling coherent noise and enhance vertically traveling signal. We investigate the effect of variations in the near-surface layer thickness on the performance of arrays by studying their impulse and wavelet responses. The models considered include the topographic variations and a channel in the base of weathering layer. The topographic variations include a dipping surface layer as well as a surface layer that follows a sine wave. The geological channel is assumed to be present under the entire receiver array as well as partially under a few receivers. We use Ricker wavelet and model plane wavefronts with incidence angles (90, 70, 45, 20, and 5 °) on a 12-element equally weighted array for the weathering layer models. We found that the array responses are more degraded for near-vertically traveling waves in all cases. The array responses are also found to be more degraded when channel variations are present underneath a few receivers as compared to the entire array length. We recommend designing a spatial filter or using single-sensor records after correcting for intra-array statics to preserve the reflection amplitudes in the presence of weathering layer variations.     

Influence of magnetic fabric anisotropy on seismic wave velocity in paramagnetic granites from NW Himalaya: Results from preliminary investigations

Anisotropy of Magnetic Susceptibility (AMS) and seismic wave velocity studies of some paramagnetic Himalayan granitoids show good correlation between magnetic fabric anisotropy and P wave velocity (Vp). Vp shows strong positive correlation with magnetic lineation (L) and degree of magnetic anisotropy (P′) having correlation coefficient (r) values of 0.93 and 0.89 respectively. Both Vp and Vs show positive correlation with the SiO₂ content of Proterozoic and Paleozoic granitoids. Velocity of S wave (Vs) shows negative correlation with mean magnetic susceptibility (K_m) having ‘r’ value of 0.86. The correlation between Vs-K_m, V_p-P′, V_p-L also shows >95% probability in Spearman’s rank correlation. Based on the results from the present sample size it is suggested that, in paramagnetic granites, V_p is proportional to intensity of deformation and preferred orientation of minerals as well as the mineralogy. On the other hand, Vs is more dependent on the mineralogy alone.     

Precambrian crystalline basement in southern Mongolia as revealed by SHRIMP zircon dating

Single zircon ages determined by ion microprobe (SHRIMP II) for granitoid gneisses from the southern slope of the Baga Bogd massif (Gobi-Altai, southern Mongolia) reveal several episodes of zircon growth, ranging from late Palaeoproterozoic to late Cambrian. The oldest events are documented by a zircon crystallization age for a gneiss protolith at 1519 ± 11 Ma and by a xenocrystic zircon from a dark grey augen-gneiss yielding an age of c. 1701 Ma. Discrete igneous events are recorded in granite-gneisses with protolith emplacement ages of 983 ± 6, 956 ± 3 and 954 ± 8 Ma. These ages provide the first record of early Neoproterozoic magmatic activity in this region. A much younger and discrete magmatic event is recorded by several dioritic to granitic orthogneisses which are tectonically interlayered with the older gneisses and have protolith emplacement ages between 502 and 498 Ma. These late Cambrian granitoids of calc-alkaline affinity are likely to have been emplaced along an active continental margin and suggest that the Baga Bogd Precambrian crustal fragment was either docked against the southward (present-day coordinates) growing margin of the CAOB or was a large enough crustal entity to develop an arc along its margin. We speculate that the Precambrian gneisses of this massif may be part of a crustal fragment rifted off the Tarim Craton.