著录项信息
专利名称 | 裂缝预测方法和装置 |
申请号 | CN201010205983.4 | 申请日期 | 2010-06-12 |
法律状态 | 授权 | 申报国家 | 中国 |
公开/公告日 | 2010-12-08 | 公开/公告号 | CN101907725A |
优先权 | 暂无 | 优先权号 | 暂无 |
主分类号 | G01V1/28 | IPC分类号 | G01V1/28;G01V1/36;E21B49/00查看分类表>
|
申请人 | 中国石油天然气股份有限公司 | 申请人地址 | 北京市东城区东直门北大***
变更
专利地址、主体等相关变化,请及时变更,防止失效 |
权利人 | 中国石油天然气股份有限公司 | 当前权利人 | 中国石油天然气股份有限公司 |
发明人 | 刘军迎;杨午阳;高建虎;雍学善;李三家;王恩利 |
代理机构 | 北京三友知识产权代理有限公司 | 代理人 | 陶海萍 |
摘要
本发明实施例提供一种裂缝预测方法和装置,该方法包括:获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数;根据叠前共中心点动校道集地震数据生成对应的索引数据,并根据叠后地震数据解释获取层位数据;将叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据;根据常规宽方位叠前共中心点动校道集数据获取每个地震道的相关道头信息;利用获取的道头信息计算每个地震道的方位角;以目标层为拾取中心,利用拾取时窗获取每个地震道的反射振幅;利用获取的方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度。通过本发明实施例,可提高裂缝储层预测的精度;并且,不仅可对裂缝储层实现三维预测,还可进行二维预测,提高探井和钻井的成功率。
裂缝预测方法和装置
技术领域
[0001] 本发明涉及裂缝预测技术,特别涉及一种裂缝预测方法和装置。
背景技术
[0002] 保存石油和天然气的地层叫含油气储层,它们多为具有一定孔隙度的地层,其中,含有大量裂缝或以裂缝为主的储层就是裂缝型储层如泥岩裂缝储层、砂岩裂缝储层、碳酸盐岩裂缝储层、火山岩裂缝储层等,形成的油气藏就是裂缝型油气藏。
[0003] 石油勘探的目的就是综合运用地球物理勘探方法、地球重/磁力勘探方法、地球物理测井勘探方法、地质勘探方法、地球化学勘探方法等来寻找地下储集了石油和天然气的储层,这就是广义上的储层预测。其中,实际应用最广的方法是综合运用各种地球物理勘探方法来定性或定量预测储层的方法,这就是通常意义下的储层预测。
[0004] 储层预测在石油勘探与开发工业中占有相当重要的地位,储层预测的主要内容包括储层分布的平面位置、地下深度、空间展布范围、空间展布形态、速度、密度、孔隙度、渗透率、裂缝方向与密度等。储层预测是利用测井、地震信息在地质理论指导下,对油气储层的空间展布和几何形态进行宏观描述,对储层的微观特征进行横向预测的一门技术。
[0005] 目前,进行储层裂隙预测方法,大概可分为以下几种方式:
[0006] 1)横波、转换波及多分量检测裂缝;2)多方位垂直地震剖面法(VSP,Vertical Seismic Prospection)检测裂缝;3)纵波分析检测裂缝。
[0007] 在发明人实现本发明的过程中,发明人发现上述方式的缺陷在于:采用横波、转换波及多分量检测裂缝时,横波信噪比低、频率低,多波采集和处理成本太高,转换波的处理比较麻烦;采用多方位VSP检测裂缝时,采集、处理麻烦、成本高,作用范围有限;
[0008] 而采用纵波分析检测裂缝时,由于裂缝介质是强烈的各向异性介质,即为横向各向同性(HTI,Horizontal Transverse Isotropic)介质或大角度倾斜各向同性(THI,Tilt Transverse Isotropic)介质,需要采用全方位或宽方位地震勘探技术来研究。目前的地震勘探从本质上来说是基于VTI介质的层状介质理论,主要利用的是叠后资料来研究。该叠后技术的缺点是所利用的叠后资料没有方位信息和偏移距信息,从理论上讲,常规勘探技术缺乏各向异性研究的理论依据。
发明内容
[0009] 本发明实施例提供一种裂缝预测方法和装置,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
[0010] 本发明实施例提供一种裂缝预测方法,所述方法包括:
[0011] 获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数,所述相关参数包括:纵测线面元大小、横测线面元大小和拾取时窗;
[0012] 根据所述叠前共中心点动校道集地震数据选取目标层,并生成对应的索引数据;
所述索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数和最大覆盖次数;
[0013] 根据所述叠前共中心点动校道集地震数据获得叠后地震数据,并利用所述叠后地震数据进行地震层位解释,以生成层位数据;所述层位数据包括:纵测线号、横测线号、顶界面时间和底界面时间;
[0014] 利用获取的所述纵测线面元大小、横测线面元大小以及索引数据,将所述叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据;
[0015] 根据所述常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;所述道头信息包括炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标和偏移距;
[0016] 利用获取的所述炮点横坐标、炮点纵坐标、检波点横坐标和检波点纵坐标计算每个地震道的方位角;
[0017] 以所述目标层为拾取中心,利用所述拾取时窗获取每个地震道的反射振幅;
[0018] 利用获取的所述方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度;
[0019] 其中,利用获取的所述纵测线面元大小、横测线面元大小以及索引数据,将所述常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据,具体包括:
[0020] 根据所述纵测线面元大小和横测线面元大小,取用所述常规窄方位叠前共中心点动校道集数据相邻的多个共中心点道集;
[0021] 根据所述索引数据对所述相邻的多个共中心点道集统一附以面元中心点的共中心点道集的纵测线号和横测线号,以使面元内所有数据具有统一的纵测线号和横测线号,生成所述常规宽方位叠前共中心点动校道集数据。
[0022] 本发明实施例提供一种裂缝预测装置,所述装置包括:
[0023] 数据获取单元,用于获取目标层的叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数,所述相关参数包括:纵测线面元大小、横测线面元大小和拾取时窗;
[0024] 第一生成单元,用于根据所述叠前共中心点动校道集地震数据选取目标层,并生成对应的索引数据;所述索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数和最大覆盖次数;
[0025] 第二生成单元,用于根据所述叠前共中心点动校道集地震数据获得叠后地震数据,并利用所述叠后地震数据进行地震层位解释,以生成层位数据;所述层位数据包括:纵测线号、横测线号、顶界面时间和底界面时间;
[0026] 第三生成单元,用于利用获取的所述纵测线面元大小、横测线面元大小以及索引数据,将所述叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据;
[0027] 信息获取单元,用于根据所述常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;所述道头信息包括炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标和偏移距;
[0028] 方位角计算单元,用于利用获取的所述炮点横坐标、炮点纵坐标、检波点横坐标和检波点纵坐标计算每个地震道的方位角;
[0029] 振幅获取单元,用于以所述目标层为拾取中心,利用所述拾取时窗获取每个地震道的反射振幅;
[0030] 拟合单元,用于利用获取的所述方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度;
[0031] 其中,所述第三生成单元具体包括:
[0032] 取用单元,用于根据所述纵测线面元大小和横测线面元大小,取用所述常规窄方位叠前共中心点动校道集数据相邻的多个共中心点道集;
[0033] 附加单元,用于根据所述索引数据对所述相邻的多个共中心点道集统一附以面元中心点的共中心点道集的纵测线号和横测线号,以使面元内所有数据具有统一的纵测线号和横测线号,生成所述常规宽方位叠前共中心点动校道集数据。
[0034] 本发明实施例的有益效果在于,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
附图说明
[0035] 此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
[0036] 图1是本发明实施例1的裂缝预测方法的流程图;
[0037] 图2是本发明实施例2的裂缝预测方法的流程图;
[0038] 图3是本发明实施例2的叠前共中心点动校道集地震数据的实例图;
[0039] 图4是本发明实施例2的层位数据中目标层的顶界面时间构造图;
[0040] 图5是本发明实施例2的层位数据中目标层的底界面时间构造图;
[0041] 图6是本发明实施例2的地震层位解释剖面图;
[0042] 图7是本发明实施例2的拾取时窗的示意图;
[0043] 图8是本发明实施例2的椭圆拟合的模拟结果实例图;
[0044] 图9是本发明实施例2的偏移距——方位角图的实例图;
[0045] 图10是本发明实施例2的inline线与crossline线间隔为3的裂缝密度预测效果图;
[0046] 图11是本发明实施例2的inline线与crossline线间隔均为1的裂缝密度预测效果图;
[0047] 图12是本发明实施例2的预测的裂缝密度曲线图;
[0048] 图13是本发明实施例2的预测的裂缝方向曲线图;
[0049] 图14是本发明实施例2的裂缝密度预测结果平面图;
[0050] 图15是本发明实施例3的裂缝预测装置的构成图;
[0051] 图16是本发明实施例4的裂缝预测装置的构成图;
[0052] 图17是本发明实施例4的第三生成单元的构成图;
[0053] 图18是本发明实施例4的振幅获取单元的构成图。
具体实施方式
[0054] 为使本发明的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例作进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
[0055] 实施例1
[0056] 本发明实施例提供一种裂缝预测方法,如图1所示,该方法包括:
[0057] 步骤101,获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数,相关参数包括:纵测线面元大小、横测线面元大小和拾取时窗;
[0058] 步骤102,根据叠前共中心点动校道集地震数据获取目标层,并生成对应的索引数据;索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数、最大覆盖次数;
[0059] 步骤103,根据叠前共中心点动校道集地震数据获得叠后地震数据,并利用叠后地震数据进行地震层位解释,以生成层位数据;
[0060] 步骤104,利用获取的纵测线面元大小、横测线面元大小以及索引数据,将叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据;
[0061] 步骤105,根据常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;道头信息包括炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标和偏移距;
[0062] 步骤106,利用获取的炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标计算每个地震道的方位角;
[0063] 步骤107,以目标层为拾取中心,利用拾取时窗获取每个地震道的反射振幅;
[0064] 步骤108,利用获取的方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度。
[0065] 在本实施例中,叠前共中心点动校道集(NMO-CMP,Normal MoveOut-Common Middle Point)地震数据可采用SEG-IBM格式,可包括道头信息:纵测线道头(Inline)、横测线道头(Crossline)、炮点横坐标道头(SOU_X)、炮点纵坐标道头(SOU_Y)、检波点横坐标道头(REC_X)、检波点纵坐标道头(REC_Y)、偏移距道头(Offset);但不限于此,可根据实际情况确定具体的信息。
[0066] 在本实施例中,可获取进行裂缝预测的相关参数,可根据经验值指定该相关参数,相关参数可包括:纵测线面元大小、横测线面元大小和拾取时窗;但不限于此,还可根据实际情况确定具体的参数。
[0067] 在本实施例中,可根据叠前共中心点动校道集地震数据生成对应的索引数据;索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数、最大覆盖次数,但不限于此,还可根据实际情况确定具体的数据。
[0068] 为节约磁盘空间、提高效率,可先对叠前数据进行分块处理;然后,为每一块叠前数据生成相应的索引数据,该索引数据可采用文本格式,以数据库文件的形式存放;
[0069] 例如,如下表1所示,索引数据可按列存放,可分为5列:
[0070] 第1列为linline线号;
[0071] 第2列为crossline线号,即共中心点动校道集(CDP)的序号,可简称为CDP号;
[0072] 第3列为CDP道集中首道在本数据体中的道序号;
[0073] 第4列为CDP道集中尾道在本数据体中的道序号;
[0074] 第5列为本CDP道集中地震道的总道数。
[0075] 表1
[0076]
[0077] 在本实施例中,还可根据叠前共中心点动校道集地震数据获得叠后地震数据,并利用叠后地震数据进行地震层位解释,以生成对应的层位数据;可通过分块处理后,按照每块叠前数据的inline线范围和Xline线范围准备相应的层位数据。
[0078] 例如,地震数据可分为如下2块:
[0079] cmp_L2531-2580.sgy,观测线范围为2531——2580;
[0080] cmp_L2581-2630.sgy,观测线范围为2581——2630;
[0081] 按照观测线的范围可分为相应的部分,其中,第1列为inline线号、第2列为crossline线号、第3列为顶界面时间、第4层为底界面时间;
[0082] 第1部分:线范围2531——2580的层位数据,可如下表2所示:
[0083] 表2
[0084]
[0085]
[0086] 第2部分:线范围2581——2630的层位数据,可如下表3所示:
[0087] 表3
[0088]
[0089] 在本实施例中,叠前共中心点动校道集地震数据是方位角很小或者等于零度的地震数据,为常规窄方位叠前共中心点动校道集数据;为了预测裂缝地层,可利用获取的纵测线面元大小、横测线面元大小以及索引数据将叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据。
[0090] 在本实施例中,可根据常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;道头信息包括炮点横坐标(SOU_X)、炮点纵坐标(SOU_Y)、检波点横坐标(REC_X)、检波点纵坐标(REC_Y)和偏移距(Offset)。
[0091] 在本实施例中,可利用获取的炮点横坐标道头、炮点纵坐标道头、检波点横坐标道头、检波点纵坐标道头计算每个地震道的方位角;计算公式可如下所示:
[0092] 方位角=atan[(SOU_X-REC_X)/(SOU_Y-REC_Y)]
[0093] 在本实施例中,以目标层为拾取中心,可利用拾取时窗获取每个地震道的反射振幅。
[0094] 在本实施例中,利用获取的方位角和反射振幅可进行椭圆拟合,以确定裂缝方向和裂缝密度;椭圆拟合公式可为:
[0095] F(θ,φ)=A(θ)+B(θ)cos2φ
[0096] 其中,F(θ,φ)表示纵波的反射振幅,θ是入射角,φ是相对裂缝方向的方位角,A,B是和方位角无关的系数。
[0097] 在本实施例中,通过选定偏移距,相当于固定了θ;通过代入拾取的反射振幅及其方位角,进行椭圆拟合,可求出长轴方向,得到裂缝的方向;可计算椭圆长轴与短轴的比,得到裂缝密度。
[0098] 由上述实施例可知,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
[0099] 实施例2
[0100] 本发明实施例提供一种裂缝预测方法,在实施例1的基础上对该裂缝预测方法进行进一步说明,如图2所示,该方法包括:
[0101] 步骤201,获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数,所述相关参数包括:纵测线面元大小、横测线面元大小和拾取时窗;
[0102] 在本实施例中,该相关参数还可包括:纵测线增量、横测线增量、最小偏移距、最大偏移距、偏移距面元大小、方位角面元大小和平滑参数;
[0103] 其中,纵测线增量和横测线增量可用于控制工作量,控制inline线数和crossline线数;最小偏移距和最大偏移距可用于控制所利用的叠前数据范围;偏移距面元大小和方位角面元大小可用于控制叠前数据的分组;平滑参数可用于预测结果的平滑。
通过上述相关参数,可使裂缝预测更加准确。相关参数不限于此,还可根据具体实施时的实际情况确定其他参数。
[0104] 图3为叠前共中心点动校道集地震数据的实例图,如图3所示,可对叠前共中心点动校道集地震数据进行浏览,可通过大量的观察、分析选取部分道集数据。
[0105] 步骤202,从叠前共中心点动校道集地震数据中选取部分数据后,进行迭代参数试验处理,以校正相关参数;
[0106] 在本实施例中,可对叠前共中心点动校道集地震数据进行分块,然后选定一小块数据,进行迭代处理,在预测结果与地质理论解释结果吻合时确定相关参数,以校正该相关参数。
[0107] 步骤203,根据叠前共中心点动校道集地震数据选取目标层,并生成对应的索引数据;所述索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数、最大覆盖次数;
[0108] 步骤204,根据叠前共中心点动校道集地震数据获得叠后地震数据,并利用叠后地震数据进行地震层位解释,以生成层位数据;
[0109] 在本实施例中,可对叠前共中心点动校道集地震数据进行叠加,共中心点道集数据在叠加前为来自于同一地下反射点的许多道,叠加后多道叠加到一起变为一道,由此获得叠后地震数据;
[0110] 在本实施例中,可利用叠后地震数据来解释要预测的目标层,进行叠后地震数据的解释,获得高信噪比的地震层位数据;层位数据可包括:纵测线号、横测线号、顶界面时间和底界面时间。
[0111] 例如,图4是层位数据中目标层的顶界面时间构造图,图5是层位数据中目标层的底界面时间构造图,图6为地震层位解释剖面图。以上仅对层位数据进行示意性说明,并不限于此。
[0112] 在本实施例中,还可对层位数据进行校正;一般情况下目标层的底界面在时间上不能小于顶界面,如果小于,则称为底界面刺穿了顶界面。通过校正,可使底界面位于顶界面之下,使层位数据更加准确。
[0113] 步骤205,根据纵测线面元大小和横测线面元大小,取用常规窄方位叠前共中心点动校道集数据相邻的多个共中心点道集。
[0114] 步骤206,根据索引数据对相邻的多个共中心点道集统一附以面元中心点的共中心点道集的纵测线号和横测线号,以使面元内所有数据具有统一的纵测线号和横测线号,生成常规宽方位叠前共中心点动校道集数据。
[0115] 步骤207,根据常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;道头信息包括炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标和偏移距。
[0116] 步骤208,利用获取的炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标计算每个地震道的方位角。
[0117] 步骤209,以目标层为拾取中心,分别向上、向下延拓半个新时窗;其中,新时窗为拾取时窗加上目标层的时间厚度的1/6至1/5。
[0118] 在本实施例中,目标层的时间厚度可通过目标层的顶界面时间和底界面时间获得,较佳地,新时窗可为拾取时窗加上目标层的时间厚度的1/5。
[0119] 步骤210,在新时窗内拾取每个地震道的反射振幅。
[0120] 图7为拾取时窗的示意图,如图7所示,可以层位数据中的目标层为拾取中心,在拾取时窗内拾取反射振幅。进一步地,还可在新时窗内获取反射振幅,这样,可扩大拾取范围,以防止动校不准而引起的目标层反射波拾取的遗漏。
[0121] 步骤211,利用获取的所述方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度。
[0122] 图8为椭圆拟合的模拟结果实例图,如图8所示,以正北为方位零度,可得到理论模拟结果。以上仅为示意性说明,但不限于此,可根据实际情况确定具体的模拟结果。
[0123] 在本实施例中,可进行偏移距、方位角参数的选取,通过获取的所述方位角和偏移距生成偏移距——方位角数据,如下表4所示:
[0124] 表4
[0125]
[0126] 表4为偏移距——方位角数据的一个实例,仅进行示意性说明,并不限于此。通过偏移距——方位角数据可得到偏移距——方位角图,然后可通过该偏移距——方位角图选取偏移距、方位角参数,这样,可使裂缝预测的精度更高,预测效果更好。
[0127] 图9为偏移距——方位角图的一个实例。选取参数的原则可为:各象限均衡考虑,若各象限的偏移距范围不同,可取几个平均值,所取偏移距范围尽量多地包含最宽广的方位角变化范围。例如,对图9所示的偏移距——方位角图进行分析,可选取偏移距范围为
200米至2300米,对应的方位角范围可为-80度至+20度。
[0128] 此外,还可对测线间隔进行试验,以进一步提高预测精度。例如,图10为inline线与crossline线间隔为3的裂缝密度预测效果图,图11为inline线与crossline线间隔均为1的裂缝密度预测效果图。
[0129] 在图10中,数据抽稀预测时,结果图上出现了横向条带分布现象,它与断层相交,这是不正常的;而在图11中,采用逐线逐点来预测时,预测结果没有横向条带分布现象,预测结果顺着断层分布,结果自然,预测效果较好。并且,前者的运算量是后者的三分之一,运算速度是后者的3倍,但效果不佳。因此,通过上述分析,可选择inline线与crossline线间隔为1。
[0130] 在本实施例中,可根据确定的参数进行椭圆拟合,以确定裂缝方向和裂缝密度。可根据预测结果数据,获得裂缝预测结果图、以及二维曲线图;
[0131] 例如,图12为预测的裂缝密度曲线图,图13为预测的裂缝方向曲线图,图14为裂缝密度预测结果平面图。如图14所示,沿着断层走向,预测出的裂缝密度最强,远离断层的区域,裂缝密度较弱。
[0132] 由上述实施例可知,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
[0133] 实施例3
[0134] 本发明实施例提供一种裂缝预测装置,如图15所示,该装置包括:数据获取单元
1501、第一生成单元1502、第二生成单元1503、第三生成单元1504、信息获取单元1505、方位角计算单元1506、振幅获取单元1507和拟合单元1508;其中,
[0135] 数据获取单元1501用于获取目标层的叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数,相关参数包括:纵测线面元大小、横测线面元大小和拾取时窗;
[0136] 第一生成单元1502用于根据叠前共中心点动校道集地震数据获取目标层,并生成对应的索引数据和层位数据;索引数据包括纵测线号、横测线号、起始道号、结束道号、覆盖次数、最大覆盖次数;
[0137] 第二生成单元1503用于根据叠前共中心点动校道集地震数据获得叠后地震数据,并利用叠后地震数据进行地震层位解释,以生成层位数据;层位数据包括:纵测线号、横测线号、顶界面时间和底界面时间;
[0138] 第三生成单元1504利用获取的纵测线面元大小、横测线面元大小以及索引数据,将叠前共中心点动校道集地震数据由常规窄方位叠前共中心点动校道集数据生成常规宽方位叠前共中心点动校道集数据;
[0139] 信息获取单元1505用于根据常规宽方位叠前共中心点动校道集数据获取每个地震道的道头信息;道头信息包括炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标和偏移距;
[0140] 方位角计算单元1506用于利用获取的炮点横坐标、炮点纵坐标、检波点横坐标、检波点纵坐标计算每个地震道的方位角;
[0141] 振幅获取单元1507用于以目标层为拾取中心,利用拾取时窗获取每个地震道的反射振幅;
[0142] 拟合单元1508用于利用获取的所述方位角和反射振幅进行椭圆拟合,以确定裂缝方向和裂缝密度。
[0143] 在本实施例中,该裂缝预测装置的工作流程可如实施例1所述,此处不再赘述。
[0144] 由上述实施例可知,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
[0145] 实施例4
[0146] 本发明实施例提供一种裂缝预测装置,如图16所示,该装置包括:数据获取单元
1601、第一生成单元1602、第二生成单元1603、第三生成单元1604、信息获取单元1605、方位角计算单元1606、振幅获取单元1607和拟合单元1608;如实施例3所述,此处不再赘述。
[0147] 在本实施例中,相关参数还可包括:纵测线增量、横测线增量、最小偏移距、最大偏移距、偏移距面元大小、方位角面元大小和平滑参数;但不限于此,还可根据实际情况确定具体的相关参数。
[0148] 如图16所示,该装置还可包括:校正单元1609;
[0149] 校正单元1609用于从叠前共中心点动校道集地震数据中选取部分数据后,进行迭代参数试验处理,以校正相关参数。
[0150] 如图17所示,第三生成单元1604具体可包括:取用单元1701和附加单元1702;
其中,
[0151] 取用单元1701用于根据纵测线面元大小和横测线面元大小,取用常规窄方位叠前共中心点动校道集数据相邻的多个共中心点道集;
[0152] 附加单元1702用于根据索引数据对相邻的多个共中心点道集统一附以面元中心点的共中心点道集的纵测线号和横测线号,以使面元内所有数据具有统一的纵测线号和横测线号,生成常规宽方位叠前共中心点动校道集数据。
[0153] 如图18所示,振幅获取单元1607具体可包括:延拓单元1801和拾取单元1802;
其中,
[0154] 延拓单元1801用于以目标层为拾取中心,分别向上、向下延拓半个新时窗;其中,新时窗为拾取时窗加上目标层的时间厚度的1/6至1/5;
[0155] 拾取单元1802用于在新时窗内拾取每个地震道的反射振幅。
[0156] 在本实施例中,该裂缝预测装置的工作流程可如实施例2所述,此处不再赘述。
[0157] 由上述实施例可知,通过获取叠前共中心点动校道集地震数据以及进行裂缝预测的相关参数后,进行椭圆拟合,可提高裂缝储层预测的精度;并且,可对裂缝储层实现三维预测的同时进行二维预测,提高探井和钻井的成功率。
[0158] 以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
引用专利(该专利引用了哪些专利)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 该专利没有引用任何外部专利数据! |
被引用专利(该专利被哪些专利引用)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 该专利没有被任何外部专利所引用! |