一种利用地质构造面曲率预测区域高压盐水层孔隙压力的 方法 【技术领域】
本发明涉及一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法。背景技术 我国大部分油气资源集中在盐下构造， 盐岩地层钻井是制约我国石油钻探的关键 技术难题。由于地层受构造挤压变形的作用， 形成异常高压盐水层， 区域分布无规律， 钻井 过程中经常与高压盐水层发生遭遇战， 当钻井液的密度不能平衡地层压力时， 并会发生盐 水将进入井筒， 污染钻井液， 产生钻井事故和复杂情况， 给钻井作业的人力、 物力带来不同 程度的损失。为此如何在钻井前事先预测高压盐水层孔隙压力是非常必要的， 若在钻井前 能预测到高压盐水层孔隙压力， 就可以为钻井设计确定现场施工时安全钻井液密度提供科 学的依据， 以有效阻止井壁失稳、 防止井下复杂情况的发生。
为此本发明中的创作人凭借其多年从事相关行业的经验与实践， 并经潜心研究与 开发， 终于创造出一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法
发明内容 本发明的目的在于提供一种利用地质构造面曲率预测区域高压盐水层孔隙压力 的方法， 利用该方法可以在钻井前预测到构造区域高压盐水层孔隙压力， 以便在钻井设计 确定现场施工时为确定安全钻井液密度提供科学依据， 以有效阻止井壁失稳、 防止井下复 杂情况的发生。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力的方法， 包括有下列 步骤 ：
1) 根据地震反射剖面和实钻盐层层段， 找出区域盐层的顶界和底界， 获得区域分 布；
2) 利用趋势面法计算构造面底部的曲率， 获得每一空间点的曲率， 计算主曲率 (ρ1， ρ2) ；
3) 根据测井资料回归统计上覆岩层压力 σv ；
4) 根据测井资料计算盐间砂岩夹层的弹性模量和泊松比 ；
5) 利用已知盐层构造面底部的主曲率、 弹性模量和泊松比， 确定所在地层的主应 力， 建立孔隙压力预测模型 ；
6) 计算孔隙压力， 根据每一空间点的孔隙压力， 作压力等值线， 获得压力分布规 律。
所述步骤 2) 中的方法是由盐层构造等值线图上的地理坐标 (xi， yi) 及相应点的高 程数据 wi 利用调和趋势面法建立 w ＝ w(x， y) 的构造图的高程方程， 由方程 w ＝ w(x， y) 计 算构造图上任意一点 (xi， yi) 的主曲率 (ρ1i， ρ2i)。
所述步骤 5) 中的方法是由计算的盐岩构造底部任一点的主曲率 (ρ1i， ρ2i) 结合
相应地层的弹性模量、 泊松比计算该点的主应力 (σ1i， σ2i)， 即获得预测地层任意一点的 三向主应力序列 (σ1i， σ2i， σvi)， 利用岩石构造变形过程中任意点岩石体积变形 Δvi 近似 为岩石孔隙变形 Δvki， 得出盐岩构造曲面上任意一点孔隙压力模型 pi ＝ f(ρ1i， ρ2i， h)。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力是通过地震资料和 实钻录井资料确定盐膏层的区域分布规律和构造等值线图， 利用调和趋势面法建立地质构 造等值线图上的区域高程方程并确定构造区域任意点的主曲率， 计算任意点的主应力， 建 立任意点的孔隙压力预测模型， 确定高压盐水层的孔隙压力， 以便在钻井设计确定现场施 工时为确定安全钻井液密度提供科学依据， 以有效阻止井壁坍塌、 防止井下复杂情况的发 生。 附图说明
图 1 是预测区域地质构造图 ；
图 2 是预测区域高压盐水层孔隙压力分布图。 具体实施方式 下面将结合附图对本发明中的具体实施例作进一步详细说明。
在漫长的地质年代里， 由于地壳的褶皱运动使盐层发生弯曲， 或由于基底隆起使 盐层上拱， 盐间封闭的盐水层受构造挤压、 抬升和盐岩塑形流动的作用形成异常高压， 构造 变形的程度与盐水层异常压力具有一定的联系。 一般来说， 构造变形越剧烈， 异常高压盐水 层孔隙压力越大， 构造变形的程度可以用构造面上任意点曲率来描述， 而曲率与岩石的力 学本构变形有着内在联系。岩石本构关系是揭示岩石受力与变形的定量化描述， 其中包含 着如地层弹性参数、 地层成分、 密度、 埋藏深度、 地质年代、 孔隙率、 构造运动等因素。 地质构 造曲率能不同程度地反映地层所受外界作用的程度， 因此利用地质构造面曲率可以进行钻 前预测区域高压盐水层孔隙压力。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力包括下列步骤 ：
1. 区域盐层分布规律的确定
在预测区域根据地震反射剖面和实钻录井资料确定区域盐层分布规律和构造等 值线分布图。
即首先确定区域地震反射剖面上盐顶和盐底的深度， 同时根据实钻录井资料获得 的盐顶和盐底的深度进行标定， 建立区域的盐层分布规律。
2. 确定盐层地质构造面的主曲率
根据划定的盐层区域的分布规律和获得的盐层构造等值线图 1， 获取地层构造等 值线图上的任意点的地理坐标 (xi， yi) 及相应点的高程数据 wi， 利用调和趋势面法建立构 造图的高程方程 w ＝ w(x， y)， 由方程 w ＝ w(x， y) 计算构造图上任意点 (xi， yi) 的主曲率 (ρ1i， ρ2i)， 具体方法如下 ：
为了计算上的方便， 设构造图高程方程满足一阶傅立叶级数趋势面方程， 具体形 式如下 ：
w ＝ a00+a10A1C0+a01A0C1+a11A1C1+b10B1C0 (2.1)
+b11B1C1+c01A0D1+c11A1D1+d11B1D1
一阶傅立叶级数趋势面方程中有 9 个特定系数 ； a00， a10， a01， a11， b10， b11， c01， c11， d11。式中 ：
为求出方程中的待定系数， 可按最小二乘法原理， 是每个待定系数对观测值与趋 势值的离差平方和
的偏导数等于 0， 即
上述 9 个方程式经整理可得到一阶傅立叶级数趋势面的正规方程组， 可写成如下 矩阵形式
由式 (2.1) 可解出式 (2.1) 中的待定系数。 根据构造图上的各点坐标 (xi， yi) 及高程数据 wi， 可建立正归方程
式中 N 为总点数， w 为海拔高程， 解出 c00， c10， c01， c20， c11…等系数， 即得 w(x， y) 的 近似表达式。
根据薄板小挠度弯曲理论， 由于 w 是微小的， 薄板中面在 x 和 y 方向的曲率及扭率 可近似的表示为 ：
利用 (2.1) 式求出构造图上任意点的三个二阶偏导数和即可计算出构造曲面上任意点的曲率 向。主方向可由式 ：
据此计算构造面上任意点的主曲率及主方确定。其中， αi 为主方向和 x 轴的夹角。构造曲面上任意点主曲率的计算公式为
即3. 根据测井资料回归统计上覆岩层压力 σv
密度测井和声波测井可以直观地反映地层压实规律， 并能获得岩石体积密度值。 由于密度测井资料容易受到井径的大小及仪器检测程度影响， 在利用密度测井数据之前， 应结合井径测井资料过滤掉非真实数据， 以获得比较可靠的密度， 并利用这些密度散点数 据， 利用式 (3.1) 计算上覆岩层压力梯度 ：
式中为任意点 i 的上覆岩层压力梯度， g/cm3 ； ρw， Hw 分别为地层水的密度及水深， g/cm3， m； ρ0， H0 分别为上部无密度测井数据段的平均密度及井深， g/cm3， m； ρi， dhi 为 3 密度测井数据及与其对应的测井层段厚度， g/cm ， m。
由测井数据计算得到的密度散点求出上覆岩层压力梯度数据以后， 通过数据回 归， 即可获得上覆岩层压力的统计规律。 考虑到幂律和二项式回归模型精度较低的缺陷， 采 用式 (3.2) 形式的三次多项式进行回归 ：
Gz ＝ a0+a1h+a2h2+a3h3 (3.2)
式中 Gz 为上覆岩层压力梯度， g/cm3 ； h 为考察点的深度， m； a 0， a1， a2， a3 为待定回归系数。即构造曲面上任意点的上覆岩层压力为 ：
σvi ＝ GZhi (3.3)
4. 根据测井资料计算预测地层的弹性模量和泊松比
地层的弹性模量和泊松比由测井数据解释获得， 具体步骤如下 ：
(1) 计算动态弹性模量和动态泊松比
式中 ：(2) 确定静态弹性模量 和静态泊松比式中 ： a1、 b1、 a2 和 b2 为系数， 视具体区域取值。
5. 建立高压盐水层孔隙压力预测模型
由计算的盐岩构造底部任一点的主曲率 (ρ1i， ρ2i) 结合相应盐岩地层的弹性模 量、 泊松比计算该点的主应力 (σ1i， σ2i)， 即获得盐层任意一点的三向主应力序列 (σ1i， σ2i， σvi)， 利用岩石构造变形过程中任意点岩石体积变形 Δvi 近似为岩石孔隙变形 Δvki， 得出盐岩构造曲面上任意一点孔隙压力模型 pi ＝ f(ρ1i， ρ2i， h)。
具体步骤如下 ：
(1) 盐岩地质构造曲面上任意点三向主应力计算
根据薄板小挠度弯曲理论， 板内一点的应变为 ：
式中 ： w ＝ w(x， y) 为薄板的挠度。 由于式 (2.3) 和 (5.1) 可得则任意点的应力分量为
式中， E、 μ 分别是弹性模量和泊松比。 由 (5.3) 中第三式， 如果 x 和 y 是一点的两个主方向， 则 τxy ＝ 0， 从而 所以 x 和 y 也是该点的两个主曲率方向， 说明薄板弹性曲面上一点的主应力方向和主曲率方 向相一致。
对于弯曲板的任一横截面， 正应力最大值出现在隆起一侧的板面上因而板面上一点的主应力为
(2) 高压盐水层孔隙压力预测模型的建立
岩石受应力作用后， 一部分由岩石孔隙中的流体承受 ( 即孔隙压力 )， 一部分由岩 石的骨架承受 ( 即有效应力 )。若岩石在封闭不排水情况下任意点受到平均应力为 σi 的 三向应力作用， 其平均有效应力 σ0i 与孔隙压力 Ppi 的关系为 ：
根据弹性理论， 在应力作用下任意点岩石的体积变化为 ：孔隙中流体发生的体积变化为 ：式中， ni 为岩石孔隙率， Kvi 为孔隙的体积压缩系数。
由于成岩以后的岩石骨架体积压缩很小， 岩石的体积变化近似等于孔隙的体积变 化， 则 ΔVi ＝ ΔVVi， 于是 ：
式中， Bi 为孔隙压力系数。 联立式 (3.3)、 (5.6)、 (5.7) 和 (5.8) 得高压盐水层孔隙压力预测模型 ：式中定义为等效曲率。6. 计算高压盐水层孔隙压力将待预测区域的地质构造面上的任意点的曲率带入步骤 5) 中的预测模型， 如某并 4464.72m 构造曲率 0.0208， 水平最大、 最小和上覆岩层地应力分别为 116.06、 83.48 和 112.94MPa， 计算出该处孔隙压力 1.86MPa/100m， 实测该点压力数值压力数值 1.92MPa/100m， 满足工程需要。根据每一空间点的孔隙压力， 作压力等值线， 获得压力分布 规律。