An improved rotated staggered grid fi nite difference scheme in coal seam

A coal seam is thin compared to the wavelength of seismic waves and usually shows strong anisotropy.It may form special geological formations,such as goafs and collapses,in coal mines.The existence of these formations may lead to instability in numerical simulations of the goaf area in a coal seam.The calculation speed of simulations is always a factor that restricts the development of simulation techniques.To improve the accuracy and effi ciency of seismic numerical simulations of goaf areas,an improved vacuum method has been incorporated into a rotated staggered grid scheme and calculations implemented by combining parallel computing and task parallelism.This ensures that the proposed numerical simulation method can be utilized in a geological model with large differences in elastic parameters among layers and improve the performance of a parallel application by enabling the full use of processor resources to expedite the calculations.We set up anisotropic coal seam models and then analyze numerically the characteristics of synthetic seismograms and snapshots of diff erent goaf areas with or without collapse.The results show that the proposed method can accurately simulate the goaf area and the calculation method can run with a high speed and parallel efficiency.The research will further advance the technology of anisotropic seismic exploration in coal fi elds,provide data for seismic inversion and build a theoretical support for coal mine disaster prediction.

Time difference correlation between earthquake lights and seismic ground accelerations

Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was possible to obtain the fi rst fi lm recording on earthquake lights.Likewise,during the earthquakes in Ecuador in 2016 and in Mexico in 2017,two fi lms of the earthquake lights were recorded.These fi lm recordings have helped in the study of earthquake lights,both for their objectivity and for their informational content.Several causal mechanisms have been proposed to explain earthquake lights:piezoelectricity,radon emanation,fluid diffusion,friction-vaporization,positive holes and dipole currents,among others.In this work a time difference correlation between earthquake lights and seismic ground accelerations was found and we use both seismic data and fi lm recordings of earthquake lights to explain its origin.In the discussion section it is suggested that fracturing of rocks manifest itself to some extent in the form of static electricity producing earthquake lights through induction The induction model proposed is new and it can explain the formation of EQL,even if the earth’s crust has layers of large electrical resistivity.The model also explains the formation of seismic lights without the need for special conditions on the earth’s surface or in the atmosphere.A better understanding of the earthquake lights generation process can improve our understanding of seismicity and help in the prediction of earthquakes.

Q estimation based on logarithmic spectral areas with different high and low frequencies

When a seismic wave propagates through subsurface viscoelastic media,the formation absorbs the high-frequency energy of the seismic wave more strongly than the low-frequency energy.As the depth and the off set increase,the diff erence between the logarithmic spectral areas with high and low frequencies gradually increases.Based on this seismic wave characteristic,we have developed a novel Q-estimation method based on logarithmic spectral area diff erence of high and low frequencies(referred to as the LSAD_LH method).In this paper,we derive the theoretical relationship between the Q value and difference of logarithmic spectral areas with high and low frequencies and prove the applicability of the LSAD_LH method using a single-layer medium numerical model.To verify the sensitivity of the LSAD_LH method to bandwidth selection and noise,we compare the LSAD_LH method with two credible methods—the logarithmic spectral ratio(LSR)and logarithmic spectral area diff erence(LSAD)methods using a synthetic model containing the random noise.The results demonstrate that the LSAD_LH method is not highly dependent on bandwidth,and in terms of noise immunity,it is signifi cantly better than the LSR method and has the same advantages as the LSAD method.To further highlight the advantages of the LSAD_LH method,we apply the LSAD_LH and LSAD methods to the vertical seismic profi ling(VSP)numerical simulation of the multilayer media and the fi eld zero-off set VSP data.The application of the two cases proves the applicability of the LSAD_LH method and the accuracy of the high Q-value estimation relative to the LSAD method.Moreover,via the transmission coeffi cient,the LSAD_LH method overcomes the weakness of the LSAD method.

Multi-temporal NDVI analysis using UAV images of tree crowns in a northern Mexican pine-oak forest

The use of unmanned aerial vehicles(UAV)for forest monitoring has grown significantly in recent years,providing information with high spatial resolution and temporal versatility.UAV with multispectral sensors allow the use of indexes such as the normalized difference vegetation index(NDVI),which determines the vigor,physiological stress and photo synthetic activity of vegetation.This study aimed to analyze the spectral responses and variations of NDVI in tree crowns,as well as their correlation with climatic factors over the course of one year.The study area encompassed a 1.6-ha site in Durango,Mexico,where Pinus cembroides,Pinus engelmannii,and Quercus grisea coexist.Multispectral images were acquired with UAV and information on meteorological variables was obtained from NASA/POWER database.An ANOVA explored possible differences in NDVI among the three species.Pearson correlation was performed to identify the linear relationship between NDVI and meteorological variables.Significant differences in NDVI values were found at the genus level(Pinus and Quercus),possibly related to the physiological features of the species and their phenology.Quercus grisea had the lowest NDVI values throughout the year which may be attributed to its sensitivity to relative humidity and temperatures.Although the use of UAV with a multispectral sensor for NDVI monitoring allowed genera differentiation,in more complex forest analyses hyperspectral and LiDAR sensors should be integrated,as well other vegetation indexes be considered.

Phylogenetic relationship of Picea mongolica with other Picea species in the same area based on chloroplast gene variations

Picea mongolica is a conifer with a limited distribution in China and its taxonomic status is controversial.In order to explore genetic diff erences between P.mongolica and other nearby Picea species and to investigate its taxonomic status,phylogenetic relationships were analyzed between P.mongolica and Picea koraiensis,Picea meyeri and Picea wilsonii by three chloroplast gene sequences matK,chlB and atpA.The length of joint chloroplast sequence is 2379 bp.The fi fteen haplotypes were identifi ed by haplotype network analysis,among which two were major haplotypes and nine were unique.In addition,the genetic diversity of the sample collection species was inferred.Based on the haplotype network and Neighbor Joining phylogenetic tree analysis,P.mongolica was located on the basal clade of the phylogenetic tree which had more primitive taxa,and the genetic diversity of P.mongolica was highest.The signifi cant diff erences between P.mongolica and these other Picea species were identifi ed by this research.

3D prestack reverse time migration of ground penetrating radar data based on the normalized correlation imaging condition

The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation.

Role of Microalloyed Y and Gd in Improving the Corrosion Resistance of Rolled Mg-3Al-1Zn Alloy

The corrosion performances of rolled AZ31 alloy with addition of Y and Gd were comparatively investigated. The corrosion rates of AZ31 alloy modified with Y and Gd were 3.3 and 3.7 mm/y immersing for 24 h in 3.5 wt% NaCl solution, respectively, which were much lower than that of AZ31 alloy(23.6 mm/y). The remarkable improvement in corrosion performance by adding Y or Gd was ascribed to preferentially form less noble Al_(2)Y or Al_(6)GdMn_(6) and more protective corrosion product layer, resulting in the reduction in H_(2) evolution rate and the enhancement of passivity.