本发明公开了一种利用录井数据的井漏裂缝宽度随钻估算方法，该方法首先建立基于漏失量的裂缝宽度随钻估算模型，模型假设裂缝面与井眼水平面的夹角为γ，γ取值范围为0≤γ≤90°；根据井眼轴线与裂缝面相交夹角γ的大小，可以代表不同的情形：当γ＝0°时，井眼轴线与裂缝面垂直相交；当0
1.一种利用录井数据的井漏裂缝宽度随钻估算方法，其特征在于，步骤如下：S1、建立基于漏失量的裂缝宽度随钻估算模型，模型假设裂缝面与井眼水平面的夹角为γ，γ取值范围为0≤γ≤90°；
当0≤γ<90°时，裂缝中的流动为径向流，无因次侵入半径rD(i)的计算公式如下：理论无因次侵入时间tD(i),T计算公式如下：
式(1)-(3)中，Vm(ti)—时间点ti的累计漏失量，L；wf—裂缝宽度，m；rwe—当量井眼半径， rw—井眼半径，m；n—钻井液流性指数；τy—钻井液的动切力，Pa；△p—井筒和漏层的压差，Pa；K—稠度系数，Pa·s；t—实测时间，s；t(i),m—第i个实测时间，s；
当γ＝90°时，裂缝中的流动为线性流，无因次侵入深度xD(i)的计算公式：理论无因次侵入时间tD(i),T计算公式如下：
式(4)-(6)中，Vm(ti)—时间点ti的累计漏失量，L；wf—裂缝宽度，mm；n—钻井液流性指数；τy—钻井液的动切力，Pa；△p—井筒和漏层的压差，MPa；K—稠度系数，Pa·s；h—裂缝高度，m；rw—井眼半径，m；t—实测时间，s；t(i),m—第i个实测时间，s；
S2、基于步骤S1的估算模型公式，采用图版曲线对比算法计算裂缝宽度，根据夹角γ的大小不同，具体的计算步骤不同：当0≤γ<90°时，按如下步骤计算：
S22、根据时间点ti的累计漏失量Vm(ti)，按照公式(1)计算无因次侵入半径rD(i)，同理，计算时间点ti+1的无因次侵入半径rD(i+1)；
S23、利用S22得到的无因次侵入半径，按照公式(2)计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T；其中，当i＝1时，取tD(i-1),T＝0；
S24、按照公式(3)计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1，S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，重复步骤S22～S25，直到满足条件Si+1≥Si为止；
S22、根据时间点ti的累计漏失量Vm(ti)，按照公式(4)计算无因次侵入深度xD(i)，同理，计算时间点ti+1的无因次侵入深度xD(i+1)；
S23、利用S22得到的无因次侵入深度，按照公式(5)计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T，其中，当i＝1时，取tD(i-1),T＝0；
S24、按照公式(6)计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1：S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，重复步骤S22～S25，直到满足条件Si+1≥Si为止。
2.如权利要求1所述的利用录井数据的井漏裂缝宽度随钻估算方法，其特征在于，步骤S21中，取裂缝宽度wf的初始值wf,j＝wf,0＝0.1mm。
3.如权利要求1所述的利用录井数据的井漏裂缝宽度随钻估算方法，其特征在于，步骤S26中，裂缝宽度增加步长Δw＝0.1mm。
4.如权利要求1所述的利用录井数据的井漏裂缝宽度随钻估算方法，其特征在于，步骤S1中，当γ＝0°时，井眼轴线与裂缝面垂直相交；当0<γ<90°时，井眼轴线与裂缝面倾斜相交；当γ＝90°时，井眼轴线与裂缝面平行。
一种利用录井数据的井漏裂缝宽度随钻估算方法
技术领域
[0001] 本发明涉及石油钻井技术领域，特别是一种利用录井数据的井漏裂缝宽度随钻估算方法。
背景技术
[0002] 在油气勘探开发的钻井过程中，常常会发生裂缝性漏失。井漏是钻井工程中常见的井内复杂情况，多数钻井过程都有不同程度的漏失，严重的井漏会导致井内压力下降，影响正常钻井、引起井壁失稳、诱发地层流体涌入井筒并井喷。
[0003] 井下裂缝具有一定的张开宽度是裂缝性漏失发生的必然条件之一。裂缝宽度不仅是影响钻井液漏失严重程度的重要因素，而且是后续进行堵漏钻井液设计的关键性基础参数。确定井下裂缝宽度特别是井壁裂缝的张开宽度对于裂缝性漏失堵漏设计至关重要，只有准确的确定井下裂缝的真实宽度，才能科学合理地设计封堵裂缝的堵漏材料。
发明内容
[0004] 本发明的目的是提供了一种利用录井数据的井漏裂缝宽度随钻估算方法。
[0005] 本发明提供的利用录井数据的井漏裂缝宽度随钻估算方法，步骤如下：
[0006] 步骤S1、建立基于漏失量的裂缝宽度随钻估算模型。为了全面反映裂缝面的走向与井眼的位置关系，模型假设裂缝面与井眼水平面的夹角为γ，γ取值范围为0≤γ≤90°。根据夹角γ的大小，可以代表不同的情形：(a)当γ＝0°时，井眼轴线与裂缝面垂直相交；
(b)当0<γ<90°时，井眼轴线与裂缝面倾斜相交；(c)当γ＝90°时，井眼轴线与裂缝面平行。
[0007] (1)当0≤γ<90°时，裂缝中的流动为径向流：
[0008] 赫巴模式本构方程：
[0009]
[0010] 式中：
[0011] z—沿裂缝宽度wf方向的垂直距离，m；
[0012] n—钻井液流性指数；
[0013] ν—在位置z的速度，m/s；
[0014] K—稠度系数，Pa·s；
[0015] τy—钻井液的动切力，Pa。
[0016] 压降梯度方程：
[0017]
[0018] 式中：
[0019] p—流体压力，Pa；
[0020] r—径向距离，m；
[0021] q—截面流量，L/s；
[0022] wf—裂缝宽度，m；
[0023] 当n＝1时，即流体为宾汉模式，压降梯度方程简化为：
[0024]
[0025] 式中：
[0026] μp—塑性粘度，Pa·s；
[0027] vm—径向流条件下钻井液质点在裂缝中的流速，m/s。
[0028] 当量井眼半径rwe：
[0029]
[0030] 式中：
[0031] γ—裂缝面与井眼水平面的夹角；
[0032] rw—井眼半径，m；
[0033] 时间点t的累计漏失量：
[0034]
[0035] 式中：rf(t)—t时刻钻井液的径向侵入距离，m；
[0036] 漏失速率：
[0037]
[0038] 流速方程：
[0039]
[0040] 整理，得
[0041]
[0042] 式中:
[0043] rf—钻井液的径向侵入距离，m；
[0044] Δp—井筒和漏层的压差，Pa
[0045] 引入无因次变量，无因次侵入半径rD和无因次侵入时间tD：
[0046]
[0047] 将公式 两边同时平方，得到，然后，两边开平方，得到
[0048]
[0049] tD＝βt (3)
[0050] 其中β：
[0051]
[0052] 进行无因次化处理以后，得：
[0053]
[0054] 即
[0055]
[0056] 其中：ΔrD＝rD(i)-rD(i-1)
[0057]
[0058] (2)当γ＝90°时，裂缝中的流动为线性流：
[0059] 时间点t的累计漏失量：
[0060] Vm(t)＝2wfh[xf(t)-rw]
[0061] 式中：
[0062] h—裂缝高度，m；
[0063] xf(t)—t时刻钻井液侵入深度，m；
[0064] wf—裂缝宽度，m；
[0065] rw—井眼半径，m。
[0066] 漏失速率：
[0067]
[0068] 即：
[0069]
[0070] 流速方程：
[0071]
[0072] 钻井液在裂缝中线性流动的压降梯度方程如下：
[0073]
[0074] 式中：
[0075] p—流体压力，Pa；
[0076] x—沿x轴方向距离，m；
[0077] n—钻井液流性指数；
[0078] K—稠度系数，Pa·s；
[0079] q—截面流量，L/s；
[0080] τy—钻井液的动切力，Pa；
[0081] 对上式两边同时积分，得
[0082]
[0083] 式中：
[0084] xf—钻井液侵入深度，m；
[0085] Δp—井筒和漏层的压差，Pa。
[0086] 整理，得单翼裂缝流量：
[0087]
[0088] 同时，单翼裂缝流量：
[0089]
[0090] 联立单翼裂缝流量公式，即得液体在裂缝内的侵入前沿方程：
[0091]
[0092] 整理，得：
[0093]
[0094] 令无因次侵入深度xD和无因次侵入时间tD：
[0095]
[0096] 进一步变换得到，
[0097] tD＝βt (6)
[0098] 式中，令
[0099] 则：
[0100]
[0101] 无因次化处理，得到：
[0102]
[0103] 因此，钻井液无因次侵入前沿：
[0104]
[0105] 即：
[0106]
[0107] 其中：
[0108]
[0109] wf—裂缝宽度，m；n—钻井液流性指数；τy—钻井液的动切力，Pa；△p—井筒和漏层的压差，MPa；h—裂缝高度，m；rw—井眼半径，m。
[0110] 上述的基于漏失量的裂缝宽度随钻估算模型利用井漏数据求得天然裂缝的宽度。而目前没有直接求取裂缝宽度的解析解，只能通过数值解得到图版曲线，然后将无因次化的实测井漏数据曲线与图版曲线对比，从而得到裂缝宽度。图版曲线对比算法步骤如下：
[0111] 步骤S2、基于步骤S1的估算模型公式，采用图版曲线对比算法计算裂缝宽度，根据夹角γ的大小不同，具体的计算步骤不同：
[0112] 当0≤γ<90°时，按如下步骤计算：
[0113] S21、取一个裂缝宽度值wf＝wf,j，例如，取裂缝宽度wf的初始值wf,j＝wf,0＝0.1mm。
[0114] S22、根据时间点ti的累计漏失量Vm(ti)，按照公式(1)计算无因次侵入半径rD(i)，同理，计算时间点ti+1的无因次侵入半径rD(i+1)；
[0115]
[0116]
[0117] S23、利用S22得到的无因次侵入半径，按照公式(2)计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T；其中，当i＝1时，取tD(i-1),T＝0；
[0118]
[0119]
[0120] S24、按照公式(3)计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
[0121]
[0122]
[0123] S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1，
[0124]
[0125]
[0126] S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，例如，取裂缝宽度增加步长Δw＝0.1mm，重复步骤S22～S25，直到满足条件Si+1≥Si为止。
[0127] 当γ＝90°时，按如下步骤计算：
[0128] S21、取一个裂缝宽度值wf＝wf,j；例如，取裂缝宽度wf的初始值wf,j＝wf,0＝0.1mm。
[0129] S22、根据时间点ti的累计漏失量Vm(ti)，按照公式(4)计算无因次侵入深度xD(i)，同理，计算时间点ti+1的无因次侵入深度xD(i+1)；
[0130]
[0131]
[0132] S23、利用S22得到的无因次侵入深度，按照公式(5)计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T，其中，当i＝1时，取tD(i-1),T＝0；
[0133]
[0134]
[0135] S24、按照公式(6)计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
[0136]
[0137]
[0138] S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1：
[0139]
[0140]
[0141] S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，例如，取裂缝宽度增加步长Δw＝0.1mm，重复步骤S22～S25，直到满足条件Si+1≥Si为止。
[0142] 与现有技术相比，本发明的有益之处在于：
[0143] 本发明利用新建立的基于漏失量的裂缝宽度随钻估算数学模型，可根据漏失速度、累积漏失量计算得出裂缝宽度，可以用于后续科学合理地设计封堵裂缝的堵漏材料。
[0144] 本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
[0145] 图1、裂缝面与井眼水平面相交几何模型图。
[0146] 图2、本发明的图版曲线对比算法流程图。
[0147] 图3、应用实例中漏失速率随时间的变化数据。
[0148] 图4、应用实例1中无因次侵入半径与无因次时间的关系图。
[0149] 图5、应用实例2中无因次侵入深度与无因次时间的关系图。
具体实施方式
[0150] 以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
[0151] 本发明的利用录井数据的井漏裂缝宽度随钻估算方法，首先是建立基于漏失量的裂缝宽度随钻估算模型。为了清楚明白地反映裂缝面的走向与井眼的位置关系，给出了图1。如图1所示，模型假设裂缝面与井眼水平面的夹角为γ，γ取值范围为0≤γ≤90°。然后基于建立的数学模型，采用图版曲线对比算法计算裂缝宽度，图2显示了图版曲线对比算法流程图。
[0152] 应用实例1
[0153] 采用本发明的裂缝宽度随钻估算方法进行具体计算：
[0154] (1)计算参数工程数据取值如下表：
[0155]
[0156] 获取的漏失速率随时间的变化数据见图3。
[0157] 由于漏层倾角(即裂缝面与井眼水平面的夹角为γ)为45°，处于0≤γ<90°范围内，采用对应的图版曲线对比算法步骤如下：
[0158] S21、取一个裂缝宽度值wf＝wf,j，取裂缝宽度wf的初始值wf,j＝wf,0＝0.1mm。
[0159] S22、根据时间点ti的累计漏失量Vm(ti)，计算无因次侵入半径rD(i)，同理，计算时间点ti+1的无因次侵入半径rD(i+1)；
[0160]
[0161]
[0162] S23、利用S22得到的无因次侵入半径，计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T；其中，当i＝1时，取tD(i-1),T＝0；
[0163]
[0164]
[0165] S24、计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
[0166]
[0167]
[0168] S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1，
[0169]
[0170]
[0171] S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，取裂缝宽度增加步长Δw＝0.1mm，重复步骤S22～S25，直到满足条件Si+1≥Si为止。
[0172] 利用上述工程数据和上述步骤，根据漏失量数据计算出每个实测时间对应的理论时间，得到图4所示曲线。当裂缝宽度为1.05mm时，满足Si+1≥Si，则取1.05mm为最接近的裂缝宽度值，该裂缝宽度可视为真实裂缝宽度。
[0173] 应用实例2
[0174] 采用本发明的裂缝宽度随钻估算方法进行具体计算：
[0175] (1)计算参数工程数据取值如下表：
[0176]当前井段井径(mm) 200 动切力τy(pa) 10
漏层倾角γ(°) 90 稠度系数K(pa·s) 0.3
漏层顶界深度(m) 1000 计算精度(mm) 0.05
漏层底界深度(m) 1001
井筒和漏层的压差(MPa) 2
流性指数n 0.9
[0177] 假设获取的漏失速率随时间的变化数据依然是图3。
[0178] 由于漏层倾角(即裂缝面与井眼水平面的夹角为γ)为90°，采用对应的图版曲线对比算法步骤如下：
[0179] S21、取一个裂缝宽度值wf＝wf,j，取裂缝宽度wf的初始值wf,j＝wf,0＝0.1mm。
[0180] S22、根据时间点ti的累计漏失量Vm(ti)，计算无因次侵入深度xD(i)，同理，计算时间点ti+1的无因次侵入深度xD(i+1)；
[0181]
[0182]
[0183] S23、利用S22得到的无因次侵入深度，计算时间点ti和ti+1对应的理论无因次时间tD(i),T和tD(i+1),T，其中，当i＝1时，取tD(i-1),T＝0；
[0184]
[0185]
[0186] S24、计算时间点ti和ti+1对应的实测无因次时间tD(i),m和tD(i+1),m；
[0187]
[0188]
[0189] S25、计算时间点ti和ti+1各自的理论无因次时间与实测无因次时间误差的平方和Si、Si+1：
[0190]
[0191]
[0192] S26、判断误差条件：若Si+1≥Si，则取wf,j为最接近的裂缝宽度值；否则，继续增加裂缝宽度值wf＝wf,j+Δw，取裂缝宽度增加步长Δw＝0.1mm，重复步骤S22～S25，直到满足条件Si+1≥Si为止。
[0193] 利用上述工程数据和上述步骤，根据漏失量数据计算出每个实测时间对应的理论时间，得到图5所示曲线。当裂缝宽度为4.65mm时，满足Si+1≥Si，则取4.65mm为最接近的裂缝宽度值，该裂缝宽度可视为真实裂缝宽度。
[0194] 以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
