Active faults pose a great threat to urban security.As the largest NW-trend active fault in Beijing area,the Nankou-Sunhe fault plays an important role in earthquake disaster and city construction.In this study,we col...Active faults pose a great threat to urban security.As the largest NW-trend active fault in Beijing area,the Nankou-Sunhe fault plays an important role in earthquake disaster and city construction.In this study,we collect continuous ambient noise data recorded by 43 temporary short-period seismograph between September 21th to October 12th 2019 to investigate the near-surface structure beneath the Nankou-Sunhe fault by using ambient noise tomography(ANT)and horizontal-to-vertical spectral ratio(HVSR)method.From ambient noise processing,fundamental-mode Rayleigh wave signals are clearly observed in the frequency band of 0.4-2.5 Hz.Then direct surface-wave inversion algorithm is applied to calculate the 3D shear-wave velocity model.Our results show that there is a sharp velocity contrast across the Nankou-Sunhe fault,with low velocities down to about 2 km on the hanging wall and high velocity on the footwall of the fault.According to the geological investigation,the low velocities are related to thicker sediments and Jurassic volcanic rock below which are the cap rock of the hydrothermal system.From the HVSR analysis,the HVSR curves of the sites near the fault shows double-peak,one less than 1 Hz and the other centered 7 Hz.After converting frequency to depth by the empirical equation,the results show that the thickness of sediments is thinned from southwest to northeast,which generally agrees well to field survey.Our results provide high-resolution near-surface structure for future study on disaster risk reduction and urban planning.展开更多
The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a m...The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.展开更多
The quality factor Q, which reflects the energy attenuation of seismic waves in subsurface media, is a diagnostic tool for hydrocarbon detection and reservoir characterization. In this paper, we propose a new Q extrac...The quality factor Q, which reflects the energy attenuation of seismic waves in subsurface media, is a diagnostic tool for hydrocarbon detection and reservoir characterization. In this paper, we propose a new Q extraction method based on the energy ratio before and after the wavelet attenuation, named the energy-ratio method（ERM）. The proposed method uses multipoint signal data in the time domain to estimate the wavelet energy without invoking the source wavelet spectrum, which is necessary in conventional Q extraction methods, and is applicable to any source wavelet spectrum; however, it requires high-precision seismic data. Forward zero-offset VSP modeling suggests that the ERM can be used for reliable Q inversion after nonintrinsic attenuation（geometric dispersion, refl ection, and transmission loss） compensation. The application to real zero-offset VSP data shows that the Q values extracted by the ERM and spectral ratio methods are identical, which proves the reliability of the new method.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
In this study,we considered the three-dimensional flow of a rotating viscous,incompressible electrically conducting nanofluid with oxytactic microorganisms and an insulated plate floating in the fluid.Three scenarios ...In this study,we considered the three-dimensional flow of a rotating viscous,incompressible electrically conducting nanofluid with oxytactic microorganisms and an insulated plate floating in the fluid.Three scenarios were considered in this study.The first case is when the fluid drags the plate,the second is when the plate drags the fluid and the third is when the plate floats on the fluid at the same velocity.The denser microorganisms create the bioconvection as they swim to the top following an oxygen gradient within the fluid.The velocity ratio parameter plays a key role in the dynamics for this flow.Varying the parameter below and above a critical value alters the dynamics of the flow.The Hartmann number,buoyancy ratio and radiation parameter have a reverse effect on the secondary velocity for values of the velocity ratio above and below the critical value.The Hall parameter on the other hand has a reverse effect on the primary velocity for values of velocity ratio above and below the critical value.The bioconvection Rayleigh number decreases the primary velocity.The secondary velocity increases with increasing values of the bioconvection Rayleigh number and is positive for velocity ratio values below 0.5.For values of the velocity ratio parameter above 0.5,the secondary velocity is negative for small values of bioconvection Rayleigh number and as the values increase,the flow is reversed and becomes positive.展开更多
We present new quantitative model describing the pressure dependence of acoustic P-and S-wave velocities.Assuming that a variety of individual mechanisms or defects(such as cracks,pore collapse and grain crushing)can ...We present new quantitative model describing the pressure dependence of acoustic P-and S-wave velocities.Assuming that a variety of individual mechanisms or defects(such as cracks,pore collapse and grain crushing)can contribute to the pressure-dependent change of the wave velocity,we order a characteristic pressure to all of them and allow a series of exponential terms in the description of the(Pand S-waves)velocity-pressure function.We estimate the parameters of the multi-exponential rock physical model in inversion procedures using laboratory measured P-and S-wave velocity data.As is known,the conventional damped least squares method gives acceptable results only when one or two individual mechanisms are assumed.Increasing the number of exponential terms leads to highly nonlinear ill-posed inverse problem.Due to this reason,we develop the spectral inversion method(SIM)in which the velocity amplitudes(the spectral lines in the characteristic pressure spectrum)are only considered as unknowns.The characteristic pressures(belonging to the velocity amplitudes)are excluded from the set of inversion unknowns,instead,they are defined in a set of fixed positions equidistantly distributed in the actual interval of the independent variable(pressure).Through this novel linear inversion method,we estimate the parameters of the multi-exponential rock physical model using laboratory measured P-and S-wave velocity data.The characteristic pressures are related to the closing pressures of cracks which are described by well-known rock mechanical relationships depending on the aspect ratio of elliptical cracks.This gives the possibility to estimate the aspect ratios in terms of the characteristic pressures.展开更多
震害现场调查和研究成果表明,场地土层结构对地震动有较大影响。为了研究场地内土层结构的放大效应,以河北省地震灾害防御与风险评价重点实验室的场地与结构观测台阵(Seismic Monitoring Array of Site and Structure,SMASS)所记录的古...震害现场调查和研究成果表明,场地土层结构对地震动有较大影响。为了研究场地内土层结构的放大效应,以河北省地震灾害防御与风险评价重点实验室的场地与结构观测台阵(Seismic Monitoring Array of Site and Structure,SMASS)所记录的古冶和滦州地震记录为基础,采用传统谱比法,将地下236 m测点处设为参考场地,把不同监测点处观测到的地震记录傅氏谱同参考场地进行处理及对比分析,得到以下结论:(1)SMASS台阵场地中浅部土层对地震动的放大效应更为明显。(2)SMASS台阵场地中土层结构对地震动的放大效应存在明显差异,在深部土层中,土层结构对水平向地震动的放大效应更为明显;在浅层土层中,土层结构对竖直向地震动的放大效应较为明显。(3)SMASS台阵场地中随着监测点所处深度不断上升,低频成分频带逐渐变宽,高频成分受到的影响较低频成分更为显著。展开更多
The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer.We present a new geometry that can be used to acquire field data by...The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer.We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation.In this study,we drilled number of receiver holes around the source hole,each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys.We also propose a novel tomographic inversion of the Q factor without the effect of the source signature,and examine its stability and reliability using synthetic data.We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield.The results show that seismic absorption in the nearsurface layer is much greater than that in the subsurface strata.Thus,it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption.In addition,we derive different Q factors from two frequency bands,which can be treated,to some extent,as evidence of a frequency-dependent Q.展开更多
In order to obtain stable interval Q factor, by analyzing the spectrum of monitoring wavelet and down-going wavelet of zero-offset VSP data and referring the spectrum expression of Ricker wavelet, we propose a new exp...In order to obtain stable interval Q factor, by analyzing the spectrum of monitoring wavelet and down-going wavelet of zero-offset VSP data and referring the spectrum expression of Ricker wavelet, we propose a new expression of source wavelet spectrum. Basing on the new expression, we present improved amplitude spectral fitting and spectral ratio methods for interval Q inversion based on zero-offset VSP data, and the sequence for processing the zero-offset VSP data. Subsequently, we apply the proposed methods to real zero-offset VSP data, and carry out prestack inverse Q filtering to zero-offset VSP data and surface seismic data for amplitude compensation with the estimated Q value.展开更多
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.展开更多
基金supported by Beijing Natural Science Foundation(No.8212041)National Natural Science Foundation of China(Nos.41874063 and U1939203)Shanghai Sheshan National Geophysical Observatory(No.2020K02).
文摘Active faults pose a great threat to urban security.As the largest NW-trend active fault in Beijing area,the Nankou-Sunhe fault plays an important role in earthquake disaster and city construction.In this study,we collect continuous ambient noise data recorded by 43 temporary short-period seismograph between September 21th to October 12th 2019 to investigate the near-surface structure beneath the Nankou-Sunhe fault by using ambient noise tomography(ANT)and horizontal-to-vertical spectral ratio(HVSR)method.From ambient noise processing,fundamental-mode Rayleigh wave signals are clearly observed in the frequency band of 0.4-2.5 Hz.Then direct surface-wave inversion algorithm is applied to calculate the 3D shear-wave velocity model.Our results show that there is a sharp velocity contrast across the Nankou-Sunhe fault,with low velocities down to about 2 km on the hanging wall and high velocity on the footwall of the fault.According to the geological investigation,the low velocities are related to thicker sediments and Jurassic volcanic rock below which are the cap rock of the hydrothermal system.From the HVSR analysis,the HVSR curves of the sites near the fault shows double-peak,one less than 1 Hz and the other centered 7 Hz.After converting frequency to depth by the empirical equation,the results show that the thickness of sediments is thinned from southwest to northeast,which generally agrees well to field survey.Our results provide high-resolution near-surface structure for future study on disaster risk reduction and urban planning.
基金supported by The National Key Research and Development Program Plane(No.2017YFC0601505)National Natural Science Foundation(No.41672325)Science&Technology Department of Sichuan Province Technology Project(No.2017GZ0393)
文摘The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.
基金supported by the Major Project of the Ministry of Science and Technology of China(No.2011ZX05024-001-01)National Nature Science Foundation of China(No.41140033)
文摘The quality factor Q, which reflects the energy attenuation of seismic waves in subsurface media, is a diagnostic tool for hydrocarbon detection and reservoir characterization. In this paper, we propose a new Q extraction method based on the energy ratio before and after the wavelet attenuation, named the energy-ratio method（ERM）. The proposed method uses multipoint signal data in the time domain to estimate the wavelet energy without invoking the source wavelet spectrum, which is necessary in conventional Q extraction methods, and is applicable to any source wavelet spectrum; however, it requires high-precision seismic data. Forward zero-offset VSP modeling suggests that the ERM can be used for reliable Q inversion after nonintrinsic attenuation（geometric dispersion, refl ection, and transmission loss） compensation. The application to real zero-offset VSP data shows that the Q values extracted by the ERM and spectral ratio methods are identical, which proves the reliability of the new method.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
文摘In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
文摘In this study,we considered the three-dimensional flow of a rotating viscous,incompressible electrically conducting nanofluid with oxytactic microorganisms and an insulated plate floating in the fluid.Three scenarios were considered in this study.The first case is when the fluid drags the plate,the second is when the plate drags the fluid and the third is when the plate floats on the fluid at the same velocity.The denser microorganisms create the bioconvection as they swim to the top following an oxygen gradient within the fluid.The velocity ratio parameter plays a key role in the dynamics for this flow.Varying the parameter below and above a critical value alters the dynamics of the flow.The Hartmann number,buoyancy ratio and radiation parameter have a reverse effect on the secondary velocity for values of the velocity ratio above and below the critical value.The Hall parameter on the other hand has a reverse effect on the primary velocity for values of velocity ratio above and below the critical value.The bioconvection Rayleigh number decreases the primary velocity.The secondary velocity increases with increasing values of the bioconvection Rayleigh number and is positive for velocity ratio values below 0.5.For values of the velocity ratio parameter above 0.5,the secondary velocity is negative for small values of bioconvection Rayleigh number and as the values increase,the flow is reversed and becomes positive.
基金supported by the European Union,co-financed by the European Social Fund and the GINOP-2.315-2016-00010"Development of enhanced engineering methods with the aim at utilization of subterranean energy resources"project in the framework of the Szechenyi 2020 Plan,funded by the European Union,co-financed by the European Structural and Investment Funds。
文摘We present new quantitative model describing the pressure dependence of acoustic P-and S-wave velocities.Assuming that a variety of individual mechanisms or defects(such as cracks,pore collapse and grain crushing)can contribute to the pressure-dependent change of the wave velocity,we order a characteristic pressure to all of them and allow a series of exponential terms in the description of the(Pand S-waves)velocity-pressure function.We estimate the parameters of the multi-exponential rock physical model in inversion procedures using laboratory measured P-and S-wave velocity data.As is known,the conventional damped least squares method gives acceptable results only when one or two individual mechanisms are assumed.Increasing the number of exponential terms leads to highly nonlinear ill-posed inverse problem.Due to this reason,we develop the spectral inversion method(SIM)in which the velocity amplitudes(the spectral lines in the characteristic pressure spectrum)are only considered as unknowns.The characteristic pressures(belonging to the velocity amplitudes)are excluded from the set of inversion unknowns,instead,they are defined in a set of fixed positions equidistantly distributed in the actual interval of the independent variable(pressure).Through this novel linear inversion method,we estimate the parameters of the multi-exponential rock physical model using laboratory measured P-and S-wave velocity data.The characteristic pressures are related to the closing pressures of cracks which are described by well-known rock mechanical relationships depending on the aspect ratio of elliptical cracks.This gives the possibility to estimate the aspect ratios in terms of the characteristic pressures.
文摘震害现场调查和研究成果表明,场地土层结构对地震动有较大影响。为了研究场地内土层结构的放大效应,以河北省地震灾害防御与风险评价重点实验室的场地与结构观测台阵(Seismic Monitoring Array of Site and Structure,SMASS)所记录的古冶和滦州地震记录为基础,采用传统谱比法,将地下236 m测点处设为参考场地,把不同监测点处观测到的地震记录傅氏谱同参考场地进行处理及对比分析,得到以下结论:(1)SMASS台阵场地中浅部土层对地震动的放大效应更为明显。(2)SMASS台阵场地中土层结构对地震动的放大效应存在明显差异,在深部土层中,土层结构对水平向地震动的放大效应更为明显;在浅层土层中,土层结构对竖直向地震动的放大效应较为明显。(3)SMASS台阵场地中随着监测点所处深度不断上升,低频成分频带逐渐变宽,高频成分受到的影响较低频成分更为显著。
基金supported by the National Natural Science Foundation of China(Grant No.41174117 and 41474109)the National Key Basic Research Development Program of China(Grant No.2013CB228606)
文摘The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer.We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation.In this study,we drilled number of receiver holes around the source hole,each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys.We also propose a novel tomographic inversion of the Q factor without the effect of the source signature,and examine its stability and reliability using synthetic data.We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield.The results show that seismic absorption in the nearsurface layer is much greater than that in the subsurface strata.Thus,it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption.In addition,we derive different Q factors from two frequency bands,which can be treated,to some extent,as evidence of a frequency-dependent Q.
基金sponsored by the National Nature Science Foundation of China(Nos.41174114 and 41274128)
文摘In order to obtain stable interval Q factor, by analyzing the spectrum of monitoring wavelet and down-going wavelet of zero-offset VSP data and referring the spectrum expression of Ricker wavelet, we propose a new expression of source wavelet spectrum. Basing on the new expression, we present improved amplitude spectral fitting and spectral ratio methods for interval Q inversion based on zero-offset VSP data, and the sequence for processing the zero-offset VSP data. Subsequently, we apply the proposed methods to real zero-offset VSP data, and carry out prestack inverse Q filtering to zero-offset VSP data and surface seismic data for amplitude compensation with the estimated Q value.
基金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.
