具体实施方式
下面将结合附图,对本发明实施例进行描述。
实施例
图1是本发明一种储层裂缝的确定方法,该方法包括以下步骤:
S101:获取储层裂缝属性及分布;
此步骤中,根据所述实际地震资料类型对所述实际地震资料进行多种方法裂缝预测分析,获取储层裂缝属性及分布;
并且在此步骤前还可进一步包括:接收微电阻率扫描成像FMI资料及同一测井旁的实际地震资料。具体可通过以下步骤实现:在同一测井上进行两类数据的采集,一类为测井资料,一类为地震资料,其中,测井资料可通过测井勘测直接获得,地震资料可通过井旁地震资料采集给予实现。
在此步骤中还进一步包括:当所述实际地震数据为叠前数据时,根据实际地震资料进行多种属性各向异性裂缝分析或基于叠前反演的应力场分析获取储层构造成因裂缝属性及分布;当所述实际地震数据为叠后数据时,根据实际地震资料及实际地质数据进行应力场分析或协同建模叠后反演裂缝预测获取储层裂缝属性及分布。
上述当所述实际地震数据为叠前数据时,根据实际地震资料进行多种裂缝各向异性分析或基于叠前反演的应力场分析获取多个储层裂缝属性及分布;通过实际地震资料进行方位各向异性椭圆拟合;以上拟合椭圆的过程具 体包括:根据实际地震资料提取道集数据;根据所述道集数据选取三个以上方位角道集;提取所述三个以上方位角道集的方位角地震属性;根据所述方位角地震属性进行方位椭圆拟合。所述道集数据可具体为CMP共中心点道集数据,从而增大数据采集的精度,其中方位角地震属性包括:振幅、层速度、时差、方位AVO梯度、方位层频率、层频率差、叠加振幅及叠加振幅方位差。所述多种属性各向异性裂缝分析还包括:振幅各向异性分析、频率各向异性分析及衰减各向异性分析。所述通过实际地震资料进行方位各向异性椭圆拟合后步骤还包括:提取所述裂缝拟合椭圆获取长轴、短轴方向矢量值及长度值信息;根据所述长轴方向矢量值确定裂缝走向,根据所述长轴及短轴长度值获取裂缝密度,并通过所述裂缝走向及裂缝密度确定裂缝属性。图2为利用叠前地震资料提取方位地震属性分析裂缝的流程图,在属性研究的基础上,进行各向异性椭圆拟合,实现对裂缝方位、裂缝密度等参数的描述,椭圆短长轴之比c/a指示裂缝密度,椭圆长轴方向指示裂缝走向,椭圆短轴方向指示裂缝法向(图3、图4)。同时根据地震扁率可以划分裂缝的级别,扁率即椭圆才长短轴之比a/c。以上过程通常可采用裂缝储层综合描述软件LD-FRSTMfracture V3.2给予实现。
同时,还可通过应力场进行分析,所述应力场分析包括:常规应力场分析。对于常规应力场分析通常根据实际地质数据几何信息、岩性信息获取主应力方向;所述根据实际地质数据几何信息、岩性信息获取主应力方向步骤还包括:根据所述采集地层的几何信息、岩性信息获取地层面的曲率张量、变形张量和应力场张量;根据所述曲率张量、变形张量和应力场张量获取主应力方向。
所述根据所述曲率张量、变形张量和应力场张量获取主应力方向步骤还包括:根据所述曲率张量、变形张量和应力场张量及以下公式获取主应力方向:
其中,
所述采集地层的几何信息具体为:所述采集地层的地表三维坐标值ux,uy,uz,其中三维坐标z向与海拔方向一致。所述岩性信息具体为:地层速度、地层密度。
首先,应力场是指地壳中或地球体内,应力状态随空间点的变化,称为应力场,或构造应力场。应力场一般随时间变化,但在一定地质阶段相对比较稳定。研究应力场,就是研究应力分布的规律性,确定地壳上某一点或某一地区,在特定地质时代和条件下,受力作用所引起的应力方向、性质、大小以及发展演化等特征。随着地质演化,一个地区常常经受多次不同方式的地壳运动,导致同一地区内,呈现出受不同时期不同形式地应力场作用所形成的各种构造极其叠加或改造的复杂景观。因此,只有最近一期地质构造事件,未经破坏或改造,才能确切地反映这个时期的应力场,同时,地质模型的建立是做好应力场模拟的先决条件,首先将储层的目的层连同上下盖层和覆盖层作为一个岩石块体的隔离体来计算,然后从地质的角度提出构造成因,构造裂缝的特征,构造应力场的宏观特征及断层发育史。在此基础上,恢复古构造剖面图,推断地质隔离体的受力方向及大小,设定边界条件并提出反演应力场及裂缝的地质标准。但是,由于地质体是一个十分复杂的地下岩石块体,其地壳中各种地质构造形态,类型,成因是在漫长历史时期的地质演化过程中形成的,这种复杂的地质演化过程不可能恢复,只能用相对静止的观点和相对简化的方法去处理构造与古应力场的问题。考虑到储油构造中主要的一种类型是背斜构造,它是由地壳受到挤压发生弯曲,或由于基底隆起使沉积地层上拱而形成。在从历史时期的连续的接受应力而形成现代的构造体系,我们可以近似的以现在的构造应力场来作为古构造应力场的一个发展模型,来模拟在这样的应力场的条件下会导致怎样的构造裂缝体系。
以弯曲薄板作为油层构造模拟的力学模型,用二维的方法来处理油层构造,这种方法计算方便,人工干预少,对于应力场模拟可以将边界作为自由边界处理,不需要另外考虑边界条件。该方法主要以背斜构造作为模型进行分析,得出可以用构造面上一点的最大曲率值作为该点裂缝发育程度的判据, 而以最小主曲率方向指示可能出现的张裂缝走向,这样就将构造裂缝的分布问题化为构造面的主曲率计算问题。
设以薄板中面为z=0的坐标面,规定按右手规则,以平行于大地坐标为X,Y坐标,以向上为正。沿X,Y正方向的位移分别为ux,uy,沿Z方向的位移为扰度w(x,y)。在直角坐标系中:
a.基本方程
直角坐标系中的变形几何方程:
uz=w
由薄板理论可知,有:
且有:
定义曲率变形分量:
因此,应变分量可写为:
εx=-zκx,εy=-zκy,γxy=-2zκxy
b.物理本构关系(广义虎克定律)
其逆关系为:
σx=2Gεx+λθ,τxy=Gγxy
σy=2Gεy+λθ,τyz=Gγyz
σz=2Gεz+λθ,τxz=Gγxz
θ=εkk
λ,G为拉梅(Lame)常数,G也就是剪切模量(Shear modulus),E杨氏模量(Yong modulus)。θ为体积应变。
将前面的式代入,得到:
因而有:
将地层厚度t=2z代入上式,得到由曲率分量表示的地层面上的应力分量:
由上式可知,当地层面向上凸时,曲率大于零,正好对应上凸地层面受拉张应力,张应力为正。
为了与地质力学符号相符,这里采用压应力为正,张应力为负的符号约定。曲率小于零,表示地层上凸。
求出该点的沿坐标的应力后,就可求出其主应力及其方向:
σmax与X轴的夹角α,σmin与X轴的夹角β:
因此,若能得到地层面的扰度方程或其面上点的曲率,就可以估算其上的应力场,进而计算由此应力产生的裂缝。
c.地层曲率的计算
由前面理论可知,若能求出地层面的曲率分量,就可以求出其上的应力场。采用趋势面拟合方法拟合地层面的趋势函数,进而计算其上点的曲率分量。采用最小二乘法拟合趋势面。设趋势面的待定系数的函数为:
w(x,y)=a0+a1x+a2y+a3x2+a4xy+a5y2
由层面散点处的坐标值(x,y,z),,建立最小二乘方程,对一个散点:
εi=zi-wi(xi,yi)
当用n个散点拟合一个趋势面时,可得到拟合方程组,解此线性方程组,就可得到趋势面函数。
趋势面的曲率计算式:
如图5所示:应力场应用流程,裂缝储层综合描述软件LD-FRSTMfractureV3.2,系统提供的应力场数值模拟,考虑了储层岩石的厚度、岩性,并且考虑了储层受构造控制的裂缝分布等因素,基于弹性薄板理论,在计算构造的曲率分布基础之上,进一步计算构造的应力、应变分布,然后根据构造的应力应变场,就可以对储层裂缝的发育程度及展布关系进行分析。构造曲率:表示构造面梯度变化的快慢;最大主应变:表示形变的大小;张应变(+值):与裂缝密度有关;压应变(-值):表示底层压实变形;最大主应力:压应力(+值):平行于裂缝方向;张应力(-值):垂直于裂缝方向。从而通过以上方式确定裂缝属性及分布。
S102:获取匹配度及吻合度;
此步骤中,根据实际地质参数进行正演,根据比对正演模拟结果及所述储层裂缝特征,获取所述储层裂缝分布的匹配度;比对FMI资料及所述储层裂缝属性,获取所述储层裂缝属性的吻合度;
所述根据实际地质参数进行正演步骤还包括:收集井旁的实际地质参数;
根据所述实际地质参数获取地质结构和储层参数;根据所述地质结构和储层参数定义地质模型参数;根据所述地质结构参数及实际地质参数获取地质模型参数;采用波动方程正演;根据所述波动方程利用所述地质模型参数进行地震波正演模拟获取储层裂缝正演结果。所述波动方程包括下属:声学介质波动方程、弹性介质波动方程、粘弹性介质波动方程、各向异性(裂隙)介质波动方程:其中:
1、声学介质波动方程
或
其中n为声波速度,r为密度,P为声压
2、弹性介质波动方程
其中l、m为拉梅常数,r为密度,ui为位移分量
3、粘弹性介质波动方程
3.1、运动方程
3.2、本构关系(应力应变关系)
3.3、记忆变量的一阶方程
其中:
当模型为3D时,n=3;为2D时,n=2。弹性常数λ、μ是松弛的或低频拉梅系数;
是未松弛的或高频拉梅系数。常数
和
对应于P波品质因子Qp的松弛时间;
和
对应于S波品质因子Qs的松弛时间。借助于SLS模型,其关系为:
4、各向异性(裂隙)介质波动方程
4.1、运动方程
4.2、本构关系(应力应变关系)
σij=cijklekl i,j,k,l=1,2,3
4.3、几何关系(应变位移关系)
其中:r为密度;cijkl为介质的弹性常数,介质为具有垂直对称轴的VTI介质时,其值为:
介质为具有水平对称轴的HTI介质时,其值为:
所述比对正演模拟结果及所述储层裂缝特征,获取所述储层裂缝特征的匹配度,步骤还包括:将所述正演结果的裂缝特征与所述储层裂缝预测结果比对,已重叠部分为匹配分布,根据匹配分布占所述储层裂缝分布的比例获取匹配度。
所述比对FMI资料及所述储层裂缝属性,获取所述储层裂缝属性的吻合度,步骤还包括:根据所述根据FMI资料获取裂缝特征,根据所述裂缝特征比对所述储层各向异性属性获取差异值,根据所述差异量值获取吻合度。
S103:确定储层裂缝;
此步骤中,根据所述匹配度及吻合度从所述多种方法的裂缝预测分析中优选一种裂缝分析方法,根据所述优选的裂缝分析方法,解释所述实际地震资料,确定储层裂缝。
所述根据所述匹配度及吻合度从所述多种方法预测裂缝的分析中优选一种裂缝分析方法,步骤还包括:根据所述匹配度从高到低得到储层多个裂缝分布特征,根据所得到的多个裂缝分布特征的吻合度选取最高吻合度的裂缝分布,并根据该分布特征从所述多种方法预测裂缝的分析中选取一种裂缝分析方法。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可以通过硬件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上公开的仅为本发明的几个具体实施例,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。