Enabling affordable omnidirectional subsurface extended image volumes via probing

Image gathers as a function of subsurface offset are an important tool for the inference of rock properties and velocity analysis in areas of complex geology. Traditionally, these gathers are thought of as multidimensional correlations of the source and receiver wavefields. The bottleneck in computing these gathers lies in the fact that one needs to store, compute, and correlate these wavefields for all shots in order to obtain the desired image gathers. Therefore, the image gathers are typically only computed for a limited number of subsurface points and for a limited range of subsurface offsets, which may cause problems in complex geological areas with large geologic dips. We overcome increasing computational and storage costs of extended image volumes by introducing a formulation that avoids explicit storage and removes the customary and expensive loop over shots found in conventional extended imaging. As a result, we end up with a matrix–vector formulation from which different image gathers can be formed and with which amplitude‐versus‐angle and wave‐equation migration velocity analysis can be performed without requiring prior information on the geologic dips. Aside from demonstrating the formation of two‐way extended image gathers for different purposes and at greatly reduced costs, we also present a new approach to conduct automatic wave‐equation‐based migration‐velocity analysis. Instead of focusing in particular offset directions and preselected subsets of subsurface points, our method focuses every subsurface point for all subsurface offset directions using a randomized probing technique. As a consequence, we obtain good velocity models at low cost for complex models without the need to provide information on the geologic dips.     

Analysis of seismic noise in the Romanian-Bulgarian cross-border region

We present the first detailed investigation of the background seismic noise recorded in the Romanian-Bulgarian cross-border region over 3 years (2012–2015). We used the power spectral densities probability density functions (PSD PDFs) to study the noise variations in the period domain (0.025–1 s) as well as in the secondary microseism band (2–10 s). Strong diurnal variations and an increase of the noise levels during working days were observed at high frequencies at all stations, thus confirming the anthropic origin of the noise at low periods. The noise variations observed at longer periods (>?1 s) are relatively small among the stations and are related to season changes. The dominant feature in the noise spectra between 2 and 10 s is the double-frequency peak (DFP) whose amplitude increases and changes during winter. For a specific interval, from 25th to 27th of January 2014, when a storm was reported in the Black Sea area, the maximum of the DFP shifted from larger periods (~?5.5 s) at stations far from the Black Sea towards smaller periods (~?1.8 s) at stations located on the coastline. The polarization analysis showed that the short period double-frequency microseisms originating from the Black Sea dominate during the winter month. Finally, we showed that site conditions vary due to noise variations related to weather conditions in the Black Sea or to changes in anthropogenic noise sources.     

2015—2018年海南省臭氧污染特征及其天气类型

根据2015—2018年海南省18个市县32个空气质量监测站O₃浓度资料,分析了区域性O₃污染（O₃-8h浓度超标市县≥3个）时空变化特征,并对造成O₃污染的天气系统进行主观分型。结果表明：2015—2018年海南省共有40 d发生了区域性O₃污染,发生概率为2.73%。其中2015年和2017年达到了13 d,发生概率为3.56%,2018年为11 d（3.01%）,2016年仅为3 d（0.82%）。发生区域性O₃污染主要有4种天气类型：冷空气偏西下型、冷空气偏东下型、变暖高压脊型和热带系统型。其中冷空气偏西下型是最主要的天气类型,共出现了14 d,占所有天数的35%,且污染较重。不同天气类型下海南省O₃污染表现出不同的分布特征。500 hPa有下沉气流、低层受东北风控制,有相对湿度低值区从中国东部向海南省延伸,地面位于冷高压底部或热带气旋西北侧,温度露点差在5 ℃以上等条件均有利于海南省区域性O₃污染天气出现。     

黏声介质平面波有限差分叠前深度偏移及在莺歌海盆地的应用

莺琼盆地诸探区中存在底劈现象,深层气源产生的气体沿底劈产生的裂隙通道向上漫溢.漫溢过程中,一是充填在遇到的砂体中形成气藏;二是弥漫在上溢通道中,使得通道中的纵波速度发生变化,进而纵波波阻抗差异变小,反射变弱.另外,通道中气体的存在,会加强地层的吸收衰减,使得地震波振幅变弱、高频成分损失导致同相轴分辨率降低.利用OBC数据进行多波地震勘探和利用黏声介质的叠前深度偏移都是改善模糊区成像质量的重要方法技术.为此,本文提出用黏声介质平面波有限差分法叠前深度偏移成像方法改善气体充填区域的成像质量.黏声介质成像目的是补偿地震波的吸收衰减;平面波偏移成像目的是适应海上单炮数量巨大,提高波动方程叠前深度偏移成像的效率;有限差分法叠前偏移的目的是适应该区浅层气分布局域性极强、Q值的空间变化大的情况.在莺歌海某探区的实际数据上的黏声介质平面波有限差分叠前深度偏移试验证明,本方法是改善模糊区成像质量的较为有效的途径.     

The dayside magnetospheric boundary layer at Mercury

Magnetic field and plasma data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft on the outbound portions of the first (M1) and second (M2) flybys of Mercury reveal a region of depressed magnetic field magnitude and enhanced proton fluxes adjacent to but within the magnetopause, which we denote as a dayside boundary layer. The layer was present during both encounters despite the contrasting dayside magnetic reconnection, which was minimal during M1 and strong during M2. The overall width of the layer is estimated to be between 1000 and 1400 km, spanning most of the distance from the dayside planetary surface to the magnetopause in the mid-morning. During both flybys the magnetic pressure decrease was ∼1.6 nPa, and the width of the inner edge was comparable to proton gyro-kinetic scales. The maximum variance in the magnetic field across the inner edge was aligned with the magnetic field vector, and the magnetic field direction did not change markedly, indicating that the change in field intensity was consistent with an outward plasma-pressure gradient perpendicular to the magnetic field. Proton pressures in the layer inferred from reduced distribution observations were 0.4 nPa during M1 and 1.0 nPa during M2, indicating either that the proton pressure estimates are low or that heavy ions contribute substantially to the boundary-layer plasma pressure. If the layer is formed by protons drifting westward from the cusp, there should be a strong morning–afternoon asymmetry that is independent of the interplanetary magnetic field (IMF) direction. Conversely, if heavy ions play a major role, the layer should be strong in the morning (afternoon) for northward (southward) IMF. Future MESSENGER observations from orbit about Mercury should distinguish between these two possibilities.     

Measurement of the radius of Mercury by radio occultation during the MESSENGER flybys

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft completed three flybys of Mercury in 2008–2009. During the first and third of those flybys, MESSENGER passed behind the planet from the perspective of Earth, occulting the radio-frequency (RF) transmissions. The occultation start and end times, recovered with 0.1 s accuracy or better by fitting edge-diffraction patterns to the RF power history, are used to estimate Mercury's radius at the tangent point of the RF path. To relate the measured radius to the planet shape, we evaluate local topography using images to identify the high-elevation feature that defines the RF path or using altimeter data to quantify surface roughness. Radius measurements are accurate to 150 m, and uncertainty in the average radius of the surrounding terrain, after adjustments are made from the local high at the tangent point of the RF path, is 350 m. The results are consistent with Mercury's equatorial shape as inferred from observations by the Mercury Laser Altimeter and ground-based radar. The three independent estimates of radius from occultation events collectively yield a mean radius for Mercury of 2439.2±0.5 km.     

Seismic and Petrophysical Relationships from UK Continental Shelf Rocks at Elevated Pressures

Abstract

There are few published seismic (P- and S-waves) properties for seafloor bedrocks. At low pressures (1 to 10 MPa), velocities and attenuations are determined mainly by open microcracks. At higher pressures, the microcracks close, and the velocities and attenuations depend primarily on the matrix porosity. We have investigated both the relationships between the acoustic, petrophysical, and geological properties of the rocks at 40 MPa pressure and the effect of microcracks on the acoustic properties at 10 MPa pressure. In this paper we discuss the former; the latter will be discussed separately.

P- and S-wave velocity and attenuation measurements were carried out on a suite of seabed sedimentary and igneous rocks at effective pressures from 10 to 40 MPa at ultrasonic frequencies. The porosities and permeabilities of the rocks ranged from 0% to 32% and 0 to 110 mDarcy, respectively. Characterization of the rocks revealed that most of the sandstones have a substantial clay content (kaolinite, illite, and chlorite) and fractures. Most of the igneous rocks are chloritized.

The seismic properties of the rocks are markedly lower than those of similar continental rocks because of the microporosity formed by the alteration of feldspars, micas, and mafic minerals to clays (e.g., chloritization of pyroxenes) and the corresponding reduction of the elastic moduli. The results of this study suggest that the values of velocities and quality factors used for ocean acoustic propagation models are lower than normally assumed.     

塔里木盆地塔北露头区的古岩溶发育模式

目前塔里木盆地海相碳酸盐岩的油气勘探已进入高峰期,为弥补钻井岩心资料不足,进一步认识井下岩溶发育规律及发育模式特征,开展塔北露头区奥陶系风化壳古岩溶缝洞系统发育特征及发育模式研究。研究表明,岩溶发育主要受控于岩性因素和不整合面,主要发育层位为一间房组地层;主要岩溶作用期为吐木休克组／一间房组沉积间断岩溶期。为井下岩溶缝洞系统的认识及岩溶地质模型的建立提供了重要理论依据。     

全球变暖背景下新疆地区近45 a来最大冻土深度变化及其突变分析

利用覆盖新疆大部分地区资料完整的93个站点资料,对1961-2005年新疆地区最大冻土深度进行了分析. 结果表明：新疆地区月最大冻土深度有明显的季节变化,低海拔区域(海拔<1 800 m)最大值出现在1月份,而高海拔区域(海拔≥1 800 m)的最大值出现在2月份,比低海拔区域要滞后. 新疆地区最大冻土深度的地理分布特征表现为北疆深于南疆,山区深于平原,且与气温的分布有很好的一致性. 全年和冬、春季最大冻土深度与气温场的空间相关系数分别为-0.795、-0.736和-0.848. 年际变化表明,近45 a来的最大冻土深度出现了较为明显的下降. 高海拔区域与低海拔区域年最大冻土深度的倾向率分别为-15.65 cm·(10a)^-1和-9.48 cm·(10a)^-1,且与气温的相关系数分别为-0.51和-0.69,均通过了0.001的信度检验. 同时发现,高海拔区域冬季下降多,而低海拔区域春季下降多. 新疆地区年最大冻土深度在近45 a有明显的突变现象,高海拔区域和低海拔区域突变发生年份分别为1996/1997年度和1978/1979年度,说明新疆地区高海拔区域的年最大冻土深度对气温变化的响应比低海拔区域要滞后. 突变年后高海拔区域与低海拔区域年最大冻土深度比突变年前的平均值分别降低了61.12 cm和26.67 cm.