An improved empirical green function method to extract microseismic signals

In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the routine methods for arrival picking are insufficiently efficient due to these strong influences. Here, we modified the empirical green function method and applied it to extract the green functions and the radiation coefficients of microseismic events with the largest energy in clusters.Multiple-channel records were used to estimate an "average" source spectrum and then the "average" source wavelet was removed from the records by deconvolution. We applied this method to the real data,and the result indicated clear improvement in extracting the dominant event of the clusters.

Source time functions of the Gonghe，China earthquake retrieved from long－period digital waveform data using empirical Green＇s function technique

An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.

Rupture process of November 6,1988,Lan-cang-Gengma,Yunnan,China,earthquake of M_s=7.6 using empirical Green’s function de-convolution method

The empirical Green's function deconvolution technique is applied to the November 6, 1988, Ms=7.6, Yunnan, China earthquake to study the rupture process of this event. An aftershock of Ms=6.3 is taken as the empirical Green's function. Recordings of these two events at six stations of China Digital Seismographic Network (CDSN) are employed. Deconvolution results show that this event is a relatively simple event. Directivity analysis indicates that the rupture was initiated at hypocenter and propagated bilaterally and nearly symmetrically towards the northwest and southeast directions with a total length of 70 km and a time duration of 19 s. The rupture velocity is estimated to be about 2.0 km/s.

A hybrid method integrating Green's function Monte Carlo and projected entangled pair states

This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function in GFMC.By leveraging PEPS’s proficiency in capturing quantum state entanglement and GFMC’s efficient parallel architecture,the hybrid method is well-suited for the accurate and efficient treatment of frustrated quantum spin systems.As a benchmark,we applied this approach to study the frustrated J_(1)–J_(2) Heisenberg model on a square lattice with periodic boundary conditions(PBCs).Compared with other numerical methods,our approach integrating PEPS and GFMC shows competitive accuracy in the performance of ground-state energy.This paper provides systematic and comprehensive discussion of the approach of our previous work[Phys.Rev.B 109235133(2024)].

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Existence of Monotone Positive Solution for a Fourth-Order Three-Point BVP with Sign-Changing Green’s Function

This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.

A NEW ONE-STEP BAND-LIMITED EXTRAPOLATION PROCEDURE USING EMPIRICAL ORTHOGONAL FUNCTIONS

A one-step band-limited extrapolation procedure is systematically developed under an a priori assumption of bandwidth. The rationale of the proposed scheme is to expand the known signal segment based on a band-limited basis function set and then to generate a set of Empirical Orthogonal Functions (EOF’s) adaptively from the sample values of the band-limited function set. Simulation results indicate that, in addi- tion to the attractive adaptive feature, this scheme also appears to guarantee a smooth result for inexact data, thus suggesting the robustness of the proposed procedure.

Singularity-free Green's function for EM sources embedded in a stratified medium

We present a method to unify the calculation of Green＇s functions for an electromagnetic（EM） transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green＇s function are derived by decomposing the partial wave solutions to Helmholtz＇s equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method（MCSEM） demonstrates its simplicity and flexibility.

The Green's function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green(s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.

纳米分子器件量子传导中的Green's Function模型

将纳米分子器件的电极和分子本体作为一个系统,考虑到入射电子在原子电极和分子本体中的传导,在电极与分子之间界面上的透射的量子状态,在满足Fisher-Lee关系式的基础上,一个新型的Green’sFunction被设计出来。这一理论模型把三个区域上的不同量子状态完整地表达在一个方程式里,符合物理学基本原理,可以系统地用来研究纳米分子器件的量子传导特性。

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

Green's functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.

NONLINEAR EVOLUTION SYSTEMS AND GREEN’S FUNCTION

In this paper, we will introduce how to apply Green＇s function method to get the pointwise estimates for the solutions of Cauchy problem of nonlinear evolution equations with dissipative structure. First of all, we introduce the pointwise estimates of the time-asymptotic shape of the solutions of the isentropic Navier-Stokes equations and show to exhibit the generalized Huygen＇s principle. Then, for other nonlinear dissipative evolution equations, we will only introduce the result and give some brief explanations. Our approach is based on the detailed analysis of the Green＇s function of the linearized system and micro-local analysis, such as frequency decomposition and so on.

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

In-plane dynamic Green's functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.

Searching for single-particle resonances with the Green’s function method

Single-particle resonances in the continuum are crucial for studies of exotic nuclei.In this study,the Green’s function approach is employed to search for single-particle resonances based on the relativistic-mean-field model.Taking^(120)Sn as an example,we identify singleparticle resonances and determine the energies and widths directly by probing the extrema of the Green’s functions.In contrast to the results found by exploring for the extremum of the density of states proposed in our recent study[Chin.Phys.C,44:084105(2020)],which has proven to be very successful,the same resonances as well as very close energies and widths are obtained.By comparing the Green’s functions plotted in different coordinate space sizes,we also found that the results very slightly depend on the space size.These findings demonstrate that the approach by exploring for the extremum of the Green’s function is also very reliable and effective for identifying resonant states,regardless of whether they are wide or narrow.

TIME-HARMONIC DYNAMIC GREEN'S FUNCTIONS FOR ONE-DIMENSIONAL HEXAGONAL QUASICRYSTALS

Quasicrystals have additional phason degrees of freedom not found in conventional crystals. In this paper, we present an exact solution for time-harmonic dynamic Green＇s function of one-dimensional hexagonal quasicrystals with the Laue classes 6/mh and 6/mhmm. Through the introduction of two new functions φ and ψ, the original problem is reduced to the determination of Green＇s functions for two independent Helmholtz equations. The explicit expressions of displacement and stress fields are presented and their asymptotic behaviors are discussed. The static Green＇s function can be obtained by letting the circular frequency approach zero.

Decoupling coefficients of dilatational wave for Biot's dynamic equation and its Green's functions in frequency domain

Green＇s functions for Blot＇s dynamic equation in the frequency domain can be a highly useful tool for the investigation of dynamic responses of a saturated porous medium. Its applications are found in soil dynamics, seismology, earthquake engineering, rock mechanics, geophysics, and acoustics. However, the mathematical work for deriving it can be daunting. Green＇s functions have been presented utilizing an analogy between the dynamic thermoelasticity and the dynamic poroelasticity in the frequency domain using the u-p formulation. In this work, a special term ＂decoupling coefficient＂ for the decomposition of the fast and slow dilatational waves is proposed and expressed to present a new methodology for deriving the poroelastodynamic Green＇s functions. The correct- ness of the solution is demonstrated by numerically comparing the current solution with Cheng＇s previous solution. The separation of the two waves in the present methodology allows the more accurate evaluation of Green＇s functions, particularly the solution of the slow dilatational wave. This can be advantageous for the numerical implementation of the boundary element method （BEM） and other applications.

Solution of pavement temperature field in“Environment-Surface”system through Green's function

In order to simplify the boundary conditions of pavement temperature field,the "Environment-Surface" system which considered the natural environment and pavement surface was established.Based on this system,the partial differential equations of the one-dimensional heat conduction in the pavement were established on the basis of the heat transfer theory.Furthermore,the function forms of the initial and boundary conditions of the equations were created through the field experiments.The general solution of the pavement one-dimensional heat conduction partial differential equations was acquired by using Green's function,and the explicit expression of pavement temperature field under specific constraint conditions was derived.For the purpose of analysis,the pavement temperatures in different seasons were calculated using the explicit expression of pavement temperature field,and the calculation accuracy was analyzed through the comparison between measured and calculated values.Then,the relationship between fitting accuracy and calculation accuracy of pavement temperatures was analyzed.The analysis results show that: the usage of "Environment-Surface" system simplifies the calculation of pavement temperature field; the relative error between calculated and measured values is generally less than 7% and is seldom influenced by seasons; there is a positive correlation between the calculation accuracy and the fitting accuracy of pavement surface temperature; high fitting accuracy would result in less error of pavement temperature prediction.