井内流动与传热的三维耦合数值模型

钻井过程中,井内的压力、流场及温度分布是钻井作业必需的重要参数。本文提出了一个描述井内流动和传热的三维数学模型,用SIMPLER方法获得了该数学模型的数值解。该数值模型考虑了温度、压力对钻井液的流变性、密度以及各种介质的热物理参数的影响。现场试验数据与数值模拟结果的对比表明：本文提出的数学模型及采用的数值计算方法是可行的。     

裂缝气液重力置换的流动实验及仿真

通过可视化模拟实验对真实裂缝下气液重力置换现象进行了模拟,得到了其气液界面现象与常规平板裂缝的差异;通过对压力、置换量等数据的处理证实了气液重力置换的机理;根据实验的液体漏失流量数据提出了钻井过程中漏失、置换、溢流区间;最后根据k-ε双方程模型进行了气液重力置换流动的仿真计算,得到三维裂缝流道中都是钻井液下部侵入,呈现一个近似的直角三角形形状。本文对现场钻井工程方案设计及安全钻井具有指导意义。     

Criticality of bifurcation in the tuned axial–torsional rotary drilling model

We study the bifurcation characteristics of a lumped-parameter model of rotary drilling with 1:1 internal resonance between the axial and the torsional modes which leads to the largest stability thresholds. For this special case, the two-degree-of-freedom model for the drill-string reduces to an effectively single-degree-of-freedom system facilitating further analysis. The regenerative effect of the cutting action due to the axial vibrations is incorporated through a delayed term in the cutting force with the delay depending on the torsional oscillations. This state dependency of the delay introduces nonlinearity in the current model. Steady drilling loses stability via a Hopf bifurcation, and the nature of the bifurcation is determined by an analytical study using the method of multiple scales. We find that both subcritical and supercritical Hopf bifurcations are present in this system depending on the choice of operating parameters. Hence, the nonlinearity due to the state-dependent delay term could both be stabilizing or destabilizing in nature, and the self-interruption nonlinearity is essential to capture the global behavior. Numerical bifurcation analysis of a global axial–torsional model of rotary drilling further confirms the analytical results from the method of multiple scales. Further exploration of the rotary drilling dynamics unravels more complex phenomena including grazing bifurcations and possibly chaotic solutions.     

Fluid-solid coupling model for studying wellbore instability in drilling of gas hydrate bearing sediments

As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.     

Modeling and experimental investigation of drilling into lunar soils

Dry drilling only with the assistance of an auger is a reliable and realistic approach to remove abundant soils from the side of a bit in the harsh, dry conditions on the Moon. Based on an elementary analysis, using Janssen's model to reflect the coupling effect among the different components of the stress, the present paper models the conveying dynamics along the helical groove and the sampling mechanism in the centering hole of the stem for an auger drilling into lunar soil simulant. Combining the two parts as well as a simple cutting model for the bit, a whole drilling model is established to investigate the complicated relation among the conveying ability of the auger, the coring rate, and drilling parameters such as the penetration and rotation speeds. The relation is revealed by the complicated transition between different sub-models with the help of the physical transition conditions. A series of experiments with constant penetration and rotation speeds are conducted to verify the model. Three aspects of characteristics of the drilling dynamics are manifested,(i) the loads on the bit are almost independent of penetration;(ii) three obvious drilling stages with respect to cut per revolution are grouped;(iii) a linear relationship is found between the coring rate and the revolution per penetration.     

旋转导向钻具组合力学分析

旋转导向系统代表了当今定向钻井的先进水平,是井眼轨道控制技术的发展方向. 为了解决旋转导向钻井轨道控制问题,进行了旋转导向钻具组合(rotary steering bottom hole assembly, RS-BHA) 力学分析. 应用纵横弯曲法建立了典型的柔性RS-BHA 三维力学分析模型,得到了钻头侧向力和导向参数间的相互关系,进而根据极限曲率法可预测旋转导向工具的造斜能力,并计算设计井眼轨道需要的导向参数. 通过力学分析,可以优化RS-BHA,了解井眼轨道控制规律,为旋转导向工具设计和定向钻井自动控制提供理论依据.     

STUDY OF CREEP PARAMETERS OF GRANITE IN STEADY STATE UNDER CONDITIONS OF HIGH TEMPERATURE AND HYDROSTATIC PRESSURE

高温高压花岗岩钻孔实验表明,温度低于500℃,静水应力低于150MPa 状态下,岩体的钻孔变形均属于稳态蠕变变形阶段. 该文选取了广义开尔文模型来反映其特征,通过拉普拉斯变换及逆变换,详细推演出了钻孔径向位移解析解,并且考虑温度-应力的耦合效应,给出了模型参数随温度及应力变化的关系式. 利用该关系进行拟合计算,说明广义开尔文模型来表达高温高压环境中的花岗岩稳态蠕变变形特性,寻求蠕变参数是合理可靠的. 该文对于高温岩体地热资源开发中的钻孔施工与维护、钻孔变形量预测等方面,具有实际的指导意义.     

Numerical investigation on dynamic characteristics of drilling shaft in deep hole drilling influenced by minimal quantity lubrication

The dynamic characteristics of drilling shaft in deep hole drilling influenced by minimal quantity lubrication (MQL) is investigated. According to the features of the compressible fluid Reynolds equation in oil/air feature, a time-dependent mathematical model is established to describe the pressure distribution of cutting fluid with nonlinearity in MQL deep hole drilling. By introducing the differential transformation approach, the time-dependent pressure equation arising from cutting fluid is solved by the use of direct integral method. The influences of the rotational speed, the transverse displacement ratio, and radial clearance on the hydrodynamic pressure distribution of cutting fluid are obtained. The advantage of this method is to overcome much of the computational cost and has its rapid convergence rate. Furthermore, the nonlinear responses of drilling shaft influenced by MQL are analyzed, and the instability rotational speeds of drilling tool are discussed while the design parameters of drilling shaft system changing.     

Equivalent boundary model of lunar soil drilling simulation by DEM

Aiming to solve the computational cost problem in the discrete element simulation for lunar soil drilling sampling, an equivalent boundary method was proposed. A high-accuracy DEM model of lunar soil was established firstly. As the novel alterable constitutive law, the accuracy of the model was verified to meet the performance of real lunar soil very much both in shear strength indices and elastic–plastic behavior. A common drill bit in the geological exploration field for sampling soil was chosen as the simulation object. In preanalysis, it was known that with the increase of drilling depth, the stress concentration area was always near the drill bit, while the affected area of the lunar soil was a cylindrical area around the drill pipe, which extended towards the drilling direction instead of extending around it. Then a big boundary drilling simulation scene was established to investigate the flow direction of lunar soil particles. The motion law of particles and the velocity field information were obtained, and a U-shape chain was described around the drill bit. Finally an equivalent boundary was set near the U-shaped chain, and the size was determined by comparing the soil stress in the fierce collision zone and around the reference boundary. This method could be a reference for other lunar soil drilling researches with other drills of different sizes.     

钻井完井液浸泡弱化页岩脆性机制

页岩脆性是页岩地层钻井、水力压裂设计的关键参数之一,目前针对钻井过程中工作液浸泡对页岩脆性的影响还未引起关注.通过开展钻井完井液浸泡前后页岩三轴力学实验,利用脆性评价模型分析了页岩脆性变化特征.结果表明,延长组页岩脆性强于龙马溪组页岩;油基和水基钻井完井液浸泡均能导致页岩脆性降低,且油基钻井完井液浸泡后的页岩脆性降低幅度更大;龙马溪组页岩浸泡后脆性减弱幅度较延长组页岩大.页岩脆性弱化机制包括:(1)由层理面胶结强度不同引起的脆性强弱差异;(2)由毛管自吸作用导致的高孔压、高应力强度因子、低临界断裂韧性;(3)由碱液侵蚀导致的页岩溶蚀孔形成及矿物颗粒碎裂;(4)由黏土矿物水化膨胀产生的膨胀应力;(5)由钻井完井液滤液润滑导致的页岩破裂面摩擦系数降低.延长组页岩层理面强度较龙马溪组页岩低,导致延长组页岩脆性强于龙马溪组页岩.其次,和水基钻井完井液相比,油基钻井完井液具有更大的自吸量、更高的pH值、更强的润滑性,因此,油基钻井完井液浸泡降低页岩脆性幅度更大.另外,由于龙马溪组页岩具有更小的润湿角、更强的毛管自吸和碱液侵蚀作用,相同浸泡条件下,龙马溪组页岩脆性降低幅度更大.本研究可为页岩地层钻井液性能优化、井壁稳定控制、水力压裂设计等提供理论指导.      

Evaluation for energy-saving effect of hybrid drilling rig system based on the logic threshold method

A simulation model of the current rotary drilling rig (RDR) system was built according to the basic configuration of the TR160D-type rotary drilling rig, and the test-rig based on the working principle of the RDR system was built. The results show that simulation model of the RDR system was accurate and reasonable by means of comparative analysis between theoretical computation and test data. The working principle and logic threshold control method of the HDR system were presented to reduce the fuel consumption and improve energy-saving efficiency via the analysis of energy flow of the RDR system, which mainly includes engine, electric motor/generator, super capacitor, AC–DC rectifier/DC–AC inverter and torque coupler, etc. Finally, the comparative analysis of the fuel economy between the RDR and the hybrid drilling rig (HDR) system was completed in the same overall typical operating conditions. The results show that energy-saving efficiency of HDR system is at 18.8% and total fuel economy efficiency of HDR reaches 15.9% relatively in a typical operating condition.     

结构力学分析方法在南海水平井钻进预测中的应用

本文探讨了结构力学分析方法-非线性有限元性在预测水平井钻进中的工程应用,在钻探水平井尤其是大位移水平井中,事先较准确地预测出钻柱的钻进扭矩及钻头前进的方向,是取得水平井成功的关键条件之一。利用非线性有限元理论,着重分析了水平井实际钻进中的钻柱大变形非线性及钻柱与井壁的接触非线性问题,并导出了一套非线性有限元钻柱力学模式,在此基础上编制了一套非线性力学计算程序软件,并在海洋石油南海西部石油公司一个钻井平台水平井钻井进行了五口井的实际验证,获得了一些有启迪性的结论。     

近钻头钻具多源动态姿态组合测量方法

在导向钻井系统的姿态测量过程中,由于近钻头强振动的影响,导致姿态参数测不准甚至不可测,为了消除有规律的干扰、振动等对测量准确性的影响,快速解算出准确的钻具姿态,提出一种新的多源动态姿态组合测量方法。采用三轴加速度计、三轴磁通门以及角速率陀螺仪等构成测量系统,建立基于四元数的姿态测量非线性模型,研究钻具运动状态与振动加速度之间的关系,根据模型及噪声特性,采用基于四元数的无迹卡尔曼滤波方法对振动干扰信号进行滤除。试验结果表明,采用提出的方法能够消除近钻头干扰对姿态参数测量的影响,井斜角在5.2°左右,工具面角误差小于10°,有效地提高了导向钻井工具姿态动态测量的准确性。     

钻井液流体特性对井下涡轮性能的影响

涡轮结构广泛应用于钻井工具,在对其进行仿真分析研究过程中通常采用清水作为介质,但该方法并不能比较确切地反映涡轮结构在井下实际的工作特征。针对该问题,本文建立了某型涡轮三维模型,采用Turbo Grid软件确定合理的网格参数,并进行了网格划分。采用Herschel-Bulkley钻井液流变模型研究了钻井液流体参数对涡轮结构的输出扭矩和效率的影响。基于响应面法,采用Box-Behnken设计方法得到了钻井液流体参数与涡轮性能指标的关系。结果表明,相较于稠度指数和流变指数而言,钻井液的密度和动切力对涡轮结构的输出扭矩和效率影响显著。因此,在探讨钻井液流变参数对涡轮性能影响时,应该优先考虑钻井液的密度和动切力特性。     

波动钻压下阶梯钻柱系统稳定性研究

处于狭长井筒中的钻柱,其动力响应受到钻具组合、内外钻井液流动以及钻井参数等因素的影响,钻柱动力失稳导致的剧烈振动是井壁坍塌和钻具失效的重要原因。考虑到钻杆和钻铤在刚度和线密度上存在很大的差别,论文将钻柱简化为单阶梯输液管柱,钻井液沿着钻柱内部向下泵入并从环空返回地面。耦合考虑钻柱自重、随时间简谐变化的波动钻压、稳定器以及钻井液的水动力和阻尼力,建立了直井中钻柱横向振动的解析模型。利用有限单元法离散为四阶常微分方程后,采用Bolotin法得到临界频率方程确定系统的不稳定区范围,研究了钻压、钻杆长度、稳定器安装位置、钻井液的流速和密度等参数对系统稳定性影响的机理。研究表明：钻压的平均值和波动幅值都是钻柱失稳的驱动因素,而系统的稳定性对处于受拉状态的钻杆的长度变化不敏感。在论文所研究的参数范围内,降低钻井液流速和密度、下移稳定器的安装位置均有助于增强系统的稳定性。     

Coupled axial-torsional dynamics in rotary drilling with state-dependent delay: stability and control

Nonlinear motions of a rotary drilling mechanism are considered, and a two degree-of-freedom model is developed to study the coupled axial-torsional dynamics of this system. In the model development, state-dependent time delay and nonlinearities that arise due to dry friction and loss of contact are considered. Stability analysis is carried out by using a semi-discretization scheme, and the results are presented in terms of stability volumes in the three-dimensional parameter space of spin speed, cutting depth, and a cutting coefficient. These stability volume plots can serve as a guide for choosing parameters for rotary drilling operations. A control strategy based on state and delayed-state feedback is presented with the goal of enlargening the stability region, and the effectiveness of this strategy to suppress stick-slip oscillations is illustrated.     

Effect of the drilling mud filtrate temperature on the resistivity of the stratum saturated by oil and gas

A mathematical model of the axisymmetric distribution of the phases in the zone of invasion of the water-based drilling mud into the productive stratum whose porous space can simultaneously contain three immiscible fluids (oil, gas, and natural water) is constructed; the model takes into account the high rate of heat transfer between the fluids and the rock matrix. It is shown that the resistivity of the invaded zone depends not only on saturation of the latter by the fluids and the concentration of salts in the water phase, but also on the drilling mud filtrate temperature. It is also shown that there is a jump in the function of stratum saturation by oil on the thermal front.     

Study on steady surge pressure for Yield-Pseudoplastic fluid in a concentric annulus

Surge pressure caused by pipe moving through the wellbore full of drilling fluid is the main factor affecting wellbore stability, on which the additional mud density is designed. Thus it is desired to predict the surge pressure with high accuracy. While this is closely related with the selection of mud rheological models. Yield-Pseudoplastic model which has three parameters and is famous for its higher accuracy has been used to describe the mud rheological properties in recent years. Based on this model the paper presents a new theoretical model for calculating surge pressure caused by mud viscosity during tripping or casing in a concentric annulus of directional wells under steady laminar situation. For convenience, the paper plots the distribution of the coefficients of surge pressure for different conditions. An example is given. These results provide the basis for controlling pressure surges and tripping velocity.     

Torsional Vibration Control and Cosserat Dynamics of a Drill-Rig Assembly

Aspects of drill-string vibrations in the context of a recently developed integrated model of a drill-rig assembly based on the Cosserat theory of rods are discussed. Computer simulations are used to compare existing rotary feedback strategies currently in use to optimise drilling performance where torsional slip-stick vibrations are a hazard. Guided by the wave nature of axial and torsional vibrations in axially symmetric drill-string configurations, we present a new control mechanism, torsional rectification, and compare its performance with existing controllers within the context of the model. The practical guidelines for the improvement of drilling rates in a wide variety of circumstances are discussed.