发明内容
针对现有技术中存在的不足,本发明的目的之一在于解决上述现有技术中存在的一个或多个问题。例如,本发明的目的之一在于克服常规窄方位叠前地震裂缝预测技术存在的反演精度不高、空间分辨率不足等缺陷,提出了一种利用宽方位偏移距矢量片地震数据进行叠前方位各向异性裂缝参数(例如,裂缝发育密度、裂缝方位)预测方法。
为了实现上述目的,本发明的提供了一种利用偏移距矢量片地震叠前道集数据预测储层裂缝的方法。所述方法包括以下步骤:
A、对于经过叠前时间偏移处理的偏移距矢量片地震叠前道集数据,根据目的储层所在时间范围,选取裂缝预测时窗。
B、对于裂缝预测时窗内每一个成像点位置,计算偏移距矢量片地震叠前道集数据的偏移距和方位角。
C、对于裂缝预测时窗内每一个成像点位置的偏移距矢量片地震叠前道集数据,构建纵波反射系数R(θ)。
其中,θ为实际地震炮检观测方位角,R(θ)为实际地震炮检观测方位角θ对应的纵波反射系数,为裂缝介质对称轴的方位角,表示炮检方向和裂缝走向之间的夹角,P为各向同性介质下的反射系数,Q为与裂缝发育密度相关的各向异性反射系数。
D、对于裂缝预测时窗内每一个成像点位置的偏移距矢量片地震叠前道集数据,构建目标函数:
其中,θi为第i个偏移距矢量片地震叠前道(简称道)的实际地震炮检观测方位角,Rreal(θi)为第i个偏移距矢量片地震叠前道的实际地震炮检观测方位角θi对应的纵波反射系数,N为该成像点位置偏移距矢量片地震叠前道集数据的最大反射道数,θN为实际最大炮检观测方位角。
在E达到最小的情况下,反演求解得到储层的弹性参数值P、Q、及
E、利用所述步骤D的反演结果,计算得到裂缝密度e和裂缝方位角
在本发明的一个示例性实施例中,所述方法还可以包括步骤F:对所有成像点偏移距矢量片地震叠前道集数据,针对目标储层时窗范围,重复以上步骤A至步骤E,以得到整个三维工区地震数据体目标储层的裂缝密度e和裂缝方位角反演结果。
在本发明的一个示例性实施例中,其中,采用基于最小二乘解的计算方法来使步骤D中的E最小,求解方程组得到储层的弹性参数值P、Q、及
所述方程组为:
在本发明的一个示例性实施例中,在所述步骤B中,根据第i个偏移距矢量片地震叠前道的横向偏移距信息OFFx、纵向偏移距信息OFFy,计算得到所述第i个偏移距矢量片地震叠前道的偏移距OFFi和观测方位角θi,
在本发明的一个示例性实施例中,优选地,所述方法还可以包括在步骤B和C之间进行地步骤G:对裂缝预测时窗内给定成像点位置的偏移距矢量片道集数据进行部分叠加处理。所述步骤G可以包括:
G1、将给定成像点位置的偏移距矢量片地震叠前道集数据的方位角划分为X个扇区,每个扇区内的道数为Mm,m=1,2,…,X。
G2、对任一扇区m内的所有道进行叠加,得到所述任一扇区m内新生成地震道的偏移距OFFm_new和方位角信息θm_new。
其中,OFFj为第m个扇区内第j道的偏移距,θj为第m个扇区内第j道的方位角,j=1,2,…,Mm;G3、重复所述步骤G2,直到所有扇区都被处理。
在本发明的一个示例性实施例中,所述X大于3且小于道集最大覆盖次数。例如,X取值为18。
在本发明的一个示例性实施例中,将采集的三维宽方位地震数据在偏移距矢量片域进行处理,得到偏移距矢量片地震叠前道集数据。
在本发明的一个示例性实施例中,所述方法还包括将根据反演的储层裂缝发育密度和裂缝方位参数进行储层裂缝发育带预测。
在本发明的一个示例性实施例中,所述步骤(A)至(D)中,成像点位置可以是预测时窗内任意一个成像点的位置,也可以是预测时窗内的一段成像点范围,也可以是全部三维工区数据体的所有成像点范围,即该技术方案适用于三维工区中任意时窗范围、任意平面范围内偏移距矢量片叠前道集数据的反演。具体反演时窗的选取是根据具体地震工区目标储层的地质层位决定的。
与现有技术相比,本发明的有益技术效果包括:
(1)本申请采用宽方位偏移距矢量片叠前道集数据进行反演计算,该道集数据保留了更全面、完整的偏移距和方位角信息,通常其道集覆盖次数比常规窄方位地震道集数据的覆盖次数更高,克服了常规窄方位叠前地震裂缝预测技术存在的反演精度不高、空间分辨率不足等缺陷,为方位各向异性研究提供了更合适的地震数据,能够更精确、更高效地研究振幅、速度等地震响应随方位角变化的特征,更有利于识别断层、裂缝和储层气水变化特征。
(2)本申请所构建的目标函数表达式清晰简洁,易于推广应用,方程组求解采用了基于最小二乘解的计算方法,保证了数值计算的稳定性和收敛性,提高了计算效率和反演精度。
具体实施方式
在下文中,将结合示例性实施例详细地描述根据本发明的利用偏移距矢量片地震叠前道集数据预测储层裂缝的方法。
本发明所阐述的偏移距矢量片地震叠前道集数据这一概念,是目前石油、天然气勘探领域宽方位地震勘探技术中广泛采用的通用概念。所谓偏移距矢量片,是指在一个共接收线炮集中按炮线距和检波线距等距离划分成许多小矩形,每一个矩形就是一个偏移距矢量片。偏移距矢量片处理技术的优势在于:偏移距矢量片道集数据保留了更完整的偏移距和方位角信息,为方位各向异性研究提供更精确的地震数据,有利于研究速度、振幅等地震响应特征随方位角的变化特征,增强了断层、裂缝和地层岩性变化的识别能力。
本发明利用同时带有方位角和偏移距信息的偏移距矢量片地震叠前道集数据直接反演出储层裂缝发育密度和裂缝方位。在一个示例性实施例中,本发明利用偏移距矢量片地震叠前道集数据预测储层裂缝的方法包括如下步骤:
(1)以本领域公知的方法进行三维宽方位地震勘探采集,以得到适合进行宽方位处理的地震数据,然后,对采集得到的三维宽方位地震数据在偏移距矢量片(OVT)域进行处理,得到能够进行储层裂缝密度和裂缝方位反演的偏移距矢量片地震叠前道集数据,根据目的储层所在时间范围,选取裂缝预测时窗。
(2)导入经过叠前时间偏移处理的偏移距矢量片地震叠前道集数据,对于裂缝预测时窗内每一个成像点位置,其最大道数为N,根据第i个偏移距矢量片地震叠前道的横向偏移距信息OFFx、纵向偏移距信息OFFy,计算该偏移距矢量片道的偏移距OFFi和观测方位角θi,i=1,2,3,...,N,其计算公式如下:
(3)基于Ruger纵波HTI介质(横向各向同性介质,Horizontal TransverseIsotropic)反射系数公式,对其进行简化,建立如下的计算公式:
对该公式进行三角函数展开,得到如下形式:
其中θ为实际地震炮检观测方位角,R(θ)为实际地震炮检观测方位角θ对应的纵波反射系数,为裂缝介质对称轴的方位角,表示炮检方向和裂缝走向之间的夹角,P为各向同性介质下的反射系数,而Q为与裂缝发育密度相关的各向异性反射系数。在上述公式中,纵波反射系数随观测方位角的变化表现为椭圆特征。
(4)对于步骤(2)输入的裂缝预测时窗内每一个成像点位置的偏移距矢量片地震叠前道集数据,构建目标函数:
其中,θi为第i个偏移距矢量片地震叠前道的实际地震炮检观测方位角,Rreal(θi)为实际地震炮检观测方位角θi对应的纵波反射系数,N为该成像点位置偏移距矢量片地震叠前道集数据的最大反射道数,θN为实际最大炮检观测方位角。
在E达到最小的情况下,反演求解得到储层的弹性参数值P、Q、及
为使E达到最小值,采用基于最小二乘解的计算方法,求解如下方程组:
计算上述方程组的解,即可获得P、Q、及的值。
(5)计算裂缝密度及裂缝方位角。在方位各向异性介质中,上述步骤(3)中纵波反射系数随观测方位角的变化表现为椭圆特征,P+Q为椭圆特征的长轴,P-Q为椭圆特征的短轴,椭圆长轴方向即为裂缝走向方位角,长轴与短轴之比指示了裂缝发育强度。
利用步骤(4)的反演结果,分别计算椭圆长轴P+Q,椭圆短轴P-Q,则可计算得到裂缝密度e和裂缝方位角如下:
(6)对于提供的所有成像点位置的偏移距矢量片地震叠前道集数据,针对目标储层时窗范围,重复以上步骤(1)至步骤(5),则可得到整个三维工区地震数据体目标储层的裂缝密度e和裂缝方位角反演结果。然后,即可根据整个地震工区目标储层的裂缝密度和裂缝方位角反演结果进行裂缝性储层预测和识别。
本发明步骤(3)至步骤(5)中所阐述的纵波反射系数随观测方位角的变化规律及裂缝密度、裂缝方位计算方法,其理论基础为各向异性理论,即地震波在各向异性介质中平行或垂直裂缝方向传播时具有不同的传播特征。当射线平面与裂缝发育带方向平行时,纵波反射振幅最大;当射线平面与裂缝发育带方向垂直时,纵波反射振幅最小;其它方向,纵波反射振幅介于二者之间,其振幅近似为周期180度的正余弦曲线。
在上述步骤(2)至步骤(5)中,成像点位置可以是预测时窗内任意一个成像点的位置,也可以是预测时窗内的一段成像点范围,也可以是全部三维工区数据体的所有成像点范围,即该技术方案既适用于三维工区中任意时窗范围、任意平面范围内偏移距矢量片叠前道集数据的反演。具体反演时窗的选取是根据具体地震工区目标储层的地质层位决定的。
以下将结合具体示例来详细说明本发明的示例性实施例。
在本发明的一个具体示例中,利用偏移距矢量片地震叠前道集数据预测储层裂缝的方法可以由以下步骤实现:
(1)宽方位纵波地震数据偏移距矢量片(OVT)域叠前时间偏移处理。
以本领域公知的方法进行三维宽方位地震勘探采集,以得到适合进行宽方位处理的地震数据,然后,对采集得到的三维宽方位地震数据在偏移距矢量片(OVT)域进行处理,得到能够进行储层裂缝密度和裂缝方位反演的偏移距矢量片地震叠前道集数据。
(2)预测时窗选取。
导入经过叠前时间偏移处理的偏移距矢量片地震叠前道集数据,根据目标储层所在时间范围,选取裂缝预测时窗[Ta,Tb]。针对不同沉积特征的目标储层,预测时窗[Ta,Tb]具有不同的取值范围。如针对浅层碎屑岩沉积储层须家河组的储层裂缝预测,[Ta,Tb]分别代表须家河组顶、底反射界面的反射时间;又如针对深层碳酸盐岩沉积储层龙王庙组的储层裂缝预测,[Ta,Tb]分别代表龙王庙组顶、底反射界面的反射时间。
(3)对于预测时窗内给定成像点位置的偏移距矢量片地震叠前道集数据,计算偏移距和方位角信息。
导入经过叠前时间偏移处理的偏移距矢量片地震叠前道集数据,其最大道数为N,根据第i个偏移距矢量片地震叠前道的横向偏移距信息OFFx、纵向偏移距信息OFFy,计算该偏移距矢量片的偏移距OFFi和观测方位角(实际地震炮检观测方位角)θi,i=1,2,3,...,N,其计算公式如下:
(4)对给定成像点位置的偏移距矢量片地震叠前道集数据进行部分叠加处理。
由于原始偏移距矢量片(OVT)叠前道集数据通常信噪比较低,一般需要对道集数据进行部分叠加处理(所谓部分叠加处理,是指将道集分为若干个扇区,分别对每个扇区内的地震道进行叠加,即分扇区叠加。相对于传统意义上的道集全部叠加,该处理过程被称为道集分部分叠加),以提高道集的信噪比,进而提高后续反演算法的稳定性和精确度。具体实施过程如下:
假设给定成像点位置的偏移距矢量片(OVT)道集数据方位角范围为0度至180度,可以将方位角范围按照10度的间隔划分为18等份:0-10度为扇区1,10-20度为扇区2,…,170-180度为扇区18,则对应每一个扇区内的数据道数为Mm,m=1,2,…,18。
对于第m个扇区,将该扇区内的地震道进行叠加处理,则该扇区内新生成地震道的偏移距OFFm_new和方位角信息θm_new由如下公式计算:
m=1,2,…,18
其中,OFFj为第m扇区内第j道的偏移距,θj为第m个扇区内第j道的方位角,j=1,2,…,Mm。
上述公式计算得到的偏移距OFFm_new和方位角θm_new,即可作为新生成的部分叠加偏移距矢量片地震叠前道集数据的更新后实际观测偏移距信息和实际观测方位角信息,并对外输出,然后将已完成部分叠加的偏移距矢量片地震叠前道集数据作为后续反演计算的输入数据。
(5)基于Ruger纵波HTI介质反射系数公式,对其进行简化,建立如下的计算公式:
对该公式进行三角函数展开,得到如下形式:
其中,θ为实际地震炮检观测方位角,R(θ)为实际地震炮检观测方位角θ对应的纵波反射系数,为裂缝介质对称轴的方位角,表示炮检方向和裂缝走向之间的夹角,P为各向同性介质下的反射系数,而Q为与裂缝发育密度相关的各向异性反射系数。在上述公式中,纵波反射系数随观测方位角的变化表现为椭圆特征。
(6)对于步骤(4)输入的裂缝预测时窗内每一个成像点位置的偏移距矢量片地震叠前道集数据,构建目标函数:
其中,θi为实际地震炮检观测方位角,Rreal(θi)为实际地震炮检观测方位角θi对应的纵波反射系数,i=1,2,3,...,N,N为该成像点位置偏移距矢量片地震叠前道集数据的最大反射道数,θN为实际最大炮检观测方位角。在E达到最小的情况下,反演求解得到储层的弹性参数值P、Q、及
为使E达到最小值,采用基于最小二乘解的计算方法,求解如下方程组:
计算上述方程组的解,即可获得P、Q、及的值。
(7)利用步骤(4)的反演结果,分别计算椭圆长轴P+Q,椭圆短轴P-Q,则可计算得到裂缝密度e和裂缝方位角如下:
(8)对于提供的所有成像点偏移距矢量片地震叠前道集数据,针对目标储层时窗范围,重复以上步骤(1)至步骤(7),则可得到整个三维工区地震数据体目标储层的裂缝密度e和裂缝方位角反演结果。
具体地讲,要达到目标储层段裂缝预测目的,需要针对目标储层段反演时窗[Ta,Tb]内的所有成像点进行反演计算。如需针对浅层碎屑岩沉积储层须家河组进行储层裂缝预测,则需要针对须家河组顶、底反射界面的反射时间[Ta,Tb]内所有成像点进行反演计算;又如需针对深层碳酸盐岩沉积储层龙王庙组进行储层裂缝预测,则需要针对龙王庙组顶、底反射界面的反射时间[Ta,Tb]内所有成像点进行反演计算。然后,即可根据整个地震工区目标储层的裂缝密度和裂缝方位角反演结果进行裂缝性储层预测和识别。
综上所述,本发明的优点包括:采用了宽方位偏移距矢量片地震叠前道集数据进行反演计算,保留了更全面、完整的偏移距和方位角信息,为方位各向异性研究提供了更合适的地震数据,克服了常规窄方位叠前地震裂缝预测技术存在的反演精度不高、空间分辨率不足等缺陷,反演计算出的储层裂缝发育密度和裂缝方位信息更加精确、分辨率更高。另外,构建的目标函数表达式清晰简洁,易于推广应用,基于最小二乘解的反演算法,保证了数值计算的稳定性和收敛性,提高了计算效率和反演精度,反演结果更有利于断层、裂缝和储层气水变化特征的有效识别。
也就是说,本发明通过利用宽方位偏移距矢量片地震叠前道集数据能够直接反演出裂缝密度和裂缝方位,进而可实现更加准确的储层裂缝发育带预测。本发明针对复杂岩性油气藏勘探中裂缝性储层识别问题具有广泛的应用前景。
尽管上面已经通过结合示例性实施例描述了本发明,但是本领域技术人员应该清楚,在不脱离权利要求所限定的精神和范围的情况下,可对本发明的示例性实施例进行各种修改和改变。