基于地震波振幅分辨塔河油田溶洞最小高度的理论探讨

塔河油田奥陶系灰岩溶洞反射特征表现为"串珠状"强振幅短同相轴反射结构,但该串珠状反射的地震属性在空间展布的视体积与其对应的实际体积差异很大,因此,需要确定地震波振幅分辨塔河油田溶洞最小高度.本文在分析塔河油田钻遇储层段储层厚度特征的基础上,分别导出洞穴内一次反射渡和含一次多次波的时域反射脉冲响应的振幅谱关于洞穴高度的解析表达式,进而对不同洞穴高度、不同类型的充填物时洞穴反射振幅谱特征以及相位关于频率的变化率的变化特征进行分析,结果表明,谱出现极大值相应的洞穴高度在相位关于频率变化率曲线上对应的洞穴高度处也出现极大,可以利用相位关于频率的变化率来检测溶洞内是否有充填物.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.     

粘弹性介质地震波传播的褶积微分算子法数值模拟研究

早期的褶积微分算子法都是基于正反傅立叶变换而实现的,其精度比四阶有限差分稍高。本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与l6阶有限差分相当。粘弹性波动方程更真实地描述了实际地下介质中弹性波的传播规律及其波场特征。本文以二维粘弹性波动方程为例,推导了粘弹性介质波动方程的离散格式,用迭积微分算子法实现了粘弹性介质的地震波场正演模拟,并对其波传播特征进行了分析。计算结果表明该算法能正确模拟粘弹性介质中的地震波,正确地反映粘弹性介质中波场的传播规律。     

基于拟微分算子的横向各向同性介质黏声波方程

由于地层的各向异性与黏滞性,地震波在传播过程中表现出方向上的各向异性和振幅值的吸收衰减特征.若简单地将地下介质视为各向同性,实际地质结构的地震波响应可能会被曲解,甚至忽略了有效信息.为了得到深部储层的高精度成像,本文将标准线性固体模型(SLS)扩展到各向异性介质中,并结合各向异性拟声波方程,推导出时间空间域的各向异性介质黏声波方程组,该方程中用伪微分算子来表述地震波的衰减特征.在逆时偏移的过程中,为了解决高频不稳定问题,文中通过引入规则化算子来构建稳定的逆时传播算子.模型的正演模拟表明该方程既能准确的描述各向异性介质下的标量波的传播规律,又能体现地层的吸收衰减效应.     

各向异性介质地震波场的优化褶积微分算子法数值模拟

在前人工作基础上,通过对窗函数参数进行优化实现了对基于Shannon奇异核理论的交错网格褶积微分算子的优化过程.应用这种优化褶积微分算子方法对各向异性介质进行了数值模拟,讨论了优化褶积微分算子法模拟的PML吸收边界条件以及稳定性条件,分析了弹性波在此类介质中的传播特征,并与高阶交错网格有限差分方法进行了对比.数值实验结果表明,该方法适用于各向异性介质中弹性波场模拟,精度高,稳定性好,是一种研究复杂介质中地震波传播的有效数值方法.     

基于广义正交多项式褶积微分算子的地震波场数值模拟方法

地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性.     

深层碳酸盐岩岩溶风化壳溶洞型油气藏地震波场的理论研究

碳酸盐岩洞缝型油气藏是我国乃至世界的重要油气藏类型之一，属典型的复杂非均匀介质。其中，岩溶洞穴型油气藏是最容易发现的。为了解决这类储层的勘探开发及预测问题，需要从理论上对它们的地震波场进行研究。我们利用随机介质模型和非均匀介质弹性波动方程计算了不同高度和宽度的充满流体藏层的合成地震剖面。结果表明当绕射波偏移为珠子状的同相轴时即使远小于四分之一波长的溶洞在常规带宽的地震剖面上也是可分辨的，绕射波是溶洞高度和宽度的函数。我们引入一个宽度-振幅因子。可以用这个宽度-振幅因子从对无限宽溶洞计算的绕射波计算有限宽度的溶洞的绕射振幅。     

油水两相渗流多孔介质弹性波传播动力学模型研究

弹性波在储层渗流场中的传播与衰减规律是研究波场强化采油动力学机理的重要基础.基于等效流体理论和饱和静态流体弹性波传播Biot理论,建立油水两非混相流体渗流条件下储层多孔介质中弹性波传播的动力学模型,通过算例求解与分析,发现含油水两相渗流储层多孔介质中同时存在着3种纵波P1、P2、P3和1种横波S;受频率和含水饱和度的影响,各波波速和品质因子呈现出不同变化规律,4种体波波速与频率、饱和度正相关,P1、P2波品质因子与饱和度正相关,P3和S波品质因子与饱和度负相关;最后,通过与传统静态弹性波模型结果对比,进一步分析了宏观渗流场对弹性波传播特征的影响规律,为揭示低频人工地震波辅助强化采油技术的动力学机理和工艺参数优化提供了重要理论依据.

     

基于NAD算子和FCT技术改进的Runge-Kutta方法及其波场模拟

地震波正演数值模拟中所产生的数值频散严重地影响波场模拟的精度和分辨率.基于声波方程,以8阶精度的近似解析离散算子来校正差分算子对空间偏导数进行离散,以三阶Runge-Kutta方法对时间导数进行离散,并结合通量校正传输技术以抑制数值频散,从而建立基于NAD算子和TCF技术改进的Runge-Kutta方法.波场模拟结果表明,该方法在抑制数值频散上明显优越于传统的高阶有限差分方法.同时,该方法对模拟地震波在复杂介质中的传播有着很强的适应能力.     

Carbon dioxide in cave atmospheres. New results in belgium and comparison with some other countries

More than 600 measurements of the carbon dioxide content of cave air in Belgium lead up to the conclusion that the main factors of its distribution are: (1) a flow originating from the biomass and diffusing in the soil and the voids of bedrock; (2) a trend, discernible in very still air only, to go down by density; (3) in some caves, draughts caused, for instance, by a swift underground stream. Results in Belgium are compared with published and unpublished data from other countries, showing that CO₂ is often less abundant in cold climate caves and in caves of semi-arid regions (influence of the biomass). Special attention is paid to human contamination during analyses: the influence of people passing through the cave nearby the operator, but also the influence of the operator himself are discussed, and the use of special precautions (including a CO₂-absorbing mask) in defined critical situations is stressed.     

碳酸盐岩孔洞型储集体地震反射特征分析

定量分析碳酸盐岩孔洞型储集体的地震响应特征,对孔洞储集体的预测和刻画研究有着重要意义.基于声波方程,在忽略多次散射的情况下,导出了深层小尺度孤立洞体的地震响应近似解析式,该解析表达式给出了孔洞尺度、充填程度与地震振幅属性之间的关系.分析表明:(1)当孔洞异常体尺度小于地震波长的1/8时,反射波振幅与异常体体积呈近似正比关系;(2)异常体横向上的尺度对反射波振幅的影响主要由菲涅耳半径控制,同时振幅随横向尺度的变化规律受到纵向尺度的影响;(3)异常体纵向上的尺度对反射振幅的影响主要表现为调谐效应,当纵向上的尺度等于调谐厚度(地震波长的1/4)时,反射振幅最强.不同尺度溶洞的数值模拟研究结果以及不同填充物的单洞物理模型试验和溶洞群物理模型试验验证了以上结论的正确性.     

窑洞外形特征对结构抗震性能影响研究

调查发现山西北部农村地区农居窑洞多是独立式土坯窑洞,结构主要受力体系是土坯砌体砌筑拱结构,其建造无统一标准,外形样式存在差异。由于拱结构的内力分布对外形尺寸变化比较敏感,为研究窑洞外形特征对结构抗震性能影响,建立相应的有限元模型进行计算分析,对比不同外形特征窑洞在地震作用下的损伤情况,并结合历史震害,指出窑洞结构的薄弱位置。研究结果表明覆土较厚的窑洞容易损伤,连拱孔数较多窑洞的两侧位置更容易发生破坏,拱的曲线形状对结构抗震性能有一定影响,合理拱曲线窑洞的抗震性能优于一般形状窑洞。     

碳酸盐岩洞穴型储层双侧向测井数值模拟和响应特征(英文)

洞穴型碳酸盐岩地层具有显著的非均质性和测井预测难度大的特点,如何利用双侧向测井对识洞穴进行准确识别与定量评价,是当前的一个研究难点与热点。利用数值模拟方法计算电测井响应,可为洞穴识别与评价提供理论依据。本文根据双侧向测井原理,采用有限元素法,研究了球状洞穴模型的测井响应特征,并考察了不同尺寸的洞穴对电测井响应的影响;然后,系统研究了不同充填程度的洞穴模型的测井响应,得出了电阻率与充填程度之间的关系;最后,研究了洞穴中充填泥质、灰质、砾岩和砂泥岩薄互层等多种充填物的测井响应特征。研究结果对碳酸盐岩洞穴型储层的识别与评价具有较高的应用价值和研究意义。     

碳酸盐岩溶洞成像要素分析与研究

碳酸盐岩油气藏储集空间以裂缝和溶洞为主,在地震剖面上主要表现为“串珠状”。为了研究溶洞成像规律,本文通过大量的数值模拟,分析溶洞的成像分辨率,探讨溶洞成像与上覆地层、信噪比(SNR)、覆盖次数、偏移速度和偏移方法之间的定量关系。研究认为:①无论是溶洞的纵向分辨率还是横向分辨率,均满足1/4波长调谐规律,可以将薄层机理转用至溶洞分辨率解释;②溶洞的存在会对上覆地层产生下拉效应,使溶洞在地震剖面上呈塌陷等“非串珠状”特征;③高覆盖和小道距(面元)可提高串珠能量,压制偏移附带的背景噪声;④叠前偏移能改善资料的SNR,提高成像精度,特别是低SNR资料;⑤偏移速度偏小,串珠边缘下拉,反之会上翘,速度偏大带来的成像质量变差的程度要大于速度偏小;⑥不同偏移方法对缝洞成像效果存在差异,对于低SNR资料,溶洞成像拟选取频率域或频波域方法。     

基于3D-FEM的洞穴型地层双侧向测井数值模拟及响应特征研究

洞穴型碳酸盐岩储层非均质性强、电测井响应复杂、测井识别和表征难度大,利用数值模拟方法明确该类储层的电测井响应特征,可为洞穴识别与评价提供理论依据.本文基于三维有限元素法（3D-FEM）,引入边界局部加密技术,实现对复杂球形洞穴的精细刻画与响应精确计算;对比分析井眼钻穿型洞穴和井旁洞穴的双侧向测井响应特征和敏感性,进而考察双侧向测井对两者敏感范围的差异.结果表明：洞穴的存在导致双侧向测井响应明显降低,受洞穴边界及仪器探测深度等影响,井眼钻穿型洞穴双侧向测井曲线复杂,而井旁洞穴曲线呈“抛物线”型;双侧向测井对井眼钻穿型洞穴的敏感性远大于对井旁洞穴的敏感性,深侧向和浅侧向测井最大可对洞穴边界距井壁0.5 m和0.3 m的井旁洞穴敏感.

     

Tectonic uplift,sea level changes and Plio‐Pleistocene evolution of a coastal karst system: the Mount Saint Paul (Palawan,Philippines)

The St. Paul karst (Palawan, Philippines) is a tropical coastal karst, consisting of towers, cones, huge depressions and large caves. This area hosts the Puerto Princesa Subterranean River (PPSR, 24 km long), whose main entrance is a large spring along the coast and which is one of the largest cave complexes in eastern Asia. A geomorphological study performed by several field surveys and a morphometric analysis of the digital terrain model (DTM) and 3D cave models, allowed formulation of a first evolutionary framework of the karst system. The DTM was extracted from maps and aerial photographs in order to find different generations of ‘relict’ landforms, through the morphometric analysis of topographic surface and karst landforms. Several features suggest a long and multi‐stage evolution of the karst, whose age ranges from Pliocene to present. The southern and northern sectors of the area differ in their altimetric distribution of caves. In the southern sector, some large caves lie between 300 and 400 m asl and were part of an ancient system that developed at the base level of a past river network. In the northern sector, some mainly vadose caves occur, with a phreatic level at 120–130 m asl. An important phase of base‐level cave development is well documented in the inactive passages of PPSR at 50–80 m asl. Morphological features, such as horizontal solution passages and terraced deposits, suggest a phase of stillstand of the base level, which is recorded in the topography as low‐relief surfaces at 40–50 m asl. The age of this phase is probably Early Pleistocene, on the basis of assumed uplift rates. The more recent caves are still active, being located at the current sea level, but they show more than one cycle of flooding and dewatering (with calcite deposition). In the PPSR, several morphologic features, such as two main water level notches at +12·4 and +7·7 m asl and terraced alluvial deposits, suggest that the lower and active level passed through more than two high‐stands of sea level and so it could have formed throughout most of the Middle‐Late Pleistocene. Copyright © 2010 John Wiley & Sons, Ltd.     

The role of GPR techniques in determining ice cave properties: Peña Castil ice cave,Picos de Europa

The structure and ice content of ice caves are poorly understood. Ground penetrating radar (GPR) can provide useful insights but has only rarely been applied to ice caves. This paper interprets GPR images (radargrams) in terms of internal structure, stratification, compaction, thickness and volume of the ice block in the Peña Castil ice cave (Central Massif of Picos de Europa, northern Spain), providing the endokarst geometry of the ice cave in GPR data reflections. Eight radargrams were obtained by applying a shielded ground‐coupled antenna with a nominal frequency of 400 MHz. Although the radargrams do not depict the ice–basal bedrock interface, they suggest that the ice block is at least 54 m deep and similarly thick. Some curved reflection signatures suggest a potential vertical displacement in the block of ice, and thus certain dynamics in the ice body. Other images show numerous interbedded clasts and thin sediment layers imaged as banded reflections. In this particular cave a direct visual inspection of the ice stratigraphy is a difficult task but GPR provides clear reflectivity patterns of some of its internal features, making GPR a suitable instrument for this and future studies to achieve a better and broader understanding of the internal behavior of ice caves.     

Speleoseismology: A critical perspective

Speleoseismology is the investigation of earthquake records in caves. Traces can be seen in broken speleothems, growth anomalies in speleothems, cave sediment deformation structures, displacements along fractures and bedding plane slip, incasion (rock fall) and co-seismic fault displacements. Where earthquake origins can be proven, these traces constitute important archives of local and even regional earthquake activity. However, other processes that can generate the same or very similar deformation features have to be excluded before cave damage can be interpreted as earthquake induced. Most sensitive and therefore most valuable for the tracing of strong earthquake shocks in caves are long and slender speleothems, such as soda straws, and deposits of well-bedded, water-saturated silty sand infillings, particularly in caves close to the earth's surface. Less easily proven is a co-seismic origin of an incasion and other forms of cave damage. The loads and creep movements of sediment and ice fillings in caves can cause severe damage to speleothems which have been frequently misinterpreted as evidence of earthquakes. For the dating of events in geological archives, it is important to demonstrate that such events happened at approximately the same time, i.e. within the error bars of the dating methods. A robust earthquake explanation for cave damage can only be achieved by the adoption of appropriate methods of direct dating of deformation events in cave archives combined with correlation of events in other geological archives outside caves, such as the deformation of lake and flood-plain deposits, locations of rock falls and active fault displacements.