法律状态公告日
法律状态信息
法律状态
2017-12-05
授权
授权
2015-12-23
实质审查的生效 IPC(主分类):E21B49/00 申请日:20150716
实质审查的生效
2015-11-25
公开
公开
技术领域
本发明涉及石油勘探测井技术领域,特别涉及一种连续判断快慢地层的方法。
背景技术
地层的横波速度大于井眼流体的速度时,此类地层被称为快地层,而当地层的横 波速度小于井眼流体的速度时,地层被称为慢地层。因为地层快慢属性不同时,评价 某些关键岩石力学参数的方法也不同,如最小水平主应力的计算等,特别是当地层具 有各向异性特征时,所以在连续深度范围内准确评价地层的快慢属性,即判断地层到 底为快地层还是慢地层,对于准确评价地层的岩石力学特征,进而为优化完井方案、 提高单井油气产量具有特别重要的意义。由于井眼流体的速度难以连续深度直接进行 测量,所以只能通过间接方法来连续深度判断地层的快慢属性。目前唯一可行的方法 是利用从阵列声波测井资料得到的横波速度径向变化剖面来实现这一目标。但这一方 法需从阵列声波测井资料中反演得到横波速度径向变化信息,而这一过程十分复杂繁 琐,很难为一般测井解释人员所掌握,并且如果每口井都采用这一方法来判断地层的 快慢属性,不仅效率低,而且成本高,不易进行工业化规模应用。因此油气勘探家们 迫切需要一种简便易行准确高效的方法来实现对地层快慢属性的判断。
发明内容
本发明实施例提供了一种连续判断快慢地层的方法,无需从阵列声波测井资料中 反演得到横波速度径向变化信息,只利用地层的泥质含量就可以判断地层的快慢属 性,使得判断地层快慢属性的过程变得简单、效率高、成本低且易进行工业化规模应 用。该方法包括:
利用采样井的横波速度径向变化剖面,连续判断所述采样井在不同地层深度的近 井壁地层的横波速度相对变化量是否为正值特征;所述近井壁地层为距离井筒中心线 小于阈值范围内的地层;
利用所述采样井的测井资料,连续计算所述采样井在不同地层深度的泥质含量;
在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层深度,统计 泥质含量的最大值;
利用所述泥质含量的最大值,连续判断所述采样井在不同地层深度的地层快慢属 性。
在一个实施例中,还包括:利用所述泥质含量的最大值和所述采样井的邻井在不 同地层深度的泥质含量,连续判断所述采样井的邻井在不同地层深度的地层快慢属 性。
在一个实施例中,所述利用采样井的横波速度径向变化剖面,连续判断所述采样 井在不同地层深度的近井壁地层的横波速度相对变化量是否为正值特征,具体按照如 下方式进行判断:
若在采样井的当前深度上每个地层径向位置上的横波速度相对变化量均大于等 于0,且随径向距离的减小而逐渐增加,则采样井的近井壁地层的横波速度相对变化 量确定为正值特征。
在一个实施例中,所述横波速度相对变化量按照如下公式计算:
其中,ΔVS为横波速度相对变化量;
VS′为近井壁地层中的横波速度,m/s;
VS为远离井壁的原状地层中的横波速度,m/s。
在一个实施例中,所述利用所述采样井的测井资料,连续计算所述采样井在不同 地层深度的泥质含量,包括:利用所述采样井的自然伽马测井资料,连续计算所述采 样井在不同地层深度的泥质含量,具体按照如下公式计算:
其中,Vsh为当前深度地层的泥质含量;
GR为当前深度地层的自然伽马测井值;
GRmin为纯砂岩段自然伽马特征值;
GRmax为纯泥岩段自然伽马特征值。
在一个实施例中,所述利用所述采样井的测井资料,连续计算所述采样井在不同 地层深度的泥质含量,包括:利用所述采样井的中子测井资料和密度测井资料,连续 计算所述采样井在不同地层深度的泥质含量,具体按照如下公式计算:
其中,Vsh为当前深度地层的泥质含量;
den_sh为干粘土的密度;
cn_sh为干粘土的中子值;
den_m为骨架矿物的密度;
cn_m为骨架矿物的中子值;
den_f为孔隙流体的密度;
cn_f为孔隙流体的中子值;
den为当前深度地层的密度;
cn当前深度地层的中子值。
在一个实施例中,所述利用所述泥质含量的最大值,连续判断所述采样井在不同 地层深度的地层快慢属性,具体按照如下方式进行判断:
若所述采样井的当前地层深度的泥质含量小于所述泥质含量的最大值,则将所述 采样井的当前地层深度的地层判断为快地层;
若所述采样井的当前地层深度的泥质含量大于所述泥质含量的最大值,则将所述 采样井的当前地层深度的地层判断为慢地层。
本发明实施例还提供了一种连续判断快慢地层的方法,无需从阵列声波测井资料 中反演得到横波速度径向变化信息,只利用地层的脆性指数就可以判断地层的快慢属 性,使得判断地层快慢属性的过程变得简单、效率高、成本低且易进行工业化规模应 用。该方法包括:
利用采样井的横波速度径向变化剖面,连续判断所述采样井在不同地层深度的近 井壁地层的横波速度相对变化量是否为正值特征;所述近井壁地层为距离井筒中心线 小于阈值范围内的地层;
利用所述采样井的密度测井资料、纵波时差资料及横波时差资料,连续计算所述 采样井在不同地层深度的脆性指数;
在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层深度,统计 脆性指数的最小值;
利用所述脆性指数的最小值,连续判断所述采样井在不同地层深度的地层快慢属 性。
在一个实施例中,还包括:利用所述脆性指数的最小值和所述采样井的邻井在不 同地层深度的脆性指数,连续判断所述采样井的邻井在不同地层深度的地层的快慢属 性。
在一个实施例中,所述利用采样井的横波速度径向变化剖面,连续判断所述采样 井在不同地层深度的近井壁地层的横波速度相对变化量是否为正值特征,具体按照如 下方式进行判断:
若在采样井的当前深度上每个地层径向位置上的横波速度相对变化量均大于等 于0,且随径向距离的减小而逐渐增加,则采样井的近井壁地层的横波速度相对变化 量确定为正值特征。
在一个实施例中,所述横波速度相对变化量按照如下公式计算:
其中,ΔVS为横波速度相对变化量;
VS′为近井壁地层中的横波速度,m/s;
VS为远离井壁的原状地层中的横波速度,m/s。
在一个实施例中,所述利用所述采样井的密度测井资料、纵波时差资料及横波时 差资料,连续计算所述采样井在不同地层深度的脆性指数,具体按照如下公式计算:
其中,Y为脆性指数;
E为杨氏模量,GPa;
Emin为杨氏模量最小值,常数;
Emax为杨氏模量最大值,常数;
ν为泊松比,小数;
νmin为泊松比最小值,常数;
νmax为泊松比最大值,常数;
按如下公式计算所述杨氏模量E:
其中,ρ为体积密度,g/cm3;
VS为远离井壁的原状地层中的横波速度,m/s;
Vp为纵波速度,m/s;
按如下公式计算所述泊松比ν:
在一个实施例中,所述利用所述脆性指数的最小值,连续判断所述采样井在不同 地层深度的地层快慢属性,具体按照如下方式进行判断:
若所述采样井的当前地层深度的脆性指数大于所述脆性指数的最小值,则将所述 采样井的当前地层深度的地层判断为快地层;
若所述采样井的当前地层深度的脆性指数小于所述脆性指数的最小值,则将所述 采样井的当前地层深度的地层判断为慢地层。
在本发明实施例中,在所述采样井的近井壁地层的横波速度相对变化量为正值特 征的地层深度内,统计泥质含量的最大值,利用所述泥质含量的最大值,连续判断所 述采样井在不同地层深度的地层快慢属性;或,在所述采样井的近井壁地层的横波速 度相对变化量为正值特征的地层深度内,统计脆性指数的最小值,利用所述脆性指数 的最小值,连续判断所述采样井在不同地层深度的地层快慢属性,与现有的方法相比, 本发明方法只需要利用地层的泥质含量或脆性指数就能判断出地层的快慢属性,使得 判断过程简单,效率高,成本低且易进行工业化规模应用。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不 构成对本发明的限定。在附图中:
图1是本发明实施例中提供的一种连续判断快慢地层征的方法流程示意图;
图2是本发明实施例中提供的一种连续判断快慢地层征的方法流程示意图;
图3是本发明实施例中提供的研究区一口资料齐全井X井的综合成果图;
图4是本发明实施例中提供的研究区中另一口井Y井的综合成果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图, 对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发 明,但并不作为对本发明的限定。
现有的判断地层的快慢属性的方法,需要从阵列声波测井资料中反演得到横波速 度径向变化信息,该过程复杂繁琐很难被一般测井解释人员掌握,且若每口井均采用 该方法来判断地层的快慢属性,效率低、成本高,不易进行工业化规模应用。如果可 以提出一种方法,利用从常规测井资料得到的泥质含量或脆性指数就可以判断地层的 快慢属性,就可以克服现有方法中存在的问题。基于此,本发明提出一种连续判断快 慢地层的方法。
图1是本发明实施例中提供的一种连续判断快慢地层的方法流程示意图;如图1 所示,该方法包括如下步骤:
步骤101:利用采样井的横波速度径向变化剖面,连续判断所述采样井在不同地 层深度的近井壁地层的横波速度相对变化量是否为正值特征;所述近井壁地层为距离 井筒中心线小于阈值范围内的地层;
步骤102:利用所述采样井的测井资料,连续计算所述采样井在不同地层深度的 泥质含量;
步骤103:在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层 深度,统计泥质含量的最大值;
步骤104:利用所述泥质含量的最大值,连续判断所述采样井在不同地层深度的 地层快慢属性。
具体实施时,在研究区中选择一口测井资料齐全的井作为采样井,其中,所说的 测井资料齐全指的是该井的测井资料至少包括常规测井资料如自然伽玛、自然电位、 井径、声波、中子、密度、深浅电阻率等,及新技术测井资料如阵列声波测井资料等。
利用从该井阵列声波测井资料中得到横波速度径向变化剖面,连续判断近井壁地 层的横波速度相对变化量是否为正值特征,其中,所说的横波速度径向变化剖面是指, 在从井筒与地层的交界处到距离井筒中心线1米远处的地层范围内,依次有2L个地 层位置(在地层的纵向剖面上,左半侧有L个,右半侧有L个),以井筒中心线为参 考点,其与这些地层位置的距离由小变大呈等间隔增加,被称为地层径向距离,这些 地层位置被称为地层径向位置;由于受钻头破坏等因素的影响,当地层径向位置从距 参考点最远处向最近处移动时,每个径向位置上对应的横波速度相对变化量逐渐增加 或减小。在从井筒中心线到井筒与地层交界处的井筒区域内,有M个空间位置,这 些位置与井筒中心线的距离也呈等间隔分布,被称为井筒径向位置。连续深度范围内, 所有地层径向位置上的横波速度相对变化量被统称为横波速度径向变化剖面。
所说的近井壁地层为距离井筒中心线小于阈值范围内的地层;此处的阈值为1 米。
所说的横波速度相对变化量是指,近井筒(井壁)附近地层中的横波速度VS′与 远离井筒的原状地层(距离井筒中心线1米范围之外的地层)中的横波速度VS有一 定差异,这种差异用横波速度相对变化量ΔVS来表示,计算公式如下:
所述利用从该井阵列声波测井资料得到的横波速度径向变化剖面,连续判断采样 井在不同地层深度的(即连续深度判断)近井壁地层的横波速度相对变化量是否为正 值特征是指,假设在当前深度(即设定的某一深度)上每个地层径向位置上的横波速 度相对变化量均大于等于0,且随径向距离的减小而逐渐增加,此时称当前深度靠近 井壁地层的横波速度相对变化量为正值特征,当前深度地层为快地层,令指示曲线 FAST在当前深度的取值为1;若当前深度上每个地层径向位置上的横波速度相对变 化量均小于等于0,且随径向距离的减小而逐渐减小,此时称当前深度靠近井壁地层 的横波速度相对变化量为负值特征,当前深度地层为慢地层,令指示曲线FAST在当 前深度的取值为0。根据这一原则来考察横波速度径向变化剖面在每个深度每个地层 径向位置上的横波速度相对变化量的数值大小及相对变化关系,连续深度对指示曲线 FAST进行赋值,形成非0即1的连续曲线,利用该曲线连续深度直观指示靠近井壁 地层的横波速度相对变化量是否为正值特征。
具体实施时,在判断完近井壁地层的横波速度相对变化量是否为正值特征后,需 要在为正值特征的深度段内统计泥质含量的最大值(上限)。地层的泥质含量的计算 方法有多种,本发明可以采用以下两种方法来计算,其中一种是利用自然伽马测井资 料按如下公式来计算地层的泥质含量:
式(2)中,Vsh为当前深度地层的泥质含量;GR为当前深度地层的自然伽马测 井值;GRmin为纯砂岩段自然伽马特征值;GRmax为纯泥岩段自然伽马特征值。
另外一种可以利用其它常规测井资料来计算地层的泥质含量,如可以利用中子、 密度测井资料按如下公式来计算:
式(3)中:Vsh为地层泥质含量;den_sh、cn_sh分别为干粘土的密度、中子 值;den_m、cn_m分别为骨架矿物的密度、中子值;den_f、cn_f分别为孔隙 流体的密度、中子值。这些值为所在地层的常数值,均为已知;den为当前深度地层 的密度;cn当前深度地层的中子值,这两个值需要测量。
按照上述两种方法中的任一种(即公式(2)或(3))计算出不同地层深度的泥 质含量之后,统计泥质含量的最大值。具体的,即在指示曲线FAST取值为1的全部 深度点内,统计采样井的泥质含量的最大值。
在找到采样井的泥质含量的最大值之后,利用泥质含量的最大值,和按照上述公 式(2)或(3)计算出的采样井在不同地层深度的泥质含量,判断采样井在不同地层 深度的地层快慢属性。具体的地层快慢属性的判断规则是:若采样井的当前深度地层 的泥质含量小于所述泥质含量的最大值,则将采样井的当前深度地层判断为快地层, 令指示曲线FASTN在当前深度的取值为1;若采样井的当前深度地层的泥质含量大于 所述泥质含量的最大值,则将采样井的邻井的当前深度地层判断为慢地层,令指示曲 线FASTN在当前深度的取值为0。其它各深度点都依此规则来对地层的快慢属性进行 判别,对指示曲线FASTN进行赋值。
同样的,还可以利用统计到的采样井的泥质含量的最大值,和按照上述公式(2) 或(3)计算的采样井的邻井在不同地层深度的泥质含量,按照上述的地层快慢属性 的判断规则,来判断采样井的邻井在不同地层深度的地层快慢属性。
图2是本发明实施例中提供的一种连续判断快慢地层的方法流程示意图;如图2 所示,该方法包括如下步骤:
步骤201:利用采样井的横波速度径向变化剖面,连续判断所述采样井在不同地 层深度的近井壁地层的横波速度相对变化量是否为正值特征;所述近井壁地层为距离 井筒中心线小于阈值范围内的地层;
步骤202:利用所述采样井的密度测井资料、纵波时差资料及横波时差资料,连 续计算所述采样井在不同地层深度的脆性指数;
步骤203:在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层 深度,统计脆性指数的最小值;
步骤204:利用所述脆性指数的最小值,连续判断所述采样井在不同地层深度的 地层快慢属性。
具体实施时,步骤201的具体实施过程与步骤101的具体实施过程相同,重复之 处不再赘述。
然后,在判断完近井壁地层的横波速度相对变化量是否为正值特征后,需要在为 正值特征的深度段内统计脆性指数的最小值(下限)。地层的脆性指数的计算方法可 以有多种,本发明方法利用密度测井资料、纵波时差资料及横波时差资料,连续计算 (连续深度计算)所述采样井在不同地层深度的脆性指数,按照如下公式来计算:
式(4)中,Y为脆性指数;
E为杨氏模量,GPa;
Emin为杨氏模量最小值,常数;
Emax为杨氏模量最大值,常数;
ν为泊松比,小数;
νmin为泊松比最小值,常数;
νmax为泊松比最大值,常数。
杨氏模量E的计算公式为:
式中,ρ为体积密度,g/cm3;
VS为远离井壁的原状地层中的横波速度,m/s;
Vp为纵波速度,m/s.
泊松比ν的计算公式为:
按照上述方法(即公式(4)、(5)和(6))计算出不同地层深度的脆性指数之后, 统计脆性指数的最小值。具体的,在指示曲线FAST取值为1的全部深度点,统计脆 性指数的最小值。
在找到脆性指数的最小值之后,利用脆性指数的最小值,和按照上述公式(4)、 (5)和(6)计算出的采样井在不同地层深度的脆性指数,判断采样井在不同地层深 度的地层快慢属性。具体的地层快慢属性的判断规则是:若采样井的当前深度地层的 脆性指数大于所述脆性指数的最小值,则将采样井的当前深度地层判断为快地层,令 指示曲线FASTN在当前深度的取值为1;若采样井的当前深度地层的脆性指数小于所 述脆性指数的最小值,则将采样井的当前深度地层判断为慢地层,令指示曲线FASTN 在当前深度的取值为0。其它各深度点都依此规则来对地层的快慢属性进行判别,对 指示曲线FASTN进行赋值。
同样的,还可以利用统计到的采样井的脆性指数的最小值,和按照上述公式(4)、 (5)和(6)计算采样井的邻井在不同地层深度的脆性指数,按照上述的地层快慢属 性的判断规则,来判断采样井的邻井在不同地层深度的地层的快慢属性。
下面结合两种判断快慢地层的方法,以实例来说明本发明两种方法的优点,同时 便于理解如何实施本发明。
以某油田某一研究区中的一口测井资料齐全的井X井及其邻井Y井为例,来具 体说明本发明的技术方案。
第一步,在X井中,利用从X井的阵列声波测井资料得到的横波速度径向变化 剖面连续深度判断近井壁地层的横波速度相对变化量是否为正值特征,判别结果如图 3中第6道所示,第6道为指示曲线FAST,若当前深度近井壁地层的横波速度相对 变化量为正值特征,则曲线值为1,若当前深度近井壁地层横波速度相对变化量为负 值特征,则曲线值为0。其中,图3中第2道为深度道,指示地层的连续深度DEPTH; 第3道为声波(AC,单位微秒/米,用US/M表示)、中子(CNL)、密度(DEN,单 位克/立方厘米,用G/C3表示)测井曲线;第5道为横波速度径向变化剖面(用VSF 表示)。
第二步,按照第一种方法计算泥质含量。由X井常规测井资料连续深度计算地 层的泥质含量,采用本步骤中采用公式(2)来计算泥质含量,即利用图3中第1道 显示的自然伽马测井曲线GR(单位GAPI),同时第1道显示的还有井径曲线(HCAL, 单位厘米CM),计算结果如图3中第9道所示(用VSH表示泥质含量曲线)。其中, 纯砂岩段自然伽马特征值GRmin为60;纯泥岩段自然伽马特征值GRmax为150。
第三步,按照第二种方法计算脆性指数。根据X井的密度测井曲线及由阵列声 波测井资料提供的纵、横波时差曲线计算脆性指数,即利用的是图3中第4道显示的 是纵波时差曲线(AC,单位微秒/米,用US/M表示)、横波时差曲线(DTSM,单位 微秒/米,用US/M表示),计算结果如图3中第8道所示(用Y表示脆性指数曲线)。
第四步,在X井中近井壁地层横波速度相对变化量为正值特征的深度段,统计 泥质含量最大值(上限)及脆性指数最小值(下限)。本步骤中,在第6道指示曲线 FAST取值为1的全部深度点上统计得到的泥质含量最大值(上限)为0.6,脆性指数 最小值(下限)为45%。
第五步,利用X井的泥质含量最大值(上限)或脆性指数最小值(下限)连续 深度判断X井的地层的快慢属性。在X井中,判别结果如图3中第7道的快慢地层 指示曲线FASTN所示,该判别结果是将脆性指数最小值(下限)45%作为判别依据 而得到的,若当前深度地层的脆性指数大于脆性指数最小值(下限)45%,则当前深 度地层为快地层,则指示曲线FASTN的值为1,若当前深度的脆性指数小于脆性指数 最小值(下限值)45%,则当前深度地层地层为慢地层,则指示曲线FASTN的值为0。
利用Y井来验证本发明方法的有效性。在判断Y井所在的连续地层的快慢属性 时,以X井中统计得到的脆性指数最小值(下限值)45%为判别依据。图4为Y井的 综合成果图。其中,第1道为自然伽马测井曲线;第2道为深度道,指示地层的连续 深度;第3道为声波(AC)、中子(CNL)、密度(DEN)测井曲线;第4道显示的 是纵、横波时差曲线;第5道为脆性指数曲线;第6道为泥质含量;第7道为快慢地 层指示曲线FASTN,显示的深度段为快地层。利用现有的从阵列声波测井资料得到横 波速度径向变化剖面来判断Y井所在的连续地层的快慢属性,如图4中第8道所示 为横波速度径向变化剖面,第9道为根据横波速度径向变化剖面得到的快慢地层指示 曲线FAST。比较第7道快慢地层指示曲线FASTN和第9道快慢地层指示曲线FAST, 可知两条曲线非常吻合。因此,与现有方法相比,利用本发明方法同样能判断地层的 快慢属性,且本发明方法过程简单。
综上所述,本发明实施例提供的技术方案,在所述采样井的近井壁地层的横波速 度相对变化量为正值特征的地层深度内,统计泥质含量的最大值,利用所述泥质含量 的最大值,连续判断所述采样井在不同地层深度的地层快慢属性;或,在所述采样井 的近井壁地层的横波速度相对变化量为正值特征的地层深度内,统计脆性指数的最小 值,利用所述脆性指数的最小值,连续判断所述采样井在不同地层深度的地层快慢属 性。与现有的方法相比,本发明方法使得判断地层快慢属性的过程变得简单、效率高 和成本低。
同样的还可以将采样井的泥质含量的最大值,和采样井邻井的连续深度地层的泥 质含量进行比较,以此来判断采样井的邻井的地层的快慢属性;或将采样井的脆性指 数的最小值,和采样井邻井的连续深度地层的脆性指数进行比较,与现有的方法相比, 本发明方法易进行工业化规模应用;该方法为提高地层快慢属性的判别效率及其工业 化应用、优选试油层位、优化完井方案奠定坚实基础。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以 用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算 装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而, 可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于 此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或 者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例 不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技 术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内, 所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
机译: 在地下地层的井眼上进行压裂操作的方法,对在地下地层中具有储层的井进行增产作业的方法,对在地下地层中具有储层的井眼进行增产作业的方法,用于对具有穿透地下地层的钻井的井场进行增产作业的系统和系统,该地下地层在其中具有不连续性
机译: 一种存在性判断装置,一种中等判别装置,存在性判断方法和一种中等判别方法
机译: 用水基钻井液通过粉笔地层稳定并在粘土地层中钻一个钻孔的井眼中稳定井眼的方法。一种在地下地层中钻井孔的钻井液,以及一种增强井眼穿透稳定性的方法孔洞黏土