;
[0062] Vp为纵波速度，m/s;
[0063] 按如下公式计算所述泊松比v :
[0064]
[0065] 在一个实施例中，所述利用所述脆性指数的最小值，连续判断所述采样井在不同 地层深度的地层快慢属性，具体按照如下方式进行判断：
[0066] 若所述采样井的当前地层深度的脆性指数大于所述脆性指数的最小值，则将所述 采样井的当前地层深度的地层判断为快地层；
[0067] 若所述采样井的当前地层深度的脆性指数小于所述脆性指数的最小值，则将所述 采样井的当前地层深度的地层判断为慢地层。
[0068] 在本发明实施例中，在所述采样井的近井壁地层的横波速度相对变化量为正值特 征的地层深度内，统计泥质含量的最大值，利用所述泥质含量的最大值，连续判断所述采样 井在不同地层深度的地层快慢属性；或，在所述采样井的近井壁地层的横波速度相对变化 量为正值特征的地层深度内，统计脆性指数的最小值，利用所述脆性指数的最小值，连续判 断所述采样井在不同地层深度的地层快慢属性，与现有的方法相比，本发明方法只需要利 用地层的泥质含量或脆性指数就能判断出地层的快慢属性，使得判断过程简单，效率高，成 本低且易进行工业化规模应用。
【附图说明】
[0069] 此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不 构成对本发明的限定。在附图中：
[0070] 图1是本发明实施例中提供的一种连续判断快慢地层征的方法流程示意图；
[0071] 图2是本发明实施例中提供的一种连续判断快慢地层征的方法流程示意图；
[0072] 图3是本发明实施例中提供的研究区一口资料齐全井X井的综合成果图；
[0073] 图4是本发明实施例中提供的研究区中另一口井Y井的综合成果图。
【具体实施方式】
[0074] 为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对 本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并 不作为对本发明的限定。
[0075] 现有的判断地层的快慢属性的方法，需要从阵列声波测井资料中反演得到横波速 度径向变化信息，该过程复杂繁琐很难被一般测井解释人员掌握，且若每口井均采用该方 法来判断地层的快慢属性，效率低、成本高，不易进行工业化规模应用。如果可以提出一种 方法，利用从常规测井资料得到的泥质含量或脆性指数就可以判断地层的快慢属性，就可 以克服现有方法中存在的问题。基于此，本发明提出一种连续判断快慢地层的方法。
[0076] 图1是本发明实施例中提供的一种连续判断快慢地层的方法流程示意图；如图1 所示，该方法包括如下步骤：
[0077] 步骤101 :利用采样井的横波速度径向变化剖面，连续判断所述采样井在不同地 层深度的近井壁地层的横波速度相对变化量是否为正值特征；所述近井壁地层为距离井筒 中心线小于阈值范围内的地层；
[0078] 步骤102 :利用所述采样井的测井资料，连续计算所述采样井在不同地层深度的 泥质含量；
[0079] 步骤103 :在所述采样井的近井壁地层的横波速度相对变化量为正值特征的地层 深度，统计泥质含量的最大值；
[0080] 步骤104 :利用所述泥质含量的最大值，连续判断所述采样井在不同地层深度的 地层快慢属性。
[0081] 具体实施时，在研究区中选择一口测井资料齐全的井作为采样井，其中，所说的测 井资料齐全指的是该井的测井资料至少包括常规测井资料如自然伽玛、自然电位、井径、声 波、中子、密度、深浅电阻率等，及新技术测井资料如阵列声波测井资料等。
[0082] 利用从该井阵列声波测井资料中得到横波速度径向变化剖面，连续判断近井壁地 层的横波速度相对变化量是否为正值特征，其中，所说的横波速度径向变化剖面是指，在从 井筒与地层的交界处到距离井筒中心线1米远处的地层范围内，依次有2L个地层位置（在 地层的纵向剖面上，左半侧有L个，右半侧有L个），以井筒中心线为参考点，其与这些地层 位置的距离由小变大呈等间隔增加，被称为地层径向距离，这些地层位置被称为地层径向 位置；由于受钻头破坏等因素的影响，当地层径向位置从距参考点最远处向最近处移动时， 每个径向位置上对应的横波速度相对变化量逐渐增加或减小。在从井筒中心线到井筒与地 层交界处的井筒区域内，有M个空间位置，这些位置与井筒中心线的距离也呈等间隔分布， 被称为井筒径向位置。连续深度范围内，所有地层径向位置上的横波速度相对变化量被统 称为横波速度径向变化剖面。
[0083] 所说的近井壁地层为距离井筒中心线小于阈值范围内的地层；此处的阈值为1 米。
[0084] 所说的横波速度相对变化量是指，近井筒（井壁）附近地层中的横波速度、'与 远离井筒的原状地层（距离井筒中心线1米范围之外的地层）中的横波速度％有一定差 异，这种差异用横波速度相对变化量AVs来表示，计算公式如下：
[0085]
⑴
[0086] 所述利用从该井阵列声波测井资料得到的横波速度径向变化剖面，连续判断采样 井在不同地层深度的（即连续深度判断）近井壁地层的横波速度相对变化量是否为正值特 征是指，假设在当前深度（即设定的某一深度）上每个地层径向位置上的横波速度相对变 化量均大于等于〇,且随径向距离的减小而逐渐增加，此时称当前深度靠近井壁地层的横波 速度相对变化量为正值特征，当前深度地层为快地层，令指示曲线FAST在当前深度的取值 为1 ;若当前深度上每个地层径向位置上的横波速度相对变化量均小于等于〇,且随径向距 离的减小而逐渐减小，此时称当前深度靠近井壁地层的横波速度相对变化量为负值特征， 当前深度地层为慢地层，令指示曲线FAST在当前深度的取值为0。根据这一原则来考察横 波速度径向变化剖面在每个深度每个地层径向位置上的横波速度相对变化量的数值大小 及相对变化关系，连续深度对指示曲线FAST进行赋值，形成非0即1的连续曲线，利用该曲 线连续深度直观指示靠近井壁地层的横波速度相对变化量是否为正值特征。
[0087] 具体实施时，在判断完近井壁地层的横波速度相对变化量是否为正值特征后，需 要在为正值特征的深度段内统计泥质含量的最大值（上限）。地层的泥质含量的计算方法 有多种，本发明可以采用以下两种方法来计算，其中一种是利用自然伽马测井资料按如下 公式来计算地层的泥质含量：
[0088]
.(.2.)
[0089] 式（2)中，Vsh为当前深度地层的泥质含量；GR为当前深度地层的自然伽马测井值； GR_S纯砂岩段自然伽马特征值；GR _为纯泥岩段自然伽马特征值。
[0090] 另外一种可以利用其它常规测井资料来计算地层的泥质含量，如可以利用中子、 密度测井资料按如下公式来计算：
[0091]
[0092] 式⑶中：Vsh为地层泥质含量；den_sh、cn_sh分别为干粘土的密度、中子值；den_ m、cn_m分别为骨架矿物的密度、中子值；den_f、cn_f分别为孔隙流体的密度、中子值。这 些值为所在地层的常数值，均为已知；den为当前深度地层的密度；cn当前深度地层的中子 值，这两个值需要测量。
[0093] 按照上述两种方法中的任一种（即公式（2)或（3))计算出不同地层深度的泥质 含量之后，统计泥质含量的最大值。具体的，即在指示曲线FAST取值为1的全部深度点内， 统计采样井的泥质含量的最大值。
[0094] 在找到采样井的泥质含量的最大值之后，利用泥质含量的最大值，和按照上述公 式（2)或（3)计算出的采样井在不同地层深度的泥质含量，判断采样井在不同地层深度的 地层快慢属性。具体的地层快慢属性的判断规则是：若采样井的当前深度地层的泥质含量 小于所述泥质含量的最大值，则将采样井的当前深度地层判断为快地层，令指示曲线FASTN 在当前深度的取值为1 ;若采样井的当前深度地层的泥质含量大于所述泥质含量的最大 值，则将采样井的邻井的当前深度