辽河油田D块中—厚层稠油油藏隔夹层的分类及识别

充分利用研究区Ng油层、XLT油层取芯井的岩芯资料和测井等资料,识别并划分出不同类型的隔夹层。研究区主要发育泥质隔夹层、物性隔夹层和砾岩隔夹层,以泥质隔夹层和物性隔夹层为主。从定性和定量两方面采用多学科、多信息对该区的隔夹层进行精细研究,建立隔夹层识别的定性、定量识别标准,为隔夹层的精细刻画奠定基础,为提高稠油热采开发效果提供地质依据     

中地壳的水和水岩相互作用实验及其地球物理涵义

本文重点报道了高温高压下流体与流体-岩石相互作用实验结果,提供了中地壳条件下流体性质和水岩反应速率数据.这些数据有助于理解中地壳的一些地球物理现象.作者进行了25℃~435℃和22~39 MPa条件下水-岩相互作用反应动力学实验.同时,研究水在近临界区至超临界区的性质.一般地说,中地壳大致位于10(15)至25 km的深度范围.各地的地壳厚度不同,但是中地壳高导-低速层的深度范围十分相似.中地壳的顶界温度处于300℃,底界大致为450℃范围,压力高达200 MPa以上.流体-岩石相互作用实验表明:硅酸盐矿物和岩石的硅最大溶解速率出现在300℃~400℃.此时,硅酸盐矿物格架解体.通常,地壳里普遍存在水、流体.地壳构造活动导致断裂空隙、减压、流体流动.这时,有可能导致中地壳处于300℃~450℃流体的压力减低,由超临界区进入临界态、亚临界态.这会引发强烈流动的水与岩石相互作用.溶解反应导致岩层的硅淋失,硅的强烈淋失又会导致硅酸盐矿物格架解体,岩石崩塌.同时,进一步促进流体的流动.实验表明300℃~400℃下的强烈水岩相互作用促进了岩石破坏,并有可能影响岩层的地球物理性质,如高导层出现.另外,实验和理论研究表明处于300℃~400℃流体具有高电导率性质.这些水岩相互作用会使中地壳出现高导-低速层.     

Dissecting large earthquakes in Japan: Role of arc magma and fluids

We synthesized information from recent high-resolution tomographic studies of large crustal earthquakes which occurred in the Japanese Islands during 1995–2008. Prominent anomalies of low-velocity and high Poisson's ratio are revealed in the crust and uppermost mantle beneath the mainshock hypocenters, which may reflect arc magma and fluids that are produced by a combination of subducting slab dehydration and corner flow in the mantle wedge. Distribution of 164 crustal earthquakes ( M 5.7–8.0) that occurred in Japan during 1885–2008 also shows a correlation with the distribution of low-velocity zones in the crust and uppermost mantle. A qualitative model is proposed to explain the geophysical observations recorded so far in Japan. We consider that the nucleation of a large earthquake is not entirely a mechanical process, but is closely related to the subduction dynamics and physical and chemical properties of materials in the crust and upper mantle; in particular, the arc magma and fluids.     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

Disruption rates for one vulnerable soil in Organ Pipe Cactus National Monument,Arizona, USA

Rates of soil disruption from hikers and vehicle traffic are poorly known, particularly for arid landscapes. We conducted an experiment in Organ Pipe Cactus National Monument (ORPI) in western Arizona, USA, on an air-dry very fine sandy loam that is considered to be vulnerable to disruption. We created variable-pass tracks using hikers, an all-terrain vehicle (ATV), and a four-wheel drive vehicle (4WD) and measured changes in cross-track topography, penetration depth, and bulk density. Hikers (one pass = 5 hikers) increased bulk density and altered penetration depth but caused minimal surface disruption up to 100 passes; a minimum of 10 passes were required to overcome surface strength of this dry soil. Both ATV and 4WD traffic significantly disrupted the soil with one pass, creating deep ruts with increasing passes that rendered the 4WD trail impassable after 20 passes. Despite considerable soil loosening (dilation), bulk density increased in the vehicle trails, and lateral displacement created berms of loosened soil. This soil type, when dry, can sustain up to 10 passes of hikers but only one vehicle pass before significant soil disruption occurs; greater disruption is expected when soils are wet. Bulk density increased logarithmically with applied pressure from hikers, ATV, and 4WD.     

Earthquake response analyses of a full‐scale five‐story steel frame equipped with two types of dampers

The seismic performance tests of a full‐scale five‐story passively controlled steel building were conducted on the E‐Defense shaking table in Japan in March 2009. Before the tests, a blind prediction contest was held to allow researchers and practitioners from all over the world to construct analytical models and predict the dynamic responses of the steel frame specimen equipped with buckling‐restrained braces (BRBs) or viscous dampers (VDs). This paper presents the details of two refined prediction models made and results obtained before the tests. When the proposed analytical modeling techniques are adopted as in the two refined prediction models, the overall prediction accuracy is about 90%. Sensitivity studies conducted after the tests are also presented in this paper. The effects of varying each modeling feature on the response simulation accuracy have been investigated. The analytical results suggest that considering concrete full‐composite actions for beam members could improve prediction accuracy by about 20% against using the simplified bare steel beam model. Adopting refined BRB stiffness computed from incorporating finite‐element gusset stiffness only improves the overall prediction accuracy by 0.9%. Considering the BRB dynamic loading test results for analytical BRB strength reduces the error by 1.9%. For the VD frame, incorporating the brace and VD stiffness could improve the overall prediction accuracy by about 15%. Copyright © 2012 John Wiley & Sons, Ltd.     

The instability of barotropic circular vortices

Abstract

The linear, normal mode instability of barotropic circular vortices with zero circulation is examined in the f-plane quasigeostrophic equations. Equivalents of Rayleigh's and Fjortoft's criteria and the semicircle theorem for parallel shear flow are given, and the energy equation shows the instability to be barotropic. A new result is that the fastest growing perturbation is often an internal instability, having a finite vertical scale, but may also be an external instability, having no vertical structure. For parallel shear flow the fastest growing perturbation is always an external instability; this is Squire's theorem. Whether the fastest growing perturbation is internal or external depends upon the profile: for mean flow streamfunction profiles which monotonically decrease with radius, the instability is internal for less steep profiles with a broad velocity extremum and external for steep profiles with a narrow velocity extremum. Finite amplitude, numerical model calculations show that this linear instability analysis is not valid very far into the finite amplitude range, and that a barotropic vortex, whose fastest growing perturbation is internal, is vertically fragmented by the instability.     

Penetrative convection in rotating fluids: A model for the base of stellar convection zones

Abstract

To model penetrative convection at the base of a stellar convection zone we consider two plane parallel, co-rotating Boussinesq layers coupled at their fluid interface. The system is such that the upper layer is unstable to convection while the lower is stable. Following the method of Kondo and Unno (1982, 1983) we calculate critical Rayleigh numbers R_c for a wide class of parameters. Here, R_c is typically much less than in the case of a single layer, although the scaling R_c～T^2/3 as T → ∞ still holds, where T is the usual Taylor number. With parameters relevant to the Sun the helicity profile is discontinuous at the interface, and dominated by a large peak in a thin boundary layer beneath the convecting region. In reality the distribution is continuous, but the sharp transition associated with a rapid decline in the effective viscosity in the overshoot region is approximated by a discontinuity here. This source of helicity and its relation to an alpha effect in a mean-field dynamo is especially relevant since it is a generally held view that the overshoot region is the location of magnetic field generation in the Sun.     

A local stability analysis of astronomical disks

Abstract

A fifth-order dispersion relation describing the local stability of a differentially rotating flow against small perturbations is derived. Finite viscosity and conductivity and both vertical (parallel to the rotation axis) and radial gradients in density, temperature and pressure are included. A general form is assumed for the equation of state, although this is not exploited in the paper. A number of special cases are studied: with negligible viscosity and conductivity, it is shown that modes can often be separated into two high frequency (modified acoustic), two intermediate frequency (combined inertial and internal waves) and a low frequency mode. In convectively unstable situations the intermediate frequency modes may be replaced by a damped/growing pair of instablities. Various criteria for mode excitation are given. It is shown that viscosity always inhibits instability at very short wavelengths, while non-zero conductivity may destabilize the flow. At intermediate wavelengths viscosity could also play a destabilizing role. A parameter study of the effects of fluctuations in the conductivity shows that it could cause mode excitation under certain circumstances.     

Large wavenumber convection in the rotating annulus

Abstract

The annulus model considers convection between concentric cylinders with sloping endwalls. It is used as a simplified model of convection in a rapidly rotating sphere. Large azimuthal wavenumbers are preferred in this problem, and this has been exploited to develop an asymptotic approach to nonlinear convection in the annulus. The problem is further reduced because the Taylor-Proudman constraint simplifies the dependence in the direction of the rotation vector, so that a nonlinear system dependent only on the radial variable and time results. As Rayleigh number is increased a sequence of bifurcations is found, from steady solutions to periodic solutions and 2-tori, typically ending in chaotic behaviour. Both the magnetic (MHD convection) and non-magnetic problem has been considered, and in the non-magnetic case our bifurcation sequence can be compared with those found by previous two-dimensional numerical simulations.