随钻电阻率测量的方案分析与实现

随钻电阻率测量的方法的研究与试验

一、课题的背景

本课题来源于胜利石油管理局，胜利石油管理局与我校钻井测控研究中心已合作多年，涉及石油生产的测井、钻井等多个领域，本课题就是在双方进一步合作的基础上，为了满足胜利石油管理局定向井开发的需要而建立的研究课题。

随钻测量（MWD—Measurement While Drilling），是一项在钻井过程中，实时对井底的各种参数进行测量的技术，MWD的最大优点在于它使得司钻和地质工作者实时的看到井下正在发生的情况，可以极大的改善决策过程。最早的随钻测量研究工作始于本世纪30年代，随着钻井技术的发展，1930年出现了最早的井场人工检测法。我国1970曾开始研制MWD系统，但由于种种原因而中断，1981年继续开展这项研究。目前有线随钻测量系统已经通过技术鉴定，井下存储MWD系统正在现场实验，该系统可以测量的参数只有方向、自然伽马和温度，已经完成电磁波传输信道可行性研究。

随钻测量技术极大的推动了钻井技术的发展，为地层评价提供了新的手段，由于可以直接观测井下工程参数，这就为钻井的进一步科学化提供了有利的条件，及时获得地层资料对于准确评价地层和进行地层对比以及油藏描述也具有重要的意义。目前随钻测量技术的研究和应用正向纵深发展。

MWD系统测量的一个十分重要的方面就是电阻率地层评价测井和地质追踪（所谓地质追踪就是用随钻地层评价数据对水平井或大角度斜井进行实时的、交互式的顺层追踪，把非垂直井眼引导到最优化的地质目的层）。1MHz和2MHz 传播工具是目前尖端的MWD电阻率测井仪器，目前Sperry-Sun Drilling Service 服务公司的多空间1~2MHz“电磁波电阻率相位测井”是工业上唯一商业化的、真正的多探测深度的电阻率测井工具。

石油需求量的不断增加和海洋钻井的发展导致了定向井技术的广泛应用，降低钻井成本的持续需求促使提高效益的新工具和新技术的产生，随钻测量技术因此备受关注，在短短的20年里，飞速发展，取得了巨大的进步。目前我国国内对于随钻电阻率测量和导向钻井方面与国外的先进技术相比，还存在着较大的差距。而我国的大型石油生产基地（如大庆油田和胜利油田）以先后将科技引入生产，定向井和大位移井的数量与日俱增，对随钻电阻率测量和导向钻井技术的发展有着迫切的需求，MWD系统具有十分广阔的应用前景。

本课题将在随钻电阻率测量的方法、理论及试验方案上作一些探讨。在随钻测量研究之前常规电阻率测量已有较大的发展，其中有许多电阻率测量方法。最早使用的电阻率测井方法叫作普通电阻率测井，经改进后，发展成为目前广泛使用的聚焦式电阻率测井（或称侧向测井），它包括三侧向、七侧向、双侧向、邻近侧向、球形聚焦等测量方法，一般用于探测深度较深的场合。对于测量深度浅的情况，通常采用微电阻率测井，其电极尺寸小，电极间的距离较近，探测深度浅，主要用于测量井眼泥浆或冲洗带的电阻率。这些都属于直流电测井的范畴。直流电测井要求井内必须充满导电的泥浆或水等井筒液体，这样才能使得测量电流进入地层，达到测量地层的目的。但是对于有些情况，为了准确的了解地层的原始含油饱和度或保持地层的原始渗透性，往往采用油基泥浆或进行空气钻井，

这就使直流电测量遇到了无法克服的困难。感应测井使用交变电磁场对地层电阻率进行测量，不会因为井内不存在导电介质而受到限制。常规感应测井使用的场频较低（多数使用1kHz ），通过测量接收机的感生电动势，对地层电阻率进行分析，它已逐步发展成为常规测井系列中一个主要的电法测井项目。常规测井方法到目前已经发展的较为成熟，取得了许多成功的应用经验。

由于随钻测量是在钻井过程中对地层参数进行测量，其难度大，技术要求高，未定因素多，缺乏足够的现场经验，由此其研究工作远不如常规测井，随钻电阻率测量的研究属于MWD 系统的一个前沿课题。随钻电阻率测量的理论依据是电磁波的传播效应，对于固定的电磁场，不同的地层由于其地质参数的差异对与电磁场的响应也各不相同，不同的钻井环境（如油基钻井液，空气钻井）同样也会影响电磁场的分布，通过检测电磁场的变化，就可以获得地层和钻井环境的有关信息，随钻的含义就是在钻井的过程中，对井下的钻井参数进行实时的测量跟踪。随钻测量的优势在于：可以获得井眼周围地层最原始的地质资料；帮助工程人员及时调整钻井方案，以求得到最佳结果。

二、课题研究的内容

课题的主要研究内容包括：

1、 研究电阻率测量的数学方法，建立电阻率测量的数学模型。 电阻率的随钻测量在MWD 中是一个较为复杂的部分，建立合理的数学模型有助于我们进行深入全面的研究，

基于 Rytov 近似的电阻率测量模型

矢量电磁场的Rytov 近似对地层电阻率的测量有着十分重要的应用价值，我们将据此导出基于电磁波传播效应测量电阻率（EWR ）的灵敏度函数。

假设辐射背景为均匀的各向同性的理想介质，在此基础上进行电磁场的分析和计算，实际上，对于其他的介质可以看作是在理想介质的电磁场中叠加一个扰动量，为了便于分析我们假设一个矩形坐标系，但其结果对其他的坐标系也是同样适用的，因为，在推导过程中，并没有矩形坐标系特有的性质。

将电磁场A 可以表示为：

z y

x i z z i y y i x x e A a e A a e A a A ψψψ000++=

（A-1）

其中0A 满足：

002

02=+∇A k A b

（A-2）

z z y y x x A a A a A a A 0000++=

0A 代表了背景介质，b k 为其波数。 同样A 满足：

022=+∇A k A

（A-3）

其中)(r k k =为场A 的波数，是一个位置的函数，且有：

ωμσμεωi k +=2

ω—辐射频率

μ—为磁导率

ε—为地层的介电常数 σ—为介质电导率

由于在矩形坐标中，矢量拉普拉斯算子与标量拉普拉斯算子等价，因此我们可以将（A-2）式和（A-3）式写成标量的形式：

002

02=+∇x b x A k A

（A-4） 0)(0202=+∇x x i x i x e A k e A ϕϕ

（A-5）

类似的0A 和A 的其他分量也可以写成如上的形式。进一步化简（A-5）式；因为有

x

i x i i x i x i x A e A e e A e A e A x x x x x 002200022)

()(∇⋅∇+∇+∇=∇⋅∇=∇ϕϕϕϕϕ （A-6）

x i i x x ie e ϕϕϕ∇=∇

（A-7） )(22x x x i i i e e x x ϕϕϕϕϕ∇⋅∇-∇=∇

（A-8）

将（A-6）、（A-7）和（A-8）式代入（A-5）式中可得：

02)(0200220=+∇⋅∇+∇+∇⋅∇-∇x x i x x x x x A k A e i A i A x ϕϕϕϕ

（A-9）

有由（A-4）式可得x b x A k A 0202-=∇，因此（A-9）式可以进一步简化得到： 0)(2)(022020=-+∇⋅∇+∇⋅∇-∇x b x i x x x x A k k A e i i A x ϕϕϕϕ

(A-10)

又

x x x x x x x x A A A A ϕϕϕϕ20002022)(∇+∇⋅∇+∇=∇

（A-11）

将（A-10）式两边同乘以i ，利用（A-11）式可以消除x x A 0∇⋅∇ϕ项，可以得到：

x x x b x x b x x A k k i A k A 0220202)()(ϕϕϕϕ∇⋅∇--=+∇

（A-12）

在（A-12）式中由于x x ϕϕ∇⋅∇的存在，使得表达式是非线性的，当满足条件：

22b x x k k -〈〈∇⋅∇ϕϕ

（A-13）