饱和流体多孔介质中AVA特征分析及储层参数反演研究

AVA(Amplitude Versus Angle)技术是基于常规介质模型(均匀各向同性介质模型)发展起来的,由于忽略了储层的孔隙结构和充填流体的影响,造成AVA特征中是部分对地震波能量的吸收和衰减作用反映不足.从振幅特征方程的实际应用出发,建立起各参数与常规岩性参数之间的关系.以Gassmann方程与Biot理论为基础,推导出振幅特征方程中弹性参数与常规岩性、储层参数(如:纵波速度、横波速度、密度、孔隙度、流体参数Kf)之间的转化关系式,详细研究饱和流体多孔介质模型中的AVA特征,比较该介质模型与常规介质模型在AVA特征上的差异,并将饱和流体多孔介质AVA技术应用于川西凹陷深层须家河组储层预测,通过多波AVA储层参数的反演研究,为直接利用地震资料进行储层识别并进一步识别其流体特征提供了一种有力手段.     

Absorbing boundary conditions for dynamic analysis of fluid-saturated porous media

Based on Biot's two-phase mixture theory and the paraxial approximation, the absorbing boundary conditions in the time domain for u-w, u-U and u-p formulations are presented for the dynamic analysis of fluid-saturated porous media. These absorbing boundaries are equivalent to the viscous boundaries in the fundamental mode. The expressions for the energy ratio and reflection coefficient between the reflected and incident waves along the absorbing boundary are given. The numerical results show that the proposed absorbing boundaries can greatly suppress the spuriously reflected wave and model the far field of the system. These results also dynamic analysis of infinite fluid-saturated porous media. for the transient dynamic analysis of infinite fluid-saturated porous media.     

二维非均匀饱和土体的地震反应分析

将饱和土体视为由弹性骨架和不可压缩流体组成的两相孔隙介质,取固相位移,液相位移,孔隙水压力为场变量。从两相孔介质的动力耦合微分方程出发,利用伽辽金原理和Ｗｉｌｓｏｎ－θ法导出了三场有限元积分方法,该方法可进行二维非均匀饱和土体的地震反应分析,文中通过算例讨论了非均匀饱和土体地震反应的一些特征。     

层状饱和场地中圆柱形基础的阻抗函数

本文利用间接边界元法,建立层状饱和场地中三维土-基础动力相互作用模型,研究了孔隙中流体的存在、场地动力特性和圆柱形基础的高度对基础阻抗函数的影响。数值结果表明,孔隙中流体的存在对基础阻抗函数有明显的影响,尤其是水平和竖向分量;场地动力特性也对基础阻抗函数有明显的影响,层状饱和场地中阻抗函数围绕相应均匀半空间结果振荡,且基岩与土层刚度比越大或土层厚度越小,振荡的幅度越大,土层厚度越小,振荡的周期越小;另外,基础高度也对阻抗函数有明显的影响,随基础高度增大,阻抗函数逐渐增大。     

流体饱和两相多孔介质动力反应计算分析

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种情况下计算得到的两相介质弹性动力反应结果的差异进行对比研究,从而揭示孔隙流体渗流对两相介质动力反应性质的影响。计算结果表明：两相介质弹性动力反应时程的波形与入射地震波的波形相同,且弹性动力反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻;孔隙流体的渗流将对两相介质的弹性动力反应性质产生显著的影响。数值计算同时表明,时域显式有限元方法是进行流体饱和两相多孔介质弹性动力反应计算分析的一种有效的方法。     

Hazards from pyroclastic density currents at Mt. Etna (Italy)

Despite the recent recognition of Mount Etna as a periodically violently explosive volcano, the hazards from various types of pyroclastic density currents (PDCs) have until now received virtually no attention at this volcano. Large-scale pyroclastic flows last occurred during the caldera-forming Ellittico eruptions, 15–16 ka ago, and the risk of them occurring in the near future is negligible. However, minor PDCs can affect much of the summit area and portions of the upper flanks of the volcano. During the past ~ 20 years, small pyroclastic flows or base-surge-like vapor and ash clouds have occurred in at least 8 cases during summit eruptions of Etna. Four different mechanisms of PDC generation have been identified during these events: (1) collapse of pyroclastic fountains (as in 2000 and possibly in 1986); (2) phreatomagmatic explosions resulting from mixing of lava with wet rock (2006); (3) phreatomagmatic explosions resulting from mixing of lava with thick snow (2007); (4) disintegration of the unstable flanks of a lava dome-like structure growing over the rim of one of the summit craters (1999). All of these recent PDCs were of a rather minor extent (maximum runout lengths were about 1.5 km in November 2006 and March 2007) and thus they represented no threat for populated areas and human property around the volcano. Yet, events of this type pose a significant threat to the lives of people visiting the summit area of Etna, and areas in a radius of 2 km from the summit craters should be off-limits anytime an event capable of producing similar PDCs occurs. The most likely source of further PDCs in the near future is the Southeast Crater, the youngest, most active and most unstable of the four summit craters of Etna, where 6 of the 8 documented recent PDCs originated. It is likely that similar hazards exist in a number of volcanic settings elsewhere, especially at snow- or glacier-covered volcanoes and on volcano slopes strongly affected by hydrothermal alteration.     

In-plane soil–structure interaction in layered,fluid-saturated,poroelastic half-space I: Structural response

Linear in-plane soil–structure interaction in two dimensions (2D) is studied in fluid-saturated, poroelastic, layered half-space using the Indirect Boundary Element Method (IBEM). The structure is a shear wall supported by a rigid embedded foundation. Exact stiffness matrices for the soil layer and half-space, and Green׳s functions of uniformly distributed loads and pore pressure on an inclined line are derived. Results of the system response in the frequency domain are presented for the special case of single soil layer over bedrock, semi-circular foundation and zero seepage force. The effects of water saturation, soil porosity, depth of soil layer, rigidity contrast between layer and bedrock are investigated in the frequency domain for incident plane P- and SV waves. The results suggest that water saturation may cause increase of the system frequency by more than 10%.     

震电效应在油气勘探开发中的应用

震电效应是流体饱和孔隙介质中客观存在的一种物理现象。开展震电勘探研究是当今地球物理领域的前沿研究课题。笔者基于当前国内外大量室内和野外试验资料，对震电效应及其产生机理、探测方式、勘探应用潜力进行了分析探讨，指出震电效应应用于油气勘探开发，特别是在寻找“剩余油”方面具有广阔的应用前景。     

基于流体饱和多孔隙介质波动方程的时间推移地震反演

针对二维流体饱和多孔隙介质波动方程,设计了自适应同伦法和小波-自适应同伦法,在此基础上,将局域化思想引入时间推移地震反演,并使用小波-自适应同伦法和自适应同伦法反演两次监测数据,提出了一种适用于时间推移地震勘探的局域化反演算法.为了验证新方法是可行的,文中把它应用于同一地点不同时间的观测模型,通过数值实验表明了这种方法的有效性.     

PML Absorbing Boundary Condition for Seismic Numerical Modeling by Convolutional Differentiator in Fluid-Saturated Porous Media

The perfectly matched layer(PML) was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator(FPCD) algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot's equations is gi...