3D near-surface P-wave velocity structure imaging with Distributed Acoustic Sensing and electric hammer source

Distributed Acoustic Sensing(DAS) is an emerging technique for ultra-dense seismic observation, which provides a new method for high-resolution sub-surface seismic imaging. Recently a large number of linear DAS arrays have been used for two-dimensional S-wave near-surface imaging in urban areas. In order to explore the feasibility of three-dimensional(3D) structure imaging using a DAS array, we carried out an active source experiment at the Beijing National Earth Observatory. We deployed a 1 km optical cable in a rectangular shape, and the optical cable was recast into 250 sensors with a channel spacing of 4 m. The DAS array clearly recorded the P, S and surface waves generated by a hammer source. The first-arrival P wave travel times were first picked with a ShortTerm Average/Long-Term Average(STA/LTA) method and further manually checked. The P-wave signals recorded by the DAS are consistent with those recorded by the horizontal components of short-period seismometers. At shorter source-receiver distances, the picked P-wave arrivals from the DAS recording are consistent with vertical component recordings of seismometers, but they clearly lag behind the latter at greater distances.This is likely due to a combination of the signal-to-noise ratio and the polarization of the incoming wave. Then,we used the Tomo DD software to invert the 3D P-wave velocity structure for the uppermost 50 m with a resolution of 10 m. The inverted P-wave velocity structures agree well with the S-wave velocity structure previously obtained through ambient noise tomography. Our study indicates the feasibility of 3D near-surface imaging with the active source and DAS array. However, the inverted absolute velocity values at large depths may be biased due to potential time shifts between the DAS recording and seismometer at large source-receiver distances.

Weighted-elastic-wave interferometric imaging of microseismic source location

Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging method.

Influence of rock fractures on the amplitude of dipole-source refl ected shear wave

To identify refl ector fractures near borehole by using dipole-source refl ected-shearwave logging, we need to understand the relation between the amplitude of the refl ected shear wave and the source radiation, borehole conditions, and attenuation owing to the surrounding formations. To assess the effect of these factors on the amplitude of the refl ected waves, we first studied the radiation performance and radiation direction of the dipole source in fast, medium, and slow formations by using the asymptotic solution in the far fi eld of the borehole. Then, the relation between the fracture parameters, and the refl ected-shear-wave amplitude as well as the ratio of the refl ected-shear-wave amplitude to the direct-wave amplitude (relative amplitude, RA) was evaluated by the three-dimensional fi nite-difference (3D FDTD) method. Finally, the fracture detection capability of the dipole reflected-shear-wave logging tool in different formations was analyzed by using the RA. The results suggest that the radiation amplitude of the SH-wave in the slow formation is weaker than those in the fast and medium formations, and the amplitude of the refl ected shear wave is lower. However, the RA in the slow formation is close to or even greater than in the fast and medium formations, which means that dipole-source shear-wave logging has the same or even better fracture detection capability in the slow formation as in the fast and medium formations. In addition, when RA is small, there is a good correlation between the RA and the various fracture parameters in the different types of formation, which can be used in determining the lower limit of the fracture parameters identifi ed by refl ection logging.

An analytical solution to the scattering of cylindrical SH waves by a partially filled semi-circular alluvial valley: near-source site effects

The earth’s surface irregularities can substantially affect seismic waves and induce amplifi cations of ground motions. This study investigates whether and how the source characteristics affect the site amplifi cation effects. An analytical model of a line source of cylindrical waves impinging on an alluvial valley is proposed to link the source and site. The analytical solution to this problem proves one aspect of the strong effect of source on site amplifi cation, i.e., the wave curvature effect. It is found that the site amplifi cation depends on the source location, especially under conditions of a small source-to-site distance. Whether the displacement is amplifi ed or reduced and the size of the amplifi cation or reduction may be determined by the location of the source. It is suggested that traditional studies of site responses, which usually ignore the source effect, should be further improved by combining the source with site effects.

Characteristics of the Seismic Waves from a New Active Source Based on Methane Gaseous Detonation

Active seismic sources are critical for obtaining high resolution images of the subsurface.For active imaging in urban areas,environment friendly and green seismic sources are required.In present work,we introduce a new type of green active source based on the gaseous detonation of methane and oxygen.When fired in a closed container,the chemical reaction,i.e.gaseous detonation,will produce high pressure air over 150 MPa.Seismic waves are produced when high pressure air is quickly released to impact the surroundings.The first field experiment of this active source was carried out in December,2017 in Jingdezhen,Jiangxi Province,where a series of active sources were excited to explore their potential in mine exploration.In current work,we analyzed the seismic waves recorded by near-field accelerators and a dense short-period seismic array and compared them with those from a mobile airgun source,another kind of active source by releasing high pressure air into water.The results demonstrate that it can be used for high resolution near surface imaging.Firstly,the gaseous detonation productions are harmless CO2 and water,making it a green explosive source.Secondly,the dominant seismic frequencies are 10-80 Hz and a single shot can be recorded up to 15 km,making it suitable for local structure investigations.Thirdly,it can be excited in vertical wells,similar to traditional powder explosive sources.It can also act as an additional on-land active source to airgun sources,which requires a suitable water body as intermediate media to generate repeating signals.Moreover,the short duration and high frequency signature of the source signals make it safe with no damage to nearby buildings.These make it convenient to excite in urban areas.As a new explosive source,the excitation equipment and conditions,such as gas ratio,sink depth and air-releasing directions,need further investigation to improve seismic wave generation efficiency.

Varying Rossby Wave Trains from the Developing to Decaying Period of the Upper Atmospheric Heat Source over the Tibetan Plateau in Boreal Summer

This study demonstrates the two different Rossby wave train（RWT） patterns related to the developing/decaying upper atmospheric heat source over the Tibetan Plateau（TPUHS） in boreal summer. The results show that the summer TPUHS is dominated by quasi-biweekly variability, particularly from late July to mid-August when the subtropical jet steadily stays to the north of the TP. During the developing period of TPUHS events, the intensifying TPUHS corresponds to an anomalous upper-tropospheric high over the TP, which acts as the main source of a RWT that extends northeastward, via North China, the central Pacific and Alaska, to the northeastern Pacific region. This RWT breaks up while the anomalous high is temporarily replaced by an anomalous low due to the further deepened convective heating around the TPUHS peak. However, this anomalous low, though existing for only three to four days due to the counteracting dynamical effects of the persisting upper/lower divergence/convergence over the TP, acts as a new wave source to connect to an anomalous dynamical high over the Baikal region. Whilst the anomalous low is diminishing rapidly, this Baikal high becomes the main source of a new RWT, which develops eastward over the North Pacific region till around eight days after the TPUHS peak. Nevertheless, the anomaly centers along this decaying-TPUHS-related RWT mostly appear much weaker than those along the previous RWT.Therefore, their impacts on circulation and weather differ considerably from the developing to the decaying period of TPUHS events.

Study of the Change in Wave Velocity Ratio before and after Two Strong Earthquakes Using Airgun Source Data

Using the signals excited by the large-volume airgun source at the Binchuan transmitting seismic station from January to June,2016,arrival-time data was acquired at four stations near the epicenter of the Eryuan MS4.5 and MS4.0 earthquakes on February 8,2016,as well as the epicenter of the Yunlong MS5.0 and Eryuan MS4.6 earthquakes on May 18,2016 through the waveform cross-correlation technique.The wave velocity ratio of the four stations was calculated using the single-station method.At the same time,the b-value and the focal mechanism consistency parameters of the study area were also calculated.The results show that:(1)the wave velocity ratio of each station experienced a process of decline-recovery-fast rise before the two strong earthquakes,and a significant quasi-synchronous rapid rise occurred within 3-12 days before the earthquake;(2)the timing of the rapid rise of the wave velocity ratio of the four stations before the Yunlong MS5.0 and Eryuan MS4.6 earthquakes were related to the epicentral distance.The station which observed the earliest increase in rapid rise is the farthest one from the epicenter,and the station where the rapid rise appeared in the latest is closest to the epicenter;(3)the form of change of the wave velocity ratio before the earthquake was different between stations located at different directions in the epicenter area;(4)the b-value and the focal mechanism consistency parameter which is commonly used to characterize the stress level both showed a downward trend before the two strong earthquakes,and were consistent with the change in the wave velocity ratio.According to the preliminary analysis,the wave velocity ratio obtained by using airgun source can better reflect the change in the stress state of the underground medium.

Decomposition of the Unsteady Wave Patterns for Bessho form Translating-Pulsating Source Green Function

In order to interpret the physical feature of Bessho form translating-pulsating source Green function, the phase function is extracted from the integral representation and stationary-phase analysis is carried out in this paper. The complex characteristics of the integral variable and segmentation of the integral intervals are discussed in m complex plane. In θ space, the interval [-π/2+φ,-π/2+φ-iε] is dominant in the near-field flow, and there is a one-to-one correspondence between the real intervals in m space and the unsteady wave patterns in far field. If 4τ>1(τ is the Brard number), there are three kinds of propagation wave patterns such as ring-fan wave pattern, fan wave pattern and inner V wave pattern, and if 0<4τ<1, a ring wave pattern, an outer V and inner V wave pattern are presented in far field. The ring-fan or ring wave pattern corresponds to the interval [-π+α,-π/2+φ] for integral terms about k2, and the fan or outer V wave pattern and inner V wave pattern correspond to [-π+α,-π/2) and(-π/2,-π/2+φ] respectively for terms about k1. Numerical result shows that it is beneficial to decompose the unsteady wave patterns under the condition of τ≠0 by converting the integral variable θ to m. In addition, the constant-phase curve equations are derived when the source is performing only pulsating or translating.

Numerical study on wave loads and motions of two ships advancing in waves by using three-dimensional translating-pulsating source

A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.

SOURCE RADIATION AND RESPONSES OF WAVE PROPAGATION

Recordings of seismic waves propagating from earthquake source to a station at the earth’s surface are a system response function. The convolution operator in time domain can be simplified as a multiplication operator in frequency domain. We discuss in frequency domain the separation of source, path and site effects for global scaling of earthquake source radiation. Also discussed are source scaling model, faulting mechanism, and the H/V inversion problems with crustal and near surface structures. Gross features of apparent source spectra appear to be not much region dependent although there may be difference between tectonic styles within a region of tectonic mixture for which we need further study as data accumulate. Vertical spectra may be a better approach to approximate source radiation, as it has less crustal amplification effects than horizontal spectra. The H/V ratio is evidently a comprehensive indicator of amplification effects from near surface to deep structure. This gives it potential as an inversion tool to deduce site crustal structure.

Development of a seven-cell S-band standing-wave RF-deflecting cavity for Tsinghua Thomson scattering X-ray source

A 2856-MHz,π-mode,seven-cell standingwave deflecting cavity was designed and fabricated for bunch length measurement in Tsinghua Thomson scattering X-ray source(TTX)facility.This cavity was installed in the TTX and provided a deflecting voltage of 4.2 MV with an input power of 2.5 MW.Bunch length diagnoses of electron beams with energies up to 39 MeV have been performed.In this article,the RF design of the cavity using HFSS,fabrication,and RF test processes are reviewed.High-power operation with accelerated beams and calibration of the deflecting voltage are also presented.

Simulation and analysis of point-source surface wave fields

The complexity of near surface intensifies the diversity of seismic wave fields, which makes study on near surface wavefields important in many aspects. The strong absorption of low velocity layer can affect the resolution of seismic data, and free boundary can cause surface wave. Considering the above problems, we focus on the Rayleigh wavefields simulation using finite-difference wave equation of higher-order staggered grids and PML boundary conditions. Free boundary, buried source and overlying low velocity layer are taken into consideration and point explosion source is adopted. Through some numerical simulation with different parameters, we quantitatively analyze relationship between wave intensity and source depth, as well as the energy variation with propagation and obtain some practical knowledge and conclusions.

SOLITARY WAVE AND PERIODIC WAVE SOLUTIONS FOR THE NON-NEWTONIAN FILTRATION EQUATIONS WITH NONLINEAR SOURCES AND A TIME-VARYING DELAY

This paper is concerned with the non-Newtonian filtration equations with non- linear sources and a time-varying delay. By an extension of Mawhin＇s continuation theorem and some analysis methods, several sufficient conditions ensuring the existence of solitary wave and periodic wave solutions are obtained. Some corresponding results in the literature are improved and extended. An example is given to illustrate the effectiveness of our results.

Source-independent wave-equation based microseismic source location using traveltime inversion

A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.

Application of passive source surface-wave method in site engineering seismic survey

Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m under- ground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.

A Homogenized Function to Recover Wave Source by Solving a Small Scale Linear System of Differencing Equations

In order to recover unknown space-dependent function G(x)or unknown time-dependent function H(t)in the wave source F(x;t)=G(x)H(t),we develop a technique of homogenized function and differencing equations,which can significantly reduce the difficulty in the inverse wave source recovery problem,only needing to solve a few equations in the problem domain,since the initial condition/boundary conditions and a supplementary final time condition are satisfied automatically.As a consequence,the eigenfunctions are used to expand the trial solutions,and then a small scale linear system is solved to determine the expansion coefficients from the differencing equations.Because the ill-posedness of the inverse wave source problem is greatly reduced,the present method is accurate and stable against a large noise up to 50%,of which the numerical tests confirm the observation.

Seismoelectric waves in a borehole excited by an external explosive source

The conversion of energy between seismic and electromagnetic wave fields has been described by Pride's coupled equations in porous media. In this paper, the seismoelectric field excited by the explosive point source located at the outside of the borehole is studied. The scattering fields inside and outside a borehole are analyzed and deduced under the boundary conditions at the interface between fluid and porous media. The influences of the distance of the point source, multipole components of the eccentric explosive source, and the receiving position along the axis of vertical borehole, on the converted waves inside the borehole are all investigated. When the distance from the acoustic source to the axis of a borehole is far enough, the longitudinal and coseismic longitudinal wave packets dominate the acoustic and electric field, respectively. The three components of both electric field and magnetic field can be detected, and the radial electric field is mainly excited and converted by the dipole component. Owing to the existence of borehole, the electric fields and magnetic fields in the borehole are azimuthal. The distance from the point where the maximum amplitude of the axial components of electric field is recorded, to the origin of coordinate indicates the horizontal distance from the explosive source to the axis of vertical borehole.

Numerical Evaluation of Longitudinal Motions of Wigley Hulls Advancing in Waves by Using Bessho Form Translating-Pulsating Source Green's Function

In the framework of 3D potential flow theory, Bessho form translating-pulsating source Green's function in frequency domain is chosen as the integral kernel in this study and hybrid source-and-dipole distribution model of the boundary element method is applied to directly solve the velocity potential for advancing ship in regular waves. Numerical characteristics of the Green function show that the contribution of local-flow components to velocity potential is concentrated at the nearby source point area and the wave component dominates the magnitude of velocity potential in the far field. Two kinds of mathematical models, with or without local-flow components taken into account, are adopted to numerically calculate the longitudinal motions of Wigley hulls, which demonstrates the applicability of translating-pulsating source Green's function method for various ship forms. In addition, the mesh analysis of discrete surface is carried out from the perspective of ship-form characteristics. The study shows that the longitudinal motion results by the simplified model are somewhat greater than the experimental data in the resonant zone, and the model can be used as an effective tool to predict ship seakeeping properties. However, translating-pulsating source Green function method is only appropriate for the qualitative analysis of motion response in waves if the ship geometrical shape fails to satisfy the slender-body assumption.

The Rossby wave train patterns forced by shallower and deeper Tibetan Plateau atmospheric heat-source in summer in a linear baroclinic model

By using a linear baroclinic model(LBM),this study investigates the different Rossby wave train(RWT)patterns associated with the Tibetan Plateau(TP)upper-atmospheric heat source(TPUHS)that is anomalously shallower and deeper in boreal summer.Observational results indicate the different RWT patterns between the developing and decaying periods of synoptic TPUHS events,when the anomalous TPUHS develops from a relatively shallower to a deeper TP heat source.Based on the different vertical heating profiles between these two periods in observation,this study forces the LBM with prescribed TPUHS profiles to mimic a shallower and deeper summer TP heat source.The results show that the atmospheric responses to a shallower and deeper TPUHS do exhibit different RWT patterns that largely resemble those in observation.Namely,corresponding RWT pattern to a shallower TPUHS stretches from the TP to the west coast of America,while that to a deeper TPUHS extends from the TP region to Alaska.

The Effect of Spatial Structure Character of Heat Source on the Ray Path and the Evolution of Wave Energy of Meridional Wave Train

This paper studies correlations between the spatial structure character of thermal forcing and deformation and the amplitude of rays of meridional wave train. It is shown that if thermal forcing appears a meridional linear variation the rays of quasi-stationary planetary wave may propagate along oblique lines and if the meridional variability of heat source has second order term the rays show distinct deformation as a great circular route. Additionally, the inhomogeneous distribution may cause lower frequency oscillations in mid- and low-latitudes. The combination of zonal and meridional wave numbers and distributive character of heat source may form an inverse mechanism of variational trend of generized wave energy, reflecting in some degree the physical process of transition between meridional and zonal flow patterns.