Scattering and diffraction of plane P-waves in a 2-D elastic half-space Ⅱ：shallow arbitrary shaped canyon

Scattering and Diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P-and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P-and SV-scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.

Intrinsic and attenuative dispersion characteristics of direct P-waves in and near the source area of the 1999 MW7.6 Chi-Chi,Taiwan,earth-quake before and after the mainshock

Based on the measurement of the arrival time of maxima magnitude from band-pass filtering signals which were determined using a new Morlet wavelet multiple-filter method, we develop a method for measuring intrinsic and attenuative dispersion of the first cycle direct P-wave. We determine relative group delays of spectral components of direct P-waves for 984 ray paths from SML and ALS stations of the Taiwan Central Weather Bureau Seismic Network （CWBSN）. Using continuous relaxation model, we deduce a new transfer function that relates intrinsic dispersion to attenuation. Based on the genetic algorithm （GA）, we put forward a new inversion procedure for determining which is defined the flat part of quality factor Q（ω） spectrum, τ1 and τ2 parameters. The results indicate that ① The distribution of Om values versus epicentral distance and depth show that Qm values linearly increase with increasing of epicentral distance and depth, and Qm values is clearly independent of earthquakes magnitude; ② In the different depth ranges, Qm residual show no correlation with variations in epicentral distance. Some significant changes of Qm residual with time is likely caused by pre-seismic stress accumulation, and associated with fluid-filled higher density fractures rock volume in the source area of 1999 Chi-Chi Taiwan earthquake. We confirm that Qm residual with time anomaly appears about 2.5 years before the Chi-Chi earthquake; ③ A comparison of Qm residual for different depth range between SML and ALS stations show that the level of stress has vertical and lateral difference; ④ The area near observation station with both anomalously increasing and decreasing averaged Qm residual is likely an unstable environment for future strong earthquake occurrence. This study demonstrates the capability of direct P-waves dispersion for monitoring attenuation characteristics and its state changes of anelastic medium of the Earth at short propagation distance using seismograms recorded from very small events.

Surface Motion of Alluvial Valley in Layered Half-Space for Incident Plane P-Waves

The indirect boundary element method (IBEM) is used to study the surface motion of an alluvial valley in layered half-space for incident plane P-waves based on Wolf’s theory. Firstly, the free field response can be solved by the direct stiffness method, and the scattering wave response is calculated by Green’s functions of distributed loads acting on inclined lines in a layered half-space. The method is verified by comparing its results with literature and numerical analyses are performed by taking the amplification of incident plane P-waves by an alluvial valley in one soil layer resting on bedrock as an example. The results show that there exist distinct differences between the wave amplification by an alluvial valley embedded in layered half-space and that in homogeneous half-space and there is interaction between the valley and the soil layer. The amplitudes are relatively large when incident frequencies are close to the soil layer’s resonant frequencies.

The computation of a finite-frequency travel time sensitive kernel for P-waves in the AK135 earth model

Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for computing the traveltime sensitive kernel many times.It is a time-consuming task.It is easy and fast to get the travel time from analytic equations in a simple model such as a homogenous or linear velocity media.However,most of the earth models are layered.It is cumbersome to get the travel time from analytic equations.In order to enhance the computation efficiency,we used the table look-up method to compute the finite-frequency travel time sensitive kernel for P-waves in a layered structure model.We chose the AK135 earth model for the velocity model.The table look-up method saved about 50% of the computation time.We enhanced the computation speed by using the table lookup method in the same velocity model,which was very useful for enhancing the computation efficiency for the finite-frequency travel time tomography.

Remnant Paleoproterozoic Subduction or Lithospheric Drip Initiation at the Yilgarn Craton Margin:Constraints from P-wave Tomography

The Paleoproterozoic was a critical time in whether modern-style plate tectonics had become globally dominant(e.g.,Wan et al.,2020).The Capricorn Orogen witnessed the assembly of the Pilbara and Yilgarn Cratons and an exotic microcontinent,the Glenburgh Terrane,to form the West Australia Craton(WAC)through two collisional orogenic events,the 2215–2145 Ma Ophthalmian and 2005–1950 Ma Glenburgh Orogenies(Johnson et al.,2013;Fig.1).Compared to other Proterozoic orogenic belts in Australia,the Capricorn Orogen preserves‘complete'opposing continental margin successions,together with intervening arc fragments associated with oceanic closure and foreland basins associated with collisional loading(Cawood et al.,2009).

Effect of the mixing of s-wave and chiral p-wave pairings on electrical shot noise properties of normal metal/superconductor tunnel junctions

We study theoretically the electrical shot noise properties of tunnel junctions between a normal metal and a superconductor with the mixture of singlet s-wave and chiral triplet p-wave pairing due to broken inversion symmetry. We investigate how the shot noise properties vary as the relative amplitude between the two parity components in the pairing potential is changed. It is demonstrated that some characteristics of the electrical shot noise properties of such tunnel junctions may depend sensitively on the relative amplitude between the two parity components in the pairing potential, and some significant changes may occur in the electrical shot noise properties when the relative amplitude between the two parity components is varied from the singlet s-wave pairing dominated regime to the chiral triplet p-wave pairing dominated regime. In the chiral triplet p-wave pairing dominated regime, the ratio of noise power to electric current is close to 2e both in the in-gap and in the out-gap region. In the singlet s-wave pairing dominated regime, the value of this ratio is close to 4e in the inner gap region but may reduce to about 2e in the outer gap region as the relative amplitude of the chiral triplet pairing component is increased. The variations of the differential shot noise with the bias voltage also exhibit some significantly different features in different regimes. Such different features can serve as useful diagnostic tools for the determination of the relative magnitude of the two parity components in the pairing potential.

Impacts of anisotropy coefficient and porosity on the thermal conductivity and P-wave velocity of calcarenites used as building materials of historical monuments in Morocco

It is essential to study the porosity,thermal conductivity,and P-wave velocity of calcarenites,as well as the anisotropy coefficients of the thermal conductivity and P-wave velocity,for civil engineering,and conservation and restoration of historical monuments.This study focuses on measuring the thermal conductivity using the thermal conductivity scanning(TCS)technique and measuring the P-wave ve-locity using portable equipment.This was applied for some dry and saturated calcarenite samples in the horizontal and vertical directions(parallel and perpendicular to the bedding plane,respectively).The calcarenites were selected from some historical monuments in Morocco.These physical properties were measured in the laboratory to find a reliable relationship between all of these properties.As a result of the statistical analysis of the obtained data,excellent linear relationships were observed between the porosity and both the thermal conductivity and porosity.These relationships are characterized by relatively high coefficients of determination for the horizontal and vertical samples.Based on the thermal conductivity and P-wave velocity values in these two directions,the anisotropy coefficients of these two properties were calculated.The internal structure and the pore fabric of the calcarenite samples were delineated using scanning electron microscopy(SEM),while their chemical and mineral compositions were studied using the energy dispersive X-ray analysis(EDXA)and X-ray diffraction(XRD)techniques.

P-wave velocity structure beneath reservoirs and surrounding areas in the lower Jinsha River

The lower reaches of the Jinsha River are rich in hydropower resources because of the high mountains,deep valleys,and swift currents in this area.This region also features complex tectonic structures and frequent earthquakes.After the impoundment of the reservoirs,seismic activity increased significantly.Therefore,it is necessary to study the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds,thus providing seismological support for subsequent earthquake prevention and disaster reduction work in reservoir areas.In this study,we selected the data of 7.670 seismic events recorded by the seismic networks in Sichuan.Yunnan,and Chongqing and the temporary seismic arrays deployed nearby.We then applied the double-difference tomography method to this data,to obtain the P-wave velocity structure and earthquake locations in the lower reaches of the Jinsha River and surrounds.The results showed that the Jinsha River basin has a complex lateral P-wave velocity structure.Seismic events are mainly distributed in the transition zones between high-and low-velocity anomalies,and seismic events are particularly intense in the Xiluodu and Baihetan reservoir areas.Vertical cross-sections through the Xiangjiaba and Xiluodu reservoir areas revealed an apparent high-velocity anomaly at approximately 6 km depth:this high-velocity anomaly plays a role in stress accumulation,with few earthquakes distributed inside the high-velocity body.After the impoundment of the Baihetan reservoir,the number of earthquakes in the reservoir area increased significantly.The seismic events in the reservoir area north of 27°N were related to the enhanced activity of nearby faults after impoundment:the earthquakes in the reservoir area south of 27°N were probably induced by additional loads(or regional stress changes),and the multiple microseismic events may have been caused by rock rupture near the main faults under high pore pressure.

Development and prospect of acoustic reflection imaging logging processing and interpretation method

Acoustic reflection imaging logging technology can detect and evaluate the development of reflection anomalies,such as fractures,caves and faults,within a range of tens of meters from the wellbore,greatly expanding the application scope of well logging technology.This article reviews the development history of the technology and focuses on introducing key methods,software,and on-site applications of acoustic reflection imaging logging technology.Based on the analyses of major challenges faced by existing technologies,and in conjunction with the practical production requirements of oilfields,the further development directions of acoustic reflection imaging logging are proposed.Following the current approach that utilizes the reflection coefficients,derived from the computation of acoustic slowness and density,to perform seismic inversion constrained by well logging,the next frontier is to directly establish the forward and inverse relationships between the downhole measured reflection waves and the surface seismic reflection waves.It is essential to advance research in imaging of fractures within shale reservoirs,the assessment of hydraulic fracturing effectiveness,the study of geosteering while drilling,and the innovation in instruments of acoustic reflection imaging logging technology.

Real-time arrival picking of rock microfracture signals based on convolutional-recurrent neural network and its engineering application

Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.

Effect of crack aperture on P-wave velocity and dispersion

We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equal-radius penny-shaped crack inclusions in each layer. The P-wave velocity and amplitude were measured parallel and perpendicular to the layers of cracks at frequencies of 0.1 MHz to 1 MHz. The experiments show that as the crack aperture increases from 0.l mm to 0.34 mm, the amplitude of the P-waves parallel to the crack layers decreases linearly with increasing frequency and the P-wave velocity dispersion varies from 1.5% to 2.1%, whereas the amplitude of the P-wave perpendicular to the crack layers decreases quadratically with increasing frequency and the velocity dispersion varies from 1.9% to 4.7%. The variation in the velocity dispersion parallel and perpendicular to the cracks intensifies the anisotropy dispersion of the P-waves in the crack models （6.7% to 83%）. The P-wave dispersion strongly depends on the scattering characteristics of the crack apertures.

Permeability Estimation of Shale Oil Reservoir with Laboratory-derived Data: A Case Study of the Chang 7 Member in Ordos Basin

The shale oil reservoir within the Yanchang Formations of Ordos Basin harbors substantial oil and gas resources and has recently emerged as the primary focus of unconventional oil and gas exploration and development.Due to its complex pore and throat structure,pronounced heterogeneity,and tight reservoir characteristics,the techniques for conventional oil and gas exploration and production face challenges in comprehensive implementation,also indicating that as a vital parameter for evaluating the physical properties of a reservoir,permeability cannot be effectively estimated.This study selects 21 tight sandstone samples from the Q area within the shale oil formations of Ordos Basin.We systematically conduct the experiments to measure porosity,permeability,ultrasonic wave velocities,and resistivity at varying confining pressures.Results reveal that these measurements exhibit nonlinear changes in response to effective pressure.By using these experimental data and effective medium model,empirical relationships between P-and S-wave velocities,permeability and resistivity and effective pressure are established at logging and seismic scales.Furthermore,relationships between P-wave impedance and permeability,and resistivity and permeability are determined.A comparison between the predicted permeability and logging data demonstrates that the impedance–permeability relationship yields better results in contrast to those of resistivity–permeability relationship.These relationships are further applied to the seismic interpretation of shale oil reservoir in the target layer,enabling the permeability profile predictions based on inverse P-wave impedance.The predicted results are evaluated with actual production data,revealing a better agreement between predicted results and logging data and productivity.

Major methods of seismic anisotropy

Seismic anisotropy reveals that seismic wave velocity, amplitude, and other physical properties show variations in different directions, which can be divided into lattice-preferred orientation(LPO) and shape-preferred orientation(SPO) according to its physical mechanisms. The main methods for studying seismic anisotropy include shearwave splitting analysis, P-wave travel time inversion and surface-wave tomography, etc. There are some differences and correlations among these methods. Seismic anisotropy is an important way to reveal the dynamic processes of crust-mantle evolution, and it is significant for monitoring crustal stress changes and improve seismic exploration studies. With the help of long-term observation, the application of machine learning techniques and combining inversion based on multiple phases would become potential developments in seismic anisotropy studies. This may improve the understanding of complex seismic anisotropic models, such as multiple layers anisotropy with an oblique axis of symmetry.

CREDIT-X1local:A reference dataset for machine learning seismology from ChinArray in Southwest China

High-quality datasets are critical for the development of advanced machine-learning algorithms in seismology.Here,we present an earthquake dataset based on the ChinArray Phase I records(X1).ChinArray Phase I was deployed in the southern north-south seismic zone(20°N-32°N,95°E-110°E)in 2011-2013 using 355 portable broadband seismic stations.CREDIT-X1local,the first release of the ChinArray Reference Earthquake Dataset for Innovative Techniques(CREDIT),includes comprehensive information for the 105,455 local events that occurred in the southern north-south seismic zone during array observation,incorporating them into a single HDF5 file.Original 100-Hz sampled three-component waveforms are organized by event for stations within epicenter distances of 1,000 km,and records of≥200 s are included for each waveform.Two types of phase labels are provided.The first includes manually picked labels for 5,999 events with magnitudes≥2.0,providing 66,507 Pg,42,310 Sg,12,823 Pn,and 546 Sn phases.The second contains automatically labeled phases for 105,442 events with magnitudes of−1.6 to 7.6.These phases were picked using a recurrent neural network phase picker and screened using the corresponding travel time curves,resulting in 1,179,808 Pg,884,281 Sg,176,089 Pn,and 22,986 Sn phases.Additionally,first-motion polarities are included for 31,273 Pg phases.The event and station locations are provided,so that deep learning networks for both conventional phase picking and phase association can be trained and validated.The CREDIT-X1local dataset is the first million-scale dataset constructed from a dense seismic array,which is designed to support various multi-station deep-learning methods,high-precision focal mechanism inversion,and seismic tomography studies.Additionally,owing to the high seismicity in the southern north-south seismic zone in China,this dataset has great potential for future scientific discoveries.

The Magnetic Longitudinal (P-) Wave’s Propagation and Energy Models Underlying the Mechanisms of Its Capacity to Absorb Free Energy

The longitudinal wave term within Faraday’s law of electromagnetic induction (Faraday’s law) underwent recovery to ensure its suitability for theoretical derivation of the equation governing longitudinal electromagnetic (LEM) waves. The revised Maxwell’s equations include the crucial parameters being the attenuation time constants of magnetic vortex potential and electric vortex potential generated by external electromagnetic field within the propagation medium. Specific expressions for them are obtained through theoretical analysis. Subsequently, a model for propagating magnetic P-wave generated by the superposition of a left-handed photo and a right-handed photon in a vacuum is formulated based on reevaluated total current law and revised Faraday’s law, covering wave equations, energy equation, as well as propagation mode involving mutual induction and conversion between scalar magnetic field and vortex electric field. Furthermore, through theoretical derivations centered around magnetic P-wave, evidence was presented regarding its ability to absorb huge free energy through the entangled interaction between zero-point vacuum energy field and the torsion field produced by the vortex electric field.

Component azimuths of the CEArray stations estimated from P-wave particle motion

The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.

Classification and assessment of rock mass parameters in Choghart iron mine using P-wave velocity

Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The basic data required may be obtained from visual observation and laboratory or field tests.However,owing to the discontinuous and variable nature of rock masses,it is difficult for rock engineers to directly obtain the specific design parameters needed.As an alternative,the use of geophysical methods in geomechanics such as seismography may largely address this problem.In this study,25 seismic profiles with the total length of 543 m have been scanned to determine the geomechanical properties of the rock mass in blocks Ⅰ,Ⅲ and Ⅳ-2 of the Choghart iron mine.Moreover,rock joint measurements and sampling for laboratory tests were conducted.The results show that the rock mass rating（RMR） and Q values have a close relation with P-wave velocity parameters,including P-wave velocity in field(V;).P-wave velocity in the laboratory(V;) and the ratio of V;V;(i.e.K;= V;/V;.However,Q value,totally,has greater correlation coefficient and less error than the RMR,In addition,rock mass parameters including rock quality designation（RQD）,uniaxial compressive strength（UCS）,joint roughness coefficient（JRC） and Schmidt number（RN） show close relationship with P-wave velocity.An equation based on these parameters was obtained to estimate the P-wave velocity in the rock mass with a correlation coefficient of 91%.The velocities in two orthogonal directions and the results of joint study show that the wave velocity anisotropy in rock mass may be used as an efficient tool to assess the strong and weak directions in rock mass.

Inconsistency of changes in uniaxial compressive strength and P-wave velocity of sandstone after temperature treatments

It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and systematic analysis of the microstructure variation of rocks with temperature rising and corresponding propagation mechanism of elastic wave,the results show that(1)There are three different trends for the changes of UCS and P-wave velocity of sandstone when heated from room temperature(20C or 25C)to 800C:(i)Both the UCS and P-wave velocity decrease simultaneously;(ii)The UCS increases initially and then decreases,while the P-wave velocity decreases continuously;and(iii)The UCS increases initially and then fluctuates,while the P-wave velocity continuously decreases.(2)The UCS changes at room temperaturee400C,400Ce600C,and 600Ce800C are mainly attributed to the discrepancy of microstructure characteristics and quartz content,the transformation plasticity of clay minerals,and the balance between the thermal cementation and thermal damage,respectively.(3)The inconsistency in the trends of UCS and P-wave velocity changes is caused by the change of quartz content,phase transition of water and certain minerals.

Deep learning for P-wave arrival picking in earthquake early warning

Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up noise,missing P-waves and inaccurate P-wave arrival estimation.To address these issues,an automatic algorithm based on the convolution neural network(DPick)was developed,and trained with a moderate number of data sets of 17,717 accelerograms.Compared to the widely used approach of the short-term average/long-term average of signal characteristic function(STA/LTA),DPick is 1.6 times less likely to detect noise as a P-wave,and 76 times less likely to miss P-waves.In terms of estimating P-wave arrival time,when the detection task is completed within 1 s,DPick′s detection occurrence is 7.4 times that of STA/LTA in the 0.05 s error band,and 1.6 times when the error band is 0.10 s.This verified that the proposed method has the potential for wide applications in EEW.

A method to model the effect of pre-existing cracks on P-wave velocity in rocks

Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.