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.展开更多
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.展开更多
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.展开更多
在城市中应用微动H/V谱比方法面对大量且复杂的人文噪声干扰,需要对噪声强度较大的微动数据进行去噪处理或信号分析。本文针对现有方法难以处理干扰较大的微动数据以及信号提取过程繁琐的问题,提出基于XGBoost(extreme gradient boosti...在城市中应用微动H/V谱比方法面对大量且复杂的人文噪声干扰,需要对噪声强度较大的微动数据进行去噪处理或信号分析。本文针对现有方法难以处理干扰较大的微动数据以及信号提取过程繁琐的问题,提出基于XGBoost(extreme gradient boosting)的多重加权谱比降噪方法。首先对采集的微动数据进行幅值和频率分析,建立幅值加权谱比、频率加权谱比和多重加权谱比;然后根据建立的多重加权谱比,通过XGBoost方法获得降噪后的谱比曲线。将本文方法与传统STA/LTA(short time average/long time average)方法进行实际高噪声数据对比分析,结果表明相比于STA/LTA方法,本文方法对高噪声数据提取效果更好。展开更多
基金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.
基金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.
基金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.
文摘在城市中应用微动H/V谱比方法面对大量且复杂的人文噪声干扰,需要对噪声强度较大的微动数据进行去噪处理或信号分析。本文针对现有方法难以处理干扰较大的微动数据以及信号提取过程繁琐的问题,提出基于XGBoost(extreme gradient boosting)的多重加权谱比降噪方法。首先对采集的微动数据进行幅值和频率分析,建立幅值加权谱比、频率加权谱比和多重加权谱比;然后根据建立的多重加权谱比,通过XGBoost方法获得降噪后的谱比曲线。将本文方法与传统STA/LTA(short time average/long time average)方法进行实际高噪声数据对比分析,结果表明相比于STA/LTA方法,本文方法对高噪声数据提取效果更好。