区域地层可钻性场数值模拟计算方法研究.pdfVIP

维普资讯 石 油 地 质 与 工 程 2008年 5月 PETROLEUM GEOLOGYANDENGINEERING 第 22卷 第 3期 文章编号 ：1673—8217(2008)03—0082—03 区域地层可钻性场数值模拟计算方法研究 杨明合 ，翟应虎 ，夏宏南 (1．长江大学石油工程学 院，湖北荆州 434023；2．中国石油大学(北京)石油天然气工程学 院) 摘要 ：岩石可钻性是评价岩石破碎难易程度的综合指标 ，是钻头选型的重要依据 ，而 目前所能给 出的可钻性剖面一 般仅为一条 曲线或某个曲面 ，应用范围小 ，难 以满足实际工程的需要 。在分析滑动最小二乘法 的理论基础上，针对 三维空间可钻性场分布规律计算这一问题进行深入地分析和讨论 ，给 出了计算区域可钻性场的理论模型和计算方 法。实例表 明，应用该数学模型计算精度完全满足工程要求 ，该方法切实可行 。 关键词 ：可钻性场 ；滑动最小二乘法 ；权 函数 ；基函数；空 间域 中图分类号：TE21 文献标识码 ：A 岩石可钻性是评价岩石被钻难易程度的综合指 则使用滑动最小二乘法可 以确定 K (x)的一 标 ，是钻头选型 的重要依据之一 。而 目前所能给出 个近似函数为 ： 的可钻性剖面仅仅是根据一 口井的测井和录井数据 r门 Kd(x)一 Pj(x)n (x)一 P (x)n(x) 得出的，依此来模拟区域可钻性剖面实际上就是一 J = l 条按井眼分布的一维 曲线，其直接的结果应用范围 (2) 小，难以满足实际工程的需要。而根据多 口井相关 式中：(x)为 m维基 向量 。为保证收敛性，P(x) 数据得到的区域可钻性横向变化规律图也是仅仅局 应取 自完全多项式 ，三维情况可以选取如下形式： 限于某个给定的曲面，计算精度较低，很难得到三维 线性基 ： (x)一 [1，z，Y，]，m=4； 空间的可钻性场分布规律 。本文针对三维空间可钻 二次基 ： (x)一 [1，z，Y，，xy， ，yz，z。，Y。， 性场分布规律计算这一问题进行深入地分析和讨 ] ， 一 lO； 论，给出了计算区域可钻性场的具体的计算理论和 三次基 ： (x)一 [1，z，Y，，z。，xy，Y。，yz，。，Y-SO， 方法 ，并到满意的结果。 ，z Y，xy。，Y。，Y ，yz。， ， 。z ，Y-SO，xy ] ， 一 20 l 滑动最 ／1,,--乘法基本理论[1．2] 显然 ，由于地层岩石可钻性分布的非线性，用线 一 个三维物理变量场的计算一般方法有 ：①根 性基函数拟合精度较低；二次、三次基 函数精度较 据已知的场函数进行计算；②如果场函数未知 ，则根 高 ，但是三次基函数计算量几乎是二次的二倍 ，这里 据该场中已知的测量数据进行拟合计算。显然地层 选用二次基函数进行计算。 区域可钻性场的计算属于后者。最小二乘法就是最 n(x)为 m维系数 向量 ，它是待求的。与传统 为常用拟合方法之一 。而经典的最小二乘法只在某 的最小二乘拟合方法不 同，这里的待