简介：本文通过人工变换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模型更接近实际地层情况。采用区分填充模式的等效介质模型，能够更好地实现对祁连山冻土区水合物粉砂岩地层的模拟。

简介：天然气水合物沉积层在常规叠后地震剖面上，一般表现为似海底反射（BSR）特征，常作为水合物识别的重要标志，但由于地震数据分辨率的限制与多解性的存在，类似的反射特征不一定是水合物的表现。本文先从南中国海神狐海域地质条件分析人手，综述研究区域天然气水合物发育的构造、沉积与运移环境；其次以加权方程与Zoeppritz公式的正演模拟为基础，讨论天然气水合物与BSR特征的关系；然后选取经过三口已钻获天然气水合物站位的地震测线，提取角度域共成像道集用于BSR的AVA特征研究；最后依据叠前地震弹性参数同步反演思路，对水合物重点赋存区域开展精细预测研究，获取了研究区域含水合物沉积层稳定带的分布特征与饱和度数。该方法可作为采用多道地震数据估算天然气水合物含量的一种有效技术。

简介：油层注入聚合物溶液后,电阻率测井响应发生与水驱不同的变化特征,严重影响水淹层测井评价的效果。本文首先分析了聚合物溶液的导电性质,其次按照聚合物驱油层不同施工流程开展岩电实验研究,分析不同驱替阶段的岩石电阻率的变化特征。并对比分析聚合物驱与水驱油层电阻率测井响应的差异。结果表明,不同聚合物驱替饱含油岩心过程,岩心电阻率随饱和度变化呈现"w"型或"S"型和单调下降的形态;清水聚合物驱油层电阻率测井响应呈现弱、未水淹的特征;污水聚合物驱油层电阻率测井响应呈现与污水水驱相似的变化特征。

简介：根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。

简介：海洋天然气水合物是一种潜在的巨大能源,而地球物理探测技术是勘探天然气水合物资源的重要方法。本文介绍了利用可控源电磁法探测海洋天然气水合物的可行性,在参考了大洋钻探计划(ODP)164航次对天然气水合物的采样数据的基础上建立了几种不同水合物含量的海底天然气水合物一维地电模型,并且利用建立的模型讨论了天然气水合物的频率域电磁响应。文章就天然气水合物的电场振幅值和相位随着收发距和频率等参数的变化关系等做了相应研究,同时就水合物厚度不同时其相应的电磁响应特征做了初步探讨,为海洋天然气水合物的勘探和资源评价提供了参考。

简介：本文利用取自南海神狐海域钻井岩芯进行了核磁共振测量,研究了岩石中所含水合物对地层渗透率的影响,该岩芯样品为一＂孔隙充填型＂水合物样品。核磁共振测量结果表明岩石中所含水合物对地层渗透率影响很大,随着含水合物饱和度的增加,地层渗透率迅速降低。通过对测量结果的进一步分析可知＂孔隙充填型＂水合物对地层渗透率的影响-存在多方面的因素,包括水合物对地层孔隙度、孔隙迂曲度的影响以及对水合物对孔隙喉道的堵塞等。最后本文利用Masuda下降模型建立了适用于南海神狐海域的相对渗透率与含水合物饱和度的经验公式,下降指数N=7.9718。

简介：这篇文章重点研究改进的Gabor小波（improvedGaborwavelet，IGW）变）并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform，andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性：1）IGWT把时域信号映射到时间一频率域，而传统Gabor小波变换把时间信号映射到时间一尺度域；2）IGWT可用于信号分频，通过固定变换的主频参数dominantfrequency，并变换能提取相应的子带信号，且其主频部分的信息与原信号相应频率部分的信息一致，通过调节变换的分辨率因子，变换能有效控制子带信号的带宽；3）用IGWT和IGWIT构建的滤波器有良好的时一频局部性，在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性，并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。

简介：地震数据的精确解释和分析非常依赖于所使用算法的稳定性。我们着重于盐丘地震勘探的稳定检测。我们讨论一个基于地震成像中最佳结构特征属性排序分类的盐丘探测新模型。该算法克服了现有的基于结构属性技术的局限性，因为该技术非常依赖于盐丘固有的地质属性与盐丘检测所用的属性数量。该算法综合了灰度共生矩阵（GLCM）属性，Gabor滤波器，以及含有使用属性特征排序信息协方差矩阵的本征结构等属性。将排序前列的属性组合起来形成一组最优的特征集，以保证算法即使在沿盐丘边界没有强反射层的情况下也能有效。与现有的盐丘检测技术相比，本文的算法稳定和计算高效，并能处理小尺度特征集。我用荷兰F3地块评价该算法的性能。实验结果表明，本文提出的基于信息理论的工作流程用于检测盐丘，其精度优于现有盐丘检测技术。

简介：逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求，井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异，本文将相应的变网格算法应用到偏移中。另外，也将Lanczos滤波引入到波场外推计算中，从而可以很好地消除变网格引起的人为反射。试算表明，对比地面、井中和井间的偏移，可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。

简介：随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明：该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。

简介：在井下微震监测中，各向异性模型用来获取事件位置和给出对介质实际描是十分有用的。在各向异性介质中，因为横波的分裂和奇异点的出现，使得波场的结构变得复杂。我们的研究展示了对强烈各向异性VTI介质利用速度模型校准和双偶极震源进行常规处理和运动学反演的结果。因小走时的拟合差，常用的质量评价标准，不是总正确的，人工加入各向同性介质层以降低拟合差，也可能产生非物理模型，负泊松比，事件位置的位移，

简介：本文描述了在全波实验中所见到的三种现象,即:1.纵横波初至的振幅(绝对值)具有相同的变化规律;2.纵横波初至的相位彼此总是反向的;3.纵横波初至的振幅变化周期为2π。经过充分的分析和解释指出这些现象应该出现在全波测井条件下,从而为利用相位差从全波中提取有用的信息奠定了基础。更多还原

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限，一是实时性差，二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处，作者基于反射式强度调制原理，研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器，实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃～5℃温度范围内冰生长过程中及5℃～-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明，与传统电阻应变式压力膜盒传感器相比，该传感器工作稳定、分辨率达0．02kPa、灵敏度达2．74×10-4／l（Pa，并可实现自动连续测量。

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限,一是实时性差,二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处,作者基于反射式强度调制原理,研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器,实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃~5℃温度范围内冰生长过程中及5℃~-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明,与传统电阻应变式压力膜盒传感器相比,该传感器工作稳定、分辨率达0.02kPa、灵敏度达2.74×10-4/kPa,并可实现自动连续测量。更多还原

简介：研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解，由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件，有效地消除了低波数干扰，同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析，在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像，利用校正因子对各角度的像进行校正，进而实现对成像结果的校正，从而实现了一种高效的倾角域幅度校正方法。最后通过SEG／EAGE模型进行数值计算验证了文中所述的计算方法。

简介：现阶段的叠前地震反演技术中用于描述反射系数与纵、横波速度和密度之间的关系几乎完全是Zoeppritz方程的近似式，由于这些近似公式在大角度和弹性参数变化剧烈时误差较大，这不仅降低了反演解的精度，而且增加了叠前反演的多解性。本文探索了直接利用Zoeppritz方程求解精确反射系数的理论方法，并基于广义线性反演理论详细推导了基于叠前大角度地震资料的纵、横波速度和密度三参数同步反演算法，同时在反演过程引入正则化约束阻尼因子和共轭梯度算法，有效降低了反演的不适定性和提高了反演收敛性。理论模型试算和实际工区应用表明，本文提出的反演方法能够有效利用大角度（一般入射角〉30°）的叠前地震数据，获得更精确的地震弹性参数反演结果，并且反演结果忠实于地震资料，与井吻合较好。

简介：本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。

简介：由于CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，从而更充分地利用了多次覆盖反射数据的信息。就目前的地震资料处理技术而言，它是最佳的零偏移距成像方式。本论文利用改进型的参数优化技术，得到高质量的CRS运动学参数剖面，并利用参数剖面计算出叠加孔径，实现了基于最优孔径的CRS叠加，使CRS参数的用途得到了充分利用。模型数据和实际资料的试算表明，基于最优孔径的CRS叠加的成像剖面与传统CRS叠加剖面相比，有着较高的信噪比和同相轴的连续性。

简介：阵列声波信号是典型的非线性、非平稳信号，Hilbert～Huang变换（HHT）是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解（EMD）和对瞬时频率的求解，可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解，任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征，在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程，并列举了一些声波测井波列实例，说明了该方法的有效性。