A new energy density calculation method for reservoir delineation using seismic low-frequency signal

In comparison to high-frequency signals, low-frequency seismic signals suffer less from scattering and intrinsic attenuation during wave propagation, penetrate deeper strata and thus can provide more energy information related to the hydrocarbon reservoirs. Based on the asymptotic representation for the frequency-dependent reflections in the fluid-saturated pore-elastic media, we first derive a novel equation of the reservoir energy density and present an efficient workflow to calculate the reservoir energy density using low-frequency seismic data. Then, within a low-frequency range (from 1 to 30 Hz), we construct an objective function to determine the optimal frequency, using the energy densities calculated from the post-stack seismic traces close to the wells. Next, we can calculate the reservoir energy density using the instantaneous spectra of optimal frequency at the low-frequency end of the seismic spectrum. Tests on examples for synthetic and field data demonstrate that the proposed reservoir energy density can produce high-quality images for the fluid-saturated reservoirs, and it produces less background artefacts caused by elastic layers. This method provides a new way to detect the location of hydrocarbon reservoirs and characterize their spatial distribution.     

Ambient noise levels and detection threshold in Norway

Ambient seismic noise is caused by a number of sources in specific frequency bands. The quantification of ambient noise makes it possible to evaluate station and network performance. We evaluate noise levels in Norway from the 2013 data set of the Norwegian National Seismic Network as well as two temporary deployments. Apart from the station performance, we studied the geographical and temporal variations, and developed a local noise model for Norway. The microseism peaks related to the ocean are significant in Norway. We, therefore, investigated the relationship between oceanic weather conditions and noise levels. We find a correlation of low-frequency noise (0.125–0.25 Hz) with wave heights up to 900 km offshore. High (2–10 Hz) and intermediate (0.5–5 Hz) frequency noise correlates only up to 450 km offshore with wave heights. From a geographic perspective, stations in southern Norway show lower noise levels for low frequencies due to a larger distance to the dominant noise sources in the North Atlantic. Finally, we studied the influence of high-frequency noise levels on earthquake detectability and found that a noise level increase of 10 dB decreases the detectability by 0.5 magnitude units. This method provides a practical way to consider noise variations in detection maps.     

河南平顶山平煤矿区天然地震、爆破、塌陷时频特征分析

选取2019年3月—8月河南平顶山市宝丰县平煤矿区发生的ML 2.0—2.9天然地震、爆破、塌陷等9次震动事件,在区域地质构造背景和波形特征分析基础上,采用短时傅里叶变换(STFT)方法开展时频波谱分析,提取不同类型事件的时频特征。结果显示:(1)天然地震频率成分丰富,且高、低频分布均匀,P波在约3 Hz和8 Hz处存在2个峰值,S波存在多个峰值;(2)爆破事件的时频谱相对集中,以低频为主,P波频率峰值约5Hz,信号主频随时间变化,大致呈线性降低至1—2 Hz;(3)塌陷事件频率成分以4 Hz以下的低频为主,P波无明显峰值且频率成分单一,主频出现在2 Hz左右的面波。本文结果可为今后该矿区震动事件类型判断提供参考依据。     

基于相似性重构低频数据的金属矿频域全波形反演

由于深部金属矿埋深和自身的复杂性,利用重、磁、电方法和一般的地震方法很难有效地对其进行高精度定位.全波形反演通过最小化模拟数据与观测数据的差异使深部金属矿的高精度探测成为可能,但全波形反演是一个局部优化过程,需要准确的低频数据作为起始,而这在一般的地震数据采集中难以做到.本文先在频域中使用伴随状态震源函数反演方法,通过震源附近的直达波能精确地反演出震源函数的形态.然后利用得到的高精度震源子波结合褶积与反褶积思想及相似性现象重构含有低频成分的自激自收数据.将该数据应用到全波形反演中,有效缓减了反演过程中出现的周波跳跃现象,并提高了模型反演的正确性.Marmousi模型和金属矿模型的数值模拟实验证明了新方法改善了在没有低频数据时的全波形反演结果,并有较好抗噪性.     

基于弹性波全波形反演的主被动源多分量混采地震数据速度建模

在常规地震采集中,被动源地震波场往往被视为噪声而去除,这就造成了部分有用信息的丢失.在目标区进行主动源和被动源弹性波地震数据的多分量混合采集,并对两种数据进行联合应用,使其在照明和频带上优势互补,能显著提高成像和反演的质量.本文针对两种不同类型的主被动源混采地震数据,分别提出了相应的联合全波形反演方法.首先,针对主动源与瞬态被动源弹性波混采地震数据,为充分利用被动源对深部照明的优势,同时有效压制被动震源点附近的成像异常值,提出了基于动态随机组合的弹性波被动源照明补偿反演策略.然后,针对低频缺失主动源与背景噪声型被动源弹性波混采地震数据,为充分利用被动源波场携带的低频信息,并避免对被动源的定位和子波估计,提出了基于地震干涉与不依赖子波算法的弹性波主被动源串联反演策略.最后,分别将两种方法在Marmousi模型上进行反演测试.结果说明,综合利用主动源和被动源弹性波混采地震数据,不仅能增强深部弹性参数反演效果,还能更好地构建弹性参数模型的宏观结构,并有助于缓解常规弹性波全波形反演的跳周问题.     

地震低频伴影的数值模拟与应用

Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.     

Frequency dependent polarization analysis of ambient seismic noise recorded at a broadband seismometer in the central United States

We present a new approach to polarization analysis of seismic noise recorded by three-component seismometers.It is based on statistical analysis of frequency-dependent particle motion properties determined from a large number of time windows via eigenanalysis of the 3-by-3,Hermitian,spectral covariance matrix.We applied the algorithm to continuous data recorded in 2009 by the seismic station SLM,located in central North America.A rich variety of noise sources was observed.At low frequencies (0.05 Hz) we observed a tilt-related signal that showed some elliptical motion in the horizontal plane.In the microseism band of 0.05-0.25 Hz,we observed Rayleigh energy arriving from the northeast,but with three distinct peaks instead of the classic single and double frequency peaks.At intermediate frequencies of 0.5-2.0 Hz,the noise was dominated by non-fundamental-mode Rayleigh energy,most likely P and Lg waves.At the highest frequencies (3 Hz),Rayleigh-type energy was again dominant in the form of Rg waves created by nearby cultural activities.Analysis of the time dependence of noise power shows that a frequency range of at least 0.02-1.0 Hz (much larger than the microseism band) is sensitive to annual,meteorologically induced sources of noise.     

汶川与昆仑山强烈地震激发的地球自由振荡频谱的对比分析

强烈地震激发的地震波信息非常丰富,高频信息能量部分会在小范围短时间内迅速衰减,而低频信息特别是较低频信息能量部分则衰减较慢,会绕地球很多圈后才会消耗殆尽,从而激发地球的自由振荡.虽然地球自由振荡的固有频率不变,但不同震源机制的地震激发的地球自由振荡频率成分会随着震级、震源机制、方位等的不同而有显著差别.本文利用2008年5月12日汶川(M_w7.9)地震与2001年11月14日昆仑山(M_w7.8)地震这两个能量震级相当但震源机制不同的地震,基于地震波传播理论,考虑横向各向同性PREM地球模型并融合Crust2.0地壳模型、S20RTS地幔模型及PREM径向衰减模型,利用谱元法结合高性能并行计算,对这两次强烈地震激发的全球地震波传播进行了数值模拟,计算它们激发的地球自由振荡频谱特征差异.计算结果显示了汶川地震的逆冲-走滑特性及昆仑山地震的水平左旋走滑特征.利用理论计算得到的地震波动记录进行功率谱密度估计对比分析,研究结果显示这两次强烈地震释放的低频信息能量大小相当,但总体上看昆仑山地震释放的能量略小于汶川地震的.对比分析表明,不同地震激发的地球自由振荡模式在同一台站处的发育程度不同;同一地震激发的地球自由振荡在不同台站处记录到的振型频率成分也有很大差异.对于不同地震激发的地球自由振荡,有一些频率成分共同出现,表现为地球自由振荡功率谱能量均较大;也有一些振型由于相长干涉,仅在某些地震中被激发而得以良好发育,表现为功率谱能量明显较大;而一些振型由于相消干涉,在一些地震激发方式下被压制,表现为功率谱能量明显很小,甚至接近于0.计算结果清楚地显示了不同地震激发方式下地球振荡模式丰富而复杂的特性,需要对不同地震激发的地球振荡信息进行综合深入的分析,才能对地球内部结构尤其横向不均匀特性给予全面的科学解释.     

A real <Emphasis Type="Italic">P</Emphasis>-wave and its dependence on the presence of gas

Attenuation of seismic compression waves leads to the real existence of a fast P₁ wave in rocks which are fully saturated with dropping fluid and a slow P₂ wave in the rocks containing gas in their pores. This accounts for the seismic blanking zones below the gas horizons for the P₁ waves. Oscillations of gaseous inclusions ensure the energy transfer to the dominant frequencies which are different for the cases of passive seismic (few Hz) and active source seismic (10–20 Hz). The intervals of dominant frequencies are determined from the negative attenuation of these low-frequency waves. According to the observations and the suggested equation, random noise amplifies the signal at these frequencies. Thus, the P₂ waves at the dominant frequency of the active source seismics are applicable for elaborating on the details of the saturation of the production layer by hydrocarbons. The relation to the AVO method (Amplitude Variation with Offset) and dilatancy effect during the preparation of an earthquake is noted.     

Locating wind farms by seismic interferometry and migration

We present a case study on the detection and quantification of seismic signals induced by operating wind turbines (WTs). We spatially locate the sources of such signals in data which were recorded at 11 seismic stations in 2011 and 2012 during the TIMO project (Deep Structure of the Central Upper Rhine Graben). During this time period, four wind farms with altogether 12 WTs were in operation near the town of Landau, Southwest Germany. We locate WTs as sources of continuous seismic signals by application of seismic interferometry and migration of the energy found in cross-correlograms. A clear increase of emitted seismic energy with rotor speed confirms that the observed signal is induced by WTs. We can clearly distinguish wind farms consisting of different types of WTs (different hub height and rotor diameter) corresponding to different stable frequency bands (1.3–1.6 Hz, 1.75–1.95 Hz and 2.0–2.2 Hz) which do not depend on wind speed. The peak frequency apparently is controlled by the elastic eigenmodes of the structure rather than the passing of blades at the tower. From this we conclude that vibrations are coupled into the ground at the foundation and propagate as Rayleigh waves (and not as infrasound). The migration velocity of 320 m/s corresponds to their group velocity. The applied migration method can contribute to the assessment of local sources of seismic noise. This topic gets growing attention in the seismological community. In particular, the recent boost of newly installed wind farms is a threat to seismological observatories such as the Black Forest Observatory (BFO) and the Gräfenberg array (GRF) or gravitational wave observatories (e.g. LIGO, VIRGO) in terms of a sensitivity degradation of such observatories.     

宽频带数字资料低频波在大震前的短临前兆信息研究

通过对高台、兰州、成都、安西地震台宽频带数字资料在2001年昆仑山口西8.1级地震,2008年汶川8.0级地震以及2003年山丹、民乐6.1级地震前低频波频谱变化的研究,总结了数字资料在频率域上的短临前兆异常特征,发现0.1~0.14 Hz频幅在震前几天成逐渐增强的趋势。并结合震源区构造特征对前兆机理进行了探讨。认为低频波的产生与静地震和慢地震有关。     

结晶岩地区深地震数据采集关键技术与方法

深反射地震是了解深部地质结构的主要手段,获取高品质的数据是给出合理地质解释的基础.在结晶岩地区,由于地层成层性差、非均质性严重,地震散射效应明显,导致地震波场复杂,同时结晶岩为非理想弹性体,不利于地震波能量转换,有效反射能量弱,信噪比低.加之深反射地震的目标层较深且受环境噪声干扰严重,高频信号衰减快,地震资料主频较低.这些因素使得在结晶岩区难以获得高品质的地震资料,为探索提高该类地区深反射地震资料质量的方法和技术,本文依托长江中下游成矿带2009—2014年深反射地震数据采集工作,在精细设计、严格施工的基础上,从激发和接收入手,开展了"轴向不耦合激发"、"宽频接收"和"宽线观测"等技术方法试验研究.结果表明,这些方法技术措施提高了下传弹性波能量,展宽了接收地震信号的频带,提高了覆盖次数和信噪比,有效改善了地震原始资料的品质和成像效果.研究结果对今后结晶岩地区深地震反射数据采集工作具有重要的实用价值和参考意义.     

新冠肺炎疫情对佘山与大洋山地震台背景噪声的影响

人类的生产生活产生的振动会以高频地震波的形式被地震台站所记录。2020年1月,新冠肺炎疫情爆发,为了应对此次疫情,各地政府分别启动应急响应,国内地震记录出现最长、最突出的人为地震降噪期。各地震台站背景噪声显著下降,在人口稠密及工业发达地区尤为明显。同时,静噪期为探测地下地震源的微弱信号提供了可能。静噪期内,佘山地震台2 Hz频点背景噪声功率谱密度值比平时降低10 dB,而大洋山地震台10 Hz频点背景噪声功率谱密度值较平时降低约5 dB;佘山地震台2—10 Hz以及大洋山地震台10 Hz以上频率的背景噪声受静噪期影响明显。结合地震台站所处位置分析,疫情期间佘山地震台附近人口流出较多,2—10 Hz频率的背景噪声变化明显;大洋山地震台背景噪声则受工厂、轮渡、高速汽车等影响较大,f ≥10 Hz的背景噪声变化显著,而频率在2—10 Hz则无明显变化,表明该台人口总数趋于平稳。地震噪声和人类活动之间的相关性表明,地震学研究可以提供实时人口动态估计。     

西北太平洋海岛地区地震背景噪声特征及海洋学解释

地震背景噪声特性及噪声源的分布研究逐渐成为深化背景噪声层析成像的关键问题.海岛地区由于特殊的地理位置,其背景噪声具有相对独特的特征.地脉动（约0.003~1 Hz）是地震背景噪声中能量最强的分量,其激发与特性被认为与海浪运动和固体地球之间的相互作用有关,但海岛地区地脉动特征与海洋波浪场之间的关系尚未被充分研究.本文利用西北太平洋海岛地震台站的连续记录数据、波浪浮标的实测数据以及WAVEWATCH-Ⅲ海浪模式的数值模拟结果,通过地震学和海洋学的交叉,分析海岛地区地脉动信号的时频特性及其与海洋波浪场之间的相关性,从海洋学角度对地脉动信号的特征及激发进行探讨与解释.结果表明,海岛地区地脉动信号相对于内陆地区更强,并具有明显且稳定的季节性变化特征：高频地脉动信号（0.12~0.32 Hz）在夏秋季节（5月-10月）相对较弱,而在冬春季节（11月-次年4月）相对较强,与北半球海洋活动季节性变化相一致.此外,海岛地区地脉动主要受周边海域波浪场影响,与周边海域波浪能功率密度及实测和数值模拟所得的有效波高均具有很好的互相关性.该研究结果同时表明可进一步发展利用地脉动观测数据反演海表波浪场的可能,为海洋科学研究中海表波浪场连续观测数据的获取提供地震学上的支持.     

Feasibility of borehole ambient noise interferometry for permanent reservoir monitoring

The analysis of seismic ambient noise acquired during temporary or permanent microseismic monitoring campaigns (e.g., improved/enhanced oil recovery monitoring, surveillance of induced seismicity) is potentially well suited for time‐lapse studies based on seismic interferometry. No additional data acquisition required, ambient noise processing can be automatized to a high degree, and seismic interferometry is very sensitive to small medium changes. Thus there is an opportunity for detection and monitoring of velocity variations in a reservoir at negligible additional cost and effort. Data and results are presented from an ambient noise interferometry study applied to two wells in a producing oil field in Romania. Borehole microseismic monitoring on three component geophones was performed for four weeks, concurrent with a water‐flooding phase for improved oil recovery from a reservoir in ca. 1 km depth. Both low‐frequency (2 Hz–50 Hz) P‐ and S‐waves propagating through the vertical borehole arrays were reconstructed from ambient noise by the virtual source method. The obtained interferograms clearly indicate an origin of the ambient seismic energy from above the arrays, thus suggesting surface activities as sources. It is shown that ambient noise from time periods as short as 30 seconds is sufficient to obtain robust interferograms. Sonic log data confirm that the vertical and horizontal components comprise first arrivals of P‐wave and S‐waves, respectively. The consistency and high quality of the interferograms throughout the entire observation period further indicate that the high‐frequency part (up to 100 Hz) represents the scattered wave field. The temporal variation of apparent velocities based on first‐arrival times partly correlates with the water injection rate and occurrence of microseismic events. It is concluded that borehole ambient noise interferometry in production settings is a potentially useful method for permanent reservoir monitoring due to its high sensitivity and robustness.     

应用低频信号提高高速玄武岩下的成像质量

针对高速玄武岩屏蔽层下深层成像困难的实际问题,采用波动方程波场数值模拟技术,根据玄武岩地层的特点,设计三个相应的简单高速玄武岩模型,通过对深层反射地震信号能量的分析,说明了低频地震信号既具有较强的穿透薄高速玄武岩屏蔽层的能力,也具有减弱因粗糙表面所产生的绕射噪音的能力。一个完整的2D玄武岩模型的模拟试验证明了利用低频信号可以提高高速玄武岩屏蔽层下深层成像的质量,实际资料的低通滤波处理也取得了预期的效果。     

2011年日本MW9.0地震引发的海啸对地震背景噪声的影响

基于北美沿岸和内陆地震台站的连续地震波形记录,并结合沿岸台站附近布设的DART系统记录的海底压力数据以及预测潮汐数据,利用时频分析和极化分析方法对2011年3月11日日本东北部海域MW9.0大地震所激发的海啸对地震背景噪声所产生的影响予以深入分析.结果显示:海啸对高频噪声(1.3—1.5 Hz)以及短周期双频微地动噪声...     

Acoustic noise and temperature monitoring of the Crater Lake of Mount Ruapehu Volcano

Acoustic signals in Ruapehu Crater Lake, which are now being telemetered via a satellite transmission system, show promise as a possible precursor of increased volcanic activity from Ruapehu. The start of a recent period of rapid heating of Crater Lake was preceded by low-frequency (2 Hz) acoustic signals. These accompanied similar frequency seismic signals, but seemed to be produced independently. Audio-frequency (350–3000 Hz) acoustic noise also showed a very clear peak shortly before the lake temperature started to rise.     

汶川大震前数字地震仪位移地脉动低频异常现象研究

选取了2008年4月25日至5月12日间的四川台网14个台宽频带地震仪记录到的地脉动波形资料,把速度波形做适当的数据校正和滤波处理后,经积分转换成位移波形.频谱分析结果表明:地震发生前存在低频异常颤动,地震前2～3天开始0.1～0.3Hz频段傅里叶谱幅值逐渐增大;各个台的位移振幅值逐渐增强.其中靠近龙门山断裂带的台站记录到的地脉动变化较明显.     

饱和砂岩的黏弹行为的实验研究

通过Metravib热机械分析仪,用正弦波进行加载,实验时固定静载为100 N,正弦波动载荷恒为60 N,将总载荷控制在屈服点以下;在温度为-50～125℃,升温速率为1℃/min,频率为5～1000 Hz的条件下,对泵油饱和长石砂岩、彭山砂岩样品进行单轴循环加载实验,求取泵油饱和长石砂岩和彭山砂岩的衰减、耗散角、模量...