【技术实现步骤摘要】
本专利技术涉及石油天然气地球物理测井，具体的，涉及一种套变井井径测井居中校正方法、一种套变井井径测井居中校正系统、一种计算机设备及一种存储有计算机程序的计算机可读存储介质。
技术介绍
1、页岩油气开发过程中，套管变形现象比较严重。一般情况下，套变后可通过井径测井表征井眼变化特征，定量评价套管变形程度。现阶段，井径测井表征套变程度技术已比较成熟并取得不错效果，但是当前工程技术发展对井径测井评价技术提出了新要求，如何采用合理的居中校正方法，准确建立套变井眼三维模型并计算常用工具串通过能力成为新的技术难题。
2、目前，常用的套变井眼居中校正方法是几何重心法和最小二乘法。申请号为“cn201511001098.3”、名称为“一种生产测井多臂井径斜井校正处理的方法”的中国专利和文献《超声成像测井井径的椭圆拟合及偏心时间图像的校正》、《基于最小二乘法的成像测井中井径偏心图像校正方法研究》、《基于最小二乘椭圆拟合的井径测量面偏角修正方法》等均是采用最小二乘法进行椭圆形井眼居中校正。但是在实际应用中，这种方式存在以下问题：①套变后井眼形状复杂，不一定是椭圆形，因此基于最小二乘法的椭圆拟合方法不适用；②对于近似椭圆的井眼截面，采用最小二乘椭圆法计算得到的通径一般偏大，而几何重心法计算得到的通径又偏小，严重影响工具通过能力计算。因此，提供套变井井径测井居中校正方法、系统、计算机设备和介质具有重要的意义。
技术实现思路
1、针对现有技术中存在的不足，本专利技术的目的在于解决上述现有技术中存在
2、为了实现上述目的，本专利技术一方面提供了一种套变井井径测井居中校正方法，所述居中校正方法可包括以下步骤：s1、预设井径测量数据向量表征，对于某一深度点k，井径数量为n，井径长度向量表示为rk，井径方位角向量表示为θk；s2、极坐标系向笛卡尔坐标系转换，各井径端点坐标表示为pk；s3、以pk围成的多边形区域的外接矩形an为搜索区域，将搜索区域网格化，网格内部网格点表示为gnij；s4、遍历每个网格点gnij，以其为圆心，计算坐标点pk到各圆心的距离，并将距离最小值作为网格点的内切圆半径，确定每个网格点处最小内切圆，记所有最小内切圆其中的内切圆半径最大值为rn，且记其内切圆圆心为on；s5、将gnij作为新搜索区域的顶点，建立新搜索区域an+1，将新搜索区域网格化，建立新网格点g(n+1),ij，确定新网格点处更新的最小内切圆，确定其内切圆半径最大值为rn+1，确定其内切圆圆心为on+1；s6、设定计算精度error，判断更新后内切圆半径与更新前内切圆半径的误差是否满足误差条件，若不满足，则转入步骤s5继续更新内切圆，若满足，则进行步骤s7；s7、更新后内切圆的圆心为套变井眼居中校正后的中心，更新井径长度r′k和各井径端点坐标p′k。
3、根据本专利技术一方面的一个或多个示例性实施例，所述将搜索区域网格化可包括：将所述搜索区域an离散成二维网格，搜索区域an由m行和n列组成，an的内部网格点gnij的二维坐标表示为{gn,i,j,gn,i,j},(1≤i≤n-1,1≤j≤m-1)，其中，n代表迭代次数，i和j代表二维网格点。
4、根据本专利技术一方面的一个或多个示例性实施例，所述网格点处最小内切圆的半径可表示为：记{inner_rnij}中的最大半径所对应的内切圆的圆心为on，并记半径最大值为rn。
5、根据本专利技术一方面的一个或多个示例性实施例，所述误差条件可包括|rn+1-rn|≤error。
6、根据本专利技术一方面的一个或多个示例性实施例，所述井径长度向量rk可表示为rk＝[r0,r1,…,rn-1]；设一号极板高端方位角为井径方位角向量θk可表示为
7、根据本专利技术一方面的一个或多个示例性实施例，所述极坐标系向笛卡尔坐标系转换可包括：设第k个深度点各井径端点的坐标为pk，各井径端点坐标pk表示为pk＝[x,y]；记所有井径端点x轴坐标为x＝[x0,x1,…,xn-1]；记所有井径端点y轴坐标为y＝[y0,y1,…,yn-1]；根据三角函数公式，井径端点x轴坐标为x＝sin(θk)·rk，井径端点y轴坐标为y＝co s(θk)·rk。
8、根据本专利技术一方面的一个或多个示例性实施例，确定满足所述误差条件的第n+1次迭代后的圆心on+1可为第k个深度点套变井眼居中校正后的中心ok，同时更新井径端点坐标，记更新后的井径端点坐标为p′k；对于第k个深度点的井径端点坐标pk，计算坐标pk与中心ok的距离，作为更新后的井径长度r′k，笛卡尔坐标系xy平面内r′k为：
9、根据本专利技术一方面的一个或多个示例性实施例，所述更新后的井径端点坐标p′k可表示为：p′k＝[x′,y′]；其x轴坐标为：x′＝sin(θ′k)*r′k，y轴坐标为：y′＝cos(θ′k)*r′k；更新后的井径方位角向量θ′k可表示为：坐标转换为当x0>0，则当x0<0，则其中，(x0，y0)为当前一号极板端点笛卡尔坐标。
10、本专利技术的另一方面提供了一种套变井井径测井居中校正系统，所述系统可包括依次连接的数据输入单元、数据处理单元、居中校正单元及绘图显示单元，其中，数据输入单元被配置为输入井轨迹二维坐标和套变后井径二维空间坐标；数据处理单元被配置为整合、关联井轨迹数据和套变后井径测井数据；居中校正单元被配置为可执行上述任意一项所述的套变井井径测井居中校正方法，求取特定深度套变井眼居中校正后的中心；绘图显示单元被配置为显示居中校正单元的计算结果。
11、本专利技术的再一方面提供了一种计算机设备，所述设备包括：处理器；和存储器，存储有计算机程序，当所述计算机程序被处理器执行时，可实现上述任意一项所述的套变井井径测井居中校正方法。
12、本专利技术的又一方面提供了一种存储有计算机程序的计算机可读存储介质，当所述计算机程序在被处理器执行时可实现上述任意一项所述的套变井井径测井居中校正方法。
13、与现有技术相比，本专利技术的有益效果包括以下内容中至少一项：
14、(1)本专利技术提出的套变井井径测井居中校正方法能够适用于非椭圆形井眼截面，计算精度高。
15、(2)本专利技术提出的套变井井径测井居中校正方法为符合套变井眼特征的居中校正方法，为套变后三维井眼建模和工具通过能力计算提供基础数据。
本文档来自技高网...【技术保护点】
1.一种套变井井径测井居中校正方法，其特征在于，所述居中校正方法包括以下步骤：
2.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述将搜索区域网格化包括：将所述搜索区域An离散成二维网格，搜索区域An由M行和N列组成，An的内部网格点Gnij的二维坐标表示为{gn，i，j，gn，i，j}，(1≤i≤N-1，1≤j≤M-1)，其中，n代表迭代次数，i和j代表二维网格点。
3.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述网格点处最小内切圆的半径表示为：
4.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述误差条件包括|Rn+1-Rn|≤error。
5.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述井径长度向量rk表示为rk＝[r0，r1，…，rn-1]；
6.根据权利要求5所述的套变井井径测井居中校正方法，其特征在于，所述极坐标系向笛卡尔坐标系转换包括：
7.根据权利要求6所述的套变井井径测井居中校正方法，其特征在于，确定满足所述误差条件的第
8.根据权利要求7所述的套变井井径测井居中校正方法，其特征在于，所述更新的井径端点坐标P′k表示为：P′k＝[x′，y′]；其x轴坐标为：x′＝sin(θ′k)*r′k，y轴坐标为：y′＝cos(θ′k)*r′k；
9.一种套变井井径测井居中校正系统，其特征在于，所述系统包括依次连接的数据输入单元、数据处理单元、居中校正单元及绘图显示单元，其中，
10.一种计算机设备，其特征在于，所述设备包括：
11.一种存储有计算机程序的计算机可读存储介质，其特征在于，当所述计算机程序在被处理器执行时实现如权利要求1～8中任意一项所述的套变井井径测井居中校正方法。
...【技术特征摘要】
1.一种套变井井径测井居中校正方法，其特征在于，所述居中校正方法包括以下步骤：
2.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述将搜索区域网格化包括：将所述搜索区域an离散成二维网格，搜索区域an由m行和n列组成，an的内部网格点gnij的二维坐标表示为{gn，i，j，gn，i，j}，(1≤i≤n-1，1≤j≤m-1)，其中，n代表迭代次数，i和j代表二维网格点。
3.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述网格点处最小内切圆的半径表示为：
4.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述误差条件包括|rn+1-rn|≤error。
5.根据权利要求1所述的套变井井径测井居中校正方法，其特征在于，所述井径长度向量rk表示为rk＝[r0，r1，…，rn-1]；
6.根据权利要求5所述的套变井井径测井居中校正方法，其特征在于，所述极坐...
【专利技术属性】
技术研发人员：周文高，刘子平，刘伟，管彬，方泽本，刘春亭，屈玲，师耀利，朱炬辉，王素兵，
申请(专利权)人：中国石油天然气集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。