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

本发明专利技术提供了一种用于定向井井眼轨迹的误差分析方法，用于井眼检测、石油测斜领域。本方法首先根据测量仪的仪器模型确定误差项，误差项都设置有误差量级和权重函数；根据测点处测得的井眼位置参数，确定误差项作用在测点处的误差向量；将第k个测点处的各误差项累加，获取位置不确定性矩阵Cnev；将地理坐标系下的矩阵Cnev转换到井眼位置坐标系下矩阵Chla；根据井眼位置误差的概率密度函数确定测点处误差椭球的各半轴长度，根据结构进行井眼位置不确定性分析。本发明专利技术选取适用光纤陀螺测斜仪模块意义上的独立误差源，适用于华北地区油井，可为轨迹解释部门提供高质量的数据以及对促进救险井相交和丛式井防碰提供可靠依据。

【技术实现步骤摘要】

本专利技术设及井眼检测、石油测斜领域，具体设及一种用于定向井井眼轨迹误差分 析的方法。
技术介绍
井眼测点由S个基本测量参数组成，包括井深D、井斜角I和方位角A。在井眼轨 迹测量过程中存在诸多影响其精度的误差因素，该些误差源W不同方式影响测量值，导致 井眼测量结果不精确。建立误差模型的目的就是综合考虑各种导致测量误差的物理因素， 最终确定在某一测量点处的=维位置误差楠球。[000引1999年，化曲Williamson在SPE年度技术大会上提出随钻测量井眼不确定性模 型。随后，该模型不断更新并在随钻测量领域得到广泛应用，但并不完全适用定向井分析过 程，误差源及误差参数都不适用于定向井的分析中。 国内现有定向井轨迹误差都是采用W抓误差分析模型来分析，该模型由Wolff和 deWart于1981年提出。但该模型只考虑系统误差而忽略随机误差的假设不适合于S维井 眼轨迹。且现代测量仪器多利用惯性测量原理，不能再简单按性能好坏而确定误差权重函 数。
技术实现思路
本专利技术的目的是提供，它克服了现有技 术在定向井测量领域的不足，可为轨迹解释部口提供高质量的数据W及对促进救险井相交 和丛式井防碰提供可靠依据。 本专利技术提供的用于定向井井眼轨迹的误差分析方法，包括如下步骤1~步骤6 : 步骤1、采集测斜仪上行过程中测得的井眼位置参数，根据测量仪的仪器模型确定 误差项。其中，井眼位置参数包括井斜角I、方位角AW及井深D;误差源从四方面考虑：深 度误差、不对中误差、传感器误差和仪器精度误差，从每个误差源确定一个或两个W上的误 差项，每个误差项至少包含有误差量级和权重函数。[000引步骤2、确定误差项作用在测点处的误差向量。 设第i个误差项Ei作用在测段1内的第k个测点处的误差向量为e^k;设L表 示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，则第i个误差项ei作 用在最终要计算的测点处的误差向量为。 步骤3、将第k个测点处的各误差项累加，获取位置不确定性矩阵C。。、。 步骤4、将地理坐标系下的位置不确定性矩阵C。^转换到井眼位置坐标系下，得到 井眼位置坐标系下的位置不确定性矩阵Chh。[001引其中，Chia= TTCnevT，转换矩阵T为；[001 引 Ik和Ak分别表示在当前测段内第k个测点处测量得到的井斜角和方位角。 步骤5、根据井眼位置误差的概率密度函数确定测点处误差楠球的各半轴长度。 所述的井眼位置误差的概率密度函数为/(Ar'y)，表示如下： 其中，芯g为考虑误差后的井底位置向量增量，扣，V，W]T为立>"经过正交变换后的 向量，则有： 其中，m为放大系数，A1、A2和^ 3为正交变换矩阵的特征值；[002U 误差楠球的各半轴长度的平方分别为m2A1、m2A2和m2入3。 步骤6、重复步骤2~步骤5,循环计算下一个测点处的误差楠球；同时，可根据获 得的测点处的误差楠球进行井眼位置不确定性分析。 本专利技术的优点与积极效果在于；（1)现有随钻测量中记录的误差源项目繁多，不 适用于定向井的使用，本专利技术针对光纤巧螺测斜仪该类测量仪器，确定了适用的模块意义 上的独立误差源，对该类仪器测量误差源选取具有指导意义；（2)可应用于W光纤巧螺为 惯性器件的测斜仪测量定向井井眼轨迹过程中，为后续井眼轨迹解释部口提供可靠依据； (3)本专利技术适用于华北地区油井，克服了井眼轨迹误差分析只能在随钻测量领域应用的不 足。【附图说明】 图1是本专利技术的定向井井眼轨迹误差分析方法的整体流程图。【具体实施方式】 下面将结合附图和实施例对本专利技术作进一步的详细说明。 本专利技术对井眼轨迹精度指导委员会（ISCWSA)公布的随钻测量模型做出优化，首 次应用于定向井测量领域。首先根据仪器类型确定误差类型及参数，然后剔除数据中粗差， 最后对测斜数据经过一系列计算得出测点处的位置误差楠球，作为判断该点处测量误差对 井眼轨迹测量数据影响程度的准则。 结合图1对本专利技术提供的用于定向井井眼轨迹的误差分析方法进行详细说明，包 括步骤1~步骤6。[002引步骤1、采集测斜仪上行过程中测得的井眼位置参数，并根据测量仪的仪器模型确 定所有的误差源。井眼位置参数包括井斜角I、方位角AW及电缆处测量读出的井深D，= 个参数组成井眼位置向量。 根据所用测量仪的类型，依据对应的仪器模型确定误差源。误差源主要从四个方 面进行考虑，分别为：深度误差、不对中误差、传感器误差和仪器精度误差。从每个误差源确 定一个或多个误差项，要根据具体仪器模型来确定误差项。误差项主要通过误差量级、权重 函数W及传播特性等来描述。对误差项的确定是关键，关系到最终模型实施效果的优劣。 本专利技术根据研制的存储式光纤巧螺测斜仪W及该测斜仪测量环境特性分析，确定 了表1所示的11种有关该类测量仪器的模块意义上的独立误差项，尤其适用于华北地区油 井。 表1误差项信息表 R表示误差传播特性为随机性传播，S表示误差传播特性为系统性传播。AD表示 两个测点间测量深度的差值。 步骤2、对每个误差项确定测段1内第k个测点处的误差向量。 误差对某测点k处的作用为前一测段和后一测段作用之和，因此，第i个误差项 ei作用在测段1内第k个测点的误差向量ei,i,k为：… 其中，0U表示测段1内第i个误差项的误差量级，本专利技术实施例中该值根据表1 确定；Pk表示井眼测量向量，即仪器进行下一测点测量相对k点处的向量方向；rk表示井眼 位置向量，即井眼到测点的向量；Ark表示测段长度，即；第k-1个测点到第k个测点间的 长度；Arw表示第k个测点到第k+1个测点间的长度；^表示井眼位置向量对误差项的 导数；^为第i个误差项的权重函数，用来表示第i个误差项在井眼坐标系上对井眼位置 C也 参数的作用大小。中的整数。一般井眼轨 迹测量长度为3000m，公式（1)适用于连测模式下的测斜仪，为避免数据处理冗杂，选取测 段长度为5m，即每5m进入下一测点。 使用平衡正切曲线计算井眼轨迹模型，可将式（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均为正整数。误差的累加在最终要计算的测点处终止。最终要计算的测点处的误差向量 如下：(4) 表明测量误差对该点的影响仅限于前一测段内，相反，在测段的中间点处， 测量误差作用在前后测段的测量过程中。0 表示测段L内第i个误差项的误差量级，已 由步骤1给出。 根据步骤1中确定的误差项及井眼位置参数，可计算每个误差项在最终要计算的 测点K处的误差向量。 步骤3、将第k个测点处各误差项累加，获取位置不确定性矩阵C。。、。 获取位置不确定性矩阵C。。、的具体过程是： 在测量点处的位置不确定性矩阵由步骤2中计算的误差向量累加获得，根据步骤 1中确定的误差项，用数学式表示误差源累积过程分为W下两种情形： (1)第i个误差项本文档来自技高网...

【技术保护点】
一种用于定向井井眼轨迹的误差分析方法，其特征在于，包括如下步骤：步骤1，采集测斜仪上行过程中测得的井眼位置参数，根据测量仪的仪器模型确定误差项；所述的井眼位置参数包括井斜角I、方位角A以及井深D；误差源从四方面考虑：深度误差、不对中误差、传感器误差和仪器精度误差；从每个误差源确定一个或两个以上的误差项，每个误差项至少设置有误差量级和权重函数；步骤2，确定误差项作用在测点处的误差向量；设第i个误差项εi作用在测段l内的第k个测点处的误差向量为ei,l,k；设L表示到达测量终点前的测段，K表示测段L最终要计算的测点处位置，则第i个误差项εi作用在最终要计算的测点处的误差向量为步骤3，将第k个测点处的各误差项累加，获取位置不确定性矩阵Cnev；步骤4，将地理坐标系下的位置不确定性矩阵Cnev转换到井眼位置坐标系下，得到井眼位置坐标系下的位置不确定性矩阵Chla；Chla＝TTCnevT其中，转换矩阵T为：T=cosIkcosAk-sinAksinIkcosAkcosIksinAkcosAksinIksinAk-sinIk0cosIk]]>Ik和Ak分别表示在当前测段内第k个测点处测量得到的井斜角和方位角；步骤5，根据井眼位置误差的概率密度函数确定测点处误差椭球的各半轴长度；所述的井眼位置误差的概率密度函数为表示如下：f(Δr→g)=1(2π)3/2|Chla|3/2e-12(Δr→g)Chla-1(Δr→g)]]>其中，为考虑误差后的井底位置向量增量，[U,V,W]T为经过正交变换后的向量，则有：U2λ1+V2λ2+W2λ3=m2]]>其中，m为放大系数，λ1、λ2和λ3为正交变换矩阵的特征值；误差椭球的各半轴长度的平方分别为m2λ1、m2λ2和m2λ3；步骤6，重复步骤2～步骤5，计算下一个测点处的误差椭球；同时，可根据获得的测点处的误差椭球进行井眼位置不确定性分析。...

【技术特征摘要】

【专利技术属性】
技术研发人员：高爽，焦焕静，孟照魁，李慧鹏，林铁，李先慕，王璐，蔡晓雯，白琳，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11