简介：分数阶S变换（FRST）具有较强的时频聚集性。利用FRST处理地震数据，通过合适的分数阶参数将频率轴旋转到适当位置，即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同，因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体，并利用共频率数据体进行了低频伴影分析，然后提出FRST和盲分离结合的方法，不需要对地震数据的最优分数阶参数进行估计，即可提取识别有效地质特征信息的独立频谱，提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据，并与已知井信息进行比对，验证了其有效性。

简介：分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。

简介：在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.

简介：特征值相干属性是一种高效、成熟的裂缝预测技术。但即使在连续、稳定的目的层段,部分裂缝隐没于相近高灰度值区域而难以识别;在破碎带发育区域,由于成像能量难以完美聚焦,中、小尺度的裂缝区域呈现相近低灰度值的云雾状模糊而无法区分。本文提出一种以直方图均衡化处理为核心技术的新型裂缝增强方法。该方法能够强化相干图像中不连续性信息与背景的差异性,突出裂缝的线性结构:采用相干图像阈值逐级调节的方式,实现对不同尺度裂缝的预测。该方法同时提升了显性、隐性裂缝识别能力和精度。更多还原

简介：边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原

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

简介：针对有限差分数值模拟的频散问题，本文将交错网格技术和紧致差分格式相结合，推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式；对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项，并利用Fourier误差分析方法分析了上述三种差分格式的近似精度；在此基础上，分别采用上述三种差分格式进行了波场数值模拟。结果表明，当差分方程阶数相同时，紧致交错网格差分格式截断误差最小，数值频散最弱，差分精度最高，证实了该方法的有效性。

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

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

简介：目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量（InSAR）中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性：新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。

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

简介：本文以中观孔隙结构的White模型为基础，构建了部分饱和孔隙介质模型，利用Biot方程的建立思路和Johnson推导的体变模量，推导了部分饱和孔隙介质中的纵波方程，并以平面波为例，求取了方程的衰减系数，分析了地震勘探频带范围内地震波的衰减特性。结果表明：在部分饱和孔隙介质中，地震波在低频段也会发生明显的衰减和频散现象，频率越大，衰减越大；且第二纵波的衰减比第一纵波更为明显；这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理，主要表现在在地震勘探频带范围内，波的衰减随孔隙度的增大而增大，随含油气饱和度的增大而减小，当孔隙内径尺寸小于二分之一外径尺寸时，波的衰减随内径尺寸的增大而增大，当内径尺寸大于二分之一外径尺寸时，波的衰减随内径尺寸增大而减小。

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

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

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

简介：中国多年来在岩石圈物理学研究中已取得了一系列的重要成果,并促进了中国地球科学的发展。然而在这21世纪之际,基于国家战略需求和自主创新的方针对岩石圈物理学提出了新的挑战。本文在这样的思维前提下对岩石圈物理学的科学内涵和发展导向进行了较全面的分析、研究和思考,并明确指出当今在这一学科领域应该做些什么,核心问题是什么,又存在哪些关键性的科学问题。研究结果提出,高精度的地球物理场观测与岩石圈内壳、幔精细结构(2维和3维)的刻划;在地球内部力系作用下,深部物质和能量的交换;深部物质运移的物理-力学-化学作用过程及深层动力学响应乃是深化对地球本体的认识和揭示成山、成盆、成岩、成矿和成灾的根本机理所在。为此,在本世纪的上、中叶,在地球科学领域中地球物理学研究必为先导。本文最后对岩石圈物理研究中尚存在的一些问题和困惑进行了分析和讨论。

简介：相干属性是检测地下介质不连续性特征的主要参数之一，但倾斜地层会对相干属性产生干扰，故在计算相干属性时需要消除地层倾角的影响。传统的相干属性计算一般使用时域倾角扫描法来消除地层倾角的影响，但这种方法的精度受到采样间隔等因素的限制，并不能很好地消除地层倾角的影响。根据振幅谱具有不受地层倾角影响的特点，本文提出将预先定义的子数据体内的地震数据变换为振幅谱，再利用振幅谱构造协方差矩阵，协方差矩阵最大特征根与矩阵迹的比值作为子数据体中心点的相干值。利用振幅谱计算相干属性的方法能够较好地消除地层倾角的影响，而且可以使用不同频段的振幅谱计算相干体，具备多尺度相干体算法的特性。模型数据和实际地震数据的应用结果表明，振幅谱相干属性可以更好地消除地层倾角对相干属性的影响。

简介：相干体是一种能够有效反映断层和裂缝等地质上非连续性的地震资料解释工具。然而由于受到附近地层,尤其是具有较强反射系数的地层的影响,微小断层和裂缝在相干体的时间切片和纵剖面上仍然难以辨别。本文提出了一种基于局部直方图规定化的相干体增强方法,实现相干体中的微小断层和裂缝的增强。与传统直方图规定化方法不同,本文方法处理三维相干体时不需要对数据进行离散化以统计其分布情况。在将相干体划分为若干子块后,以相干体整体的分布函数作为每个子块内的目标分布函数进行直方图规定化。另外,相邻子块的部分区域需要重叠,以克服分块产生的边界效应。对实际数据的处理结果表明,该方法可以提供更多的微小断层和裂缝的细节信息。

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

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