CN102536201A - 多臂成像数据校正方法 - Google Patents

Abstract

本发明提供一种多臂成像数据校正方法。该数据校正方法，在测井仪器中对应于各个探测臂均设置有位移传感器，测井仪器中还设置有用以获取该测井仪器实时方位及倾角信息的辅助测量传感器，其中，所述方位是由水平旋转的绝对或相对角度表达；所述数据校正方法包括两方面的数据校正：(1)对由于仪器在井下自然旋转而造成各个臂发生位置移动的方位校正；(2)对由于井斜等原因造成的多臂测量仪器不居中而使探测臂的测量数据产生偏差的倾斜校正。本发明针对多臂测量仪器与数值校正算法相结合，根据仪器测量各探测臂数值及辅助校正数值进行各探测臂数值的校正处理，得到更精确的探测臂数值，从而为分析、解释套管情况提供正确的数值。

Description

多臂成像数据校正方法

技术领域

本发明涉及一种测量套管形变、腐蚀、沾污、错段等情况的仪器的数据校正方法，具体涉及一种多臂成像数据校正方法。 

背景技术

在测井类型中，有一种工程测井用来测量井下套管情况。目前使用比较多的为具有多个探测臂的多臂井径测量仪器。在测井过程中，张开的各探测臂紧贴在套管内壁，仪器缓慢上提，通过记录各个臂在套管内壁上的轨迹数值，可反映出套管的变化情况，例如：腐蚀、形变、结蜡、沾污、错段等。现在使用的多臂仪器多为36臂、40臂、60臂等。大量的测量数据，可以很方便的进行套管成像处理。然而，直接测量的各探测臂的数据存在许多的偏差，对套管情况分析会存在一些误解。由于仪器的旋转，造成各探测臂位置的不断变化，而数值记录中只是按固定的探测臂顺序进行数据记录，在成像显示中，对于一竖直的套管变形，会造成是旋转变形的误解；另外，由于井斜，达到一定的倾斜度，仪器的扶正臂力量不能完全承受的情况下，会造成仪器中心不居中，各探测臂测量数据异常，会产生一部分探测臂数值偏大，一部分探测臂数值偏小，这些情况，都会给套管的具体状况解释带来偏差。 

发明内容

本发明的目的在于提供一种多臂成像数据校正方法，以尽可能减小测量数据的偏差及错误解释，为分析、解释套管情况提供正确的数值。 

本发明的技术解决方案如下： 

一种多臂成像数据校正方法，在测井仪器中对应于各个探测臂均设置有位移传感器，测井仪器中还设置有用以获取该测井仪器实时方位及倾角信息的辅助测量传感器，其中，所述方位是由水平旋转的绝对或相对角度表达；所述数据校正方法包括两方面的数据校正： 

(1)对由于仪器在井下自然旋转而造成各个臂发生位置移动的方位校正： 

在经过的每个深度测量点上，通过辅助测量传感器得到实时方位信息，将该实时方位信息与第一个深度测量点的方位信息比较，若方位变化净值超过设定的角度阈值，则对测井仪器的各探测臂的测量数据进行方位校正，否则，不予校正； 

对测井仪器的各探测臂的测量数据进行方位校正的计算如下： 

Ang＝Ang₁-Ang₀

M_x＝f_x(a_i，a_i+1) 

R_i＝M_{((i*n)-Ang+360)mod 360}

其中：i为探测臂的序号； 

x为每个角度，数值为1°，2°，3°……，360°； 

Ang₀为仪器初始角度； 

Ang₁为仪器当前角度； 

Ang为角度差，精度为1度； 

M_x为位于某两个测量角度之间的插值数值，每度插一个值；插值完成后，可得到沿圆周360个数值； 

R_i为校正后某探测臂数值； 

n为相邻两臂之间的夹角； 

a_i是序号为i的探测臂的测量数据； 

a_i+1是序号为i+1的探测臂的测量数据； 

mod表示求余； 

f_x为拉格朗日插值公式； 

R_i＝M_{((i*n)-Ang+360)mod 360}表示某个探测臂在圆周上旋转角度Ang后所对应的数值； 

(2)对由于井斜等原因造成的多臂测量仪器不居中而使探测臂的测量数据产生偏差的倾斜校正： 

在经过的每个深度测量点上，通过辅助测量传感器得到实时倾角信息； 

假设将仪器的各探测臂和上、下扶正器水平投影于xy坐标系，其中仪器的套管中心为坐标原点，第一探测臂作为x轴正方向，上扶正器的中心点为A (acosα，asinα)，下扶正器的中心点为B(bcosθ，bsinθ)，探测臂距离上扶正器距离S_A，距离下扶正器距离S_B；计算得到 

$c=\frac{S_B}{S_A+S_B}\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}$

$\cos \mathrm{\β}=\frac{\frac{S_A}{S_B}a\cos \mathrm{\α}+b\cos \mathrm{\θ}}{\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}}$

设测井仪器的倾斜角为 

仪器弯曲引起的仪器侧向形变量ΔL计算式为： 

其中G为仪器的重量； 

设仪器直径为Φ，芯杆长度为S，仪器弯曲后内层芯杆和外层芯杆变形量为ΔS，两个扶正器之间距离为s，探测臂的测量值与距离传感器自身位移量的比例是n，那么由于重力形变引起的探测臂测量值偏差D为： 

偏心方位角为γ； 

根据以上参数计算探测臂中心点P的偏心距离和偏心坐标为 

$p=\sqrt{c^2+{(\mathrm{\ΔL}+D)}^2+2c(\mathrm{\ΔL}+D)(\cos \mathrm{\β}\cos \mathrm{\γ}+\sin \mathrm{\β}\sin \mathrm{\γ})};$

$\cos \mathrm{\η}=\frac{c\cos \mathrm{\β}+(\mathrm{\ΔL}+D)\cos \mathrm{\γ}}{p};$

即得到仪器总体偏心量 然后根据各探测臂的测量数据来进行偏心校正；r_i的模为探测臂测量值，λ_i为每个探测臂相对于第一个探测臂的角度，i为探测臂序号； 

则根据以上参数计算探测臂测量值的校正值为 

$d_i=\sqrt{p^2+r_i^2+2p{r}_i\cos (\mathrm{\η}-{\mathrm{\λ}}_i)}.$

上述述设定的角度阈值为5度。 

上述设定的探测臂的测量数据方位补偿角为相对方位。 

本发明具有以下技术效果： 

本发明针对多臂测量仪器与数值校正算法相结合，根据仪器测量各探测 臂数值及辅助校正数值进行各探测臂数值的校正处理，得到更精确的探测臂数值，从而为分析、解释套管情况提供正确的数值。 

附图说明

图1为多臂井径仪的示意图； 

图2是将仪器的各探测臂和上、下扶正器水平投影于xy坐标系的示意图； 

图3为仪器自重的受力分析示意图； 

图4为仪器受到侧面的压力弯曲时的示意图； 

图5为仪器受到侧面的压力弯曲时形变的计算模型； 

图6为xy坐标系中P点的偏移量示意图。 

具体实施方式

本发明需要测量的主要参数值有：各探测臂测量值、方位值及倾角值。以40臂仪器为例，在仪器上，围绕360度圆周均匀分布着40个独立的探测臂，即每9度一个测量探测臂。每个探测臂均对应一个非接触式位移传感器，在仪器移动的过程中，实时采集探测臂的变化数据，并发送到地面接收器，由地面处理软件对数据进行运算处理。各探测臂的数据为仪器采集的主要数据，但为了提高仪器的精确度，在仪器上也设置了辅助测量传感器：方位及倾角传感器。同时在每一个测量深度点上，方位与倾角的数据也传送到地面处理软件中，对多臂数据进行必要的校正处理。 

方位校正的是仪器在测量过程中，由于旋转造成各探测臂位置的变化而引起的数值误解。在每个深度测量点上，可得到的相对(或绝对)方位，通过与第一个深度测量点上相对(或绝对)方位相比较，可得到仪器旋转后方位变化的差值，若方位变化净值超过设定的最小度数(例如5度)，则对多臂仪器的各探测臂的测量数据进行方位校正，否则，不予校正。 

方位校正方法算法描述： 

Ang＝Ang₁-Ang₀

M_x＝f_x(a_i，a_i+1) 

R_i＝M_{((i*n)-Ang+360)mod 360}

其中：i为探测臂的序号； 

x为每个角度，数值为1°，2°，3°……，360°； 

Ang₀为仪器初始角度； 

Ang₁为仪器当前角度； 

Ang为角度差(精度到1度)； 

M_x为位于某两个测量角度之间的插值数值，每度插一个值；插值完成后，可得到沿圆周360个数值； 

R_i为校正后某探测臂数值； 

n为相邻两臂之间的夹角； 

a_i是序号为i的探测臂的测量数据； 

a_i+1是序号为i+1的探测臂的测量数据； 

mod表示求余； 

f_x为拉格朗日插值公式； 

R_i＝M_{((i*n)-Ang+360)mod 360}表示某个探测臂在圆周上旋转角度Ang后所对应的数值。 

现在以一组数据说明方位算法对仪器探测臂的校正： 

测点1数据：方位为0，各探测臂对应数值为 

209.85，208.57，206.81，208.44，208.40，209.16，208.59，208.53，206.95，209.03，210.21，213.22，215.56，219.71，224.13，227.05，231.26，234.94，236.31，239.60，242.91，245.70，247.82，247.30，251.46，250.55，248.44，247.36，243.99，243.47，240.27，237.93，234.90，232.55，230.28，227.59，221.59，218.13，214.87，212.72。 

在仪器旋转了266度的情况下，未校正前各探测臂的数值为： 

244.18，241.38，237.11，234.96，231.42，226.69，222.61，219.18，214.82，212.24，208.92，208.20，204.95，204.37，205.07，205.95，206.18，207.89，207.36，208.54，209.60，212.55，216.03，218.78，224.64，227.39，230.82，233.32，237.17，240.21，242.40，244.48，247.01，248.78，251.79，251.83，250.01，248.77，248.46，245.77。 

经过方位校正算法校正后的各探测臂的数值为： 

208.62，206.84，204.71，204.66，205.44，206.05，206.89，207.67，207.85，208.98，210.83，214.01，217.18，221.23，225.79，228.82，231.87，234.93，238.44，241.12，243.27，245.53，247.75，250.03，251.81，251.07，249.49，248.64，247.34，245.11，243.02，239.60，236.21，233.48，229.44，224.99，221.18，217.36，213.74，210.86。 

对比第一个测点数据及旋转校正后的各探测臂的数值，明显地看到方位校正对各探测臂数值的修正功能。 

倾角校正方法算法描述及说明：(以40臂仪器为例) 

多臂井径仪在斜井中，由于仪器的自身重量，会产生一定程度的形变。因此造成在测量时，引起仪器测量偏心问题。 

(1)仪器机械加工导致的偏心 

如图1所示，上下扶正器主要使仪器处于套管的正中心，但是由于机械加工的误差等原因，都会造成扶正器本身的偏心。 

建立如图2所示的xy坐标系，假设仪器的第一个探测臂为相对方位原点，仪器从上往下看，逆时针为方位的正方向。 

假定∠α为60度，∠θ为-150度。令套管的中心为图2的坐标原点，上扶正器的中心点为A，下扶正器的中心点为B，A离原点距离为a，B离原点距离为b，上扶正器往α角方向偏距离a，下扶正器往θ角方向偏距离b。如果探测臂距离上扶正器距离S_A，距离下扶正器距离S_B。探测臂的中心点为T，偏离角度为β，偏移距离为c。那么就有： 

A点的坐标为(acosα，asinα)，B点的坐标为(bcosθ，bsinθ)，T点的坐标为(ccosβ，csinβ，)由于仪器是直线的，所以，T点就应该在AB点的连接线上，如图2： 

有： $\frac{\mathrm{TA}}{\mathrm{BT}}=\frac{S_A}{S_B}$

那么： $\frac{a\cos \mathrm{\α}-c\cos \mathrm{\β}}{c\cos \mathrm{\β}-b\cos \mathrm{\θ}}=\frac{a\sin \mathrm{\α}-c\sin \mathrm{\β}}{c\sin \mathrm{\β}-b\sin \mathrm{\θ}}=\frac{S_A}{S_B}$

于是得出： $c=\frac{S_B}{S_A+S_B}\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}$

$\cos \mathrm{\β}=\frac{\frac{S_A}{S_B}a\cos \mathrm{\α}+b\cos \mathrm{\θ}}{\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}}$

可见，仪器的机械结构相对固定，那么(由于机械加工造成的)扶正器的偏心造成探测臂偏心就是固定的。 

(2)斜井中仪器弯曲引起的偏心 

另外，在斜井中，由于仪器的重量，仪器会产生弯曲形变，以上下的扶正器为支点，仪器向下弯曲的程度与仪器的侧面受力情况有关，如图3所示： 

仪器的重量为G，侧面的压力为N，仪器侧向下的拉力为F，仪器的倾斜角为 

有： 

可以认为，仪器的侧向形变量与N的大小成正比，系数为k，系数k与仪器的本身有关系。那么侧向形变量ΔL计算式如下： 

偏心的方向始终指向地心的方向，当仪器在转动的时候，偏心的方向就会跟着转动，这个角度为γ。 

(3)传感器芯杆弯曲导致传感器测量误差 

芯杆长度为S，仪器弯曲后内层芯杆和外层芯杆变形量为ΔS。在仪器受到侧面的压力弯曲时，仪器的形变会使传感器的位移也发生变化，如图4。 

仪器侧向形变位移为ΔL，两个扶正器之间距离为s。在仪器受到较小的形变时，形变后的形状可以近似看作是一个弧，弧对应的弧度为ω，如图5。 

则： $\sin \frac{\mathrm{\ω}}{4}=\frac{\mathrm{\ΔL}}{\sqrt{\frac{s^2}{4}+\mathrm{\Δ}{L}^2}}$

由于： $\frac{s}{2}>>\mathrm{\ΔL}$

那么： $\sin \frac{\mathrm{\ω}}{4}=\frac{\mathrm{\ΔL}}{\sqrt{\frac{s^2}{4}+\mathrm{\Δ}{L}^2}}=\frac{2\mathrm{\ΔL}}{s}$

可以得出传感器芯杆部分弯曲的弧度ε为： 

当ω→0时：sinω＝ω 

那么： $\mathrm{\ϵ}=\frac{4S}{s}\·4\·\frac{\mathrm{\ω}}{4}=\frac{4S}{s}\·4\·\sin \frac{\mathrm{\ω}}{4}=\frac{4S}{s}\·4\·\frac{2\mathrm{\ΔL}}{s}=\frac{32\mathrm{S\ΔL}}{s^2}$

仪器直径为Φ，那么 $\mathrm{\ΔS}=\frac{\mathrm{\ϵ\Φ}}{2}=\frac{16\mathrm{S\ΔL\Φ}}{s^2},$

如果探测臂的测量值和传感器的位移之间的比例是n，那么由于形变引起的偏差D为： 

可以看出，D的偏离方向也是朝着地心的方向，所以它的偏心方位角为γ。 

现在将仪器加工导致的偏心，斜井中仪器弯曲导致的偏心，传感器芯杆弯曲导致传感器测量误差三个因素合并，仪器总体偏心情况如下，偏心点为P，其偏心距离为p，偏心方位角为η，P点的坐标为(pcosη，psinη)，那么： 

(pcosη，psinη)＝(ccosβ，c sinβ)+(ΔLcosγ，ΔLsinγ)+(Dcosγ，Dsinγ) 

可以得出： 

$p=\sqrt{c^2+{(\mathrm{\ΔL}+D)}^2+2c(\mathrm{\ΔL}+D)(\cos \mathrm{\β}\cos \mathrm{\γ}+\sin \mathrm{\β}\sin \mathrm{\γ})}$

$\cos \mathrm{\η}=\frac{c\cos \mathrm{\β}+(\mathrm{\ΔL}+D)\cos \mathrm{\γ}}{p}$

得到仪器总体偏心量 后，可以根据探测臂测量数据来进行偏心校正。 

如图6： 

为探测臂测量值的向量， 

的模为探测臂测量值， 

的方向为探测臂相对于第一个探测臂的角度∠λ_i， 

为校正后探测臂的测量值向量， 的模为校正后探测臂测量值， 

的方向为校正后的向量方向∠ξ_i，i为探测臂序号，那么 ${\stackrel{\→}{r}}_i=(r_i{\cos \mathrm{\λ}}_i,r_i{\sin \mathrm{\λ}}_i),$ ${\stackrel{\→}{d}}_i=(d_i{\cos \mathrm{\ξ}}_i,d_i,{\sin \mathrm{\ξ}}_i),$

${\stackrel{\→}{d}}_i=\stackrel{\→}{p}+{\stackrel{\→}{r}}_i$

那么： $d_i=\sqrt{p^2+r_i^2+2p{r}_i\cos (\mathrm{\η}-{\mathrm{\λ}}_i)}$

一个圆周有40个探测臂，那么每个探测臂之间的弧度为π/20，那么 

λ₁＝0， ${\mathrm{\λ}}_2=\frac{1}{20}\mathrm{\π},$ ${\mathrm{\λ}}_3=\frac{2}{20}\mathrm{\π},$ ${\mathrm{\λ}}_4=\frac{3}{20}\mathrm{\π}......$ ${\mathrm{\λ}}_{20}=\frac{19}{20}\mathrm{\π},$ λ₂₁＝π 

由于第22臂到第40臂已经在第三、四象限，所以 

${\mathrm{\λ}}_{22}=-\frac{19}{20}\mathrm{\π},$ ${\mathrm{\λ}}_{23}=-\frac{18}{20}\mathrm{\π},$ ${\mathrm{\λ}}_{24}=-\frac{17}{20}\mathrm{\π}......$ ${\mathrm{\λ}}_{39}=-\frac{2}{20}\mathrm{\π},$ ${\mathrm{\λ}}_{40}=-\frac{1}{20}\mathrm{\π}$

分别取每个臂的测量值r_i值，根据λ_i的值可得到每个臂的校正值d_i： 

$d_1=\sqrt{p^2+{r_1}^2+2p{r}_1\cos (\mathrm{\η}-{\mathrm{\λ}}_1)}=\sqrt{p^2+{r_1}^2+2p{r}_1\cos (\mathrm{\η}-0)}$

$d_2=\sqrt{p^2+{r_2}^2+2p{r}_2\cos (\mathrm{\η}-{\mathrm{\λ}}_2)}=\sqrt{p^2+{r_2}^2+2p{r}_2\cos (\mathrm{\η}-\frac{1}{20}\mathrm{\π})}$

$d_3=\sqrt{p^2+{r_3}^2+2p{r}_3\cos (\mathrm{\η}-{\mathrm{\λ}}_3)}=\sqrt{p^2+{r_3}^2+2p{r}_3\cos (\mathrm{\η}-\frac{2}{20}\mathrm{\π})}$

$......$

$d_{40}=\sqrt{p^2+{r_{40}}^2+2p{r}_{40}\cos (\mathrm{\η}-{\mathrm{\λ}}_{40})}=\sqrt{p^2+{r_{40}}^2+2p{r}_{40}\cos (\mathrm{\η}+\frac{1}{20}\mathrm{\π})}$

下面是倾角校正前的40臂数据： 

123.18 124.69 126.27 127.46 128.71 130.15 130.39 131.6 132.03 132.07 131.43 131.36 131.11 130.22 128.84 128.07 126.6 124.41 123.48 121.3 119.7 118.17 115.95 114.83 113.33 111.97 111.25 110.53 109.77 109.63 110.21 110.64 111.17 112.46 113.6 114.44 116.36 118.24 119.37 121.33 

仪器出厂测试仪器加工导致的偏心量为2，偏心角度为-0.25π，斜井中仪器引起的偏心系数k为0.044，仪器重量为36kg，测量时仪器相对方位为-0.6π，井斜为40°，根据这些数据校正后的数据如下： 

122 121.79 121.71 121.34 121.19 121.42 120.68 121.16 121.14 121.02 120.52 120.88 121.34 121.42 121.24 121.86 121.94 121.4 122.19 121.75 121.87 121.99 121.33 121.62 121.37 121.06 121.18 121.08 120.7 120.69 121.15 121.22 121.15 121.61 121.71 121.31 121.83 122.16 121.65 121.89 

对比校正前的数据，明显看到校正后的数据没有明显的偏心，充分体现了井斜偏心校正的功能。 

Claims (2)

1.一种多臂成像数据校正方法，在测井仪器中对应于各个探测臂均设置有位移传感器，其特征在于：测井仪器中还设置有用以获取该测井仪器实时方位及倾角信息的辅助测量传感器，其中，所述方位是由水平旋转的绝对或相对角度表达；所述数据校正方法包括两方面的数据校正：

(1)对由于仪器在井下自然旋转而造成各个臂发生位置移动的方位校正：

在经过的每个深度测量点上，通过辅助测量传感器得到实时方位信息，将该实时方位信息与第一个深度测量点的方位信息比较，若方位变化净值超过设定的角度阈值，则对测井仪器的各探测臂的测量数据进行方位校正，否则，不予校正；

对测井仪器的各探测臂的测量数据进行方位校正的计算如下：

Ang＝Ang₁-Ang₀

M_x＝f_x(a_i，a_i+1)

R_i＝M_{((i*n)-Ang+360)mod 360}

其中：i为探测臂的序号；

x为每个角度，数值为1°，2°，3°……，360°；

Ang₀为仪器初始角度；

Ang₁为仪器当前角度；

Ang为角度差，精度为1度；

M_x为位于某两个测量角度之间的插值数值，每度插一个值；插值完成后，可得到沿圆周360个数值；

R_i为校正后某探测臂数值；

n为相邻两臂之间的夹角；

a_i是序号为i的探测臂的测量数据；

a_i+1是序号为i+1的探测臂的测量数据；

mod表示求余；

f_x为拉格朗日插值公式；

R_i＝M_{((i*n)-Ang+360)mod 360}表示某个探测臂在圆周上旋转角度

Ang后所对应的数值；

(2)对由于井斜等原因造成的多臂测量仪器不居中而使探测臂的测量数据产生偏差的倾斜校正：

在经过的每个深度测量点上，通过辅助测量传感器得到实时倾角信息；

假设将仪器的各探测臂和上、下扶正器水平投影于xy坐标系，其中仪器的套管中心为坐标原点，第一探测臂作为x轴正方向，上扶正器的中心点为A(acosα，asinα)，下扶正器的中心点为B(bcosθ，bsinθ)，探测臂距离上扶正器距离S_A，距离下扶正器距离S_B；计算得到

$c=\frac{S_B}{S_A+S_B}\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}$

$\cos \mathrm{\β}=\frac{\frac{S_A}{S_B}a\cos \mathrm{\α}+b\cos \mathrm{\θ}}{\sqrt{\frac{{S_B}^2}{{S_A}^2}{a}^2+b^2+2\frac{S_B}{S_A}\mathrm{ab}\cos (\mathrm{\α}-\mathrm{\θ})}}$

设测井仪器的倾斜角为仪器弯曲引起的仪器侧向形变量ΔL计算式为：其中G为仪器的重量；

设仪器直径为Φ，芯杆长度为S，仪器弯曲后内层芯杆和外层芯杆变形量为ΔS，两个扶正器之间距离为s，探测臂的测量值与距离传感器自身位移量的比例是n，那么由于重力形变引起的探测臂测量值偏差D为：

偏心方位角为γ；

根据以上参数计算探测臂中心点P的偏心距离和偏心坐标为

$p=\sqrt{c^2+{(\mathrm{\ΔL}+D)}^2+2c(\mathrm{\ΔL}+D)(\cos \mathrm{\β}\cos \mathrm{\γ}+\sin \mathrm{\β}\sin \mathrm{\γ})};$

$\cos \mathrm{\η}=\frac{c\cos \mathrm{\β}+(\mathrm{\ΔL}+D)\cos \mathrm{\γ}}{p};$

即得到仪器总体偏心量

然后根据各探测臂的测量数据来进行偏心校正；r_i的模为探测臂测量值，λ_i为每个探测臂相对于第一个探测臂的角度，i为探测臂序号；

则根据以上参数计算探测臂测量值的校正值为

$d_i=\sqrt{p^2+r_i^2+2p{r}_i\cos (\mathrm{\η}-{\mathrm{\λ}}_i)}.$

2.根据权利要求1所述的多臂成像数据校正方法，其特征在于：所述设定的角度阈值为5度。

Images