基于多井约束优化方法的低渗油藏应力场反演与裂缝预测技术及应用.pdfVIP

基于多井约束优化方法的低渗油藏 应力场反演与裂缝预测技术及应用 闫相祯1 杨秀娟1 王建军1 刘 红2 张建良2 (1中国石油大学(华东)，山东东营2中石化江苏油田地质研究院，江苏扬州) 摘要本文针对钠气储集层地应力场反演与裂缝预测问题，运用有限元法和最优化方法建 立了地应力场与裂缝分析的有限元优化反演方法。对苏北盆地高邮凹陷W8断块的地应力场分析 表明．计算结果与文测值吻合较好，说明这种方法具有良好的实用价值，可为地应力场的反演与裂 缝分析提供一定参考。该方法克服了回归方法及边界调整法等一般反分析的不足，具有较好的实 用价值，能够方便地反演区域地应力场。 关键词低渗油藏地应力场裂缝有限元反演 引 言 在低渗油藏的开发中，迫切需要知道地应力场的分布状态。本文从弹性力学基本方程出 发，针对油气储集层地应力场反问题，结合对区域构造应力场认识，建立地应力场反演的基本 理论，采用混合罚函数法，建立应力场反演的约束优化模型，提出了多井约束优化方法和有限 元法相结合的数值反演方法，根据研究区少量地应力实测资料进行最优化计算，使计算应力场 与实测应力达到最优拟合，以弥补目前的回归反演方法和边界载荷调整法的不足，从而提高计 算精度。 地应力场反演的基本理论 在油气勘探开发中，所知道的仅是某一区域几个点位的地应力测量资料，而要知道整个区 域的应力场状态，就需要根据获得的舰测结果进行反演计算求解应力场。这类反问题在数学 上称为微分方程反分析。由于微分方程反问题的解算困难，要获得其解析解通常是不可能的， 因此须采用数值计算方法。 1．1地应力场的基本理论 设含有断层和边界S的空问域V，根据三维边值问题的有限元方程可以得到节点应力与 边界载荷满足如下关系： 口=TF(X) (1) T：DLK。1 (2) edu(“ 和科研』=作。Email：y“xz@hdpu ——142—— 式中一——节点应力； F(X)——节点力列阵； x—一描述储集层所承受未知边界载荷的参数向量； 卜传递矩阵； D——为弹性矩阵； L——微分算子矩眸； x——总体刚度矩阵。 在式(1)中，|_{{jr待研究区块储集层中只有部分油水井的应力值已知，需要求解边界载 荷，是求解微分方程的反问题。利用有限元法将所讨论的边值问题转化为一组由系统节点应 力表示的方程组，给定介质的力学参数以及边界条件即可求解节点的应力场。 1．2地应力场优化反演模型 储集层应力场反演就是在一定的准则下寻找能够最优拟合实际观测数据的模型或参数。 利用获取的观测数据，对优化准则或参数加以约束。储集层应力场获取的观测数据通常为区 域部分井的压裂、井壁崩落或岩心测试等的应力。 为方便优化方法应用于有限元模拟之中，利用序列无约束极小化技术，通过添加惩罚函数 的方法把约束问题转化为无约束问题，推导出与有限元计算相关联的优化数值计算模型。 考虑约束最优化问题： min f=“x) J．．t g。(x)≤莓， i=1，2，3，…，m1 壹。≤h。(z)， i=1，2，3，…，m2 (3) 型。≤埘。(x)≤叵， i=1，2，3，…，m3 z∈XCR8， x=(z1，z2，x3，·一，z．)’ 式中产一目标函数； x——设计变量； g。、h。、W；——状态变量，可为主应力方向和由应力表示的函数。 对约束问题进行无约束演绎，无约束日标函数如下