Research on Array Lateral Logging Real-Time Inversions Based on Hybrid Simulated Annealing Algorithms
-
摘要:
测井资料现场实时反演是判断测井质量及进行精细评价的基础。为在保证反演速度的同时,进一步提高反演精度,结合模拟退火算法和马奎特算法,提出了混合模拟退火算法,用以对阵列侧向测井资料进行3参数(地层电阻率、冲洗带半径和冲洗带电阻率)混合反演。研究发现,初始值的选择对混合反演速度具有明显的影响,引入基于电阻率幅度差信息的冲洗带半径初始值选取策略,能够避免反演结果陷入局部最小值,加快反演速度。实例计算结果表明,混合反演算法的反演速度满足实时反演要求,且反演精度相较传统马奎特算法有明显提高;反演结果与试油结果一致,验证了混合反演算法的适用性。混合模拟退火算法为现场阵列侧向测井资料的反演处理提供了新的技术。
Abstract:The real-time inversion of logging data is the basis for judging the field logging quality on-site and the subsequent detailed evaluations. By using the simulated annealing algorithm and Marquette algorithm, the hybrid inversion of three parameters (formation resistivity, radius of flushing zone and flushing zone resistivity) was carried out with the array lateral logging data, so as to further improve the inversion accuracy while ensuring the inversion speed of the conventional real-time inversion methods. The research showed that the selection of initial value had a clear effect on the results of hybrid inversion. A strategy for selecting the initial value of the radius of flushing zone can avoid the inversion result falling into the local minimum based on the resistivity amplitude difference. Further, it can speed up the inversion. The example application results indicated that the inversion speed of hybrid inversion algorithm could meet the real-time inversion requirements. The research also found that the inversion accuracy was higher than that of the traditional Marquette algorithm significantly. Remarkably, the inversion results were consistent with the test results, which verified the applicability of hybrid inversion algorithm. The hybrid simulated annealing algorithm has provided new idea and technical support for the inversion processing of field array lateral logging data.
-
深水钻井与陆地钻井最大的区别是,钻具和钻井液要通过几百甚至数千米长的隔水管系统到达海底实现钻进[1]。隔水管系统上端连接钻井平台,下端通过海底防喷器与井筒相连。钻井过程中,隔水管和井筒与钻柱不断发生碰撞、摩擦等相互作用。因此,对这一管中管系统进行准确的动力学分析和优化设计是确保深水钻井安全的关键,具有重要的工程实际意义。文献调研发现:相关研究大多侧重于对隔水管本身的力学分析[2],一般采用时域分析[3]或频域分析[4–6]来展开。在建立物理模型时,将隔水管系统视为一个上端自由、下端铰接的梁[7],然后用Morison方程[2]来描述其外部所受的海浪、海流等作用。部分模型考虑了隔水管内部钻井液流动的影响[8],有些模型也考虑了其内部钻柱振动的影响[9]。建立的模型大多采用有限差分和有限元法进行求解,并主要分析了顶张力、海水深度、平台漂移、浮力块配置和海况条件等因素对管柱系统横向弯曲变形的影响[10–11]。总体而言,现有研究对深水钻井管柱系统内外管柱相互作用的接触非线性的影响考虑较少,对设计参数组合优化以及优化结果可靠性的分析也不够。为此,笔者建立了管中管结构模型,分析了深水钻井管柱系统的动力学特性,并在此基础上结合可靠度分析进行多目标优化,提出了深水钻井管柱系统的优化设计方法。
1. 深水钻井管柱系统动力学模型
深水钻井管柱系统动力学模型主要包括管中管模型和内外管柱相互作用模型(见图1,图中:μ为摩擦系数;K为碰撞刚度,N/m;N-i和N-i+1为管柱分段编号),前者主要用于描述隔水管和钻柱的弯曲变形以及隔水管外部受到的海洋环境载荷,后者主要用于描述内外管柱之间的相互作用。这2部分统筹起来形成的整体模型包含2个独立坐标系(即O-xyz和O′-x′y′z′),分别用于描述外管和内管系统。
在对深水钻井管柱系统进行动力学建模时,首先建立隔水管和钻柱微元段的运动微分方程,然后将离散的微元段向整体坐标系转化,经拼装得到整个深水钻井管柱系统的有限元模型。就隔水管微元段而言,需对其外部的海洋环境载荷进行描述,采用Ekman漂流理论[1]计算海流流速,然后采用Airy波理论[1]计算海浪的速度和加速度,根据计算的海流和海浪参数,采用Morison方程[2]计算隔水管受到的由海洋环境载荷引起的拖曳力和惯性力。对于隔水管的弯曲变形,采用Euler-Bernoulli梁[2]进行描述。假设管柱系统在y方向发生弯曲,考虑弯曲截面剪切力、张紧力和重力三者在y方向的分量,同时根据达朗伯原理[2]引入惯性力。此外,为进一步考虑内外管柱相互作用的影响,当内外管柱发生接触时,它们之间的相互作用由沿管柱系统引入的N个弹簧–摩擦单元来模拟,并通过设置Heaviside方程来判断是否触发接触位置的弹性力和摩擦力。由此得到隔水管微元段的运动微分方程为:
12CdDρ(uw+uc−∂y∂t)|uw+uc−∂y∂t|+14D2πρ[Cmaw−(Cm−1)∂2y∂t2]+∂2∂z2[EI(z)∂2y∂z2]+p(z)∂2y∂z2+(1−Bf)m(z)g∂y∂z+m(z)∂2y∂t2+H(y−y′−G)[K(z)(y−y′−G)+μ(z)K(z)(y−y′−G)]=0 (1) 其中
H(y−y′−G)={0(y−y′<G)1(y−y′⩾ (2) 式中:Cd为拖曳力系数;D为隔水管外径,m;ρ为海水密度,kg/m3;uw为海浪速度,m/s;uc为海流流速,m/s;Cm为惯性系数;aw为海浪加速度,m/s2; EI为隔水管抗弯刚度,N·m2;p为张紧力,N;m为隔水管的线质量,kg/m;Bf为浮力块浮力系数;g为重力加速度,m/s2;H为Heaviside方程;G为初始间距,m。
同理可得钻柱微元段的运动微分方程为:
\begin{array}{l} \dfrac{{{\partial ^2}}}{{\partial z{{\rm{'}}^2}}}\left[\! {EI{\rm{'}}\left(\! {z{\rm{'}}} \!\right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial z{{\rm{'}}^2}}}} \right] \!+\! p{\rm{'}}\left(\! {z{\rm{'}}}\! \right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial z{{\rm{'}}^2}}} \!+\! m{\rm{'}}\left(\! {z{\rm{'}}}\!\right)g\dfrac{{\partial y{\rm{'}}}}{{\partial z{\rm{'}}}} + \\ m{\rm{'}}\left( {z{\rm{'}}} \right)\dfrac{{{\partial ^2}y{\rm{'}}}}{{\partial {t^2}}} + H\left( {y - y' - G} \right) \left[ {K\left( z \right)}\right. \left( {y' - y - G} \right) +\\\quad\qquad\qquad\left. {\mu \left( z \right)K\left( z \right)\left( {y' - y - G} \right)} \right] = 0 \end{array} (3) 式中:
EI' 为钻柱抗弯刚度,N·m2;m {\rm{'}} 为钻柱的线质量,kg/m;p' 为钻柱悬重,N。设置初始条件时,隔水管与钻柱的初始位移与速度均设为0,即:
\left\{ {\begin{array}{l} \!\!\!{y\left( {z,0} \right) = 0}\\ \!\!\!{\dfrac{{\partial y\left( {z,0} \right)}}{{\partial t}} = 0}\\ {\begin{array}{l} \!\!\!\!\!\!\!{y'\left( {z',0} \right) = 0}\\ \!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {z',0} \right)}}{{\partial t}} = 0} \end{array}} \end{array}} \right. (4) 设置边界条件时,隔水管上端设置为自由、下端固定于井口,上下两端由于存在挠性接头均可旋转;钻柱设置为上端随隔水管偏移、下端受岩石反作用力影响进行轴向振动,上下两端均不能旋转,即:
\left\{ {\begin{array}{l} {y\left( {0,t} \right) = 0}\\ {EI\left( 0 \right)\dfrac{{{\partial ^2}y\left( {0,t} \right)}}{{\partial {z^2}}} = {k_{\rm{b}}}\dfrac{{\partial y\left( {0,t} \right)}}{{\partial z}}}\\ {\begin{array}{l} \!\!\!\!{y\left( {L,t} \right) = S}\\ \!\!\!\!\!{EI\left( L \right)\dfrac{{{\partial ^2}y\left( {L,t} \right)}}{{\partial {z^2}}} = {k_{\rm{t}}}\dfrac{{\partial y\left( {L,t} \right)}}{{\partial z}}}\\ {\begin{array}{l} \!\!\!\!\!\!\!\!\!{y'\left( {0,t} \right) = A\sin \left( {\varOmega t} \right)}\\ \!\!\!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {0,t} \right)}}{{\partial z'}} = 0}\\ {\begin{array}{l} \!\!\!\!\!\!\!\!\!\!\!\!\!{y'\left( {L,t} \right) = S}\\ \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{\dfrac{{\partial y'\left( {L,t} \right)}}{{\partial z'}} = 0} \end{array}} \end{array}} \end{array}} \end{array}} \right. (5) 式中:A为钻头振幅,m;Ω为钻头振动角速度,rad/s;S为隔水管顶部偏移,m;kt为上挠性接头旋转刚度,(N·m)/rad;kb为下挠性接头旋转刚度,(N·m)/rad。
2. 动力学模型的有限元分析
由于建立的动力学模型考虑了内外管柱之间的相互作用,模型的求解难度和计算时间显著增加。为了更加高效,利用有限元软件Abaqus进行动态模拟分析。在Abaqus中,采用B31型梁单元描述内外管柱的弯曲变形。将井筒简化为一个与隔水管外径相等且厚度为0.10 m的梁,认为当其与钻柱接触时发生弹塑性变形。内外管柱之间的相互作用使用general-contact模块中的edge-to-edge模型来设定。隔水管顶部和底部分别与钻井平台和井口铰接,采用Abaqus中的CONN3D2型铰链模拟顶部和底部的铰链,铰链的连接方式选用Joint-Rotation组合。隔水管外部受到的海洋环境载荷利用Aqua模块施加,Aqua模块算法的基础是Morison方程[2]。Aqua模块中的海流为分层运动,南海海洋监测数据分析表明,从海面到海底的海流流速从1.35 m/s到0 m/s呈指数式递减,符合Ekman漂流理论[2]。建立有限元模型后,首先进行静态分析,然后进行时长100个海浪周期的动态分析,以确保管柱系统达到动态稳定状态。
有限元模拟用到的主要参数:隔水管和钻柱的管材密度均为7 850 kg/m3,弹性模量均为210 GPa;钻柱的外径、壁厚分别为149.20和9.17 mm,钻柱的悬重比为0.9;钻头的振幅为0.02 m,隔水管的外径、壁厚分别为533.4和19.1 mm,隔水管的张力比为1.2;内外管柱间的碰撞刚度为18 000 kN/m,摩擦系数为0.05;上、下挠性接头的旋转刚度分别为5 000和7 300 kN·m/rad,浮力块的浮力系数为0.9;海水密度为1 050 kg/m3,钻井液密度为1 200 kg/m3;拖曳力系数为0.7,惯性力系数为2.0;海浪振幅为3.0 m,海浪周期为9.0 s,海面风速为5.0 m/s,风速因子为0.03。基于以上参数,使用Abaqus软件进行数值模拟,结果见图2。其中,图2(a)所示为深水钻井管柱系统达到动态稳定时的空间姿态,为便于观察,管柱横向变形放大了10倍,图中红色部分是管柱横向偏移较大的位置,其主要位于管柱中上部的1/3处。图2(b)为图2(a)的局部截图,展示了内外管柱之间发生碰撞、摩擦时的情况。由于在建立模型时考虑了内外管柱间的相互作用,因此首先比较了这种相互作用对泥线以上管柱整体偏移的影响,以2 000 m水深、1 000 m钻进深度为例进行数值模拟,分别模拟管中管系统和不考虑钻柱的隔水管系统的动态响应,结果见图2(c)。
由图2(c)可知,相比于单独考虑隔水管的情况,将内外管柱的碰撞和摩擦等相互作用纳入模型考虑时,泥线以上管柱系统的整体偏移量受到了明显抑制。这说明对于存在管中管结构的钻井管柱系统,在给定的安全偏移范围内能够承受更大的海洋环境载荷。
3. 参数敏感性分析
由上文所建有限元模型可知,主要设计参数包括隔水管外径与壁厚、顶部偏移量、钻压和顶张力。为了分析这些设计参数对钻井管柱系统动态响应的影响,将图2所示模拟结果作为标准状态,每次只改变1个设计参数重复进行数值模拟,当管柱系统达到动态稳定时,导出计算结果进行分析。5个模拟实例的计算结果与标准状态的对比如图3所示。
由图3可知,隔水管外径与壁厚、顶部偏移、钻压和顶张力等参数对管柱系统的动态响应都有不同程度的影响。增大隔水管壁厚、减小隔水管顶部偏移量、钻压及增大顶张力,都可以减小管柱系统的整体弯曲,并能同时减小隔水管及上、下挠性接头的转角。特别是,由于隔水管外部受到较强海洋环境载荷的影响,其应力分布呈现明显的弯曲效应;而钻柱除了在井口处由于与井壁接触产生应力集中之外,其余主要体现为轴向应力,因此呈现出近似线性分布的特征。
除了前面提到的4项主要设计参数,上、下挠性接头的旋转刚度也是深水钻井管柱系统配置的关键。因此,需分析上、下挠性接头旋转刚度组合对钻井管柱系统动态响应的影响。旋转刚度组合从(10 kN·m/rad,150 kN·m/rad)到(5 000 kN·m/rad,7 300 kN·m/rad)选取4组,模拟其对钻井管柱系统动态响应的影响,结果如图4所示(为了简便,各图例中省去了数值的单位kN·m/rad)。
由图4(a)可知,上、下挠性接头的旋转刚度对管柱系统整体偏移的影响很小,可以忽略;但是,其不同组合对隔水管顶、底两端的受力状态有明显影响(见图4(b)和图4(c));此外,随着上、下挠性接头的旋转刚度增大,其转角在较小范围内波动(见图4(h)和图4(i)),这有助于在钻井过程中将挠性接头的转角控制在允许的最大转角范围之内。因此,在设计深水钻井管柱系统时,需要将上、下挠性接头的旋转刚度组合跟前述主要设计参数一起优化。
4. 基于可靠度分析的多目标优化设计
根据参数敏感性分析结果,探讨隔水管壁厚、顶部偏移量、钻柱悬重、顶张力及上、下挠性接头旋转刚度的优化设计问题。
4.1 优化设计方法
根据深水钻井工程实践,进行管柱设计时应考虑的主要约束条件有:1)管柱的最大Mises应力小于其材料屈服应力的67%,如X80钢屈服应力的67%为370 MPa;2)管柱的最大偏移量小于海水深度的2%,如2 000 m水深下管柱的偏移极限为40 m;3)上挠性接头的最大转角小于5°(即0.087 rad),下挠性接头的最大转角小于2°(即0.035 rad)。
为了提高效率,采用Isight优化分析软件进行多目标优化设计。Isight软件集成了多项成熟的优化分析算法,并且可以将在Abaqus中建立的有限元模型直接嵌入其优化分析的流程之中,并通过识别Abaqus模型中的设计参数和动态响应设置参数范围和动态响应约束条件。基于此,建立了基于可靠度分析的深水钻井管柱系统多目标优化设计方法,设计流程如图5所示。
该设计方法的设计思路是:1)选取邻域培植遗传算法(NCGA)作为多目标优化算法;2)海洋深水环境具有极不稳定的气候特性,需要分析由于系统参数波动所引起的优化方案的可靠性降低问题,因此选取Monte Carlo算法进行可靠度分析;3)为了将多目标优化设计和可靠度分析贯穿起来,采用可靠度优化算法Six Sigma进行算法整合。具体设计步骤为:用Six Sigma法驱动NCGA算法进行多目标优化设计,并将得到的非劣解传递给Monte Carlo算法;然后,用Monte Carlo算法计算该非劣解在参数波动范围内的可靠度,如果该可靠度满足Six Sigma算法设定的最低可靠度要求,则输出该非劣解;反之,则将可靠度不满足要求的情况反馈给NCGA算法,以帮助它调整并优化策略。
4.2 可靠度验证
采用图5所示流程进行了优化设计测试,设计参数取值:隔水管壁厚9.0~28.0 mm,顶部偏移量-20~20 m,隔水管顶张力400~1 500 kN,钻柱悬重600~1 000 kN,上挠性接头旋转刚度100~7 000 kN·m/rad,下挠性接头旋转刚度100~10 000 kN·m/rad。动态响应约束条件:最大偏移量40 m,最大Mises应力370 MPa,上挠性接头最大转角5.0°,下挠性接头最大转角2.0°。
采用Six Sigma法进行可靠度优化运算,将各项动态响应的可靠度最低要求都设为98%,得到的计算结果见表1(不收敛或违反约束的算例标为红色,寻优过程中间解标为白色,非劣解标为蓝色,推荐的最优解标为绿色)。
表 1 基于可靠度分析的多目标优化设计结果Table 1. Results of the multi-objective optimization design based on reliability analysis算例
序号旋转刚度/(kN·m·rad–1) 隔水管
壁厚/mm顶部
偏移/m隔水管
顶张力/kN钻柱
悬重/kN最大
偏移量/m最大Mises
应力/MPa最大转角/rad 下挠性接头 上挠性接头 下挠性接头 上挠性接头 1 5 861. 9 3 005.0 18.8 –8.18 930.1 977.9 15.40 242.4 –0.011 3 0.077 4 2 4 708.8 214.0 16.9 6.01 848.0 914.6 27.94 236.9 –0.024 8 0.090 1 3 4 383. 4 3 542.9 25.0 –12.54 1 391.7 960.8 5.90 226.6 –0.003 1 0.055 4 4 7 194.9 5 758.1 11.2 18.56 1 059.6 805.1 32.06 199.4 –0.018 2 0.059 5 5 2 360.0 1 198.2 10.6 6.25 1 339.5 850.5 21.27 212.5 –0.017 1 0.069 5 6 541.6 1 195.5 23.7 9.35 486.7 966.4 34.55 251.8 –0.038 6 0.096 4 7 10 009.6 4 037.1 15.1 12.03 1 109.4 798.1 27.70 198.4 –0.012 8 0.065 6 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 887 10 314.4 6 971.1 25.1 –16.75 1 468.0 735.7 2.48 163.1 0.009 7 0.048 3 ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ ┇ 995 2 170.9 6 401.2 25.1 –17.35 1 468.4 917.9 0.75 207.8 0.000 4 0.044 1 996 735.7 6 407.1 25.1 –16.76 1 468.7 918.9 1.15 208.0 –0.001 7 0.043 9 997 7 560.4 7 063.5 9.7 –16.75 1 471.6 983.8 0.05 227.9 0.011 5 0.043 4 998 9 299.3 6 880.8 9.7 –16.33 1 474.8 834.2 1.69 191.8 0.012 1 0.049 7 999 7 446.5 6 401.2 27.2 –17.35 1 468.4 754.2 2.18 167.7 0.005 5 0.049 5 1 000 9 375.0 6 554.6 9.7 –16.75 1 471.6 983.8 0.25 229.1 0.013 3 0.045 4 为了验证表1中各个非劣解在设计参数波动情况下的可靠性,对推荐的最优解(第887号算例的解)进行了Monte Carlo可靠度分析,允许各项设计参数在±10%范围内扰动,并将样本点设置为正态分布,最终得到的可靠度分析结果如图6所示。
由图6可以看出,最大偏移量的变化范围为8.35~16.85 m,最大Mises应力的变化范围为148~240 MPa,上挠性接头转角的变化范围为3.14°~4.50°(即0.054 8~0.078 5 rad),下挠性接头转角的变化范围为–0.43°~0.10°(即–0.007 5~0.001 7 rad)。与动态响应约束条件相比,各项动态响应值均符合设置的约束条件,即设计参数在±10%范围内扰动情况下,该设计的可靠度仍达到了100%。由此可知,设计的该参数组合满足要求。
进一步研究可知,多目标优化设计方法可以避免由于设计结果靠近约束边界而在参数波动情况下可靠度达不到要求的问题,因此所得到的非劣解集能够同时满足工程可行性和安全可靠性2方面的要求。
5. 结论与建议
1)相比于仅考虑隔水管的模型,根据管中管模型计算得到的管柱整体弯曲程度更小,即内外管柱的相互作用对管柱系统的整体偏移有抑制作用。这说明对于实际存在管中管结构的深水钻井管柱系统而言,在给定的管柱安全偏移范围内能够承受更大的海洋环境载荷。
2)建立了基于可靠度分析的多目标优化设计方法,采用该设计方法得到的非劣解集能够同时满足工程可行性和安全可靠性2方面的要求。
3)建立管中管模型时,没有考虑钻井液的影响,而钻井液可以等效为一个贯穿于管中管模型的弹塑性减振装置,其会在一定程度上减轻内外管柱之间的相互作用,从而影响管柱系统的整体动态响应,因此需要在接下来的研究中予以考虑。
-
表 1 3种算法的性能对比
Table 1 Comparison of the performances with three algorithms
算法 寻优成功率,% 最优适应度 收敛平均代数 混合模拟退火算法 90 0.004 526 45 模拟退火算法 95 0.003 443 104 马奎特算法 56 0.010 742 10 表 2 混合模拟退火算法反演时长与仪器测量时长对比
Table 2 Comparison of the inversion time of hybrid simulated annealing algorithm with the instrument measurement time
井号 处理层段/
m反演时长/
s单点反演时长/
s仪器单点测量
时长/sW1P–3 1 980~2 010 60 0.200 0 0.500 W1P–7 780~820 82 0.205 0 W1P–9 980~1 020 87 0.217 5 -
[1] 谭茂金, 高杰, 邹友龙, 等. 盐水泥浆条件下定向井双侧向测井环境校正方法研究[J]. 地球物理学报, 2012, 55(4): 1422–1432. doi: 10.6038/j.issn.0001-5733.2012.04.038 TAN Maojin, GAO Jie, ZOU Youlong, et al. Environment correction method of dual laterolog in directional well[J]. Chinese Journal of Geophysics, 2012, 55(4): 1422–1432. doi: 10.6038/j.issn.0001-5733.2012.04.038
[2] 杨震, 杨锦舟, 韩来聚. 随钻电磁波电阻率测井实时井眼影响校正[J]. 石油勘探与开发, 2013, 40(5): 625–629. doi: 10.11698/PED.2013.05.19 YANG Zhen, YANG Jinzhou, HAN Laiju. A real-time borehole correction of electromagnetic wave resistivity logging while drilling[J]. Petroleum Exploration and Development, 2013, 40(5): 625–629. doi: 10.11698/PED.2013.05.19
[3] 朱天竹, 杨守文, 白彦, 等. 利用2.5维数值模式匹配算法高效高精度建立垂直井眼中多分量阵列感应井眼校正库[J]. 地球物理学报, 2017, 60(3): 1221–1233. doi: 10.6038/cjg20170332 ZHU Tianzhu, YANG Shouwen, BAI Yan, et al. Efficient and high-precision establishment of borehole correction database for multicomponent array induction logging in vertical boreholes by a 2.5D NMM algorithm[J]. Chinese Journal of Geophysics, 2017, 60(3): 1221–1233. doi: 10.6038/cjg20170332
[4] 邓少贵, 李竹强, 李智强. 水平井双侧向测井响应及层厚/围岩影响快速校正[J]. 石油勘探与开发, 2009, 36(6): 725–729. doi: 10.3321/j.issn:1000-0747.2009.06.007 DENG Shaogui, LI Zhuqiang, LI Zhiqiang. Response of dual laterolog and fast correction for layer thickness and shoulder bed in horizontal wells[J]. Petroleum Exploration and Development, 2009, 36(6): 725–729. doi: 10.3321/j.issn:1000-0747.2009.06.007
[5] 张冲, 毛志强, 孙中春, 等. 玛河气田盐水泥浆侵入条件下双侧向测井电阻率校正方法[J]. 石油地球物理勘探, 2010, 45(5): 757–761. ZHANG Chong, MAO Zhiqiang, SUN Zhongchun, et al. Dual laterolog resistivity correction method under the conditions of salt mud invasion for Mahe Gas Field[J]. Oil Geophysical Prospecting, 2010, 45(5): 757–761.
[6] 范翔宇, 夏宏泉, 陈平, 等. 测井计算钻井泥浆侵入深度的新方法研究[J]. 天然气工业, 2004, 24(5): 68–70. doi: 10.3321/j.issn:1000-0976.2004.05.022 FAN Xiangyu, XIA Hongquan,CHEN Ping, et al. New method of calculating drilling mud invasion depth by logs[J]. Natural Gas Industry, 2004, 24(5): 68–70. doi: 10.3321/j.issn:1000-0976.2004.05.022
[7] FRENKEL M A. Real-time interpretation technology for new multi-laterolog array logging tool[R]. SPE 102772, 2006.
[8] 王才经. 三侧向测井电阻率反演[J]. 石油大学学报(自然科学版), 1989, 13(4): 75–81. WANG Caijing. Formation resistivity reversion for laterolog 3[J]. Journal of the University of Petroleum, China(Natural Science Edition), 1989, 13(4): 75–81.
[9] 倪小威, 冯加明, 徐观佑, 等. 仪器偏心情况下阵列侧向测井井眼校正图版研制[J]. 断块油气田, 2018, 25(4): 469–472. NI Xiaowei, FENG Jiaming, XU Guanyou, et al. Development of borehole correction plates for array lateral logging under instrumental eccentricity condition[J]. Fault-Block Oil & Gas Field, 2018, 25(4): 469–472.
[10] 冯周, 李宁, 武宏亮, 等. 缝洞储集层测井最优化处理[J]. 石油勘探与开发, 2014, 41(2): 176–181. doi: 10.11698/PED.2014.02.06 FENG Zhou, LI Ning, WU Hongliang, et al. Logging optimization processing method for fractured-vuggy reservoirs[J]. Petroleum Exploration and Development, 2014, 41(2): 176–181. doi: 10.11698/PED.2014.02.06
[11] 顿月芹, 袁建生. 阵列侧向电法测井的快速反演[J]. 清华大学学报(自然科学版), 2009, 49(11): 1871–1875. doi: 10.3321/j.issn:1000-0054.2009.11.032 DUN Yueqin, YUAN Jiansheng. Fast inversion of array lateral electric-logging[J]. Journal of Tsinghua University(Science and Technology), 2009, 49(11): 1871–1875. doi: 10.3321/j.issn:1000-0054.2009.11.032
[12] 蔡军, 张恒荣, 曾少军, 等. 随钻电磁波电阻率测井联合反演方法及其应用[J]. 石油学报, 2016, 37(3): 371–381. doi: 10.7623/syxb201603009 CAI Jun, ZHANG Hengrong, ZENG Shaojun, et al. Joint inversion method of electromagnetic wave resistivity logging while drilling and its application[J]. Acta Petrolei Sinica, 2016, 37(3): 371–381. doi: 10.7623/syxb201603009
[13] SMITS J W, DUBOURG I, LULING M G, et al. Improved resistivity interpretation utilizing a new array laterolog tool and associated inversion processing[R]. SPE 49328, 1998.
[14] 冯加明, 刘迪仁, 倪小威, 等. 定向井阵列侧向测井层厚-围岩影响分析及快速校正[J]. 断块油气田, 2018, 25(5): 593–597, 625. FENG Jiaming, LIU Diren, NI Xiaowei, et al. Impact analysis and rapid correction of layer thickness-shoulder bed for response of array lateral logging in directional wells[J]. Fault-Block Oil & Gas Field, 2018, 25(5): 593–597, 625.
[15] 李智强, 范宜仁, 邓少贵, 等. 基于改进差分进化算法的阵列侧向测井反演[J]. 吉林大学学报(地球科学版), 2010, 40(5): 1199–1204. LI Zhiqiang, FAN Yiren, DENG Shaogui, et al. Inversion of array laterolog by improved difference evolution[J]. Journal of Jilin University(Earth Science Edition), 2010, 40(5): 1199–1204.
[16] 潘克家, 汤井田, 杜华坤, 等. 轴对称地层中高分辨率阵列侧向测井信赖域反演法[J]. 地球物理学报, 2016, 59(8): 3110–3120. doi: 10.6038/cjg20160833 PAN Kejia, TANG Jingtian, DU Huakun, et al. Trust region inversion algorithm of high-resolution array lateral logging in axisymmetric formation[J]. Chinese Journal of Geophysics, 2016, 59(8): 3110–3120. doi: 10.6038/cjg20160833
[17] 王才志, 石广仁, 张丽君. 基于极快速模拟退火算法的地层横波各向异性反演[J]. 石油学报, 2007, 28(1): 67–72. doi: 10.3321/j.issn:0253-2697.2007.01.013 WANG Caizhi, SHI Guangren, ZHANG Lijun. Computing method of shear wave anisotropy of formation based on very fast simulated reannealing[J]. Acta Petrolei Sinica, 2007, 28(1): 67–72. doi: 10.3321/j.issn:0253-2697.2007.01.013
[18] 刘财, 符伟, 郭智奇, 等. 基于贝叶斯框架的各向异性页岩储层岩石物理反演技术[J]. 地球物理学报, 2018, 61(6): 2589–2600. doi: 10.6038/cjg2018L0176 LIU Cai, FU Wei, GUO Zhiqi, et al. Rock physics inversion for anisotropic shale reservoirs based on Bayesian scheme[J]. Chinese Journal of Geophysics, 2018, 61(6): 2589–2600. doi: 10.6038/cjg2018L0176
[19] 吴大奎, 李亚林, 伍志明, 等. 地震、测井资料联合反演的神经网络算法研究[J]. 天然气工业, 2004, 24(3): 55–57. doi: 10.3321/j.issn:1000-0976.2004.03.016 WU Dakui, LI Yalin, WU Zhiming, et al. Research on the neural network algorithm of joint inversion of seismic and log data[J]. Natural Gas Industry, 2004, 24(3): 55–57. doi: 10.3321/j.issn:1000-0976.2004.03.016
[20] 刘丽丽, 邢光龙. 基于支持向量机的电磁波传播测井反演参数约束条件提取方法[J]. 测井技术, 2008, 32(3): 215–219. doi: 10.3969/j.issn.1004-1338.2008.03.005 LIU Lili, XING Guanglong. Constraint extraction of inversion parameters based on support machine vector for electromagnetic propagation logging[J]. Well Logging Technology, 2008, 32(3): 215–219. doi: 10.3969/j.issn.1004-1338.2008.03.005
[21] 汪宏年, 陶宏根, 其木苏荣, 等. 水平层状介质中双侧向资料的全参数正则化迭代反演与应用[J]. 地球物理学报, 2002, 45(增刊1): 387–399. WANG Hongnian, TAO Honggen, CHIMED Surong, et al. Regularized entire-parameter iterative inversion of dual-laterolog in horizontally stratified media and its application[J]. Chinese Journal of Geophysics, 2002, 45(supplement 1):387–399.
[22] 倪小威, 徐观佑, 敖旋峰, 等. 阵列侧向测井曲线极化角影响因素研究[J]. 石油钻探技术, 2018, 46(2): 120–126. NI Xiaowei, XU Guanyou, AO Xuanfeng, et al. The influencing factors on the polarizing angle of array laterolog curves[J]. Petroleum Drilling Techniques, 2018, 46(2): 120–126.
[23] 倪小威, 徐思慧, 别康, 等. 不同井眼偏心距下水平井阵列侧向测井围岩校正研究[J]. 石油钻探技术, 2018, 46(4): 121–126. NI Xiaowei, XU Sihui, BIE Kang, et al. Surrounding rock influence correction for array laterolog responses with borehole eccentricities in horizontal wells[J]. Petroleum Drilling Techniques, 2018, 46(4): 121–126.
[24] 余昭平, 燕善俊, 张文政, 等. 基于模拟退火算法的DES最佳线性逼近[J]. 计算机工程, 2006, 32(16): 158–159. doi: 10.3969/j.issn.1000-3428.2006.16.060 YU Zhaoping, YAN Shanjun, ZHANG Wenzheng, et al. Best linear approximations of DES based on simulate annealing algorithms[J]. Computer Engineering, 2006, 32(16): 158–159. doi: 10.3969/j.issn.1000-3428.2006.16.060
[25] 王杜娟, 刘锋, 王延章. 恶化效应下加工时间可控的新工件到达干扰管理[J]. 系统管理学报, 2016, 25(5): 895–906, 913. WANG Dujuan, LIU Feng, WANG Yanzhang. Disruption management for new jobs arrivals with deteriorating effect and controllable processing times[J]. Journal of Systems & Management, 2016, 25(5): 895–906, 913.
[26] 王璞, 吴国忱, 李伟. 基于改进模拟退火算法的横波速度求取[J]. 断块油气田, 2015, 22(3): 330–333. WANG Pu, WU Guochen, LI Wei. Shear-wave velocity prediction based on revised and simulated annealing algorithm[J]. Fault-Block Oil & Gas Field, 2015, 22(3): 330–333.
[27] 顾元宪, 项宝卫, 赵国忠. 一种改进的连续变量全局优化模拟退火算法[J]. 系统工程理论与实践, 2005, 25(4): 103–109. doi: 10.3321/j.issn:1000-6788.2005.04.016 GU Yuanxian, XIANG Baowei, ZHAO Guozhong. An improved simulated annealing algorithm for global optimization problems with continuous variables[J]. Systems Engineering—Theory & Practice, 2005, 25(4): 103–109. doi: 10.3321/j.issn:1000-6788.2005.04.016
[28] 冯思臣, 王绪本, 阮帅. 一维大地电磁测深几种反演算法的比较研究[J]. 石油地球物理勘探, 2004, 39(5): 594–599. doi: 10.3321/j.issn:1000-7210.2004.05.019 FENG Sichen, WANG Xuben, RUAN Shuai. Comparison among several inversion algorithms of 1D MT[J]. Oil Geophysical Prospecting, 2004, 39(5): 594–599. doi: 10.3321/j.issn:1000-7210.2004.05.019
[29] 汪勇, 金菲, 张瑞军. 引导函数支配的进化模糊聚类算法[J]. 系统工程理论与实践, 2011, 31(2): 302–307. doi: 10.12011/1000-6788(2011)2-302 WANG Yong, JIN Fei, ZHANG Ruijun. Evolutionary fuzzy clustering algorithm dominated by guided function[J]. Systems Engineering—Theory & Practice, 2011, 31(2): 302–307. doi: 10.12011/1000-6788(2011)2-302
-
期刊类型引用(13)
1. 马永乾,赵鹏,畅元江,王仕超,张晏铭. 基于LAGRANGE方程的深水钻井隔水管–水下井口系统动力分析. 应用科技. 2024(01): 151-157 . 百度学术
2. 王文昌,徐祖凯,周星,王昭彬,陈锋. 超深井钻柱动态疲劳失效特征及参数优选. 石油钻探技术. 2024(02): 118-125 . 本站查看
3. 畅元江,王仕超,李健,刘秀全. 深水隔水管-防喷器组耦合系统动力学建模与仿真验证. 振动与冲击. 2024(11): 1-8 . 百度学术
4. 狄勤丰,骆大坤,秦垦,王文昌,陈锋. 考虑“筒中筒”耦合作用的深水隔水管柱与钻柱碰摩运动规律. 石油钻采工艺. 2024(01): 25-32 . 百度学术
5. 王高伟,廖茂林,蒋宏伟,郜志英,房超. 基于钻井液上返流动作用下的隔水管动力学行为研究. 船舶力学. 2023(05): 719-730 . 百度学术
6. 孙巧雷,张辉易,姚力萍,李乐勤,张圆圆,冯定. 海上管中管结构力学研究现状与发展. 科学技术与工程. 2023(21): 8907-8915 . 百度学术
7. 曹峰,胡波,王安伟,刘永亮. 压裂泵的数值分析及仿真标定验证. 机械设计与研究. 2023(06): 205-210 . 百度学术
8. 李旭,朱志强,冯颖韬,赵琥,温达洋,山永林,林黎明. 海洋钻井隔水管系统动力学行为研究. 钻采工艺. 2022(03): 31-36 . 百度学术
9. 张昊,雷雨辰,高云飞,张维维,闫启方,刘林超. 分数黏弹性可压缩油气井管杆准静态响应分析. 信阳师范学院学报(自然科学版). 2021(02): 331-338+344 . 百度学术
10. 李泉新,褚志伟,许超,杨冬冬. 煤矿井下双动力复合定向钻进轨迹调控研究. 工矿自动化. 2021(12): 25-31 . 百度学术
11. 王江帅,李军,柳贡慧,罗晓坤. 气侵条件下新型双梯度钻井环空出口流量变化规律研究. 石油钻探技术. 2020(04): 43-49 . 本站查看
12. 张光伟,乔阳,田帆,高嗣土,尹福来. 全旋转内置式导向钻井工具动力学分析及仿真. 机械强度. 2020(06): 1417-1423 . 百度学术
13. 董学莲,杨芯萍,赵雄,张聪林,郑亮,费凌. 液缸式隔水管张紧系统建模及性能分析. 液压与气动. 2019(08): 39-48 . 百度学术
其他类型引用(9)