简介：随着煤矿的深部开采，多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁（简称TEM）法因场结构与地层的耦合关系，对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间，而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演，估算了探测目的层所需要的时间长度，根据噪声对观测数据造成的影响，给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线，表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中，利用电测井资料作为初始参数进行定量反演，约束了煤系薄层的等值性。所推断的三层和两层积水采空区，为钻孔所证实。研究结果表明，当观测时间有足够的长度、下伏地层的异常显示大于观测误差，中心回线TEM法探测多层采空区积水是可行性的。

简介：岩石物理学研究中Gassmann方程被广泛应用于预测岩石中的地震波速度,由于输入的基质矿物体积模量参数不准确,极大的影响预测结果的可靠性,特别是复杂基质矿物组合的碳酸盐岩储层。因此本文结合Russell流体因子和Gassmann-Boit-Geertsma方程计算式,通过引入干岩石骨架泊松比,提出了一种基质矿物体积模量提取方法,能够自适应反演岩石基质矿物等效体积模量,提高流体影响预测的可靠性,通过实际资料流体替换验证,该方法的预测结果是可靠的,并且计算效率高、适应性强。

简介：为了分析沁水盆地太原组煤层岩石物理性质、镜质组反射率与微观结构之间的关系，采集了三个不同煤矿井下的11个煤岩样品，利用MTS岩石物理参数测试系统对煤样进行弹性参数测试，通过薄片鉴定和扫描电镜对11个煤样的有机显微组分、镜质组最大反射率和微观结构特征进行了分析。对上述实验结果进行交会分析发现：一些弹性参数之间具有线性关系；纵波速度与横波速度均具有各向异性特征，其各向异性特征与镜质组最大反射率有良好的负相关性；纵横波速度、杨氏模量与微孔含量呈负相关关系。获得的煤岩岩石物理参数之间的关系构成了煤岩地震岩石物理建模和煤层气地震反演技术的基础。

简介：礁滩储层是我国海相碳酸盐岩油气勘探的重要目标，也是我国油气产能的重要接替领域之一。由于礁滩储层储集空间类型多样，非均质性极强，依靠常规测井判别其沉积相和岩性非常困难。成像测井能够清晰地反映礁滩储层的结构组分和沉积构造，为沉积相和岩性的判别提供了可靠的依据。在对大量礁滩储层的成像和岩心对比观察的基础上，我们提出了9种典型的成像解释模式，建立了利用成像解释模式自动判别礁滩储层沉积相和岩性的方法并研制了相应的处理软件。该方法在塔中和川东北地区礁滩储层的实际应用中取得了良好的效果。

简介：常规地震观测系统设计方法基于地下水平层状介质的假设,通常不能适应复杂构造情况。我们从控制照明的思想出发,提出了一种面向目标成像的地震观测系统设计方法,该方法需要一个由初步地震解释得到的速度模型。利用单程傅立叶有限差分波场传播算子将目标层的平面源延拓到地表,通过分析从目标层延拓到地表的波场能量的分布,可以确定目标层成像所需要的炮点或者检波点的位置。利用SEG-EAGE盐丘模型数值试算结果表明,该方法用于设计面向目标成像的特定地震采集系统。

简介：探索地震前兆是一件具有挑战性的研究。本文采用高通滤波的方法，处理分析了来自世界不同位置4个超导重力台站5次大地震前的秒级固体潮重力数据，发现在这些地震发生之前数天内，能够接收到稳定的重力高频扰动信号，其中2008年5月12日的汶川地震前出现一个主频为0.1-0.3Hz的重力扰动，而其他4次地震前出现两个主频为0.12-0.17Hz和0.06-0.085Hz的重力扰动。另一个发现是当地震发生在大陆与海洋板块断裂带时这种重力扰动信号具有远场效应，即重力扰动异常出现在距离震中较远的超导重力仪上，而较近的台站异常较小或没有异常。以上研究结果表明，这种重力扰动信号与地震发生具有较密切的相关性，可作为一种有潜力的临震预报信号；重力扰动信号的远场效应也许能够揭示某些地震发生前地球内部板块之间的相互作用机制。但由于地球上重力固体潮台站分布的非常不均匀，上述结果还需要今后研究更多的样本数据加以证实。

简介：本文提出了基于测井、VSP和地震数据拟合的子波估计方法,从输入输出都包含随机噪声的统计模型出发,采用相关性拟合技术来提取子波。拟合度和误差分析为整个过程提供了定量的质量控制手段,可以评估数据拟合和子波估计的可靠性。实际数据试算表明,该方法在含有噪声的实际数据中稳定而有效,在地震频带内的子波估计和数据拟合是可靠的。该方法无需对子波相位和振幅谱进行任何假设,其主要优点在于确定相位的能力。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

简介：本文利用x射线CT获取反映岩心微观结构的三维数字岩心，利用数学形态学中的开运算模拟了岩石的油驱水排驱过程中，不同含水饱和度下油和地层水在孔隙空间中的分布。利用有限元方法计算了岩石电阻率，进而得到岩石地层因素和电阻率指数，并考查了岩石润湿性对岩石电阻率指数的影响。数值模拟结果表明：基于数字岩心的水湿岩石地层因素和电阻率指数数值模拟结果与实验结果一致，拓展了岩石电阻率实验的能力；岩石润湿性对岩石电性有重要影响，在相同含水饱和度下，油湿岩石电阻率高于水湿岩石电阻率，油湿岩石饱和度指数远大于水湿岩石饱和度指数。

简介：目前观测系统设计中接收道数的选择基本上是基于水平层状介质假设，这显然不能满足复杂地区采集的需求。本文在朱金平研究基础上，通过野外试验法获得最优接收道数问题，进行数值模拟实现，在固定道间距情况下，建立接收总能量与平均效率能量之间的数学模型，获得最优接收道数，并利用一个复杂工区的地震数据对方法的正确性进行验证，模型数据计算的结果与实际数据处理的结果表明，本文方法与实际处理结论一致，说明该方法具有较高的科学性，可以为复杂地区观测系统设计提供新的思路。

简介：渗透率评价是储层评价中的一个复杂问题，传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响，因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型，然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系，以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析，发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算，结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好，证明该方法具有一定的可行性。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

简介：通过常规测井曲线小波变换来评价裂缝参数是一个较新的研究课题。本文针对火山岩裂缝储层,对常规测井曲线做coif5、bior4.4和db5基小波的小波变换,通过变换后的分解信号与成像测井裂缝密度的对比研究,寻找出与裂缝密度相匹配的基小波,再将常规测井曲线在这个基小波下的分解信号与曲线变化率法相结合建立裂缝指示曲线来识别裂缝发育段,最后通过裂缝指示曲线与成像测井的裂缝密度之间的量化关系来比较准确的评价裂缝密度。将裂缝指示曲线法在松辽盆地南部进行了应用,通过裂缝指示曲线求取的裂缝密度与成像测井的裂缝密度之间的相对误差为0.53。

简介：地震纹理属性与地震相和储层特征密切相关,广泛地应用在地震资料的解释中。传统的地震纹理属性基于叠后数据提取,受叠加作用的影响,易造成地层特征信息的损失,不利于复杂储层的描述。本文提出叠前纹理属性,其不仅可以精细地刻画不同反射点波形的横向连续性,也能体现AVO、各向异性和介质的均质性。基于叠前纹理属性丰富的地层特征表达能力,结合SOM聚类算法,形成了利用叠前数据进行地震相分析的方法。该方法应用于中国某工区宽方位地震资料,通过对比证实了叠前纹理属性描述地层横向变化的优越性,并能揭示各向异性特征及非均质性特征,基于叠前纹理的分类结果能有效区分不同的地震反射模式,为地震相分析提供了可靠的依据。

简介：为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区（近区、过渡区和远区）资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。

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

简介：本文给出一种既能有效衰减地震噪音又可保护地层及构造的不连续性的新方法。构造约束保边平滑技术需要已知反射局部方位和边界信息，通常这些信息由全频率地震资料估算获得，但在资料信噪比很低的情况下，噪音往往会降低估算的可靠度。对于信噪比极低的地震资料，其主频成分相对非主频成分信噪比高，所以由主频资料获取的方位和边界信息比由其它频率成分获取的更可靠。方位和边界信息通常用倾角和相干值差异来描述。由于不同频率所引起的倾角和相干值差异的变化均比地震记录的变化缓慢，所以由主频资料获取的倾角及边界信息能够近似代表所有频率成分的倾角及边界信息。Ricker子波广泛用于地震勘探，Marr小波与Ricker子波在时间和频率域均具有相同的形态，所以选用Marrl小波变换将地震数据按照倍频程分为几个分频体。扫描主频分频体，用不等权二次曲面拟合并求解极大值来获取视倾角，通过比较9个滑动窗口的相干值来确定反射边界。将这些信息用构造约束保边平滑技术可选择性地（selectively）对主频、低频、高频分频体做平滑处理，最后将平滑后的各频段地震记录合成为滤波去噪后的地震记录。理论模型和实际资料处理效果表明该方法能有效压制噪音，保护边界，保护同相轴的连续性，且灵活地保留地震记录中的有用信息。

简介：在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现：能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。

简介：常规的时间一空间域和频率一空间域预测滤波方法假设地震记录由地震信号和随机噪声两部分构成，即所谓的加噪声模型，但是，在对随机噪声进行估算时，又假设随机噪声可以通过预测误差滤波器由地震记录中进行预测，即所谓的源噪声模型。这种前后不一致的噪声模型降低了该类方法的去噪能力和保幅性能。为此，本文提出了一种基于反演的时空域随机噪声衰减方法。它首先从地震数据中估算预测滤波算子，该算子表征了地震信号的可预测性，自适应地描述了地震信号的空间结构。在得到预测误差算子之后，将该算子作为正则化约束引入到地震信号反演系统，由含有随机噪声的地震数据直接反演地震信号。不同于常规随机噪声衰减方法，该方法将随机噪声衰减问题归结为正则化约束下的地震信号反演问题，克服了常规方法噪声模型的不一致性问题。我们采用模型数据和实际数据进行了实验分析，并与常规方法进行了效果对比。实验结果表明：与常规方法相比，本文方法在噪声压制的同时，没有对有效信号产生明显伤害，具有更好的振幅保持能力。

简介：地下断层深度的估算是重力解释难题之一，我们试利用支持向量分类（SvC）法进行计算。使用正演和非线性反演技术，通过相关误错使检测地下断层深度成为可能。但必要有一个深度初始猜测值，而且这猜测值通常不是由重力资料得。本文我们介绍以SVC作为利用重力数据估算断层深度的一种手段。在这项研究中，我们假设一种地下断层深度可归为一种类型，SVC作为一个分类算法。为了有效地利用此SVC算法，我们基于一个正确的特征选择算法去选择正确的深度特征。本次研究中我们建立了一套基于不同深度地下断层的合成重力剖面训练集，用以训练用于计算实际的地下断层深度的SVC代码。然后用其它合成重力剖面训练集测试我们训练的SVC代码，同时也用实际资料验证了我们的训练SVC代码。