一种时频连续小波变换的井震标定方法技术

本发明专利技术公开了一种时频连续小波变换的井震标定方法，包括下列步骤：1)利用测井资料制作合成地震记录；2)采用时频连续小波变换方法分别对合成地震记录和井旁地震记录做时频变换，得到二者的时频分布图；3)调整时频分布图使二者的时频轴对齐，同时对合成地震记录和/或井旁地震记录进行同步调整，然后进行井震标定；4)计算合成地震记录和井旁地震记录的相关分析系数；如该相关分析系数≥0.9，则标定结果可靠；如该相关分析系数<0.9，则标定结果不可靠，返回步骤1)重新标定。通过本发明专利技术提供的井震标定方法，借助源自测井和地震的数据，能够更高效、准确地确定井震之间的时深关系，标定结果准确、可靠。

【技术实现步骤摘要】

本专利技术属于油气地球物理勘探
，具体涉及一种时频连续小波变换的井震标定方法。
技术介绍
众所周知，测井资料的垂向深度是以钻井深度表示，而地震资料的垂向深度是以时间表示，因此，井震标定是联系测井、地震资料的纽带和首要环节。常规的井震标定方法主要利用测井资料的声波、密度测井曲线制作合成地震记录，再与井旁地震道进行对比分析，从而进行时深标定。该方法对于构造勘探，以及大套地层的对比是可靠的。但是对于岩性油气藏的目标层位精细标定，由于测井中存在漂移、地层复杂等因素，难以达到准确标定层位。虽然现在已经发展了多元信息融合的井震标定，钻井、地质联合标定，井震分频标定等方法，也取得一些成功应用实例，但这些方法的成功应用受到多种因素制约，特别是对于一般工程技术人员对工区地质、测井和地震资料的熟悉程度和技术人员本身的经验影响较大，且标定结果的精度可靠性不高。因此，有必要专利技术一种受人为因素影响较小，而又能够简单、快速、准确地标定岩性油气藏目标层的井震标定方法。
技术实现思路
本专利技术的目的是提供一种时频连续小波变换的井震标定方法，借助时频连续小波变换对测井资料的合成地震记录和井旁地震记录的时频谱对比分析，确定目标层的层位的精细标定，提高井震标定的精度。为了实现以上目的，本专利技术所采用的技术方案是：一种时频连续小波变换的井震标定方法，包括下列步骤：1)利用测井资料制作合成地震记录；2)采用时频连续小波变换方法分别对步骤1)所得的合成地震记录和井旁地震记录做时频变换，得到二者的时频分布图；3)调整步骤2)所得时频分布图使二者的时频轴对齐，同时对合成地震记录和/或井旁地震记录进行同步调整，然后进行井震标定；4)计算井震标定结果中合成地震记录和井旁地震记录的相关分析系数；如该相关分析系数≥0.9，则标定结果可靠；如该相关分析系数<0.9，则标定结果不可靠，返回步骤1)重新标定。步骤1)中，将所得合成地震记录与井旁地震记录进行对比分析，进行初步井震标定，确定勘探目标层段；步骤2)中，在对时频分布图进行调整时，主要对勘探目标层段进行对比，使二者的时频轴对齐。在调整二者的时频轴时，应尽量使其对齐；调整的方法可以采用拉伸或压缩的方式，具体根据二者的时频轴对应关系调整。步骤2)中，所述时频连续小波变换的计算公式为：式中，f(w)为时频变换结果；w为圆频率；WTf(a,τ)为小波变换结果；是实数域中的一个连续函数；a为小波尺度；为小波变换系数。该时频连续小波变换具有较高时间、频率分辨率。步骤3)中，采用拉伸或压缩的方式调整时频分布图，同时对合成地震记录和/或井旁地震记录采用同样的方式进行同步调整。步骤4)中，所述相关分析系数的计算公式为：rXY=Σi=1N(Xi-X‾)(Yi-Y‾)Σi=1N(Xi-X‾)2Σi=1N(Yi-Y‾)2]]>式中，N为计算时间段内合成地震记录及井旁道地震数据的采样点数；Xi为第i点合成地震记录振幅值；为N点合成地震记录振幅平均值；Yi为第i点井旁道地震记录振幅值；为N点井旁道地震记录振幅平均值。步骤4)中，可以进行标定结果交互验证；结合已有钻井、测井、地质等资料验证标定结果的准确性，使主要勘探目标层的勘探结果与已钻井成果能够吻合。本专利技术的时频连续小波变换的井震标定方法，利用时频连续小波变换具有较高时间、频率分辨率的特点，结合井旁地震道的地震记录的时频特征与测井资料合成地震记录的时频特征之间对比和相关函数分析法，优选能够适合二者最佳的时频特征谱，进而实现测井资料与地震资料的精细标定。通过本专利技术提供的井震标定方法，借助源自测井和地震的数据，能够更高效、准确地确定井震之间的时深关系，标定结果准确、可靠。附图说明图1为本专利技术的时频连续小波变换的井震标定方法的实施流程图；图2为实施例1中常规合成地震记录与井旁地震记录对比图；图3为实施例1中利用时频连续小波变换获得的合成地震记录时频分布图和井旁地震记录时频分布图；图4为实施例1的方法处理后合成地震记录与井旁地震记录对比图。具体实施方式下面结合具体实施方式对本专利技术作进一步的说明。本专利技术是利用时频连续小波变换进行井震标定，其核心问题是计算具有时间、频率分辨率的时频连续小波变换。具体实施方式中，其时频连续小波变换具体实现原理如下：(a)连续小波变换：小波变换与传统傅里叶变换相比，具有变时域和变频域特性等优点，因此在地震信号时频分析和谱分解中得到了广泛应用。f(t)的连续小波变换(也称为积分小波变换)定义为：WTf(a,b)=|a|-1/2∫-∞∞f(t)ψ(t-ba)‾dt,a≠0---(1)]]>或用内积形式：WTf(a,b)＝<f,ψa,b>(2)式中，要使逆变换存在，ψ(t)要满足允许性条件：Cψ=∫-∞∞|ψ^(ω)|2|ω|dω<∞---(3)]]>式中，是ψ(t)的傅里叶变换。这时，逆变换为：f(t)=Cψ-1∫-∞∞∫-∞∞ψa,b(t)WTf(a,b)dbda|a|2---(4)]]>Cψ这个常数限制了能作为“基小波(或母小波)”的属于L2(R)的函数ψ的类，尤其是若还要求ψ是一个窗函数，那么ψ还必须属于L1(R)，即：∫-∞∞|ψ(t)|dt<∞---(5)]]>故是R中的一个连续函数。由式(4)可得在原点必定为零，即：ψ^(0)=∫-∞∞ψ(t)dt=0---(6)]]>从式(5)可以发现小波函数必然具有振荡性。(b)时频连续小波变换对式f(t)做傅里叶变换：f(w)=Cψ-1∫-∞∞∫-∞+∞∫-∞∞1a5/2ψa,b(t-τa)WTf(a,τ)dadτdt---(7)]]>根据傅里叶变换的性质：式中，是对作傅里叶变换的结果。结合式(7)和(8)，经过整理得到：式(9)就是时频连续小波变换的计算公式。实施例1本实施例的时频连续小波变换的井震标定方法，是对春光油田某井进行井震精细标定，其实施流程如图1所示，包括如下具体步骤：1)利用春光油田某井的声波测井速度、密度测井资料制作合成地震记录，并与井旁地震记录进行对比分析，进行初步井震标定，其结果如图2所示。图2中，P-wave_corr为声波测井速度曲线；Density为密度测井曲线；syn为合成地震记录；Inline98441为井所在地震道的位置，以此为中心左右各显示5个地震道，N1S2为沙湾组二段底界面的地震解释层位，C为石炭系顶界面的地震解释层位；2)采用具有较高时间、频率分辨率的时频连续小波变本文档来自技高网...

【技术保护点】
一种时频连续小波变换的井震标定方法，其特征在于：包括下列步骤：1)利用测井资料制作合成地震记录；2)采用时频连续小波变换方法分别对步骤1)所得的合成地震记录和井旁地震记录做时频变换，得到二者的时频分布图；3)调整步骤2)所得时频分布图使二者的时频轴对齐，同时对合成地震记录和/或井旁地震记录进行同步调整，然后进行井震标定；4)计算井震标定结果中合成地震记录和井旁地震记录的相关分析系数；如该相关分析系数≥0.9，则标定结果可靠；如该相关分析系数<0.9，则标定结果不可靠，返回步骤1)重新标定。

【技术特征摘要】
1.一种时频连续小波变换的井震标定方法，其特征在于：包括下列步骤：1)利用测井资料制作合成地震记录；2)采用时频连续小波变换方法分别对步骤1)所得的合成地震记录和井旁地震记录做时频变换，得到二者的时频分布图；3)调整步骤2)所得时频分布图使二者的时频轴对齐，同时对合成地震记录和/或井旁地震记录进行同步调整，然后进行井震标定；4)计算井震标定结果中合成地震记录和井旁地震记录的相关分析系数；如该相关分析系数≥0.9，则标定结果可靠；如该相关分析系数<0.9，则标定结果不可靠，返回步骤1)重新标定。2.根据权利要求1所述的时频连续小波变换的井震标定方法，其特征在于：步骤1)中，将所得合成地震记录与井旁地震记录进行对比分析，进行初步井震标定，确定勘探目标层段；步骤2)中，在对时频分布图进行调整时，主要对勘探目标层段进行对比，使二者的时频轴对齐。3.根据权利要求1或2所述的时频连续小波变换的井震标定方法，其特征在于：步骤2)...

【专利技术属性】
技术研发人员：肖学，田仁飞，刘焕，王鹏飞，邓继新，李秋菊，邢丽娟，黎腾，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司河南油田分公司物探研究院，
类型：发明
国别省市：北京;11