一种用于定向井井眼轨迹的误差分析方法

一种用于定向井井眼轨迹的误差分析方法
【技术领域】
[0001] 本发明设及井眼检测、石油测斜领域，具体设及一种用于定向井井眼轨迹误差分 析的方法。
【背景技术】
[0002] 井眼测点由S个基本测量参数组成，包括井深D、井斜角I和方位角A。在井眼轨 迹测量过程中存在诸多影响其精度的误差因素，该些误差源W不同方式影响测量值，导致 井眼测量结果不精确。建立误差模型的目的就是综合考虑各种导致测量误差的物理因素， 最终确定在某一测量点处的=维位置误差楠球。
[000引1999年，化曲Williamson在SPE年度技术大会上提出随钻测量井眼不确定性模 型。随后，该模型不断更新并在随钻测量领域得到广泛应用，但并不完全适用定向井分析过 程，误差源及误差参数都不适用于定向井的分析中。
[0004] 国内现有定向井轨迹误差都是采用W抓误差分析模型来分析，该模型由Wolff和 deWart于1981年提出。但该模型只考虑系统误差而忽略随机误差的假设不适合于S维井 眼轨迹。且现代测量仪器多利用惯性测量原理，不能再简单按性能好坏而确定误差权重函 数。

【发明内容】

[0005] 本发明的目的是提供一种用于定向井井眼轨迹的误差分析方法，它克服了现有技 术在定向井测量领域的不足，可为轨迹解释部口提供高质量的数据W及对促进救险井相交 和丛式井防碰提供可靠依据。
[0006] 本发明提供的用于定向井井眼轨迹的误差分析方法，包括如下步骤1~步骤6 :
[0007] 步骤1、采集测斜仪上行过程中测得的井眼位置参数，根据测量仪的仪器模型确定 误差项。其中，井眼位置参数包括井斜角I、方位角AW及井深D;误差源从四方面考虑：深 度误差、不对中误差、传感器误差和仪器精度误差，从每个误差源确定一个或两个W上的误 差项，每个误差项至少包含有误差量级和权重函数。
[000引步骤2、确定误差项作用在测点处的误差向量。
[0009] 设第i个误差项Ei作用在测段1内的第k个测点处的误差向量为e^k;设L表 示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，则第i个误差项ei作 用在最终要计算的测点处的误差向量为。
[0010] 步骤3、将第k个测点处的各误差项累加，获取位置不确定性矩阵C。。、。
[0011] 步骤4、将地理坐标系下的位置不确定性矩阵C。^转换到井眼位置坐标系下，得到 井眼位置坐标系下的位置不确定性矩阵Chh。
[001引其中，Chia= TTCnevT，转换矩阵T为；
[001 引
[0014] Ik和Ak分别表示在当前测段内第k个测点处测量得到的井斜角和方位角。
[0015] 步骤5、根据井眼位置误差的概率密度函数确定测点处误差楠球的各半轴长度。
[0016] 所述的井眼位置误差的概率密度函数为/(Ar'y)，表示如下：
[0017]
[0018] 其中，芯g为考虑误差后的井底位置向量增量，扣，V，W]T为立>"经过正交变换后的 向量，则有：
[0019]
[0020] 其中，m为放大系数，A1、A2和^ 3为正交变换矩阵的特征值；
[002U 误差楠球的各半轴长度的平方分别为m2A1、m2A2和m2入3。
[0022] 步骤6、重复步骤2~步骤5,循环计算下一个测点处的误差楠球；同时，可根据获 得的测点处的误差楠球进行井眼位置不确定性分析。
[0023] 本发明的优点与积极效果在于；（1)现有随钻测量中记录的误差源项目繁多，不 适用于定向井的使用，本发明针对光纤巧螺测斜仪该类测量仪器，确定了适用的模块意义 上的独立误差源，对该类仪器测量误差源选取具有指导意义；（2)可应用于W光纤巧螺为 惯性器件的测斜仪测量定向井井眼轨迹过程中，为后续井眼轨迹解释部口提供可靠依据； (3)本发明适用于华北地区油井，克服了井眼轨迹误差分析只能在随钻测量领域应用的不 足。
【附图说明】
[0024] 图1是本发明的定向井井眼轨迹误差分析方法的整体流程图。
【具体实施方式】
[0025] 下面将结合附图和实施例对本发明作进一步的详细说明。
[0026] 本发明对井眼轨迹精度指导委员会（ISCWSA)公布的随钻测量模型做出优化，首 次应用于定向井测量领域。首先根据仪器类型确定误差类型及参数，然后剔除数据中粗差， 最后对测斜数据经过一系列计算得出测点处的位置误差楠球，作为判断该点处测量误差对 井眼轨迹测量数据影响程度的准则。
[0027] 结合图1对本发明提供的用于定向井井眼轨迹的误差分析方法进行详细说明，包 括步骤1~步骤6。
[002引步骤1、采集测斜仪上行过程中测得的井眼位置参数，并根据测量仪的仪器模型确 定所有的误差源。井眼位置参数包括井斜角I、方位角AW及电缆处测量读出的井深D，= 个参数组成井眼位置向量。
[0029] 根据所用测量仪的类型，依据对应的仪器模型确定误差源。误差源主要从四个方 面进行考虑，分别为：深度误差、不对中误差、传感器误差和仪器精度误差。从每个误差源确 定一个或多个误差项，要根据具体仪器模型来确定误差项。误差项主要通过误差量级、权重 函数W及传播特性等来描述。对误差项的确定是关键，关系到最终模型实施效果的优劣。
[0030] 本发明根据研制的存储式光纤巧螺测斜仪W及该测斜仪测量环境特性分析，确定 了表1所示的11种有关该类测量仪器的模块意义上的独立误差项，尤其适用于华北地区油 井。
[0031] 表1误差项信息表
[0032]
[0033] R表示误差传播特性为随机性传播，S表示误差传播特性为系统性传播。AD表示 两个测点间测量深度的差值。
[0034] 步骤2、对每个误差项确定测段1内第k个测点处的误差向量。
[0035] 误差对某测点k处的作用为前一测段和后一测段作用之和，因此，第i个误差项 ei作用在测段1内第k个测点的误差向量ei,i,k为：
[0036]
…
[0037] 其中，0U表示测段1内第i个误差项的误差量级，本发明实施例中该值根据表1 确定；Pk表示井眼测量向量，即仪器进行下一测点测量相对k点处的向量方向；rk表示井眼 位置向量，即井眼到测点的向量；Ark表示测段长度，即；第k-1个测点到第k个测点间的 长度；Arw表示第k个测点到第k+1个测点间的长度；^表示井眼位置向量对误差项的 导数；^为第i个误差项的权重函数，用来表示第i个误差项在井眼坐标系上对井眼位置 C也 参数的作用大小。
[003引本发明实施例中，设置误差项有11个，则i取值为[1，11]中的整数。一般井眼轨 迹测量长度为3000m，公式（1)适用于连测模式下的测斜仪，为避免数据处理冗杂，选取测 段长度为5m，即每5m进入下一测点。
[0039] 使用平衡正切曲线计算井眼轨迹模型，可将式（1)推导过程表示如下：
[00创其中，j=k，k+1 和Aj.分别表示测段1内第j个测点处测量得到的井深、井 斜角和方位角；和AW分别表示在测段1内第j-1个测点处测量得到的井深、井斜 角和方位角；Ik和Ak分别表示在测段1内第k个测点处测量得到的井斜角和方位角；Ik__i 和Ak__i分别表示在测段1内第k-1个测点处测量得到的井斜角和方位角。
[00创设L表示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，L和K均为正整数。误差的累加在最终要计算的测点处终止。最终要计算的测点处的误差向量 如下：
[0044]
(4)
[0045] 表明测量误差对该点的影响仅限于前一测段内，相反，在测段的中间点处， 测量误差作用在前后测段的测量过程中。0 表示测段L内第i个误差项的误差量级，已 由步骤1给出。
[0046] 根据步骤1中确定的误差项及井眼位置参数，可计算每个误差项在最终要计算的 测点K处的误差向量。
[0047] 步骤3、将第k个测点处各误差项累加，获取位置不确定性矩阵C。。、。
[0048] 获取位置不确定性矩阵C。。、的具体过程是：
[0049] 在测量点处的位置不确定性矩阵由步骤2中计算的误差向量累加获得，根据步骤 1中确定的误差项，用数学式表示误差源累积过程分为W下两种情形：
[0化0] (1)第i个误差项为随机性参考误差时，通过均方和累加获取影响矩阵；
[0051] 在测段1中（不包括计算终点），某随机传播误差项对测点不确定性影响表示为：
[00巧
（5)< br>[005引上角标r表示随机性误差（Random)，馬表示第i个误差项属于随机误差范围时、 第i个误差项在测段内1的协方差矩阵，Ki表示测段1内的测点总数。
[0化4] 所有测段内总的贡献为：
[0化引
（6)
[0化6] Gx表示到计算终点K处、第i个误差项造成的总协方差矩阵。
[0化7] (2)第i个误差项为系统性参考误差时，通过代数和累加获取影响矩阵；
[005引对于具有系统传播特性的误差项i，在测段1中，对测点位置不确定性的贡献为：
[0059]
(7)
[0060] 其中，上角标S表示系统误差（System)，q,表示假设第i个误差项属于系统误差 范围时、第i个误差项在测段1内的协方差矩阵，Ki表示第1个测段内测点总数。
[0061] 所有的测段内总的贡献为：
[0062]
(8)
[0063] 表示到测段L的计算终点K处、第i个误差项造成的总协方差矩阵。
[0064] (3)误差项累加。通过对某测点处的各个误差项依次进行判断和计算，获取位置不 确定性矩阵C。。、。
[0065]总的位置协方差阵cr为：
[0066]
(9)
[0067] 式（9)表明，到计算终点K处，总的位置不确定性矩阵为所有误差项（系统误差和 随机误差）造成的协方差矩阵之和。其中，ieR表示第i个误差项属于随机性参考误差， iGS表示第i个误差项属于系统性参考误差。
[0068] 所要获取的位置不确定性矩阵C。。、为推导出的总的位置协方差阵CT。
[0069] 步骤4、将地理坐标系下的位置不确定性矩阵C。。、转换到井眼位置坐标系下C hh。
[0070] 步骤3中公式都是基于地理位置坐标系下的，地理位置坐标系中S轴分别指向 北、东和天，用下角标nev来标记地理位置坐标系下的值。协方差矩阵需要转换到更直观的 井眼位置坐标系中才方便分析，井眼位置坐标系中=轴分别指向高边、侧边和沿井眼方向， 用下角标hla来标记井眼参考坐标系下的值。其转换方法为：
[007UChia=TTCnevT(10)
[0072] 其中，转换矩阵T为；
[0073]
(11)
[0074] Ik和A k分别表示在当前测段内第k个测点处测量得到的井斜角和方位角。
[0075] 步骤5、根据公式（12)的概率密度函数计算该测点处误差楠球的各半轴长度。
[0076] 根据正态分布，井眼位置误差的分布概率密度函数/'(方为：
[0077]
(12)
[007引其中，击:g为考虑误差后的井底位置向量增量，其等概率密度面为：
[0079]
(13)
[0080] 其中，m为常数。因为Chh为实对称矩阵，设楠球族主轴为OU、OV和OW，其中 扣，V，W]T为友:g经过正交变换后的向量。A1、入2和A3为正交变换矩阵的特征值。将上式 经过正交变换可改写为：
[0084] 该式为误差楠球方程，其半轴长的平方分别为a2=m2A1，b2=m2 ^2,c2=m2入3。 m为放大系数，可根据给定的概率确定。放大系数m的取值参考表2。
[0085] 表2放大系数的取值
[0086]
[0087] 例如当选取标准偏差为±2时，m的取值为73. 85%。
[008引结果表明，轨迹上各点是W轨迹计算点为中屯、的一定概率下的误差楠球域。即步 骤1中所述误差源对定向井井眼轨迹某计算终点处的总的作用可描述为一个楠球域，该测 点W-定概率落在所求误差楠球面内。
[0089] 步骤6、步骤2至步骤5介绍了到计算终点K处的总的位置不确定性矩阵计算方 法，该里的K可代表井眼轨迹上的任一测点。推而广之，可循环计算第K+1个测点。重复步 骤2至步骤5即可计算第K+1个测点处井眼轨迹在一定概率下的误差楠球。
[0090] 利用步骤1至步骤6即可计算井眼轨迹任一测点处的位置误差楠球，从而可参照 该计算结果进行井眼位置不确定性分析。例如可分析其与邻井误差楠球交碰概率大小，或 可计算其与救险井相交的概率大小，用于指导钻井施工。
【主权项】
1. 一种用于定向井井眼轨迹的误差分析方法，其特征在于，包括如下步骤： 步骤1，采集测斜仪上行过程中测得的井眼位置参数，根据测量仪的仪器模型确定误差 项；所述的井眼位置参数包括井斜角I、方位角A以及井深D ; 误差源从四方面考虑：深度误差、不对中误差、传感器误差和仪器精度误差；从每个误 差源确定一个或两个以上的误差项，每个误差项至少设置有误差量级和权重函数； 步骤2,确定误差项作用在测点处的误差向量； 设第i个误差项ε i作用在测段1内的第k个测点处的误差向量为e u，k; 设L表示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，则第i个误 差项Si作用在最终要计算的测点处的误差向量为 步骤3,将第k个测点处的各误差项累加，获取位置不确定性矩阵Cmv; 步骤4,将地理坐标系下的位置不确定性矩阵Cnev转换到井眼位置坐标系下，得到井眼 位置坐标系下的位置不确定性矩阵Chla; Chla= T 1CnevT 其中，转换矩阵T为：Ik和A k分别表示在当前测段内第k个测点处测量得到的井斜角和方位角； 步骤5,根据井眼位置误差的概率密度函数确定测点处误差椭球的各半轴长度； 所述的井眼位置误差的概率密度函数为/(Δ,)，表示如下：其中，Δ，为考虑误差后的井底位置向量增量，[u，v，w]TS 经过正交变换后的向 量，则有：其中，m为放大系数，λ i、λ 2和λ 3为正交变换矩阵的特征值； 误差椭球的各半轴长度的平方分别为m2 λ pm2 λ 2和m 2 λ 3; 步骤6,重复步骤2~步骤5,计算下一个测点处的误差椭球；同时，可根据获得的测点 处的误差椭球进行井眼位置不确定性分析。2. 根据权利要求1所述的一种用于定向井井眼轨迹的误差分析方法，其特征在于，所 述的步骤1中，确定了如下11个误差项：其中，R表示随机性传播，S表示系统性传播，AD表示两个测点间测量深度的差值。3. 根据权利要求1或2所述的一种用于定向井井眼轨迹的误差分析方法，其特征在于， 所述的步骤2的具体实现过程是： 第i个误差项ε i作用在测段1内的第k个测点的误差向量e ia，k为：其中，表示测段1内第i个误差项的误差量级，Pk表示井眼测量向量，rk表示井眼 位置向量，Ark表示第k-Ι个测点到第k个测点间的长度，Ark+1表示第k个测点到第k+1 个测点间的长度； 设L表示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，则第i个误 差项ε i作用在最终要计算的测点处的误差向量<^￡为：其中，σ u表示测段L内第i个误差项的误差量级。4. 根据权利要求1所述的一种用于定向井井眼轨迹的误差分析方法，其特征在于，所 述的步骤3具体实现方法是： (1)当第i个误差项为随机性参考误差时，第i个误差项在测段内1的协方差矩阵CG 为：其中，1表示测段1内的测点总数； 则第i个误差项到计算终点K处所造成的总协方差矩阵?,Α·为：(2) 当第i个误差项为系统性参考误差时，第i个误差项在测段1的协方差矩阵为：则第i个误差项到计算终点κ处所造成的总协方差矩阵为：(3) 获取位置不确定性矩阵Cmv:其中，i e R表示第i个误差项属于随机性参考误差，i e S表示第i个误差项属于系 统性参考误差。
【专利摘要】本发明提供了一种用于定向井井眼轨迹的误差分析方法，用于井眼检测、石油测斜领域。本方法首先根据测量仪的仪器模型确定误差项，误差项都设置有误差量级和权重函数；根据测点处测得的井眼位置参数，确定误差项作用在测点处的误差向量；将第k个测点处的各误差项累加，获取位置不确定性矩阵Cnev；将地理坐标系下的矩阵Cnev转换到井眼位置坐标系下矩阵Chla；根据井眼位置误差的概率密度函数确定测点处误差椭球的各半轴长度，根据结构进行井眼位置不确定性分析。本发明选取适用光纤陀螺测斜仪模块意义上的独立误差源，适用于华北地区油井，可为轨迹解释部门提供高质量的数据以及对促进救险井相交和丛式井防碰提供可靠依据。
【IPC分类】G06F17/50
【公开号】CN104881547
【申请号】CN201510303420
【发明人】高爽, 焦焕静, 孟照魁, 李慧鹏, 林铁, 李先慕, 王璐, 蔡晓雯, 白琳
【申请人】北京航空航天大学
【公开日】2015年9月2日
【申请日】2015年6月5日