Precise source location of AE doublets by spectral matrix analysis of triaxial hodogram

Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 36-45 ◽  
Author(s):  
Hirokazu Moriya ◽  
Koji Nagano ◽  
Hiroaki Niitsuma

We have developed a precise relative source location technique using acoustic emission doublets (AE doublets) in the triaxial hodogram method to evaluate the direction and distance of subsurface extension cracks. An AE doublet is a pair of acoustic emissions with similar waveforms and adjacent locations on the same crack but which occur at different times. The relative source location is estimated by an analysis in the frequency domain. The relative distance between two AE sources is determined from the difference of P-S arrival time delays by cross‐spectrum analysis. The relative direction is derived using a spectral matrix from the difference in P‐wave polarization directions. We also propose a method to optimize the estimated relative location by using a group of AE doublets. The accuracy of the estimated source location was confirmed by performing field experiments. The relative locations of artificial wave sources about 150 m from a triaxial detector can be estimated with distance errors of less than 1 m, and direction errors of less than 3.8 degrees in both azimuth and inclination. Results of the application of this analysis on AE doublets in a geothermal field demonstrate its ability to evaluate deeper subsurface fractures.

2020 ◽  
Vol 110 (6) ◽  
pp. 3115-3129
Author(s):  
Yusuke Mukuhira ◽  
Oleg V. Poliannikov ◽  
Michael C. Fehler ◽  
Hirokazu Moriya

ABSTRACT Detection and analysis of small magnitude events is valuable for better characterization and understanding of reservoirs in addition to developing strategies for mitigating induced seismicity. Three-component (3C) receivers, which are now widely used, are commonly deployed in boreholes to provide continuous seismic data amenable to novel and powerful analysis. Using multicomponent continuous records of ground motion, we utilize two principal features of the direct P-wave arrival: (1) linearity and (2) polarization in the direction along the ray path to the source region to detect small magnitude events undetectable by conventional methods. We evaluate the linearity of polarization and direction of arrival in the time and frequency domains by introducing the spectral matrix analysis method and combine them into a scalar characteristic function that is thresholded for event detection purposes. We boost the signal-to-noise ratio by stacking characteristic functions obtained at different 3C receivers along an empirical moveout of a master event known to have occurred in an area of interest. This allows us to detect smaller events and spatially tie them to a relatively small area around the large event. We apply our method to field data recorded at the Groningen gas field in the Netherlands. Our method detects all catalog events as well as several previously undetected events.


Geophysics ◽  
1989 ◽  
Vol 54 (4) ◽  
pp. 508-513 ◽  
Author(s):  
K. Nagano ◽  
H. Niitsuma ◽  
N. Chubachi

An automatic acoustic emission (AE) source location algorithm has been developed for downhole AE measurement of subsurface cracks by using the triaxial hodogram method. The P-wave arrival time is detected by analyzing crosscorrelation coefficients among three components of AE signal energy; the P-wave direction is determined by the method of least squares. For detection of S-wave arrival time, a maximum‐likelihood method analyzes a distribution of instantaneous values of the SH-wave component amplitude. This algorithm can locate an AE source as accurately as human analysis. For field measurements, it takes less than 4 s to locate an AE source using a 16-bit personal computer with a program in C language. Automatic AE source location by the triaxial hodogram method has been realized with this algorithm.


Geophysics ◽  
1996 ◽  
Vol 61 (5) ◽  
pp. 1453-1466 ◽  
Author(s):  
Hirokazu Moriya ◽  
Hiroaki Niitsuma

We have developed a signal processing technique for three‐component microseismic data that allows the precise determination of P‐wave arrival times. The method is based on a time‐frequency representation of the signal that allows the evaluation of the 3-D particle motion from seismic waves in both time and frequency domains. A spectral matrix is constructed using the time‐frequency distributions. A crosscorrelation coefficient for the three‐component signal is derived through eigenvalue analysis of the spectral matrix. The P‐wave arrival time is determined through a statistical test of hypotheses using the crosscorrelation coefficient. This signal processing method is evaluated using a synthetic signal and it is compared to the local stationary autoregressive method for determining P‐wave arrival times. We also show that the proposed method is capable of determining the arrival time of a synthetic P‐wave to within 1 ms (five points in the discrete time series) in the presence of a signal‐to‐noise ratio of −5dB. The method can detect the arrival time of different frequency components of the P‐wave, which is a possibility for the evaluation of velocity dispersion of the seismic wave. We demonstrate the feasibility of the method further by applying it to microseismic data from a geothermal field.


Geophysics ◽  
2004 ◽  
Vol 69 (3) ◽  
pp. 681-689 ◽  
Author(s):  
J. Helen Isaac ◽  
Don C. Lawton

The location of any event imaged by P‐wave reflection seismic data beneath a tilted transversely isotropic (TTI) overburden is shifted laterally if isotropic velocities are used during data processing. The magnitude of the shift depends on five independent parameters: overburden thickness, angle of tilt, symmetry‐axis velocity, and the Thomsen anisotropy parameters and δ. The shift also varies with source–receiver offset. We have developed a procedure to estimate these five parameters when the tilt of the symmetry axis from the vertical is equal to the dip of the TTI layer (except in the special cases transverse isotropy with vertical or horizontal axis of symmetry). We observe three attributes of seismic data processed using isotropic velocities: the zero‐offset arrival time of a selected reflection, the difference in arrival time between a near‐offset and a far‐offset arrival, and the difference in imaged location (smear) of this target event between the same offsets. We then perform a cascaded scan of the five parameters to determine those combinations of the five that result in calculated attributes equivalent to the observed attributes. The multiple solutions are averaged to give the parameter estimates. Application of this method to synthetic and physical model reflection data results in multiple solutions, which are constrained and averaged to obtain the effective imaging parameters. These effective parameters are close estimates of the true model parameters in both cases. For field seismic data this procedure requires that there be a suitable observable event below the TTI overburden and assumes that the measured times and shifts are reasonably accurate.


Author(s):  
V. Dumych ◽  

The purpose of research: to improve the technology of growing flax in the Western region of Ukraine on the basis of the introduction of systems for minimizing tillage, which will increase the yield of trusts and seeds. Research methods: field, laboratory, visual and comparative calculation method. Research results: Field experiments included the study of three tillage systems (traditional, canning and mulching) and determining their impact on growth and development and yields of trusts and flax seeds. The traditional tillage system included the following operations: plowing with a reversible plow to a depth of 27 cm, cultivation with simultaneous harrowing and pre-sowing tillage. The conservation system is based on deep shelfless loosening of the soil and provided for chiseling to a depth of 40 cm, disking to a depth of 15 cm, cultivation with simultaneous harrowing, pre-sowing tillage. During the implementation of the mulching system, disking to a depth of 15 cm, cultivation with simultaneous harrowing and pre-sowing tillage with a combined unit was carried out. Tillage implements and machines were used to perform tillage operations: disc harrow BDVP-3,6, reversible plow PON-5/4, chisel PCh-3, cultivator KPSP-4, pre-sowing tillage unit LK-4. The SZ-3,6 ASTPA grain seeder was used for sowing long flax of the Kamenyar variety. Simultaneously with the sowing of flax seeds, local application of mineral fertilizers (nitroammophoska 2 c/ha) was carried out. The application of conservation tillage allows to obtain the yield of flax trust at the level of 3,5 t/ha, which is 0,4 t/ha (12.9 %) more than from the area of traditional tillage and 0,7 t/ha (25 %) in comparison with mulching. In the area with canning treatment, the seed yield was the highest and amounted to 0,64 t/ha. The difference between this option and traditional and mulching tillage reaches 0,06 t/ha (10,3 %) and 0.10 t/ha (18.5 %), respectively. Conclusions. Preservation tillage, which is based on shelf-free tillage to a depth of 40 cm and disking to a depth of 15 cm has a positive effect on plant growth and development, yield and quality of flax.


Water ◽  
2020 ◽  
Vol 13 (1) ◽  
pp. 52
Author(s):  
Xiang He ◽  
Sijia Wang ◽  
Bingjian Zhang

Dew is a common but important phenomenon. Though water is previously considered to be a threat to earthen heritage sites, artificial dew is showing potential in relic preservation. A model of dew prediction on earthen sites will be essential for developing preventive protection methods, but studies of dew formation processes on relics are limited. In this study, a two parameter model is proposed. It makes approximations according to the features of earthen heritage sites, assuming that a thin and steady air layer exists close to the air–solid interface. This semi-theoretical model was based on calculations of the mass transfer process in the air layer, and was validated by simulations of laboratory experiments (R > 0.9) as well as field experiments. Additionally, a numerical simulation, performed by the commercial software COMSOL, confirmed that the difference between fitting parameter δ and the thickness of assumed mass transfer field was not significant. This model will be helpful in developing automatic environmental control systems for stabilizing water and soluble salts, thus enhancing preventive protection of earthen heritage sites.


Water ◽  
2020 ◽  
Vol 12 (1) ◽  
pp. 292
Author(s):  
Kai Zhang ◽  
Meijian Bai ◽  
Yinong Li ◽  
Shaohui Zhang ◽  
Di Xu

The broadcast fertilization method is widely used under basin irrigation in China. A reasonable broadcast fertilization method can effectively improve application performance of fertilization and reduce pollution from non-point agricultural sources. In this study, firstly, a non-uniform broadcast fertilization method and a non-uniform application coefficient were proposed. The value of non-uniform application coefficient is defined in this paper. It represents the ratio of the difference between the maximum and the average fertilization amount of fertilizer applied on the basin surface to the average fertilization amount of fertilizer applied on the basin surface. Secondly, field experiments were conducted to study the movement characteristics of fertilizer under non-uniform broadcast fertilization for basin irrigation. Field experiment results showed that under the condition of basin irrigation, the non-uniform broadcast fertilization method could weaken the non-uniform distribution of fertilizer due to erosion and transport capacity of solid fertilizer by irrigation water flow, which could significantly improve the uniformity of soil solute content. Thirdly, the solute transport model for broadcast fertilization was corroborated by the field experiment results. The variation rule of fertilization performance with non-uniform application coefficient under different basin length and inflow rate was achieved by simulation. The simulation results showed that fertilization uniformity and fertilization storage efficiency increased first and then decreased with the increase of non-uniform application coefficient. In order to be practically applicable, suitable irrigation programs of non-uniform application coefficient under different basin length and inflow rate conditions were proposed by numerical simulation.


Soil Research ◽  
2000 ◽  
Vol 38 (2) ◽  
pp. 345 ◽  
Author(s):  
G. D. Schwenke ◽  
D. R. Mulligan ◽  
L. C. Bell

At Weipa, in Queensland, Australia, sown tree and shrub species sometimes fail to establish on bauxite-mined land, possibly because surface-soil organic matter declines during soil stripping and replacement. We devised 2 field experiments to investigate the links between soil rehabilitation operations, organic matter decline, and revegetation failure. Experiment 1 compared two routinely practiced operations, dual-strip (DS) and stockpile soil, with double-pass (DP), an alternative method, and subsoil only, an occasional result of the DS operation. Other treatments included variations in stripping-time, ripping-time, fertiliser rate, and cultivation. Dilution of topsoil with subsoil, low-grade bauxite, and ironstone accounted for the 46% decline of surface-soil (0–10 cm) organic C in DS compared with pre-strip soil. In contrast, organic C in the surface-soil (0–10 cm) of DP plots (25.0 t/ha) closely resembled the pre-strip area (28.6 t/ha). However, profile (0–60 cm) organic C did not differ between DS (91.5 t/ha), DP (107 t/ha), and pre-strip soil (89.9 t/ha). Eighteen months after plots were sown with native vegetation, surface-soil (0–10 cm) organic C had declined by an average of 9% across all plots. In Experiment 2, we measured the potential for post-rehabilitation decline of organic matter in hand-stripped and replaced soil columns that simulated the DS operation. Soils were incubated in situ without organic inputs. After 1 year’s incubation, organic C had declined by up to 26% and microbial biomass C by up to 61%. The difference in organic C decline between vegetated replaced soils (Expt 1) and bare replaced soils (Expt 2) showed that organic inputs affect levels of organic matter more than soil disturbance. Where topsoil was replaced at the top of the profile (DP) and not ploughed, inputs from volunteer native grasses balanced oxidation losses and organic C levels did not decline.


Sign in / Sign up

Export Citation Format

Share Document