CN107831542B - Ddw高精度深度域井震匹配方法 - Google Patents
Ddw高精度深度域井震匹配方法 Download PDFInfo
- Publication number
- CN107831542B CN107831542B CN201710608976.0A CN201710608976A CN107831542B CN 107831542 B CN107831542 B CN 107831542B CN 201710608976 A CN201710608976 A CN 201710608976A CN 107831542 B CN107831542 B CN 107831542B
- Authority
- CN
- China
- Prior art keywords
- depth
- point
- well
- seismic
- domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000008569 process Effects 0.000 title claims abstract description 31
- 239000002131 composite material Substances 0.000 claims abstract description 24
- 230000001186 cumulative effect Effects 0.000 claims abstract description 19
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000005553 drilling Methods 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 6
- 230000008901 benefit Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002203 pretreatment Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 230000000153 supplemental effect Effects 0.000 claims description 3
- 238000003325 tomography Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 description 7
- 230000005012 migration Effects 0.000 description 6
- 238000013508 migration Methods 0.000 description 6
- 230000006835 compression Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000155 melt Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 210000002706 plastid Anatomy 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/59—Other corrections
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种DDW高精度深度域井震匹配方法,该方法包括:步骤1,输入地震数据与地震地质层位,并对数据进行预处理;步骤2,选择目标井的井旁地震道X,进行深度域合成记录制作得到Y,参照地质分层进行合成记录样点选取;步骤3,在地质分层约束下扫描井旁道;步骤4,逐点进行距离计算,寻找与合成记录样点距离最短的点;步骤5,求取累计距离,规划最短路径;步骤6,建立逐点对应的井震匹配关系,计算井震误差校正量;步骤7,建立偏差约束体校正地震数据。该DDW高精度深度域井震匹配方法以井震深度误差影响因素量版为约束,进行深度域井旁地震道与合成记录道之间的动态深度规划相似性分析,首次实现了全井段的井震高效匹配。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种DDW高精度深度域井震匹配方法。
背景技术
目前深度域地震资料在油田勘探开发中的应用并不广泛,其主要原因之一是地震深度与实际深度之间存在一定的误差,从浅层到深层误差变化较大,且有正有负,规律性较难把握。如何确定该井震深度误差并对其合理校正、实现井震高精度匹配是发挥深度域地震优势、实现构造特征与地质体分布直接预测的关键环节。
利用叠前深度偏移方法得到的深度域地震资料具有其独特的优势,无需时深转换就可以与钻井深度直接对应,直接进行构造解释与构造成图。叠前深度偏移技术的研究是近10多年来全球油气地球物理勘探领域的热点,其理论体系不断完善,也是目前国际上公认的解决复杂构造成像问题的有效途径。深度域成像技术将地震成像、模型模拟和解释融于一体,成像结果为深度剖面,特别是随着逆时偏移方法的发展,其成像点与地质模型绕射点位置完全一致,成像位置准确,因此它具有时间域成像不可比拟的优势。
但深度域地震资料在油田勘探开发中的应用并不广泛,其主要原因之一是地震深度与实际深度之间存在一定的误差,从浅层到深层误差变化较大,且有正有负,该误差究竟与哪些因素有关?其分布特征如何?如何将其用于深度域地震资料的校正?目前国内外对深度域地震资料处理及速度建模方面的研究较多,而深度域资料解释应用方面的较少,深度域井震匹配、误差校正等技术尚不成熟,也是造成深度域地震资料应用具有较大局限性的重要原因之一。近年来,为更大程度上减少井震误差,针对高成熟探区地震数据处理技术上不断探索各向异性地震数据处理方法,如各向异性高斯束叠前深度偏移、各向异性逆时偏移等的研究逐渐成为目前的热点问题,但由于速度拾取不准确等问题造成仍存在偏移深度与实际地层深度不符的现象。解释方面对深度域层位标定、误差校正、属性、反演等新的方法技术不断出现,对井震标定的适用性、深度域的深度校正问题进行了探讨。林伯香等,何惺华等,张雪健等,王永刚等对深度域的子波渡、褶积和Fourier变换等基本问题进行探讨并提出了深度域合成地震记录制作方法,胡中平等提出了伪深度变换的方法,解决了深度域中子波随深度变化的问题,全紫荆以此为基础探讨了深度域地质统计学反演方法,刘斌采用多尺度分析的方法对合成地震道与井旁道进行分解后,利用最小平方法求取匹配滤波因子实现时深两域的井震匹配。但是,如何确定井震深度误差并对其合理校正、实现井震高精度匹配是当前发挥深度域地震资料优势、实现构造特征与地质体分布直接预测的关键。
目前深度域地震资料部分只用于构造解释,在实际应用过程中亦主要是按照时间域地震资料解释的工作流程来进行解释。具体是通过深度域叠前偏移速度场将深度域资料转换到时间域,在时间域进行层位标定、构造解释及储层研究。这样虽然能够解决深度域解释中的一些问题,但增加了多套数据间的转换带来的误差,并且增加了研究人员多套资料解释的工作量。要进一步推动叠前深度偏移资料的应用就必须首先要解决构造解释中存在的井震深度不匹配和地震地质层位标定等一系列问题,如果所有工作(构造解释、储层预测等)能够在深度域完成,就能使得井震数据实现真正的统一,其结果能够真实的反映地下地质特征。为此我们发明了一种新的DDW高精度深度域井震匹配方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种能够刻画信号的局部拉伸或压缩关系,将深度域地震信号和深度域合成记录信号按同一间隔均匀采样,通过寻找地震信号与测井合成记录信号之间的最优匹配,即可得到二者之间的深度误差校正量,实现对深度域地震信号的深度校正的DDW高精度深度域井震匹配方法。
本发明的目的可通过如下技术措施来实现:DDW高精度深度域井震匹配方法,该DDW高精度深度域井震匹配方法包括:步骤1,输入地震数据与地震地质层位,并对数据进行预处理;步骤2,选择目标井的井旁地震道X,进行深度域合成记录制作得到Y,参照地质分层进行合成记录样点选取;步骤3,在地质分层约束下扫描井旁道;步骤4,逐点进行距离计算,寻找与合成记录样点距离最短的点;步骤5,求取累计距离,规划最短路径;步骤6,建立逐点对应的井震匹配关系,计算井震误差校正量;步骤7,建立偏差约束体校正地震数据。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,输入的数据包括:深度域地震数据、测井数据、地震解释层位、地质分层数据。
在步骤1中,进行深度域地震数据预处理,对深度域地震数据开展高阶累计量时频分析和频谱分析,获取相应的频谱类参数数据体,提取地震子波。
在步骤1中,进行深度域测井数据预处理,提取声波时差、密度这些测井曲线,完成拼接、测井环境校正等曲线预处理。
在步骤2中,根据地震数据和测井数据的位置信息自动将测井位置定位在相应的地震位置上,将测井的分层信息直接标在深度域地震剖面上,统计井点处地质分层与地震解释的深度偏差表,建立各个解释层位的深度偏差量版。
在步骤2中,对地震子波按照声波时差及其对应深度关系开展深变处理,利用褶积原理完成深度域合成记录制作。
在步骤4中,首先计算地震信号与合成记录信号这两个序列中每个采样点之间的距离,即相似度,振幅值差越小,距离越小,二者越相似;定义D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离,表示为D(i,j)=|R(i)-S(j)|(式1)
R(i)与S(j)分别为两列信号各个元素的振幅值,这样D(i,j)构成一个ixj的矩阵。
在步骤4中,进行的逐点距离计算包括波形幅度、频率、时频分布。
在步骤4中,寻找与合成记录样点距离最短的点时,把地震信号序列作为X坐标轴,即X轴坐标为1,2,…,N;把合成地震记录信号作为Y坐标轴,即Y轴坐标为1,2,…,M,每个地震深度点有一个相应的测井深度点相对应,定义第一个对应点为w1,第二个点为w2,直到最后一个深度点wk,这样{w1,w2,…,wk}构成一条路径,该路径反映了地震信号在每个深度点对应的测井深度点。
在步骤4中,在未找到与合成记录样点距离最短的点时,流程返回到步骤3。
在步骤5中,采用基于递推的动态规划的方法求出该累计距离,从后向前递推,若(N,M)是一个匹配点,即地震深度N对应测井深度M,则该路径点前一个格点可能是(N-1,M),(N-1,M-1),(N,M-1),找出这三个点中距离最小者即为其上一最优路径点,其累计距离公式为
TD(i,j)=D(i,j)+min{D(i-1,j),D(i-1,j-1),D(i,j-1)} (式2)
其中,TD(i,j)为某一路径上的累计距离,D(i,j)表示两记录中相应采样点i与j的距离;从最后一点或某个确定点出发搜索,反复递推,直到到达地震深度的起始点,构建出一条累计距离最小的路径,该路径即为两信号匹配的最优路径。
在步骤5中,以某一个采样点为中心的一定时窗内的信号频率信息对于确定两者的相似度是一个非常有效的参数,基于此提出在距离计算上定义为
D(i,j)=A1*|R(i)-S(j)|+A2*|F(R(i))-F(S(j))|+A3*|T(R(i))-T(S(j))| (式3)
其中,A1、A2、A3为加权因子,F(i)为主频,T(i)为时频,R(i)与S(j)分别为两列信号各个元素的振幅值,D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离。
在步骤5中,在距离归一化的基础之上,将确定好的地震深度与测井深度对应的数据点振幅值修改为相同的并与其它数据点对有明显区别的数值,使得优选的路径必须通过这一确定的深度点,实现了已知点层位的控制。
在步骤6中,将解释层位、断层按照地层的接触关系包括平行、上超、剥蚀建立三维框架模型。
在步骤7中,利用DDW高精度井震匹配得到的井震对应关系即可得到井点处全井段井旁地震道的误差校正量,以此为基础,应用克里金三维插值方法将上述偏差校正量进行空间插值,建立偏差约束体,与地震数据进行运算即可得到与实钻井深度高度一致的客观合理的深度域地震体。
本发明中的DDW高精度深度域井震匹配方法,对深度域井震资料预处理,从深度域地震资料中提取子波,根据声波速度得到深变子波,制作深度域合成记录,建立地质分层与深度域地震解释层位关联;对动态深度规划算法及其改进,提出一种基于动态深度规划(Dynamic Depth Warping,DDW)算法的高精度井震匹配方法,以误差量版为约束,进行深度域井旁地震道与合成记录道之间的相似度分析,利用递推规划算法综合计算样点距离、波形幅度、波数等特征确定二者之间的对应关系,并求取各小层内地震同相轴与实际深度的深度偏差,得到井旁地震道数据的偏差校正量。建立逐点对应的井震匹配关系,计算三维井震偏差校正量,实现高精度井震匹配。
附图说明
图1为本发明的DDW高精度深度域井震匹配方法的一具体实施例的流程图;
图2为本发明的一具体实施例中DDW最优匹配路径方法原理图;
图3为本发明的一具体实施例中高精度井震匹配结果图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的DDW高精度深度域井震匹配方法的流程图。
步骤101,输入深度域地震数据、测井数据、地震解释层位、地质分层数据,并对数据进行预处理。包括:深度域地震数据预处理,对深度域地震数据开展高阶累计量时频分析和频谱分析,获取相应的频谱类参数数据体,提取地震子波;深度域测井数据预处理,提取声波时差、密度等测井曲线,完成拼接、测井环境校正等曲线预处理。
步骤102,选择目标井的井旁地震道X。
根据地震数据和测井数据的位置信息自动将测井位置定位在相应的地震位置上,将测井的分层信息直接标在深度域地震剖面上,统计井点处地质分层与地震解释的深度偏差表,建立各个解释层位的深度偏差量版;将深度域地震信号和深度域合成记录信号按同一间隔均匀采样,其中地震信号序列定义为R(i)(i=1,2,…,N),R(i)为每个采样点的振幅值,该序列长度为N,测井合成地震信号定义为S(j)(j=1,2,…,M),S(j)为每个采样点的振幅值,其长度为M。
步骤103,进行深度域合成记录制作得到Y。对地震子波按照声波时差及其对应深度关系开展深变处理,利用褶积原理完成深度域合成记录制作。
步骤104,参照地质分层进行合成记录样点选取。
步骤105,在地质分层约束下扫描井旁道。
步骤106,逐点进行距离计算(波形幅度、频率、时频分布等)。
首先计算地震信号与合成记录信号这两个序列中每个采样点之间的距离,即相似度,振幅值差越小,距离越小,二者越相似。定义D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离,表示为D(i,j)=|R(i)-S(j)| (式1)
R(i)与S(j)分别为两列信号各个元素的振幅值,这样D(i,j)构成一个ixj的矩阵。
步骤107,寻找与合成记录样点距离最短的点。如图2所示,把地震信号序列作为X坐标轴,即X轴坐标为1,2,…,N;把合成地震记录信号作为Y坐标轴,即Y轴坐标为1,2,…,M。这样构成一个二维网格,二维网格中的交汇点为地震信号某个深度与测井(实际)某个深度的交汇点,显然地震深度与测井深度不是一一对应的,存在位置偏差,地震深度局部存在拉伸或压缩现象。每个地震深度点应有一个相应的测井深度点相对应,定义第一个对应点为w1,第二个点为w2,直到最后一个深度点wk,这样{w1,w2,…,wk}构成一条路径,该路径反映了地震信号在每个深度点对应的测井(实际)深度点。
步骤108,求取累计距离,规划最短路径。
如图2所示,采用基于递推的动态规划的方法求出该累计距离,从后向前递推,若(N,M)是一个匹配点,即地震深度N对应测井深度M,则该路径点前一个格点可能是(N-1,M),(N-1,M-1),(N,M-1),找出这三个点中距离最小者即为其上一最优路径点,其累计距离公式为
TD(i,j)=D(i,j)+min{D(i-1,j),D(i-1,j-1),D(i,j-1)} (式2)
其中,TD(i,j)为某一路径上的累计距离,D(i,j)表示两记录中相应采样点i与j的距离。从最后一点或某个确定点出发搜索,反复递推,直到到达地震深度的起始点,构建出一条累计距离最小的路径,该路径即为两信号匹配的最优路径。由于地震信号与测井合成记录信号是多套地层的褶积叠合反映,其距离(相似度)的计算虽以振幅值为主体,但显然以某一个采样点为中心的一定时窗内的信号频率等信息对于确定两者的相似度也是一个非常有效的参数,基于此提出在距离计算上定义为
D(i,j)=A1*|R(i)-S(j)|+A2*|F(R(i))-F(S(j))|+A3*|T(R(i))-T(S(j))| (式3)
其中,A1、A2、A3等为加权因子,F(i)为主频,T(i)为时频。
已知的确定性地震深度与测井深度的对应关系应用到路径动态规划选取过程中很有必要,在路径最优选取过程中能够起到关键节点约束的作用。在距离归一化的基础之上,将确定好的地震深度与测井深度对应的数据点振幅值修改为相同的并与其它数据点对有明显区别的数值,这样在动态规划路径选取中由于距离最小原则,这一对振幅值的的距离为零,而与其它任何一点的距离都会非常巨大,从而使得优选的路径必须通过这一确定的深度点,实现了已知点层位的控制。
相应地将前述得到的深度偏差量版应用到最优路径选取中,即两列信号在某个深度段上下位置差在一定的范围内(偏差量版)。
距离参数的优选组合,地震信号与测井合成记录信号是多套地层的褶积叠合反映,其距离(相似度)的计算虽以振幅值为主体,以某一个采样点为中心的一定时窗内的信号频率等信息对于确定两者的相似度也是一个非常有效的参数。在最优匹配方面,可以加入如主频、时频等信息,并进行优化组合。如距离计算上也可以定义为
D(i,j)=w1*|R(i)-S(j)|+w2*|F(R(i))-F(S(j))|+w3*|T(R(i))-T(S(j))|
其中,w1、w2、w3等为加权因子,F(i)为主频,T(i)为时频。
步骤109,建立逐点对应的井震匹配关系,得到井点处全井段井旁地震道的井震误差校正量,实现高精度井震匹配。将解释层位、断层等按照地层的接触关系(平行、上超、剥蚀等)建立三维框架模型。如图3所示,左图为井震匹配关系图,应用步骤108中(图2)得到的两列信号的位置匹配关系建立相应的井震对应关系,二者相减得到深度域地震深度与实际深度(钻井深度)的偏差量。
步骤110,利用克里金三维插值方法将步骤109求得的偏差校正量进行空间插值,建立偏差约束体。
利用DDW高精度井震匹配得到的井震对应关系即可得到井点处全井段井旁地震道的误差校正量,以此为基础,应用克里金三维插值方法将上述偏差校正量进行空间插值,建立偏差约束体,与地震数据进行运算即可得到与实钻井深度高度一致的客观合理的深度域地震体。流程结束。
在应用本发明的一具体实施例中,如图2所示,将深度域地震信号和深度域合成记录信号按同一间隔均匀采样,其中地震信号序列定义为R(i)(i=1,2,…,N),R(i)为每个采样点的振幅值,该序列长度为N,测井合成地震信号定义为S(j)(j=1,2,…,M),S(j)为每个采样点的振幅值,其长度为M。定义D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离,表示为
D(i,j)=|R(i)-S(j)| (式1)
把地震信号序列作为X坐标轴,即X轴坐标为1,2,…,N;把合成地震记录信号作为Y坐标轴,即Y轴坐标为1,2,…,M。这样构成一个二维网格,二维网格中的交汇点为地震信号某个深度与测井(实际)某个深度的交汇点,显然地震深度与测井深度不是一一对应的,存在位置偏差,地震深度局部存在拉伸或压缩现象。每个地震深度点应有一个相应的测井深度点相对应,定义第一个对应点为w1,第二个点为w2,直到最后一个深度点wk,这样{w1,w2,…,wk}构成一条路径,该路径反映了地震信号在每个深度点对应的测井(实际)深度点。
首先计算地震信号与合成记录信号这两个序列中每个采样点之间的距离,即相似度,振幅值差越小,距离越小,二者越相似。采用基于递推的动态规划的方法求出该累计距离,从后向前递推,若(N,M)是一个匹配点,即地震深度N对应测井深度M,则该路径点前一个格点可能是(N-1,M),(N-1,M-1),(N,M-1),找出这三个点中距离最小者即为其上一最优路径点,其累计距离公式为
TD(i,j)=D(i,j)+min{D(i-1,j),D(i-1,j-1),D(i,j-1)} (式2)
其中,TD(i,j)为某一路径上的累计距离,D(i,j)表示两记录中相应采样点i与j的距离。从最后一点或某个确定点出发搜索,反复递推,直到到达地震深度的起始点,构建出一条累计距离最小的路径,该路径即为两信号匹配的最优路径。
图3为本发明的一具体实施例中高精度井震匹配结果图。其中,左图为实际合成道,右图为低频速度下的合成道。左图可以看出,匹配前由于错误的低频速度进行地震处理得到的地层深度位置与实际模型存在一定的偏差,在高速层信号存在压缩现象,而在低速层信号存在拉伸现象,井旁道与合成记录深度与波形之间均有偏差。右图为利用发明方法匹配之后的结果,地震层位与测井层位对应关系较好,各个深度点得到归位,两信号的相关系数从归位匹配前的44%提升到96%,由此实现对低频速度模型得到的合成地震道进行深度校正,地震资料的深度误差得到较好的修正。
本发明的DDW高精度深度域井震匹配方法,针对陆相断陷盆地典型地质体开展实钻井的深度域井震误差研究,分析井震匹配度及误差产生的原因,定量分析误差与各地质要素之间的关系,建立井震误差成因量版;为更好的将地震层位与测井深度相对应,以此误差成因量版为约束,提出一种DDW(动态深度规整,dynamic depth warping)的高精度井震匹配算法,该算法能够刻画信号的局部拉伸或压缩关系。通过计算地震信号与合成记录信号这两个序列中每个采样点之间的距离,即相似度,振幅值差越小,距离越小,则二者越相似。采用基于递推的动态规划的方法求出该累计距离,两个信号累计距离和最小的路径就是找到的最优路径,从而寻找到地震信号与测井合成记录信号之间的最优匹配;同时,该动态深度规整算法针对地震深度信号匹配上有其特殊性,首先深度不能往回转,由于地震信号是多个地震反射界面的褶积而成,这些反射界面的真实深度与频率相关,因此在距离计算方面需要考虑频率等参数信息。解释层位或确定性的地震深度在匹配路径规划中也需要特殊考虑,同时不同类型地质体的深度误差量版在一定程度上可以提高路径规划的效率和精度。本发明的核心是克服了现有时间域地震不能与钻井资料直接对应及时间域地震资料对复杂构造及特殊地质体成像存在误差的现状,以井震深度误差影响因素量版为约束,进行深度域井旁地震道与合成记录道之间的动态深度规划相似性分析,首次实现了全井段的井震高效匹配。
Claims (13)
1.DDW高精度深度域井震匹配方法,其特征在于,该DDW高精度深度域井震匹配方法包括:
步骤1,输入地震数据与地震地质层位,并对数据进行预处理;
步骤2,选择目标井的井旁地震道X,进行深度域合成记录制作得到Y,参照地质分层进行合成记录样点选取;
步骤3,在地质分层约束下扫描井旁道;
步骤4,逐点进行距离计算,寻找与合成记录样点距离最短的点;
步骤5,求取累计距离,规划最短路径;
步骤6,建立逐点对应的井震匹配关系,计算井震误差校正量;
步骤7,建立偏差约束体校正地震数据;
在步骤5中,采用基于递推的动态规划的方法求出该累计距离,从后向前递推,若(N,M)是一个匹配点,即地震深度N对应测井深度M,则该路径点前一个格点可能是(N-1,M),(N-1,M-1),(N,M-1),找出这三个点中距离最小者即为其上一最优路径点,其累计距离公式为
TD(i,j)=D(i,j)+min{D(i-1,j),D(i-1,j-1),D(i,j-1)} (式2)
其中,TD(i,j)为某一路径上的累计距离,D(i,j)表示两记录中相应采样点i与j的距离;从最后一点或某个确定点出发搜索,反复递推,直到到达地震深度的起始点,构建出一条累计距离最小的路径,该路径即为两信号匹配的最优路径;在距离计算上定义为
D(i,j)=A1*|R(i)-S(j)|+A2*|F(R(i))-F(S(j))|+A3*|T(R(i))-T(S(j))| (式3)
其中,A1、A2、A3为加权因子,F(i)为主频,T(i)为时频,R(i)与S(j)分别为两列信号各个元素的振幅值,D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离。
2.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤1中,输入的数据包括:深度域地震数据、测井数据、地震解释层位、地质分层数据。
3.根据权利要求2所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤1中,进行深度域地震数据预处理,对深度域地震数据开展高阶累计量时频分析和频谱分析,获取相应的频谱类参数数据体,提取地震子波。
4.根据权利要求2所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤1中,进行深度域测井数据预处理,提取声波时差、密度这些测井曲线,完成拼接、测井环境校正这些曲线预处理。
5.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤2中,根据地震数据和测井数据的位置信息自动将测井位置定位在相应的地震位置上,将测井的分层信息直接标在深度域地震剖面上,统计井点处地质分层与地震解释的深度偏差表,建立各个解释层位的深度偏差量版。
6.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤2中,对地震子波按照声波时差及其对应深度关系开展深变处理,利用褶积原理完成深度域合成记录制作。
7.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤4中,首先计算地震信号与合成记录信号这两个序列中每个采样点之间的距离,即相似度,振幅值差越小,距离越小,二者越相似;定义D(i,j)为R(i)与S(j)的振幅值之差,即两帧信号采样点之间的距离,表示为D(i,j)=|R(i)-S(j)| (式1)R(i)与S(j)分别为两列信号各个元素的振幅值,这样D(i,j)构成一个ixj的矩阵。
8.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤4中,进行的逐点距离计算包括波形幅度、频率、时频分布。
9.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤4中,寻找与合成记录样点距离最短的点时,把地震信号序列作为X坐标轴,即X轴坐标为1,2,…,N;把合成地震记录信号作为Y坐标轴,即Y轴坐标为1,2,…,M,每个地震深度点有一个相应的测井深度点相对应,定义第一个对应点为w1,第二个点为w2,直到最后一个深度点wk,这样{w1,w2,…,wk}构成一条路径,该路径反映了地震信号在每个深度点对应的测井深度点。
10.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤4中,在未找到与合成记录样点距离最短的点时,流程返回到步骤3。
11.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤5中,在距离归一化的基础之上,将确定好的地震深度与测井深度对应的数据点振幅值修改为相同的并与其它数据点对有明显区别的数值,使得优选的路径必须通过这一确定的深度点,实现了已知点层位的控制。
12.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤6中,将解释层位、断层按照地层的接触关系包括平行、上超、剥蚀建立三维框架模型。
13.根据权利要求1所述的DDW高精度深度域井震匹配方法,其特征在于,在步骤7中,利用DDW高精度井震匹配得到的井震对应关系即可得到井点处全井段井旁地震道的误差校正量,以此为基础,应用克里金三维插值方法将上述误差校正量进行空间插值,建立偏差约束体,与地震数据进行运算即可得到与实钻井深度高度一致的客观合理的深度域地震体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710608976.0A CN107831542B (zh) | 2017-07-24 | 2017-07-24 | Ddw高精度深度域井震匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710608976.0A CN107831542B (zh) | 2017-07-24 | 2017-07-24 | Ddw高精度深度域井震匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107831542A CN107831542A (zh) | 2018-03-23 |
CN107831542B true CN107831542B (zh) | 2019-05-31 |
Family
ID=61643163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710608976.0A Active CN107831542B (zh) | 2017-07-24 | 2017-07-24 | Ddw高精度深度域井震匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831542B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109740427B (zh) * | 2018-11-26 | 2020-10-23 | 浙江财经大学 | 基于蓝噪声采样的标准井筛选可视分析方法 |
CN109932749B (zh) * | 2019-03-05 | 2020-12-01 | 中国石油天然气股份有限公司 | 一种井震标定方法、装置 |
CN110045416B (zh) * | 2019-04-29 | 2021-02-12 | 吉奥索特(北京)科技有限公司 | 一种跨三维地震工区的地震剖面数据的快速提取方法 |
CN110579806B (zh) * | 2019-09-11 | 2022-03-22 | 大庆油田有限责任公司 | 一种密井网条件下快速的井震标定方法 |
CN112649876A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 建立地震偏移速度模型的方法和装置 |
CN113514904B (zh) * | 2020-04-09 | 2023-02-28 | 中国石油天然气集团有限公司 | 地层参数模型建立方法及装置 |
CN112529981B (zh) * | 2020-11-24 | 2023-05-30 | 中国煤炭地质总局地球物理勘探研究院 | 一种深度域煤层底板等高线成图的方法 |
CN112965112B (zh) * | 2021-02-25 | 2022-06-10 | 中国海洋石油集团有限公司 | 地下储层微构造三维散点的生成方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8818729B1 (en) * | 2013-06-24 | 2014-08-26 | Hunt Advanced Drilling Technologies, LLC | System and method for formation detection and evaluation |
CN104237938B (zh) * | 2014-09-04 | 2017-01-25 | 中国石油天然气集团公司 | 利用动态时间调整进行纵波和转换波数据匹配的方法 |
CN106094019A (zh) * | 2016-04-12 | 2016-11-09 | 中国石油化工股份有限公司 | 基于地质信息映射的深度域地层结构反演方法 |
CN106610505B (zh) * | 2016-12-29 | 2019-03-22 | 中国石油大学(华东) | 一种基于dtw和aba联合的井震资料匹配方法 |
-
2017
- 2017-07-24 CN CN201710608976.0A patent/CN107831542B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107831542A (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107831542B (zh) | Ddw高精度深度域井震匹配方法 | |
CN111239802B (zh) | 基于地震反射波形和速度谱的深度学习速度建模方法 | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN103492910B (zh) | 时域中的同步小波提取和反卷积 | |
CN106291677B (zh) | 一种基于匹配追踪方法的叠后声波阻抗反演方法 | |
CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
CN104570067A (zh) | 一种地球物理勘探中相控地震反演方法 | |
CN102937720A (zh) | 井控提高地震资料分辨率的方法 | |
CN108445538A (zh) | 基于反射地震资料建立深度域层q模型的方法和系统 | |
CN113050157B (zh) | 一种基于露头资料的碳酸盐岩地震储层反演方法及系统 | |
CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
CN103119472B (zh) | 利用同时和顺序源方法进行全波形反演的混合方法 | |
US20150362608A1 (en) | Combined interpolation and primary estimation | |
US6519532B2 (en) | Method and apparatus for 3D depth migration | |
CN117546051A (zh) | 使用s波速度模型和机器学习进行地震成像的方法和系统 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
WO2021127382A1 (en) | Full waveform inversion in the midpoint-offset domain | |
Pafeng et al. | Prestack waveform inversion of three-dimensional seismic data—An example from the Rock Springs Uplift, Wyoming, USA | |
CN106199694A (zh) | 基于深变子波的合成记录制作方法 | |
CA2926179C (en) | System and method for seismic adaptive optics | |
CN110515122B (zh) | 正演网格搜索定位及微地震信号识别方法及装置 | |
RU2011148308A (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления | |
CN109143335B (zh) | 一种合成地震记录的制作方法、系统、介质及设备 | |
Kazemi et al. | Block row recursive least-squares migration | |
CN107561581B (zh) | 基于相关系数法建立多井模型的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |