Full waveform inversion of SH‐ and Love‐wave data in near‐surface prospecting

We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls.     

天津港“8·12”爆炸事件记录波形分析

利用天津地震台网23个短周期地震仪记录到的2015年8月12日天津港爆炸事件波形数据,对事件进行初步分析和研究。分别使用能量包络线法和波形匹配法,共识别出4次爆炸事件,并结合波形互相关的主事件定位法对爆炸事件进行定位。分析结果表明:第一次爆炸发生的时间为2015年8月12日23时34分4秒,1.77s后发生了规模更大的二号爆炸事件;34.11s后又发生了最强的三号爆炸事件;86.33s再一次发生了能量较弱的四号爆炸事件,其中二号与三号爆炸事件相距65~70m,二号与四号事件相距160~170m。此外,通过选用爆破、矿塌、核爆和天然地震等不同类型事件的波形数据,从时间域和频率域分别对这些事件进行对比分析,发现其震相特征表现出明显差异。     

2017年6月3日阿拉善左旗5.0级地震震源深度测定

2017年6月3日18时11分内蒙古阿拉善左旗（37.99° N,103.56° E）发生5.0级地震,利用近震深度震相sPL测定法、CAP波形反演法、sPn与Pn震相走时差法,进行震源深度精确测定。研究结果表明,3种方法测得的该地震震源深度基本一致,分布范围为8.5—11 km,表明该地震发生在上地壳。     

Improving gravity anomalies over China marginal sea from retracked Geosat and ERS-1 data

The quality of altimeter data and ocean tide model is critical to the recovery of coastal gravity anomalies. In this contribution, three retracking methods (threshold, improved threshold and Beta-5) are investigated with the aim of improving the altimeter data over a shallow water area. Comparison indicates that the improved threshold is the best retracking method over China Sea. Two ocean tide models, NAO99b and CSR4.0, are analyzed. Results show that different tide models used in the processing of altimeter data may result in differences more than 10 mGal in recovered coastal gravity anomalies. Also, NAO99b is more suitable than CSR4.0 over the shallow water area of China Sea. Finally, gravity anomalies over China Sea are calculated from retracked Geosat/GM and ERS-1/GM data by least squares collocation. Comparison with shipborne gravimetry data demonstrates that gravity anomalies from retracked data are significantly superior to those from non-retracked data. Our results have the same order as the other two altimeter-derived gravity models: Sandwell&Smith(V16) and DNSC08.     

基于全波形采样的激电多信息提取方法研究与应用

激电法是地质勘查中的一种重要方法,时、频域测量各具优势,但传统的时、频域测量实现方法相对独立.根据频域和时域可相互转换的理论,本文提出一种基于全波形采样的时频激电多参数提取的地球物理数据处理方法.采取时间域激电法的观测模式,接收机高精度同步记录整周期电压-电流全部采样点的波形数据,通过该文提出的处理方法,仅一次供电和测量即可提取多种时、频域激电参数:时域激电的视电阻率、视极化率和频域激电的多频视电阻率、视相位、视频散率、去耦后的视相位等参数.在云南保山某典型矿区开展方法的试验应用,提取了地质体大量的相关激电信息,总结了目标体的激电参数响应规律.应用结果表明,该方法便捷高效,抗干扰能力强,相较于传统的时域激电测量,多种参数组合可提升对目标地质体的认识,增强激电法的应用效果.     

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 MW6.4 PISHAN EARTHQUAKE

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.     

基于三维应变格林函数反演中小震震源机制

随着复杂速度结构反演的发展和高性能计算能力的提升,基于高精度3D介质模型计算格林函数反演震源机制更具可行性。中小地震因具有更好的覆盖和近似点源效应,在区域结构成像中有着广泛的应用。基于波形类的反演方法如波动方程层析成像\,全波形反演都需要震源机制解,而传统的震源机制反演方法不能很好地应用于中小地震。采用有限差分法构建应变格林张量(Strain Green Tensor,SGT)数据库,将合成波形和实际波形按震相截窗并滤波到不同的频带范围,先通过最小化互相关走时差来进行震源重定位,再通过最小化波形残差反演震源机制。通过合成数据测试验证方法的正确性,随后将该方法应用于青藏高原东部边缘龙门山地区,反演一系列M_W3.4～5.7的中小地震震源机制。由于应变格林张量数据库可预先构建,该方法可以应用于(近)实时震源机制解反演。     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.     

基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加,地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响,提高多参数反演的精度,而截断牛顿法是一种可以比较准确地估计Hessian逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明,在VTI介质声波双参数同时反演中,截断牛顿法比有限内存BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno,L-BFGS)法能更准确地估计Hessian逆算子,进而较好地平衡两类不同参数的同时更新,得到了比较精确的反演结果.