Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations.Localization based on the relative counts of different detectors has b...Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations.Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors.There are two major methods for this count distribution localization:χ^(2)minimization method and the Bayesian method.Here we propose a modified Bayesian method that could take advantage of both the accuracy of the Bayesian method and the simplicity of the χ^(2)method.With comprehensive simulations,we find that our Bayesian method with Poisson likelihood is generally more applicable for various bursts than the χ^(2)method,especially for weak bursts.We further proposed a location-spectrum iteration approach based on the Bayesian inference,which could alleviate the problems caused by the spectral difference between the burst and location templates.Our method is very suitable for scenarios with limited computation resources or timesensitive applications,such as in-flight localization software,and low-latency localization for rapidly follow-up observations.展开更多
In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effec...In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effects of three-dimensional gas flow and acoustic damping on the acoustic attenuation characteristics of marine engine exhaust silencers, a dual reciprocity boundary element method (DRBEM) was developed. The acoustic governing equation in three-dimensional potential flow was derived first, and then the DRBEM numerical procedure is given. Compared to the conventional boundary element method (CBEM), the DRBEM considers the second order terms of flow Mach number in the acoustic governing equation, so it is suitable for the cases with higher Mach number subsonic flow. For complex exhaust silencers, it is difficult to apply the single-domain boundary element method, so a substructure approach based on the dual reciprocity boundary element method is presented. The experiments for measuring transmission loss of silencers are conducted, and the experimental setup and measurements are explained. The transmission loss of a single expansion chamber silencer with extended inlet and outlet were predicted by DRBEM and compared with the measurements. The good agreements between predictions and measurements are observed, which demonstrated that the derived acoustic governing equation and the DRBEM numerical procedure in the present study are correct.展开更多
Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burni...Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burning,a wave attenuation method is proposed to assess the nozzle damping characteristics numerically.In this method,a periodic pressure oscillation signal which frequency equals to the first acoustic mode is superimposed on a steady flow at the head end of the chamber.When the pressure oscillation is turned off,the decay rate of the pressure can be used to determine the nozzle attenuation constant.The damping characteristics of three other nozzle geometries are numerically studied with this method under the same operating condition.The results show that the convex nozzle provides more damping than the conical nozzle which in turn provides more damping than the concave nozzle.All the three nozzles have better damping effect than that of basic nozzle geometry.At last,the phase difference in the chamber is analyzed,and the numerical pressure distribution satisfies well with theoretical distribution.展开更多
Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-...Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-facility.In spite of the known medium and large scale facilities all over the world,we developed a laboratory scale instrument for component fraction measurements.It has a two-phase flow homogenizer loop with the clamp-on potential of the meters to provide a regime independent measurement.It is capable of delivering a complete homogenization by γ-ray densitometer.With an error of±5%in component fraction measurements,this instrument is appropriate for testing and calibrating other meters.展开更多
Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteri...Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteristics and attenuation coefficient.By designing an appropriate vector filter,phase velocity,attenuation coefficient and amplitude can be inverted from the waveform recorded by the receiver array.Performance analysis of this algorithm is compared with Extended Prony Method(EPM)and Forward and Backward Matrix Pencil(FBMP)method.Based on the analysis results,the proposed method is capable of achieving high resolution and precision as the parametric spectrum estimation methods.At the meantime,it also keeps high stability as the other nonparametric spectrum estimation methods.At last,applications to synthetic waveforms modeled using finite difference method and real data show its efficiency.The real data processing results show that the P-wave attenuation log is more sensitive to oil formation compared to S-wave;and the S-wave attenuation log is more sensitive to shale formation compared to P-wave.展开更多
A two-component lattice Boltzmann method(LBM) with a multiple-relaxation-time(MRT) collision operator is presented to improve the numerical stability of the single relaxation time(SRT) model. The macroscopic and the m...A two-component lattice Boltzmann method(LBM) with a multiple-relaxation-time(MRT) collision operator is presented to improve the numerical stability of the single relaxation time(SRT) model. The macroscopic and the momentum conservation equations can be retrieved through the Chapman–Enskog(C-E) expansion analysis. The equilibrium moment with the diffusion term is calculated, a diffusion phenomenon is simulated by utilizing the developed model, and the numerical stability is verified. Furthermore, the binary mixture channel model is designed to simulate the sound attenuation phenomenon, and the obtained simulation results are found to be consistent with the analytical solutions. The sound attenuation model is used to study the numerical stability and calculation accuracy of the LBM model. The simulation results show the stability and accuracy of the MRT model and the SRT model under different viscosity conditions. Finally,we study the influence of the error between the macroscopic equation of the MRT model and the standard incompressible Navier–Stokes equation on the calculation accuracy of the model to demonstrate the general applicability of the conclusions drawn by the sound attenuation model in the present study.展开更多
In this paper full-energy peak (photopeak) efficiency and photopeak attenuation coefficient of 3'' × 3'' NaI(Tl) well-type scintillation detector were calculated using gamma-rayisotropic radiating...In this paper full-energy peak (photopeak) efficiency and photopeak attenuation coefficient of 3'' × 3'' NaI(Tl) well-type scintillation detector were calculated using gamma-rayisotropic radiating point sources (with photon energy: 0.245, 0.344, 0.662, 0.779, 0.964, 1.1732, 1.333 and 1.408 MeV) placed outside the detector well. These energies were obtained from <sup>152</sup>Eu, <sup>137</sup>Cs and <sup>60</sup>Co. The relations between the full energy peak efficiency and photopeak attenuation coefficients, were plotted vs. photon energy at different sources to detector distance, and it found that the full energy peak efficiency decreased by increasing the distance between the source and the detector.展开更多
Seismic attenuation has been inherent media characteristics in which an interesting topic of research, for it reflects the seismic waves propagate. There are many factors that cause seismic wave attenuation, such as g...Seismic attenuation has been inherent media characteristics in which an interesting topic of research, for it reflects the seismic waves propagate. There are many factors that cause seismic wave attenuation, such as geometry attenuation caused by energy dissipating during propagation, friction attenuation by relative sliding among rock grains, and scattering attenuation by rock heterogeneity. In this paper we study P-wave scattering attenuation in a random elastic medium by numerical simulations from a statistical point of view. A random elastic medium model is built based on general stochastic process theory. Then a staggered-grid pseudo-spectral method is used to simulate wave propagation. Scattering attenuation is estimated by the spectral ratio method based on virtual detector records. Random elastic media numerical scatter results with various heterogeneity levels show that the higher heterogeneous levels cause greater scattering attenuation. When the scatter sizes are smaller than a wave length, the larger scatters give a greater attenuation. Finally, we propose a method to evaluate fluid-flow attenuation in porous media. The fluid- flow attenuation is derived from total attenuation and scattering attenuation in random porous media and the attenuation is estimated quantitatively. Results show that in the real seismic frequency range when the heterogeneous scale is about 10^1 meters （less than one wave length）, scattering attenuation is larger than fluid-tlow attenuation in random porous media and scattering attenuation is the main factor of seismic attenuation in real heterogeneous porous media.展开更多
Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which conside...Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which considers nonlinear interactions between waves and the arc-plate breakwater, has been constructed by using the velocity wave- generating method, the volume of fluid (VOF) method and the finite volume method. The results show that the relative width, relative height and relative submergence of the breakwater are three main influencing factors and have significant influence on wave energy dissipation of the lower arc-plate open breakwater. The transmission coefficient is found to decrease with the increasing relative width, and the minimum transmission coefficient is 0.15 when the relative width is 0.45. The reflection coefficient is found to vary slightly with the relative width, and the maximum reflection coefficient is 0.53 when the relative width is 0.45. The transmission and reflection coefficients are shown to increase with the relative wave height for approximately 85% of the experimental tests when the relative width is 0.19 0.45. The transmission coefficients at relative submergences of 0.04, 0.02 and 0 are clearly shown to be greater than those at relative submergences of 0.02 and 0.04, while the reflection coefficient exhibits the opposite relationship. After the wave interacts with the lower arc-plate breakwater, the wave energy is mainly converted into transmission, reflection and dissipation energies. The wave attenuation performance is clearly weakened for waves with greater heights and longer periods.展开更多
On the condition that the velocity ratio of compressive wave to shear wave is stable during the propagating process of seismic wave, this paper develops an estimation method for medium quality factor based on single s...On the condition that the velocity ratio of compressive wave to shear wave is stable during the propagating process of seismic wave, this paper develops an estimation method for medium quality factor based on single station＇s travel time difference between direct S and P waves and the first period signal of direct P wave. 8 774 high SNR wave data altogether recorded by 75 stations are analyzed. The results show that： （1） under the normal regional stress field, the quality factors QmsP in the stable tectonic area are higher than that in the active tectonic area around the Weifang-Jiashan section and its surrounding area in the Tanlu fault zone; （2） in the Juxian-Tancheng section the seismic wave attenuation is relatively quick, and the media is relatively broken, suggesting no tectonic stress accumulation; （3） the Xinyi-Sihong section is currently locking and in accumulating elastic strain energy stage, which has the deep environment and conditions of strong earthquake generation similar to those of the 1668 M8.5 Tancheng strong earthquake.展开更多
The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which...The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which provide important information for reservoir identification. This paper derives P-wave attenuation anisotropy in the ATI media where the symmetry axis is in the arbitrary direction theoretically and modifies the spectral ratio method to measure attenuation anisotropy in the ATI media, thus avoiding a large measurement error when applied to wide azimuth or full azimuth data. Fracture dip and azimuth can be estimated through attenuation anisotropy analysis. For small-scale fractures, fracture scale and fracture density can be determined with enhanced convergence if velocity and attenuation information are both used. We also apply the modified spectralratio method to microseismic field data from an oilfield in East China and extract the fracture dip through attenuation anisotropy analysis. The result agrees with the microseismie monitoring.展开更多
Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(...Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(RTM).We can easily implement attenuation-compensated RTM using the constant Q viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion terms.However,the nonphysical amplitude-compensation process will inevitably amplify the high-frequency noise in the wavefield in an exponential form,causing the numerical simulation to become unstable.This is due to the fact that the amplitude of the compensation grows exponentially with frequency.In order to achieve stable attenuation-compensated RTM,we modify the analytic expression of the attenuation compensation extrapolation operator and make it only compensate for amplitude loss within the effective frequency band.Based on this modified analytic formula,we then derive an explicit time-space domain attenuation compensation extrapolation operator.Finally,the implementation procedure of stable attenuation-compensated RTM is presented.In addition to being simple to implement,the newly proposed attenuation-compensated extrapolation operator is superior to the conventional low-pass filter in suppressing random noise,which will further improve the imaging resolution.We use two synthetic and one land seismic datasets to verify the stability and effectiveness of the proposed attenuationcompensated RTM in improving imaging resolution in viscous media.展开更多
The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the res...The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the resolution of reverse time migration(RTM).As an effective high-resolution imaging method,attenuation-compensated RTM(ACRTM)can eff ectively compensate for the energy loss caused by the attenuation related to media absorption under the influence of resistivity.Therefore,constructing an accurate resistivity-media model to compensate for the attenuation of electromagnetic wave energy is crucial for realizing the ACRTM imaging of GPR data.This study proposes a resistivity-constrained ACRTM imaging method for the imaging of GPR data by adding high-density resistivity detection along the GPR survey line and combining it with its resistivity inversion profile.The proposed method uses the inversion result of apparent resistivity data as the GPR RTM-resistivity model for imposing resistivity constraints.Moreover,the hybrid method involving image minimum entropy and RTM is used to estimate the medium velocity at the diff raction position,and combined with the distribution characteristics of the reflection in the GPR profile,a highly accurate velocity model is built to improve the imaging resolution of the ACRTM.The accuracy and eff ectiveness of the proposed method are verified using the ACRTM test of the GPR simulated data of a typical attenuating media model.On this basis,the GPR and apparent resistivity data were observed on a field survey line,and use the GPR resistivity-constrained ACRTM method to image the observed data.A comparison of the proposed method with the conventional ACRTM method shows that the proposed method has better imaging depth,stronger energy,and higher resolution,and the obtained results are more conducive for subsequent data analysis and interpretation.展开更多
We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yi...We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yimen, and Luquan, were selected and their attenuations of P and S waves were measured by using the extended coda-normalization method. The six stations were classified into three regions according to their location, that is, Baoshan area, Yongsheng-Lijiang-Heqing (YLH) area and Luquan-Yimen (LY) area. The values of QS-1 and QP-1 are expressed as QS-1=0.00867 f-0.86, QP-1=0.01155 f-0.93, QS-1=0.01824 f-0.92, QP-1=0.02288 f-0.92, and QS-1=0.01647 f-0.91, QP-1=0.02826 f-0.97 in Baoshan, YLH, LY areas respectively. The attenuation of YLH, LY are closer to each other, however, Baoshan area is apparently lower. Comparing attenuation in the three areas with other areas of the world using the same method, it is suggested that the attenuations of P and S waves in YLH and LY areas are close to Kanto of Japan, but much higher than southeast of South Korea. The QS-1 and QP-1 in Baoshan area are slightly higher than southeastern South Korea. Furthermore, the results indicate that our QS-1 in Yunnan area is close to others by analyzing the coda attenuation.展开更多
[ Objective] The research aimed to study application of the attenuation correction technology in C-band radar precipitation estimation. [ Method~ Based on CINRAD-CB radar data in Shaanxi, we conducted the attenuation ...[ Objective] The research aimed to study application of the attenuation correction technology in C-band radar precipitation estimation. [ Method~ Based on CINRAD-CB radar data in Shaanxi, we conducted the attenuation correction experiment by using iteration method and Kufa method respectively. Moreover, we conducted application expedment of the Kufa attenuation correction method in the quantitative precipitation esti- mation. [ Result~ Attenuation correction technology could compensate for attenuation problem of the echo at the distant range. Calculation result of the iteration method finally tended to that of the Kufa method. Moreover, iteration method spent more time. Therefore, Kufa attenuation correction technology was more suitable for business operation. When strong echo was near radar, generated attenuation was more obvious, and application value of the attenuation correction was bigger. Attenuation correction technology was used for quantitative precipitation estimation, which was favor- able for improving accuracy of the precipitation estimation. But we should conduct detailed planning on calculation scheme of the precipitation esti- mation because that different calculation schemes had great influences on accuracy of the quantitative precipitation estimation. [ Cendusien] This research provided a basis for improving accuracy of the quantitative precipitation estimation in Shaanxi. Key words Attenuation correction展开更多
The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor(GECAM)mission is designed to monitor the Gamma-Ray Bursts(GRBs)associated with gravitational waves and other high-energy transient sources...The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor(GECAM)mission is designed to monitor the Gamma-Ray Bursts(GRBs)associated with gravitational waves and other high-energy transient sources.The mission consists of two microsatellites which are planned to operate at the opposite sides of the Earth.Each GECAM satellite could detect and localize GRBs in about 8 keV-5 MeV with its 25 Gamma-Ray Detectors(GRDs).In this work,we report the in-flight energy calibration of GRDs using the characteristic gamma-ray lines in the background spectra,and show their performance evolution during the commissioning phase.Besides,a preliminary cross-calibration of energy response with Fermi GBM data is also presented,validating the energy response of GRDs.展开更多
Accurate acquisition and prediction of acoustic parameters of seabed sediments are crucial in marine sound propagation research.While the relationship between sound velocity and physical properties of sediment has bee...Accurate acquisition and prediction of acoustic parameters of seabed sediments are crucial in marine sound propagation research.While the relationship between sound velocity and physical properties of sediment has been extensively studied,there is still no consensus on the correlation between acoustic attenuation coefficient and sediment physical properties.Predicting the acoustic attenuation coefficient remains a challenging issue in sedimentary acoustic research.In this study,we propose a prediction method for the acoustic attenuation coefficient using machine learning algorithms,specifically the random forest(RF),support vector machine(SVR),and convolutional neural network(CNN)algorithms.We utilized the acoustic attenuation coefficient and sediment particle size data from 52 stations as training parameters,with the particle size parameters as the input feature matrix,and measured acoustic attenuation as the training label to validate the attenuation prediction model.Our results indicate that the error of the attenuation prediction model is small.Among the three models,the RF model exhibited the lowest prediction error,with a mean squared error of 0.8232,mean absolute error of 0.6613,and root mean squared error of 0.9073.Additionally,when we applied the models to predict the data collected at different times in the same region,we found that the models developed in this study also demonstrated a certain level of reliability in real prediction scenarios.Our approach demonstrates that constructing a sediment acoustic characteristics model based on machine learning is feasible to a certain extent and offers a novel perspective for studying sediment acoustic properties.展开更多
基金supported by the National Key R&D Program of China(2021YFA0718500)support from the Strategic Priority Research Program on Space Science,the Chinese Academy of Sciences(grant Nos.XDA15360102,XDA15360300,XDA15052700 and E02212A02S)+1 种基金the National Natural Science Foundation of China(grant Nos.12173038 and U2038106)the National HEP Data Center(grant No.E029S2S1)。
文摘Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations.Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors.There are two major methods for this count distribution localization:χ^(2)minimization method and the Bayesian method.Here we propose a modified Bayesian method that could take advantage of both the accuracy of the Bayesian method and the simplicity of the χ^(2)method.With comprehensive simulations,we find that our Bayesian method with Poisson likelihood is generally more applicable for various bursts than the χ^(2)method,especially for weak bursts.We further proposed a location-spectrum iteration approach based on the Bayesian inference,which could alleviate the problems caused by the spectral difference between the burst and location templates.Our method is very suitable for scenarios with limited computation resources or timesensitive applications,such as in-flight localization software,and low-latency localization for rapidly follow-up observations.
基金the National Natural Science Foundation of China under Grant No.10474016.
文摘In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effects of three-dimensional gas flow and acoustic damping on the acoustic attenuation characteristics of marine engine exhaust silencers, a dual reciprocity boundary element method (DRBEM) was developed. The acoustic governing equation in three-dimensional potential flow was derived first, and then the DRBEM numerical procedure is given. Compared to the conventional boundary element method (CBEM), the DRBEM considers the second order terms of flow Mach number in the acoustic governing equation, so it is suitable for the cases with higher Mach number subsonic flow. For complex exhaust silencers, it is difficult to apply the single-domain boundary element method, so a substructure approach based on the dual reciprocity boundary element method is presented. The experiments for measuring transmission loss of silencers are conducted, and the experimental setup and measurements are explained. The transmission loss of a single expansion chamber silencer with extended inlet and outlet were predicted by DRBEM and compared with the measurements. The good agreements between predictions and measurements are observed, which demonstrated that the derived acoustic governing equation and the DRBEM numerical procedure in the present study are correct.
文摘Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burning,a wave attenuation method is proposed to assess the nozzle damping characteristics numerically.In this method,a periodic pressure oscillation signal which frequency equals to the first acoustic mode is superimposed on a steady flow at the head end of the chamber.When the pressure oscillation is turned off,the decay rate of the pressure can be used to determine the nozzle attenuation constant.The damping characteristics of three other nozzle geometries are numerically studied with this method under the same operating condition.The results show that the convex nozzle provides more damping than the conical nozzle which in turn provides more damping than the concave nozzle.All the three nozzles have better damping effect than that of basic nozzle geometry.At last,the phase difference in the chamber is analyzed,and the numerical pressure distribution satisfies well with theoretical distribution.
基金financially supported by Radiation Application Research School as a part of Nuclear Science and Technology Research Institute
文摘Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-facility.In spite of the known medium and large scale facilities all over the world,we developed a laboratory scale instrument for component fraction measurements.It has a two-phase flow homogenizer loop with the clamp-on potential of the meters to provide a regime independent measurement.It is capable of delivering a complete homogenization by γ-ray densitometer.With an error of±5%in component fraction measurements,this instrument is appropriate for testing and calibrating other meters.
基金This research was supported by the National Natural Science Foundation of China(No.42274141)Science Foundation of China University of Petroleum,Beijing(No.2462020YXZZ007).
文摘Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteristics and attenuation coefficient.By designing an appropriate vector filter,phase velocity,attenuation coefficient and amplitude can be inverted from the waveform recorded by the receiver array.Performance analysis of this algorithm is compared with Extended Prony Method(EPM)and Forward and Backward Matrix Pencil(FBMP)method.Based on the analysis results,the proposed method is capable of achieving high resolution and precision as the parametric spectrum estimation methods.At the meantime,it also keeps high stability as the other nonparametric spectrum estimation methods.At last,applications to synthetic waveforms modeled using finite difference method and real data show its efficiency.The real data processing results show that the P-wave attenuation log is more sensitive to oil formation compared to S-wave;and the S-wave attenuation log is more sensitive to shale formation compared to P-wave.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 12174085, 11874140, and 11574072)the State Key Laboratory of Acoustics, Chinese Academy of Sciences (Grant No. SKLA201913)the Postgraduate Research and Practice Innovation Program of Jiangsu Province, China (Grant No. KYCX21 0478)。
文摘A two-component lattice Boltzmann method(LBM) with a multiple-relaxation-time(MRT) collision operator is presented to improve the numerical stability of the single relaxation time(SRT) model. The macroscopic and the momentum conservation equations can be retrieved through the Chapman–Enskog(C-E) expansion analysis. The equilibrium moment with the diffusion term is calculated, a diffusion phenomenon is simulated by utilizing the developed model, and the numerical stability is verified. Furthermore, the binary mixture channel model is designed to simulate the sound attenuation phenomenon, and the obtained simulation results are found to be consistent with the analytical solutions. The sound attenuation model is used to study the numerical stability and calculation accuracy of the LBM model. The simulation results show the stability and accuracy of the MRT model and the SRT model under different viscosity conditions. Finally,we study the influence of the error between the macroscopic equation of the MRT model and the standard incompressible Navier–Stokes equation on the calculation accuracy of the model to demonstrate the general applicability of the conclusions drawn by the sound attenuation model in the present study.
文摘In this paper full-energy peak (photopeak) efficiency and photopeak attenuation coefficient of 3'' × 3'' NaI(Tl) well-type scintillation detector were calculated using gamma-rayisotropic radiating point sources (with photon energy: 0.245, 0.344, 0.662, 0.779, 0.964, 1.1732, 1.333 and 1.408 MeV) placed outside the detector well. These energies were obtained from <sup>152</sup>Eu, <sup>137</sup>Cs and <sup>60</sup>Co. The relations between the full energy peak efficiency and photopeak attenuation coefficients, were plotted vs. photon energy at different sources to detector distance, and it found that the full energy peak efficiency decreased by increasing the distance between the source and the detector.
基金supported by the State Key Program of National Natural Science of China (Grant No. 40839901)
文摘Seismic attenuation has been inherent media characteristics in which an interesting topic of research, for it reflects the seismic waves propagate. There are many factors that cause seismic wave attenuation, such as geometry attenuation caused by energy dissipating during propagation, friction attenuation by relative sliding among rock grains, and scattering attenuation by rock heterogeneity. In this paper we study P-wave scattering attenuation in a random elastic medium by numerical simulations from a statistical point of view. A random elastic medium model is built based on general stochastic process theory. Then a staggered-grid pseudo-spectral method is used to simulate wave propagation. Scattering attenuation is estimated by the spectral ratio method based on virtual detector records. Random elastic media numerical scatter results with various heterogeneity levels show that the higher heterogeneous levels cause greater scattering attenuation. When the scatter sizes are smaller than a wave length, the larger scatters give a greater attenuation. Finally, we propose a method to evaluate fluid-flow attenuation in porous media. The fluid- flow attenuation is derived from total attenuation and scattering attenuation in random porous media and the attenuation is estimated quantitatively. Results show that in the real seismic frequency range when the heterogeneous scale is about 10^1 meters （less than one wave length）, scattering attenuation is larger than fluid-tlow attenuation in random porous media and scattering attenuation is the main factor of seismic attenuation in real heterogeneous porous media.
基金supported by the NSFC-Shandong Joint Fund(Grant Nos.U1706220 and U1806227)the National Natural Science Foundation of China(Grant Nos.51709140 and 51879019)the Key Laboratory of Coastal Disasters and Defence of Ministry of Education(Grant No.201703)
文摘Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which considers nonlinear interactions between waves and the arc-plate breakwater, has been constructed by using the velocity wave- generating method, the volume of fluid (VOF) method and the finite volume method. The results show that the relative width, relative height and relative submergence of the breakwater are three main influencing factors and have significant influence on wave energy dissipation of the lower arc-plate open breakwater. The transmission coefficient is found to decrease with the increasing relative width, and the minimum transmission coefficient is 0.15 when the relative width is 0.45. The reflection coefficient is found to vary slightly with the relative width, and the maximum reflection coefficient is 0.53 when the relative width is 0.45. The transmission and reflection coefficients are shown to increase with the relative wave height for approximately 85% of the experimental tests when the relative width is 0.19 0.45. The transmission coefficients at relative submergences of 0.04, 0.02 and 0 are clearly shown to be greater than those at relative submergences of 0.02 and 0.04, while the reflection coefficient exhibits the opposite relationship. After the wave interacts with the lower arc-plate breakwater, the wave energy is mainly converted into transmission, reflection and dissipation energies. The wave attenuation performance is clearly weakened for waves with greater heights and longer periods.
基金supported by National Natural Science Foundation of China(grant No.40974031)Shandong Provincial Natural Science Foundation(grant No.Y2007E09)Jiangsu Provincial Social Development Project(grant No.BE2009691)
文摘On the condition that the velocity ratio of compressive wave to shear wave is stable during the propagating process of seismic wave, this paper develops an estimation method for medium quality factor based on single station＇s travel time difference between direct S and P waves and the first period signal of direct P wave. 8 774 high SNR wave data altogether recorded by 75 stations are analyzed. The results show that： （1） under the normal regional stress field, the quality factors QmsP in the stable tectonic area are higher than that in the active tectonic area around the Weifang-Jiashan section and its surrounding area in the Tanlu fault zone; （2） in the Juxian-Tancheng section the seismic wave attenuation is relatively quick, and the media is relatively broken, suggesting no tectonic stress accumulation; （3） the Xinyi-Sihong section is currently locking and in accumulating elastic strain energy stage, which has the deep environment and conditions of strong earthquake generation similar to those of the 1668 M8.5 Tancheng strong earthquake.
基金supported by 973 Program of China(No.2013CB228602)National Science and Technology Major Project of China(No.2016ZX05004003-002)863 Program of China(No.2013AA064202)
文摘The existence of aligned fractures in fluid-saturated rocks leads to obvious attenuation anisotropy and velocity anisotropy. Attenuation anisotropy analysis can be applied to estimate fracture density and scale, which provide important information for reservoir identification. This paper derives P-wave attenuation anisotropy in the ATI media where the symmetry axis is in the arbitrary direction theoretically and modifies the spectral ratio method to measure attenuation anisotropy in the ATI media, thus avoiding a large measurement error when applied to wide azimuth or full azimuth data. Fracture dip and azimuth can be estimated through attenuation anisotropy analysis. For small-scale fractures, fracture scale and fracture density can be determined with enhanced convergence if velocity and attenuation information are both used. We also apply the modified spectralratio method to microseismic field data from an oilfield in East China and extract the fracture dip through attenuation anisotropy analysis. The result agrees with the microseismie monitoring.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003).
文摘Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(RTM).We can easily implement attenuation-compensated RTM using the constant Q viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion terms.However,the nonphysical amplitude-compensation process will inevitably amplify the high-frequency noise in the wavefield in an exponential form,causing the numerical simulation to become unstable.This is due to the fact that the amplitude of the compensation grows exponentially with frequency.In order to achieve stable attenuation-compensated RTM,we modify the analytic expression of the attenuation compensation extrapolation operator and make it only compensate for amplitude loss within the effective frequency band.Based on this modified analytic formula,we then derive an explicit time-space domain attenuation compensation extrapolation operator.Finally,the implementation procedure of stable attenuation-compensated RTM is presented.In addition to being simple to implement,the newly proposed attenuation-compensated extrapolation operator is superior to the conventional low-pass filter in suppressing random noise,which will further improve the imaging resolution.We use two synthetic and one land seismic datasets to verify the stability and effectiveness of the proposed attenuationcompensated RTM in improving imaging resolution in viscous media.
基金supported by the National Natural Science Foundation of China (No.41604102)the Guangxi Natural Science Foundation project (No.2020GXNSFAA159121).
文摘The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the resolution of reverse time migration(RTM).As an effective high-resolution imaging method,attenuation-compensated RTM(ACRTM)can eff ectively compensate for the energy loss caused by the attenuation related to media absorption under the influence of resistivity.Therefore,constructing an accurate resistivity-media model to compensate for the attenuation of electromagnetic wave energy is crucial for realizing the ACRTM imaging of GPR data.This study proposes a resistivity-constrained ACRTM imaging method for the imaging of GPR data by adding high-density resistivity detection along the GPR survey line and combining it with its resistivity inversion profile.The proposed method uses the inversion result of apparent resistivity data as the GPR RTM-resistivity model for imposing resistivity constraints.Moreover,the hybrid method involving image minimum entropy and RTM is used to estimate the medium velocity at the diff raction position,and combined with the distribution characteristics of the reflection in the GPR profile,a highly accurate velocity model is built to improve the imaging resolution of the ACRTM.The accuracy and eff ectiveness of the proposed method are verified using the ACRTM test of the GPR simulated data of a typical attenuating media model.On this basis,the GPR and apparent resistivity data were observed on a field survey line,and use the GPR resistivity-constrained ACRTM method to image the observed data.A comparison of the proposed method with the conventional ACRTM method shows that the proposed method has better imaging depth,stronger energy,and higher resolution,and the obtained results are more conducive for subsequent data analysis and interpretation.
基金The specialized fund for pre-study of national key basic study (2002CCD01700)national key project of science and technology from China Ministry of Science and Technology during the 10th Five-year Plan (2004BA601B01-04-03)
文摘We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yimen, and Luquan, were selected and their attenuations of P and S waves were measured by using the extended coda-normalization method. The six stations were classified into three regions according to their location, that is, Baoshan area, Yongsheng-Lijiang-Heqing (YLH) area and Luquan-Yimen (LY) area. The values of QS-1 and QP-1 are expressed as QS-1=0.00867 f-0.86, QP-1=0.01155 f-0.93, QS-1=0.01824 f-0.92, QP-1=0.02288 f-0.92, and QS-1=0.01647 f-0.91, QP-1=0.02826 f-0.97 in Baoshan, YLH, LY areas respectively. The attenuation of YLH, LY are closer to each other, however, Baoshan area is apparently lower. Comparing attenuation in the three areas with other areas of the world using the same method, it is suggested that the attenuations of P and S waves in YLH and LY areas are close to Kanto of Japan, but much higher than southeast of South Korea. The QS-1 and QP-1 in Baoshan area are slightly higher than southeastern South Korea. Furthermore, the results indicate that our QS-1 in Yunnan area is close to others by analyzing the coda attenuation.
文摘[ Objective] The research aimed to study application of the attenuation correction technology in C-band radar precipitation estimation. [ Method~ Based on CINRAD-CB radar data in Shaanxi, we conducted the attenuation correction experiment by using iteration method and Kufa method respectively. Moreover, we conducted application expedment of the Kufa attenuation correction method in the quantitative precipitation esti- mation. [ Result~ Attenuation correction technology could compensate for attenuation problem of the echo at the distant range. Calculation result of the iteration method finally tended to that of the Kufa method. Moreover, iteration method spent more time. Therefore, Kufa attenuation correction technology was more suitable for business operation. When strong echo was near radar, generated attenuation was more obvious, and application value of the attenuation correction was bigger. Attenuation correction technology was used for quantitative precipitation estimation, which was favor- able for improving accuracy of the precipitation estimation. But we should conduct detailed planning on calculation scheme of the precipitation esti- mation because that different calculation schemes had great influences on accuracy of the quantitative precipitation estimation. [ Cendusien] This research provided a basis for improving accuracy of the quantitative precipitation estimation in Shaanxi. Key words Attenuation correction
基金supported by the Strategic Priority Program on Space Science,the Chinese Academy of Sciences,grant Nos.XDA15360102,XDA15360300 and E02212A02Sthe National Natural Science Foundation of China(Project:12061131007)。
文摘The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor(GECAM)mission is designed to monitor the Gamma-Ray Bursts(GRBs)associated with gravitational waves and other high-energy transient sources.The mission consists of two microsatellites which are planned to operate at the opposite sides of the Earth.Each GECAM satellite could detect and localize GRBs in about 8 keV-5 MeV with its 25 Gamma-Ray Detectors(GRDs).In this work,we report the in-flight energy calibration of GRDs using the characteristic gamma-ray lines in the background spectra,and show their performance evolution during the commissioning phase.Besides,a preliminary cross-calibration of energy response with Fermi GBM data is also presented,validating the energy response of GRDs.
基金funded by the Basic Scientific Fund for National Public Research Institutes of China(No.2022 S01)the National Natural Science Foundation of China(Nos.42176191,42049902,and U22A2012)+5 种基金the Shandong Provincial Natural Science Foundation,China(No.ZR2022YQ40)the National Key R&D Program of China(No.2021YFF0501202)the Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(No.SML2023 SP232)the Fundamental Research Funds for the Central Universities,Sun Yat-sen University(No.241gqb006)Data acquisition and sample collections were supported by the National Natural Science Foundation of China Open Research Cruise(Cruise No.NORC2021-02+NORC2021301)funded by the Shiptime Sharing Project of the National Natural Science Foundation of China。
文摘Accurate acquisition and prediction of acoustic parameters of seabed sediments are crucial in marine sound propagation research.While the relationship between sound velocity and physical properties of sediment has been extensively studied,there is still no consensus on the correlation between acoustic attenuation coefficient and sediment physical properties.Predicting the acoustic attenuation coefficient remains a challenging issue in sedimentary acoustic research.In this study,we propose a prediction method for the acoustic attenuation coefficient using machine learning algorithms,specifically the random forest(RF),support vector machine(SVR),and convolutional neural network(CNN)algorithms.We utilized the acoustic attenuation coefficient and sediment particle size data from 52 stations as training parameters,with the particle size parameters as the input feature matrix,and measured acoustic attenuation as the training label to validate the attenuation prediction model.Our results indicate that the error of the attenuation prediction model is small.Among the three models,the RF model exhibited the lowest prediction error,with a mean squared error of 0.8232,mean absolute error of 0.6613,and root mean squared error of 0.9073.Additionally,when we applied the models to predict the data collected at different times in the same region,we found that the models developed in this study also demonstrated a certain level of reliability in real prediction scenarios.Our approach demonstrates that constructing a sediment acoustic characteristics model based on machine learning is feasible to a certain extent and offers a novel perspective for studying sediment acoustic properties.
