地层判断为慢地层，令指示曲线FASTN在当前深度的取 值为0。其它各深度点都依此规则来对地层的快慢属性进行判别，对指示曲线FASTN进行赋 值。
[0095] 同样的，还可以利用统计到的采样井的泥质含量的最大值，和按照上述公式（2) 或（3)计算的采样井的邻井在不同地层深度的泥质含量，按照上述的地层快慢属性的判断 规则，来判断采样井的邻井在不同地层深度的地层快慢属性。
[0096] 图2是本发明实施例中提供的一种连续判断快慢地层的方法流程示意图；如图2 所示，该方法包括如下步骤：
[0097] 步骤201 :利用采样井的横波速度径向变化剖面，连续判断所述采样井在不同地 层深度的近井壁地层的横波速度相对变化量是否为正值特征；所述近井壁地层为距离井筒 中心线小于阈值范围内的地层；
[0098] 步骤202 :利用所述采样井的密度测井资料、纵波时差资料及横波时差资料，连续 计算所述采样井在不同地层深度的脆性指数；
[0099] 步骤203 :在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层 深度，统计脆性指数的最小值；
[0100] 步骤204 :利用所述脆性指数的最小值，连续判断所述采样井在不同地层深度的 地层快慢属性。
[0101] 具体实施时，步骤201的具体实施过程与步骤101的具体实施过程相同，重复之处 不再赘述。
[0102] 然后，在判断完近井壁地层的横波速度相对变化量是否为正值特征后，需要在为 正值特征的深度段内统计脆性指数的最小值（下限）。地层的脆性指数的计算方法可以有 多种，本发明方法利用密度测井资料、纵波时差资料及横波时差资料，连续计算（连续深度 计算）所述采样井在不同地层深度的脆性指数，按照如下公式来计算：
[0103] (4)
[0104] 式⑷中，Y为脆性指数；
[0105] E为杨氏模量，GPa;
[0106] E_为杨氏模量最小值，常数；
[0107] E_为杨氏模量最大值，常数；
[0108] v为泊松比，小数；
[0109] 乂_为泊松比最小值，常数；
[0110] vnax为泊松比最大值，常数。
[0111] 杨氏模量E的计算公式为：
[0112]
(5)
[0113] 式中，P为体积密度，g/cm3;
[0114] Vs为远离井壁的原状地层中的横波速度，m/s ;
[0115] Vp为纵波速度，m/s.
[0116] 泊松比v的计算公式为：
[0117]
(6)
[0118] 按照上述方法（即公式（4)、（5)和（6))计算出不同地层深度的脆性指数之后，统 计脆性指数的最小值。具体的，在指示曲线FAST取值为1的全部深度点，统计脆性指数的 最小值。
[0119] 在找到脆性指数的最小值之后，利用脆性指数的最小值，和按照上述公式（4)、（5) 和（6)计算出的采样井在不同地层深度的脆性指数，判断采样井在不同地层深度的地层快 慢属性。具体的地层快慢属性的判断规则是：若采样井的当前深度地层的脆性指数大于所 述脆性指数的最小值，则将采样井的当前深度地层判断为快地层，令指示曲线FASTN在当 前深度的取值为1 ;若采样井的当前深度地层的脆性指数小于所述脆性指数的最小值，则 将采样井的当前深度地层判断为慢地层，令指示曲线FASTN在当前深度的取值为0。其它各 深度点都依此规则来对地层的快慢属性进行判别，对指示曲线FASTN进行赋值。
[0120] 同样的，还可以利用统计到的采样井的脆性指数的最小值，和按照上述公式（4)、 (5)和（6)计算采样井的邻井在不同地层深度的脆性指数，按照上述的地层快慢属性的判 断规则，来判断采样井的邻井在不同地层深度的地层的快慢属性。
[0121] 下面结合两种判断快慢地层的方法，以实例来说明本发明两种方法的优点，同时 便于理解如何实施本发明。
[0122] 以某油田某一研究区中的一 口测井资料齐全的井X井及其邻井Y井为例，来具体 说明本发明的技术方案。
[0123] 第一步，在X井中，利用从X井的阵列声波测井资料得到的横波速度径向变化剖面 连续深度判断近井壁地层的横波速度相对变化量是否为正值特征，判别结果如图3中第6 道所示，第6道为指示曲线FAST，若当前深度近井壁地层的横波速度相对变化量为正值特 征，则曲线值为1，若当前深度近井壁地层横波速度相对变化量为负值特征，则曲线值为0。 其中，图3中第2道为深度道，指示地层的连续深度DEPTH ;第3道为声波（AC，单位微秒/ 米，用US/M表示）、中子（CNL)、密度（DEN，单位克/立方厘米，用G/C3表示）测井曲线；第 5道为横波速度径向变化剖面（用VSF表示）。
[0124] 第二步，按照第一种方法计算泥质含量。由X井常规测井资料连续深度计算地层 的泥质含量，采用本步骤中采用公式（2)来计算泥质含量，即利用图3中第1道显示的自然 伽马测井曲线GR(单位GAPI)，同时第1道显示的还有井径曲线（HCAL，单位厘米CM)，计算 结果如图3中第9道所示（用VSH表示泥质含量曲线）。其中，纯砂岩段自然伽马特征值 GR_S 60 ;纯泥岩段自然伽马特征值GR_为150。
[0125] 第三步，按照第二种方法计算脆性指数。根据X井的密度测井曲线及由阵列声波 测井资料提供的纵、横波时差曲线计算脆性指数，即利用的是图3中第4道显示的是纵波时 差曲线（AC，单位微秒/米，用US/M表示）、横波时差曲线（DTSM，单位微秒/米，用US/M表 示），计算结果如图3中第8道所示（用Y表示脆性指数曲线）。
[0126] 第四步，在X井中近井壁地层横波速度相对变化量为正值特征的深度段，统计泥 质含量最大值（上限）及脆性指数最小值（下限）。本步骤中，在第6道指示曲线FAST取 值为1的全部深度点上统计得到的泥质含量最大值（上限）为〇. 6,脆性指数最小值（下 限）为45%。
[0127] 第五步，利用X井的泥质含量最大值（上限）或脆性指数最小值（下限）连续深 度判断X井的地层的快慢属性。在X井中，判别结果如图3中第7道的快慢地层指示曲线 FASTN所示，该判别结果是将脆性指数最小值（下限）45%作为判别依据而得到的，若当前 深度地层的脆性指数大于脆性指数最小值（下限）45%，则当前深度地层为快地层，则指示 曲线FASTN的值为1，若当前深度的脆性指数小于脆性指数最小值（下限值）45%，则当前 深度地层地层为慢地层，则指示曲线FASTN的值为0。
[0128] 利用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，可知两条曲线非常吻合。因此，与现有方法相比，利用本发明方法同 样能判断地层的快慢属性，且本发明方法过程简单。
[0129] 综上所述，本发明实施例提供的技术方案，在所述采样井的近井壁地层的横波速 度相对变化量为正值特征的地层深度内，统计泥质含量的最大值，利用所述泥质含量的最 大值，连续判断所述采样井在不同地层深度的地层快慢属性；或，在所述采样井的近井壁地 层的横波速度相对变化量为正值特征的地层深度内，统计脆性指数的最小值，利用所述脆 性指数的最小值，连续判断所述采样井在不同地层深