The added value of joint PP–PS inversion for reservoir characterization: A case study using Jubarte PRM seismic data (offshore Brazil)

The added value of the joint pre-stack inversion of PP (incident P-wave and reflected P-wave) and PS (incident P-wave and reflected S-wave) seismic data for the time-lapse application is shown. We focus on the application of this technique to the time-lapse (four-dimensional) multicomponent Jubarte field permanent reservoir monitoring seismic data. The joint inversion results are less sensitive to noise in the input data and show a better match with the rock physics models calibrated for the field. Further, joint inversion improves S-impedance estimates and provides a more robust quantitative interpretation, allowing enhanced differentiation between pore pressure and fluid saturation changes, which will be extremely useful for reservoir management. Small changes in reservoir properties are expected in the short time between the time-lapse seismic acquisitions used in the Jubarte project (only 1 year apart). The attempt to recover subtle fourth-dimensional effects via elastic inversion is recurrent in reservoir characterization projects, either due to the small sensitivity of the reservoirs to fluid and pressure changes or the short interval between the acquisitions. Therefore, looking for methodologies that minimize the uncertainty of fourth-dimensional inversion outputs is of fundamental importance. Here, we also show the differences between PP only and joint PP–PS inversion workflows and parameterizations that can be applied in other projects. We show the impact of using multicomponent data as input for elastic seismic inversions in the analysis of the time-lapse differences of the elastic properties. The larger investment in the acquisition and processing of multicomponent seismic data is shown to be justified by the improved results from the fourth-dimensional joint inversion.     

北京周口店雾迷山组碳酸盐岩风暴沉积序列研究*

北京周口店地区中元古界雾迷山组以含硅质条带的白云岩为主,发育不规则硅质砾屑。这些硅质砾屑形态各异,具有强烈扰动的迹象,目前对其成因解释一直存在争议,如泄水构造、冲刷构造、震积岩、风暴岩等。文中选择北京周口店地区黄山店村恒顺厂剖面为研究对象,通过野外精细的沉积学分析,系统描述了岩性及沉积构造特征,确定了该套沉积地层发育典型的碳酸盐岩风暴沉积序列。该风暴沉积序列自下往上可分为5段: A段,为风暴前正常沉积的中厚层状泥晶白云岩;B段,具口袋构造的侵蚀冲刷面及硅质砾屑,为风暴高潮期产物;C段,灰色硅质条带白云岩,普遍发育平行层理和波状、丘状、洼状交错层理,是风暴衰减期产物;D段,灰白色薄层状泥晶白云岩,代表了风暴间歇期的正常沉积;E段,含硅质团块白云岩,为正常天气下海水局部扰动形成。该剖面上递变层理几乎不发育,结合侵蚀冲刷面特征、砾屑高含量以及杂乱堆积方式,认为该套沉积为典型的近原地风暴沉积。结合古地理资料,推测研究区在雾迷山组沉积时期可能处于热带海洋环境,热带气旋引发的频繁风暴潮导致了近原地碳酸盐岩风暴沉积的形成。上述研究成果不仅补充了燕山地区雾迷山组的风暴沉积记录,而且为该时期华北板块的古地理、古纬度和古气候研究提供了重要的沉积学证据。     

气候变暖对湖泊热力及溶解氧分层影响研究进展

气候变暖对湖泊物理、化学、生物和生态系统有着复杂而深刻的直接和间接影响,而具体影响随研究区域和水体表现不尽相同。气候变暖通过改变湖泊热力和溶解氧分层进而影响湖泊生物过程和生态系统结构与功能。从全球湖泊变暖趋势、长期缓慢气温上升、极端高温事件以及气候情景模拟等方面详细综述了气候变暖对湖泊热力及溶解氧分层影响的研究进展。研究表明,全球不同区域湖泊均存在不同程度的变暖趋势;长期缓慢气温上升和短期极端高温均会造成湖泊热力分层提前,分层结束推迟,分层时间延长,混合层和温跃层深度下降,以及热稳定性增加;相伴随的是溶解氧扩散深度和氧跃层深度明显下降,加剧了湖泊底部好氧和厌氧环境。除了这种直接影响外,气候变暖引起的流域降水、入湖物质的变化以及风速的变化也会对湖泊热力和溶解氧分层产生许多间接的影响,因此未来仍然需要更多的实验证据、经验和数值模型来验证和预测气候变暖对湖泊热力及溶解氧分层的影响。     

两次秸秆焚烧导致辽宁重污染过程的气象条件对比分析

利用环境监测站大气污染物数据、地面自动气象站观测资料、L波段加密探空资料和0.125°×0.125°的EC再分析资料,结合MODIS遥感火点监测和HYSPLIT4后向轨迹模拟结果,对比分析了2015年11月8日和2016年11月5日的两次由于东北地区秸秆焚烧导致辽宁重污染天气过程的大气边界层特征、气象扩散条件和大气污染物输送来源等。结果表明:两次过程地面PM_(2.5)浓度均出现快速上升和下降,其中2015年11月8日重污染过程的污染强度较2016年11月5日强,且持续时间更长。2015年11月8日重污染过程的混合层高度较低,其上层的中性层结转变为逆温层结,抑制混合层高度的发展。同时低层冷平流不断侵入到暖平流下方,使得大气层结稳定性增强,维持时间较2016年11月5日重污染过程更长,低层下沉运动和黑龙江西南部、吉林西部污染物的远距离输送增强使得辽宁地面污染物浓度快速累积。而2016年11月5日重污染天气过程主要受深厚冷空气影响,东北地区西部污染物的区域输送和地面风场辐合是地面污染物浓度快速上升的主要原因。     

杭州市典型雨转雪天气成因及预报模型

利用2008—2018年的NCEP(1°×1°)再分析资料、常规气象观测资料和降雪加密观测资料,选出杭州地区10次典型的雨转雪天气过程,从大尺度环流背景和动力、水汽以及热力因子等物理量场结构方面展开研究,最终得出杭州冬季典型雨转雪天气的预报模型:①大尺度环流配置需满足能为雨转雪天气的形成提供有利的水汽、动力抬升以及中低层上暖下冷的逆温或等温层结条件;②水汽和动力因子等物理量须满足产生纯雪的特定条件;③杭州温度层结须为T_(2m)≤1.5℃、T_(925)≤-4.0℃、T_(850)≤0℃、T_(700)≤-1.0℃和T_(500)≤-10.0℃。此外,进一步补充了杭州可能产生大雪甚至暴雪量级降雪的特定条件。最终选取2019年初的2次典型降水过程进行预报回报检验。     

“尤特”强降水过程大气层结分析

1311号强台风“尤特”登陆后给广东带来持续性大范围强降水,对流降水特征显著.分析了“尤特”影响期间大尺度环流背景,重点讨论了此次持续性强降水过程中大气层结问题.发现低空急流向广东输送强的暖平流,是广东大气层结不稳定得以持续维持的根本原因.进一步分析发现,低空急流本身并不是“暖”的,当“尤特”趋向陆地时,陆地上的暖气团在“尤特”环流强迫下向南传播扩散,低空急流穿越这一暖区时温度升高才具备“暖”的特性.这一事实在以前并未被关注到.通过个例反查,在许多登陆后造成连续强降水的台风过程中均发现了这一特征.因此,台风登陆引起环境温度场的演变以及与低空急流的配置需引起重视.     

The F-type asteroids with small inversion angles of polarization

UBVRI polarimetric observations carried out in 1997-2004 for the F-type Asteroids 302 Clarissa, 419 Aurelia, 704 Interamnia, and 762 Pulcova (V band only) are presented. Asteroid 419 Aurelia is characterized by a negative polarization branch which is unusual for low-albedo asteroids. Its depth is about 1%, while the inversion angle, close to 14°, reaches the smallest value ever observed for asteroids. This is the first definite example of a minor body exhibiting a considerable decrease of both the depth and width of the negative polarization branch in comparison with polarization properties of other low-albedo bodies, mimicking a behavior previously found in laboratory measurements of extremely dark surfaces. The F-type Asteroids 302 Clarissa and 704 Interamnia are also characterized by unusually small inversion angles compared to other asteroid types, while Asteroid 762 Pulcova seems to have an ordinary negative polarization branch. Laboratory measurements of low-albedo samples and computer simulations of light scattering by particles of irregular shapes were made to interpret observational data. We find that an optical homogeneity of regolith microstructure at scales of the order of visible light wavelengths may be responsible for relatively small values of the depth of the negative polarization branch and of the inversion angle. Peculiar features of the F-type asteroids compared to other taxonomic classes are discussed.     

Multi-scale Tomography for Crustal <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Velocities in Southern California

Seismic tomography is a viable tool in building depth-velocity models in the presence of strong lateral velocity variations. In this study 3-D P- and S-velocity models for the crust of southern California are constrained using more than 1,000,000 P-wave first arrivals and 130,000 S-wave arrivals from local earthquakes. To cope with the uneven distribution of raypaths, a multi-scale tomography is applied with overlapping model cells of different sizes. Within the 300 × 480 × 39 km³ model volume, the smallest cell size is 10 × 10 × 3 km³. During the iterations of velocity updating, earthquake hypocenters are determined using both P and S arrivals, and full 3-D ray tracing is implemented. Except near the edges and in the lower crust, the resultant models are robust according to various tests on the effects of reference models, resolution and signal-to-noise ratio. The tomographic velocities at shallow depths correlate very well with the regional geology of southern California. In the upper crust the P-wave and S-wave models exhibit slow velocities in major sedimentary basins and fast velocities in areas of crystalline rocks. Mid-crustal low velocity zones are present under the Coso Range, San Gabriel Mountains, and a large portion of the Mojave Desert. P- and S-velocity patterns maintain their similarity in the lower crust though the models are less reliable there.     

Processes affecting deposition of sediment in a small, morphologically complex lake

Devil Lake is morphologically complex as a result of Pleistocene glacial erosion of the Frontenac Axis of the Canadian Shield. In order to assess the processes causing highly variable sedimentation in the lake, we monitored currents, suspended sediment and temperature in the lake before and during autumn overturn in 2002. Strong summer thermal stratification (stability number to 0.11 s^–1 declining with the approach of overturn) was insufficient to prevent a dynamic response in the hypolimnion to wind forcing. Superimposed on a gradual increase in suspended sediment concentration in the last weeks of stratification from less than 2 g/l to about 30 g/l were shorter-term rises lasting up to several days. Associated with these events was an increase in particle size of the sediment from a mode of 40–50 to 150–200 m ascribed to flocculation from primary particles. These events culminated in rapid (<1 h) clearing of the water associated with strong, sustained winds over the lake, especially from the southwest. After overturn, the events were more frequent, and flocculation was unable to develop as well in the more vigorous circulation. However, currents in the hypolimnion occurred throughout the period before, during and after overturn with speed related to wind speed, but direction largely independent of wind direction. The results represent an approach to understanding the nature of sedimentary processes and thus to strengthening the use of sedimentary records as proxy in environmental and paleoenvironmental assessment.     

Observations and modelling of the travel time delay and leading negative phase of the 16 September 2015 Illapel,Chile tsunami

The systematic discrepancies in both tsunami arrival time and leading negative phase (LNP) were identified for the recent transoceanic tsunami on 16 September 2015 in Illapel, Chile by examining the wave characteristics from the tsunami records at 21 Deep-ocean Assessment and Reporting of Tsunami (DART) sites and 29 coastal tide gauge stations. The results revealed systematic travel time delay of as much as 22 min (approximately 1.7％ of the total travel time) relative to the simulated long waves from the 2015 Chilean tsunami. The delay discrepancy was found to increase with travel time. It was difficult to identify the LNP from the near-shore observation system due to the strong background noise, but the initial negative phase feature became more obvious as the tsunami propagated away from the source area in the deep ocean. We determined that the LNP for the Chilean tsunami had an average duration of 33 min, which was close to the dominant period of the tsunami source. Most of the amplitude ratios to the first elevation phase were approximately 40％, with the largest equivalent to the first positive phase amplitude. We performed numerical analyses by applying the corrected long wave model, which accounted for the effects of seawater density stratification due to compressibility, self-attraction and loading (SAL) of the earth, and wave dispersion compared with observed tsunami waveforms. We attempted to accurately calculate the arrival time and LNP, and to understand how much of a role the physical mechanism played in the discrepancies for the moderate transoceanic tsunami event. The mainly focus of the study is to quantitatively evaluate the contribution of each secondary physical effect to the systematic discrepancies using the corrected shallow water model. Taking all of these effects into consideration, our results demonstrated good agreement between the observed and simulated waveforms. We can conclude that the corrected shallow water model can reduce the tsunami propagation speed and reproduce the LNP, which is observed for tsunamis that have propagated over long distances frequently. The travel time delay between the observed and corrected simulated waveforms is reduced to <8 min and the amplitude discrepancy between them was also markedly diminished. The incorporated effects amounted to approximately 78％ of the travel time delay correction, with seawater density stratification, SAL, and Boussinesq dispersion contributing approximately 39％, 21％, and 18％, respectively. The simulated results showed that the elastic loading and Boussinesq dispersion not only affected travel time but also changed the simulated waveforms for this event. In contrast, the seawater stratification only reduced the tsunami speed, whereas the earth's elasticity loading was responsible for LNP due to the depression of the seafloor surrounding additional tsunami loading at far-field stations. This study revealed that the traditional shallow water model has inherent defects in estimating tsunami arrival, and the leading negative phase of a tsunami is a typical recognizable feature of a moderately strong transoceanic tsunami. These results also support previous theory and can help to explain the observed discrepancies.