简介:本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger'sT2公式(即SDR模型),二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。
简介:针对祁连山冻土区DK-4井孔含水合物岩层构建岩石物理模型,分别采用K-T方程模型方法和区分填充模式的等效介质模型方法(模式I和模式II)。K-T方程主要模拟地震波在两相介质中传播,基于弹性模量计算速度;而等效介质模型主要依据对水合物地层介质的两种假设:模式I是将水合物作为孔隙填充物的一部分,模式II是将水合物作为岩石骨架的一部分。首先根据粉砂岩层段的测井数据提取了水合物地层骨架的物性参数,包括纵波速度、横波速度、密度、体积模量和剪切模量。然后依据水合物地层各主要成分的物性参数,建立了基于K-T方程的岩石物理模型和区分填充模式的等效介质模型。将两类模型的速度曲线分别与实际地层数据进行了对比:由K-T方程建立的岩石物理模型,其理论计算的速度偏离实际值;而区分填充模式的等效介质模型,其理论计算的速度符合实际值,并且填充模式II模型的速度曲线比填充模式I模型更接近实际地层情况。采用区分填充模式的等效介质模型,能够更好地实现对祁连山冻土区水合物粉砂岩地层的模拟。
简介:根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震(OBS)共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。
简介:天然气水合物沉积层在常规叠后地震剖面上,一般表现为似海底反射(BSR)特征,常作为水合物识别的重要标志,但由于地震数据分辨率的限制与多解性的存在,类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手,综述研究区域天然气水合物发育的构造、沉积与运移环境;其次以加权方程与Zoeppritz公式的正演模拟为基础,讨论天然气水合物与BSR特征的关系;然后选取经过三口已钻获天然气水合物站位的地震测线,提取角度域共成像道集用于BSR的AVA特征研究;最后依据叠前地震弹性参数同步反演思路,对水合物重点赋存区域开展精细预测研究,获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。
简介:本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.
简介:本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。
简介:这篇文章重点研究改进的Gabor小波(improvedGaborwavelet,IGW)变)并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform,andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性:1)IGWT把时域信号映射到时间一频率域,而传统Gabor小波变换把时间信号映射到时间一尺度域;2)IGWT可用于信号分频,通过固定变换的主频参数dominantfrequency,并变换能提取相应的子带信号,且其主频部分的信息与原信号相应频率部分的信息一致,通过调节变换的分辨率因子,变换能有效控制子带信号的带宽;3)用IGWT和IGWIT构建的滤波器有良好的时一频局部性,在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性,并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。
简介:本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。
简介:转换波技术是解决气云区成像问题的有效手段,而二极化现象则是制约转换波构造成像精度提高的一个重要因素。解决该问题的一个有效办法是二极化校正:将转换波速度分解为基本速度和速度扰动,分两步对转换波地震数据进行叠加速度分析和叠前时间偏移速度分析。柴达木盆地三湖地区天然气资源丰富,气云广泛发育,通过对三湖地区二维三分量地震数据的研究分析,发现该地区转换波二极化现象明显,对转换波的准确成像有较大影响。通过对转换波二极化校正方法的应用,不论是在叠加成像,还是叠前偏移成像都取得了较好的效果,提高了气云区转换波构造成像的精度。二极化校正对解决气云区转换波剩余静校正有重要促进作用,二极化校正与转换波剩余静校正的结合是陆上气云区转换波高精度成像流程中必不可少的环节。
简介:现阶段的叠前地震反演技术中用于描述反射系数与纵、横波速度和密度之间的关系几乎完全是Zoeppritz方程的近似式,由于这些近似公式在大角度和弹性参数变化剧烈时误差较大,这不仅降低了反演解的精度,而且增加了叠前反演的多解性。本文探索了直接利用Zoeppritz方程求解精确反射系数的理论方法,并基于广义线性反演理论详细推导了基于叠前大角度地震资料的纵、横波速度和密度三参数同步反演算法,同时在反演过程引入正则化约束阻尼因子和共轭梯度算法,有效降低了反演的不适定性和提高了反演收敛性。理论模型试算和实际工区应用表明,本文提出的反演方法能够有效利用大角度(一般入射角〉30°)的叠前地震数据,获得更精确的地震弹性参数反演结果,并且反演结果忠实于地震资料,与井吻合较好。
简介:本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的P和S波被同一台站记录)来联合反演纵波速度(Vp)和纵、横波速度比(Vp/Vs),然后单独反演横波速度Vs,在反演过程中同时对地震参数进行定位。该方法不需要假设P和S波的射线路径一致,它是沿着P和S波射线路径计算相对慢度扰动值。该方法直接把Vp/Vs作为一个模型参数,由此能获得比采用从独立反演获得的Vp和Vs计算出Vp/Vs的方法更精确的速度比值。该新方法被应用到反演日本东北地区的壳幔速度及波速比结构的研究中,获得了较好的效果。反演结果表明,在日本东北地区,太平洋俯冲板块为一高Vp,高Vs和低Vp/Vs异常区,而在活火山下方的浅部地幔楔以及背弧深部地区为低Vp,低VS和高Vp/VS异常。虽然这些特征在前人的研究中已经报道过,但与前人的研究结果相比,本次研究所获得的Vp/Vs的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。