[发明专利]一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统

摘要

本发明涉及油气田勘探开发领域，尤其是一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统。本发明基于岩心统计、薄片观测以及测井解释，通过现今地应力分析，确定裂缝地下开度，利用古地磁定向、岩层产状法或地层倾角法确定裂缝的三维展布，进而建立真三维裂缝渗透率计算模型，推导可编程算法，编制相应的程序，准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。本发明对于真三维空间裂缝渗透率张量及井网设计、确定注水井与采油井的空间位置关系、寻找最优钻井轨迹等油田开发方案的部署具有较高的实用价值，并且预测成本低廉、可操作性强，预测结果可靠。

主权项

1.一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统，步骤如下：1)利用岩心观测、测井资料解释对裂缝的开度测量，确定裂缝的组系特征；依据成像资料、测井或者岩心观察，求取裂缝的线密度；2)在岩心裂缝观测中，对非定向岩心采用古地磁定向、岩层产状法或地层倾角法确定每条裂缝的产状；在成像资料、声波测井解释裂缝时，利用地层倾角测井或者成像测井资料对裂缝定向；3)利用水力压裂法、波速法、声发射法、应力解除法、光弹性应力法对工区现今地应力方位进行了判断，并通过测井、压裂资料以及物理实验计算了关键井现今地应力的数值，利用地震反演、测井解释、岩石力学实验确定岩石的力学属性；在此基础上，利用软件建立有限元模型，采用多井多层段实测地应力的大小、方向为约束条件，对研究区现今地应力场进行了数值模拟；利用模拟结果，确定井筒裂缝发育处的现今地应力大小、方向；利用公式(1)计算裂缝在地下的开度：公式(1)中，b为裂缝的现今地下实际开度，单位：m；b₀为裂缝的原始开度，单位：m；σ′_n为有效正应力，单位：MPa；b_res代表裂缝面承受最大有效正应力时的裂缝开度，单位：m；σ_nref为使裂缝开度减小90％的有效正应力，单位：MPa；4)如图2所示，建立真三维裂缝渗透率张量计算几何模型，设裂缝面Φ的倾角为δ，倾向为ω，以裂缝为参照物建立O‑ABC坐标系，O‑ABC坐标系中的三个坐标轴分别对应于裂缝面的法线方向、裂缝走向方向、裂缝面内垂直于裂缝走向线的方向；设任意渗流面ψ的倾角为α，倾向为β，其单位法向矢量为m；以渗流面ψ为参照物建立坐标系O‑XYZ，其中，OX轴为渗流面ψ的法线方向，OY轴、OZ轴位于渗流面ψ内，OP为渗流面ψ的走向线，HQ为渗流面ψ的倾斜线，HR为渗流面ψ倾斜线的水平投影；线OS、ON、OE、OW、HR、OH、OP、OB在水平面内，定义θ为渗流面ψ内OY轴与OP轴的夹角，通过调整θ的大小，求取裂缝在θ中不同方向的渗透率；5)建立真三维裂缝渗透率张量数学模型：在已知裂缝线密度、地下开度的基础上，单组裂缝的平行渗透率K表示为：公式(2)中，b为裂缝的现今地下实际开度，单位：m；D_lf为单组裂缝的线密度，条/m；设第i条裂缝的渗透率张量为K，则其在坐标系O‑XYZ中的表达式表示为：同理，渗透率张量K在坐标系O‑ABC中的表达式表示为：OA轴在大地坐标系中的三个分量为：OX轴在大地坐标系中的三个分量为：设裂缝的渗透率张量为K，其在O‑ABC坐标系中的表达式表示为：同理，渗透率张量K则其在O‑XYZ坐标系中的表达式表示为：公式(7)中δ为裂缝面Φ的倾角，公式(8)中e、f、g的含义如公式(6)所示；在O‑XYZ坐标系中，每个单元体内单组裂缝的渗透率张量K_ABC可表示为：由公式(7)‑(9)得到：K_XYZi＝T·K_ABC·T^T                           (10)其中，O‑ABC坐标系转换为O‑XYZ坐标系的旋转矩阵T表示为：公式(11)的旋转矩阵T中各参数分别表示为：在O‑XYZ坐标系中，单元体内每组裂缝的渗透率张量K_XYZi表示为：在复杂的地质条件下，储层裂缝经历多期构造应力场演化，每组裂缝的产状、线密度、开度往往不同，依据公式(10)、公式(11)及公式(13)，单元体内发育多组裂缝时，渗透率张量K_XYZ表示为：公式(14)中，k为单元体内裂缝的组数；b_i为第i组裂缝的开度，单位：m；T_i为第i组裂缝在O‑ABC坐标轴分量转换为O‑XYZ坐标轴分量的旋转矩阵，单位：m；D_lfi为第i组裂缝的线密度，单位：条/m；6)利用公式(1)‑(14)推导的算法编制相应的程序，通过循环迭代变量α、β、θ，所述的α为渗流面ψ的倾角，β为渗流面ψ的倾向，θ为OY轴与OP轴的夹角，其中的α范围为[0‑90°]，β的范围为[0‑360°)，θ的范围为[0‑90°)，根据计算精度需求，设置循环迭代步长⊿sp，将α、β、θ逐次累加⊿sp，利用公式(12)‑(14)，计算坐标系O‑XYZ不同的α、β、θ对应的K_XXαβθ、K_YYαβθ、K_ZZαβθ；利用公式(1)‑(14)，设置α＝0°，β为[0‑360°)的任意数值；调节θ的大小，步长⊿sp，范围为[0‑90°)，输出K_YYαβθ、K_ZZαβθ，得到水平面的渗透率变化；利用公式(1)‑(14)，设置α＝90°，β＝0°；调节θ的大小，步长⊿sp，范围为[0‑90°)，输出K_YYαβθ、K_ZZαβθ，得到东西向剖面的渗透率变化；利用公式(1)‑(14)，设置α＝90°，β＝90°；调节θ的大小，步长⊿sp，范围为[0‑90°)，输出K_YYαβθ、K_ZZαβθ，得到南北向剖面的渗透率变化；利用公式(1)‑(14)，设置α范围为[0‑90°]，β的范围为[0‑360°)，步长，⊿sp；计算不同的α、β对应的平面内——即不同的渗流面ψ内，K_YYαβθ、K_ZZαβθ，每个渗流面ψ内计算得到的K_YYαβθ、K_ZZαβθ数目为2n，2n个K_YYαβθ、K_ZZαβθ的平均值为K_aver；计算不同的α、β对应的方差函数W_αβ：公式(15)中，K_YYαβθi、K_ZZαβθi分别为特定的α、β对应的第i个Y方向渗透率主值、Z方向渗透率主值，单位：10^‑3μm²；n为每个渗流面ψ计算的Y或Z方向渗透率主值的数目，单位：个；K_aver为2n个K_YYαβθ、K_ZZαβθ的平均值，单位：10^‑3μm²；所述的K_XXαβθ为X方向渗透率主值，K_YYαβθ为Y方向渗透率主值，K_ZZαβθ为Z方向渗透率主值；通过公式(15)，计算不同的α、β对应的方差函数W_αβ，比较不同的α、β的W_αβ；其中，W_αβ最大值对应的渗流面ψ为最大差渗透率面，记录对应的倾角α，倾向β；W_αβ最小值对应的渗流面ψ为等渗透率面，记录对应的倾角α，倾向β。