沿阿穆尔板块西边界的活动断层(蒙古领土)

阿穆尔板块西部边界在蒙古境内的空间位置尚不清楚,并且活动断层构造及其沿线地壳的应力状态研究较少。本文在沿此边界的三个区域——杭爱—肯特构造鞍部、布尔古特地块(鄂尔浑—土拉交汇处)和色楞格地块(包括色楞格凹陷和布伦—努鲁隆起),利用空间图像解译、地形起伏度分析、地质构造资料以及构造压裂和沿裂缝位移资料重建构造古应力,对活动断层进行研究。研究表明,活动断裂继承了古生代和中生代古构造的非均质性。这些断层沿着板块边界并不是单一的带,而是成簇的。它们的运动取决于走向:亚纬向断层是具有一定逆分量的左旋走滑断层,北西向断层是逆断层或逆冲断层,通常具有右旋走滑分量,海底断层是右旋走滑断层,北东向断层是正断层。位于色楞格凹陷和杭爱东部的断裂构造的活动始于上新世。逆断层和走滑断层与上新世情况不符,但多与更新世地貌相符,表明其活动年代较晚,为更新世时期。利用构造断裂和沿断裂的位移,重建活动断裂带变形末阶段的应力应变状态,结果表明断裂在最大挤压轴的北北东和北东方向上以压缩和走滑为主。只有在色楞格凹陷内,以扩张和走滑类型的应力张量为主,且在最小挤压轴的北西走向尤为显著。在南部,杭爱东部(鄂尔浑地堑)内有1个以扩张机制为主的局部区域,说明蒙古中部断裂在更新世—全新世阶段的活动以及现代地震活动主要受与印度斯坦和欧亚大陆汇聚过程相关的东北方向的附加水平挤压的控制。使研究区地壳产生走滑变形、贝加尔湖裂谷发散活动以及阿穆尔板块东南运动的另一个因素是东南方向软流圈流动对岩石圈底部的影响。阿穆尔板块和蒙古地块之间的边界在构造结构上是零碎的,代表了覆盖整个蒙古西部变形带的边缘部分。     

Influence of background heterogeneity on traveltime shifts for compacting reservoirs

Compaction induced by pore‐pressure decrease inside a reservoir can be monitored by measuring traveltime shifts of reflection events on time‐lapse seismic data. Recently we introduced a perturbation‐based formalism to describe traveltime shifts caused by the 3D stress‐induced velocity field around a compacting reservoir. Application of this method to homogeneous background models showed that the offset variation of traveltime shifts is controlled primarily by the anisotropic velocity perturbations and can provide valuable information about the shear and deviatoric stresses. Here, we model and analyse traveltime shifts for compacting reservoirs whose elastic properties are different from those of the surrounding medium. For such models, the excess stress is influenced primarily by the contrast in the rigidity modulus μ across the reservoir boundaries. Synthetic examples demonstrate that a significant (25% or more) contrast in μ enhances the isotropic velocity perturbations outside the reservoir. Nevertheless, the influence of background heterogeneity is mostly confined to the reservoir and its immediate vicinity and the anisotropic velocity changes are still largely responsible for the offset dependence of traveltime shifts. If the reservoir is stiffer than the host rock, the background heterogeneity reduces anisotropic velocity perturbations inside the reservoir but increases them in the overburden. As a result, in this case, the magnitude of the offset variation of traveltime shifts is generally higher for reflections from interfaces above the reservoir. We also study compaction‐induced stress/strain and traveltime shifts for a stiff reservoir embedded in a softer layered model based on velocity profiles from the Valhall Field in the North Sea. Despite producing discontinuities in strain across medium interfaces, horizontal layering does not substantially alter the overall behaviour of traveltime shifts. The most pronounced offset variation of traveltime shifts is observed for overburden events recorded at common midpoints close to the reservoir edges. On the whole, prestack analysis of traveltime shifts should help better constrain compaction‐induced velocity perturbations in the presence of realistic background heterogeneity.     

Prestack data enhancement with phase corrections in time-frequency domain guided by local multidimensional stacking

We present a new approach to enhancing weak prestack reflection signals without sacrificing higher frequencies. As a first step, we employ known multidimensional local stacking to obtain an approximate ‘model of the signal’. Guided by phase spectra from this model, we can detect very weak signals and make them visible and coherent by ‘repairing’ corrupted phase of original data. Both presented approaches – phase substitution and phase sign corrections – show good performance on complex synthetic and field data suffering from severe near-surface scattering where conventional processing methods are rendered ineffective. The methods are mathematically formulated as a special case of time-frequency masking (common in speech processing) combined with the signal model from local stacking. This powerful combination opens the avenue for a completely new family of approaches for multi-channel seismic processing that can address seismic processing of land data with nodes and single sensors in the desert environment.     

Seismic characterization of vertical fractures described as general linear-slip interfaces

Fluid flow in many hydrocarbon reservoirs is controlled by aligned fractures which make the medium anisotropic on the scale of seismic wavelength. Applying the linear‐slip theory, we investigate seismic signatures of the effective medium produced by a single set of ‘general’ vertical fractures embedded in a purely isotropic host rock. The generality of our fracture model means the allowance for coupling between the normal (to the fracture plane) stress and the tangential jump in displacement (and vice versa). Despite its low (triclinic) symmetry, the medium is described by just nine independent effective parameters and possesses several distinct features which help to identify the physical model and estimate the fracture compliances and background velocities. For example, the polarization vector of the vertically propagating fast shear wave S₁ and the semi‐major axis of the S₁‐wave normal‐moveout (NMO) ellipse from a horizontal reflector always point in the direction of the fracture strike. Moreover, for the S₁‐wave both the vertical velocity and the NMO velocity along the fractures are equal to the shear‐wave velocity in the host rock. Analysis of seismic signatures in the limit of small fracture weaknesses allows us to select the input data needed for unambiguous fracture characterization. The fracture and background parameters can be estimated using the NMO ellipses from horizontal reflectors and vertical velocities of P‐waves and two split S‐waves, combined with a portion of the P‐wave slowness surface reconstructed from multi‐azimuth walkaway vertical seismic profiling (VSP) data. The stability of the parameter‐estimation procedure is verified by performing non‐linear inversion based on the exact equations.     

S 235 B explained: an accreting Herbig Be star surrounded by reflection nebulosity

The intent of this study is to determine the nature of the star and associated nebulosity S 235 B, which are located in a region of active star formation still heavily obscured by the parent molecular cloud. Low-resolution  ( R = 400)  long-slit spectra of the star and nebulosity, and medium-  ( R = 1800)  and high-resolution  ( R = 60 000)  spectra of the central star are presented along with the results of Fabry–Perot interferometric imaging of the entire region. Based on the long-slit and Fabry–Perot observations, the nebulosity appears to be entirely reflective in nature, with the stellar component S 235 B★ providing most of the illuminating flux. The stellar source itself is classified here as a B1V star, with emission-line profiles indicative of an accretion disc. S 235 B★ thus belongs to the relatively rare class of early-type Hebrig Be stars. Based on the intensity of the reflected component, it is concluded that the accretion disc must be viewed nearly edge-on. Estimates of the accretion rate of S 235 B★ from the width of the Hα profile at 10 per cent of maximum intensity, a method which has been used lately for T Tauri stars and Brown Dwarfs, appear to be inconsistent with the mass outflow rate and accretion rate implied from previous infrared observations by Felli et al., suggesting this empirical law does not extend to higher masses.     

Geodetic observations of ice flow velocities over the southern part of subglacial Lake Vostok, Antarctica, and their glaciological implications

In the austral summer seasons 2001/02 and 2002/03, Global Positioning System (GPS) data were collected in the vicinity of Vostok Station to determine ice flow velocities over Lake Vostok. Ten GPS sites are located within a radius of 30 km around Vostok Station on floating ice as well as on grounded ice to the east and to the west of the lake. Additionally, a local deformation network around the ice core drilling site 5G-1 was installed.
The derived ice flow velocity for Vostok Station is  2.00 m a⁻¹± 0.01 m a⁻¹  . Along the flowline of Vostok Station an extension rate of about 10⁻⁵ a⁻¹ (equivalent to 1 cm km⁻¹ a⁻¹) was determined. This significant velocity gradient results in a new estimate of 28 700 years for the transit time of an ice particle along the Vostok flowline from the bedrock ridge in the southwest of the lake to the eastern shoreline. With these lower velocities compared to earlier studies and, hence, larger transit times the basal accretion rate is estimated to be 4 mm a⁻¹ along a portion of the Vostok flowline. An assessment of the local accretion rate at Vostok Station using the observed geodetic quantities yields an accretion rate in the same order of magnitude. Furthermore, the comparison of our geodetic observations with results inferred from ice-penetrating radar data indicates that the ice flow may not have changed significantly for several thousand years.     

中国东北兴安岭地区年冻融频次的分布规律北大核心CSCD

冻融循环是影响寒区工程和环境变化的重要因素之一,而年冻融频次(冻融循环的年累计量)是可以直接标定冻融循环对寒区工程和环境造成影响大小的物理量,对寒区环境下年冻融频次的调查有重要意义。以中国东北兴安岭地区为研究对象,选取兴安岭地区18个气象站点1990—2017年地表温度数据进行整理,得出兴安岭地区各个站点的年冻融频次,并在空间分布和时间序列上分析了年冻融频次的变化情况。结果表明:年冻融频次在中国东北兴安岭地区随着纬度增加逐渐增多,随着经度的增加逐渐减少,随着海拔升高而逐渐增多,且均具有良好的线性趋势;年冻融频次在1990—2017年呈明显的减少趋势,在2004年前后发生突变。进一步分析发现,在影响年冻融频次空间分布的3个因素中,经纬度影响较大,海拔影响相对较小;在时间上年冻融频次变化与东北地区NDVI变化情况呈现良好的负相关性,表明年冻融频次在一定程度上与当地植被覆盖度有关,且随植被覆盖度的增加而减少。     

On the effects of resolution in dissipationless cosmological simulations

We present a study of numerical effects in dissipationless cosmological simulations. The numerical effects are evaluated and studied by comparing the results of a series of 64³-particle simulations of varying force resolution and number of time-steps, performed using three of the N -body techniques currently used for cosmological simulations: the Particle–Mesh (PM), the Adaptive Particle–Particle–Particle–Mesh (AP³M) and the newer Adaptive Refinement Tree (ART) codes. This study can therefore be interesting both as an analysis of numerical effects and as a systematic comparison of different codes.
We find that the AP³M and the ART codes produce similar results given that convergence is reached within the code type. We also find that numerical effects may affect the high-resolution simulations in ways that have not been discussed before. In particular, our study revealed the presence of two-body scattering, the effects of which can be greatly amplified by inaccuracies in time integration. This process appears to affect the correlation function of matter, the mass function, the inner density of dark matter haloes and other statistics at scales much larger than the force resolution, although different statistics may be affected in a different fashion. We discuss the conditions for which strong two-body scattering is possible and discuss the choice of the force resolution and integration time-step. Furthermore, we discuss recent claims that simulations with force softening smaller than the mean interparticle separation are not trustworthy and argue that this claim is incorrect in general, and applies only to the phase-sensitive statistics. Our conclusion is that, depending on the choice of mass and force resolution and the integration time-step, a force resolution as small as 0.01 of the mean interparticle separation can be justified.