基于在线监测的井筒压力控制智能决策方法、设备及可读存储介质与流程

本发明属于石油天然气勘探开发钻井领域，具体地说涉及一种基于在线监测的井筒压力控制智能决策方法。

背景技术：

1、随着石油钻井对安全作业愈加重视，控压钻井及精细控压钻井技术因其在降低钻井作业及井控风险、提高钻井效率的优势，在川渝、长庆、塔里木等部分油田逐步规模化应用。控压钻井及精细控压钻井技术其原理是通过人工或自动调节井口套压控制井筒压力，实现井筒压力与地层压力快速平衡，有效解决因安全密度窗口窄、压力系统多造成的井漏、井涌、垮塌等钻井复杂问题。

2、长庆油田油区注水开发为主体开发模式，导致高压注水区、油气叠合区储层压力异常，储层出水、出水/漏失同存等复杂问题普遍存在，钻井井控风险高、窄安全密度窗口钻进困难。同时，长庆油田为确保实现油气当量7000万吨，需规模化部署新井实现增产与稳产，年钻井超过4500口，安全钻井技术成为油田快速有效开发的基础，其中油田高压注水区、油气叠合区年部署新井近3000口井，存在异常高压出水、高气油比、窄密度窗口等技术难题，导致钻井施工难度增加、井控安全面临严峻挑战。

3、近年来，川庆钻探针对致密油注水开发区存在的钻井难题，攻关研发了具有长庆特色和成本优势的控压钻井关键技术及装备，有效降低了钻井井控风险，缩短了复杂处置时间，提高了窄安全密度窗口地层钻进能力，截止2023年8月应用1397口，出水复杂停钻处置控制在3h内，较常规钻井减少85％，单井减少出水500m3以上，钻进时间缩短20％，其中2023年截至8月已应用625口井，控压工作量逐年增多。但多为人工判断与手动控压，亟需一种基于现场参数在线监测的井筒压力控制智能决策方法及系统，准确计算钻井过程中井筒环空压力，随钻分析井筒压力与井况，实现井口控压值自动推荐控制和控压决策建议，推动井控智能化、全自动化发展，提高控压钻井的可靠性、时效性，提升窄安全密度窗口地层钻进能力。

4、现有技术中，公开号为cn114482885a，公开日为2022年5月13日的中国发明专利文献，公开了一种控压钻井智能控制系统，包括在线监测模块、智能控制模块和人为控制模块；在线监测模块分别与智能控制模块和人为控制模块连接；在线监测模块用于实时采集钻井过程中的钻井工况数据；智能控制模块用于收集钻井参数及工程技术数据进行整合分析，形成样本库，并计算实际钻井工况的数据向量，结合样本库识别出钻井作业模式，根据钻井作业模式对钻井压力进行控制调节；人为控制模块用于钻井工作人员根据钻井工况数据和钻井作业模式对钻井压力进行人为调控。本发明通过机器学习智能算法可实时得出合理的调控指令，并自动将其传输至执行机构，可在实际运行过程中，根据在线监测数据复现或模拟人为控压决策，实现智能控压。但是上述技术在使用过程中，一方面，该技术作业工况识别少，仅3种(正常钻进、接单根、起下钻)，形成的3种工况下的控制模式不能完全满足各种作业工况的控压作业，控制精准度较差；智能化较差，需要辅以人工确认选项，另一方面该技术方案为提供实现井底压力、井筒压力、控压值的高精度计算方法，仅所利用高精度钻井水力学计算系统计算井底压力；专利采用大数据的方法利用样本库进行对比进行控制指令下达，未进行理论以及实际井筒状况、复杂环境下的校正，难于满足井筒状况、复杂多变的影响，为提供多种模式下的智能决策推荐，人为控制模块用于钻井工作人员根据钻井工况数据和钻井作业模式对钻井压力进行人为调控的决策措施有限，精准度不高。

技术实现思路

1、为解决上述技术问题，针对现有控压系统在无井底压力测量工具的条件下实时分析计算能力有限，无法实现准确的控压作业实时分析计算与控制决策推荐的现状，本发明提供一种基于在线监测的井筒压力控制智能决策方法，在控压钻井过程中，与采集的钻井录井参数信息结合，准确识别水侵或油侵，确定环空流型，自动判断优选钻井液流变模式，预测地层流体侵入后井筒压力的变化，确定合理的井口回压值和推荐控压施工建议，实现全自动控压作业，提高控压钻井井口回压控制的可靠性、时效性及理论性，进一步解决长庆致密油注水开发区储层出水控制难、窄安全密度窗口地层钻进复杂等技术难题，为致密油注水开发区以及其他需要应用控压钻井解决的同类复杂区域、风险井的安全钻井提供技术支持。

2、本发明是采用如下技术方案实现的：

3、一种基于在线监测的井筒压力控制智能决策方法，包括以下步骤：

4、s1、根据现场采集的实时工程录井和控压数据，利用分支树判断方法实现钻井工况识别；

5、s2、结合现场采集的实时工程录井和控压数据，建立环空多相流模型，利用环空多相流模型计算当前钻井参数情况下的井筒压力，并反推出平衡地层压力所需要的井口回压预控值；

6、s3、根据控压钻井操作规范及行业标准，设置不同钻井工况下的控压钻井不同控压级别的控压界限值，针对不同钻井工况，修正井口回压预控值，得到不同钻井工况下的井口回压推荐值，根据井口回压推荐值实现井筒压力控制智能决策控制。

7、所述钻井工况包括控压钻进下的工况和常规钻进下的工况，所述控压钻进下的工况和常规钻进下的工况共22种，均包括起钻、下钻、旋转钻进、滑动钻进、划眼、倒划眼、循环、循环划眼、循环倒划眼、灌泥浆及悬停。

8、所述s1中，所述实时工程录井和控压数据包括日期、时间、控压值、井深、大钩载荷、钻压、转速、进口钻井排量、出口钻井排量、灌浆排量、扭矩、进口钻井液密度、出口钻井液密度、大钩高度、大钩速度、气测(全烃、轻烃、重烃)、氯根在内的测量参数。

9、所述s1中，利用分支树判断方法实现钻井工况识别的步骤为：

10、s11、判断控压值是否大于0，如果是，则后续方法识别出的各种工况属于控压钻进下的工况；如果否，则后续方法识别出的各种工况属于常规钻进下的工况；控压钻进下的工况和常规钻进下的工况均按照下述的步骤s12-s14进行；

11、s12、判断大钩速度是否小于设定的最大速度vmax，如果是，则进行后续步骤s13，如果否，则跳过步骤s13直接进行步骤s14；

12、s13、判断δt时间内的钻压wob、立管压力spp、转速rpm、进口钻井排量qi、出口钻井排量qo、灌浆排量qg这六个参数；当同时满足wob＞0，spp＞0，rpm＞0，qi＞0时，确定钻井工况为旋转钻进；当同时满足wob＞0，spp＞0，rpm＝0，qi＞0时，确定钻进工况为滑动钻进；当同时满足wob＝0，spp＞0，rpm＝0，qi＞0时，确定钻进工况为循环；当同时满足wob＝0，spp＝0，rpm＝0，qo＞0或qg＞0时，确定钻进工况为灌泥浆；如果wob、spp、rpm、qi、qo、qg的值不属于上述的任何一种情况，确定钻进工况为悬停；

13、s14、判断钻压wob、立管压力spp、转速rpm、扭矩t、进口钻井排量qi这五个参数；当同时满足wob＝0，spp＝0，rpm＝0，t＝0时，若大钩高度减小，确定钻井工况为下钻，若未减小，确定钻井工况为起钻；当同时满足wob＝0，spp＝0，rpm＞0，t＞0时，若大钩高度减小，确定钻井工况为划眼，若未减小，确定钻井工况为倒划眼；当同时满足wob＝0，spp＞0，qi＞0，rpm＞0，t＞0时，若大钩高度减小，确定钻井工况为循环划眼，若未减小，确定钻井工况为循环倒划眼。

14、所述s12中，设定的最大速度vmax应大于目标工区的最大机械钻速，并满足目标工区规定的套管内和裸眼井段避免产生激动压力的起下钻速度要求，此处推荐最大速度vmax为0.1m/s。

15、所述s2中，利用环空多相流模型计算当前钻井参数情况下的井筒压力，并反推出平衡地层压力所需要的井口回压预控值，具体包括：

16、s21、自动优选流变模式；

17、以现场实测的钻井液旋转粘度测试6档转速下的读数为输入，分别拟合出宾汉流变模式、幂率流变模式、赫巴流变模式的公式；利用三种流变模式公式分别计算6档转速下对应的剪切应力，并与测试得到的6档转速下的读数(转换为剪切应力)进行对比，计算三种流变模式公式的相关系数；选择相关系数最近接1的流变模式为优选流变模式；

18、s22、所述环空多相流模型中，液、固两相的质量守恒方程：

19、

20、式(1)中；v表示速度，m/s；下标l，s分别表示液相和固相；

21、液－固两相混合物的动量守恒方程：

22、

23、式(2)中：ff表示环空混合物摩擦压降，kg/m2/s2，下标m表示液－固混合物；θ表示井斜角；

24、环空混合物能量守恒方程：

25、

26、式中：u表示内能，m2/s2；qtotal表示环空流体与周围环境间的热量交换，kg·m/s3；aan表示环空截面积，m2；表示相单位体积的焓通量，kg·m/s3；p表示混合物的压力；

27、岩屑和液相间的滑移表示为：

28、vs＝c0vm-vsr       (4)

29、式(4)中：vs表示岩屑速度，m/s；c0表示岩屑和液相间的分布参数，vm表示环空中液相和固相混合物的速度，m/s；vsr表示岩屑沉降末速，m/s；

30、在地层流体侵入发生前，井筒中的流动为单相流，则环空任一位置处的压力、液相速度相关参数都可以确定；

31、s23、井筒地层流体侵入量与井筒压力耦合计算模型；

32、直井地层水或原油侵入量计算采用裘比公式：

33、

34、式中：pe表示地层压力，mpa；pw表示井底压力，mpa；μl表示流体粘度，mpa·s；bl表示流体体积系数，无因次；k表示有效渗透率，md；h表示有效厚度，m；re表示供给边界半径，m；rw表示井底半径，m；s表示表皮系数；

35、水平井中地层流体侵入量计算采用joshi公式：

36、

37、

38、

39、

40、式(6)、(7)中：kh表示水平渗透率，md；rp表示内部渗流场半径，m；l表示水平段长度，m；a表示椭圆长轴，m；reh表示拟圆形驱动半径，m；

41、s24、环空流动摩阻计算；

42、环空摩阻是影响井筒压力分布的一个重要因素，而单相流摩阻与流体的流动型态、流变模式和流变参数密切相关，现场应用的钻井液包括宾汉流体、幂率流体和h-b流体；

43、当宾汉流体处于层流状态时，单位长度的环空压耗表达式为：

44、

45、式(8)中：p表示压力，pa；z表示流动方向上的坐标，m；fabh表示宾汉流体环空摩阻系数，无因次；dhy代表环空水力半径，m；ρ代表流体密度，kg/m3；u代表流速，m/s；μp代表宾汉流体的塑性黏度，pa·s；τ0为宾汉流体动切力，pa；

46、当幂率流体处于层流状态时，单位长度的环空压耗表达式为：

47、

48、式(9)中：fapl表示幂率流体环空摩阻系数，无因次；k为幂率流体稠度系数，pa·s-n；n为流性指数，无因次；

49、当h-b流体处于层流状态时，单位长度的环空压耗表达式为：

50、

51、式(10)中，式中，q为钻井液排量，m3/s；θ和r为柱坐标系中的坐标轴；

52、当钻井液处于紊流状态下，其环空摩阻计算采用下式进行计算：

53、

54、对于环空范宁摩阻系数fa的求解，采用reed-pilehvari提出的公式：

55、

56、式(12)中，n为h-b流体的流性指数，无因次；reg为雷诺数，无因次；

57、s25、实验研究地层水侵入对钻井液性能的影响；

58、选取现场使用的钻井液体系进行不同含水量的性能影响评价实验，得到不同含水量条件下钻井液的密度和粘度变化；

59、发生地层水侵入后，利用实验数据对钻井液的流变参数进行修正；

60、s26、结合s21-s25，在实时读取了录井数据后，利用环空多相流模型计算当前钻井参数情况下的井筒压力，并反推出平衡地层压力所需要井口回压预控值。

61、所述步骤s2中，建立环空多相流模型的假设条件如下：

62、钻井液和固相岩屑混相，固相岩屑均匀分布在钻井液中，仅考虑固相岩屑在钻井液中质量和体积增加，计算过程视为一相；

63、混相在井筒中沿轴线方向作一维非定常流动；

64、在井筒同一过流断面上，混相具有相同的温度和压力；

65、在井筒过流断面的任一位置，混相所特有的热物性参数和流动参数均相同。

66、所述s3中，对井口回压预控值进行修正包括以下九种情况：p漏＞p孔＞pb0；p漏＝p孔＞pb0；p漏＞pb0＞p孔；pb0＞p漏＞p孔；pb0＞p漏＝p孔；p漏＝pb0＞p孔；p漏＝p孔＝pb0；p漏＞p孔＝pb0；其它修正情况。

67、所述p漏＞p孔＞pb0情况时，井口回压预控值修正如下：

68、a、p2级＝p1级

69、1)p2级≥ps上初＞ps下初，井口回压预控值修正如下：

70、①若ps上≥ps下+δp附，则ps＝ps下+δp附；

71、②若ps上＜ps下+δp附，则ps＝0.5*(ps上+ps下)；

72、2)ps上初＞p2级≥ps下初，井口回压预控值修正如下：

73、①若ps上≥ps下+δp附，则ps＝ps下+δp附；

74、②若ps上＜ps下+δp附，则ps＝0.5*(ps上+ps下)；

75、3)ps上初＞ps下初＞p2级，井口回压预控值修正如下：

76、①钻进工况

77、a、若ps下初＞pkdb，ps下初≤pklb，则ps＝pkdb；

78、b、若ps下初＞pkdb，ps下初＞pklb，则ps＝pkdb；

79、②循环工况

80、若ps下初＞pklb，则ps＝pklb；

81、③控压起下钻工况

82、若ps下初＞pktb，则ps＝pktb；

83、④需要控压的其它工况

84、若ps下初＞pkom，则ps＝pkom；

85、b、p2级＞p1级

86、1)p2级＞p1级≥ps上初＞ps下初，井口回压预控值修正如下：

87、①若ps上≥ps下+δp附，则ps＝ps下+δp附；

88、②若ps上＜ps下+δp附，则ps＝0.5*(ps上+ps下)；

89、2)p2级≥ps上初＞p1级≥ps下初，井口回压预控值修正如下：

90、①若ps上≥ps下+δp附，则ps＝ps下+δp附；

91、②若ps上＜ps下+δp附，则ps＝0.5*(ps上+ps下)；

92、3)p2级≥ps上初＞ps下初＞p1级，井口回压预控值修正如下：

93、①钻进工况

94、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

95、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

96、②循环工况

97、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

98、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

99、③控压起下钻工况

100、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

101、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

102、④需要控压的其它工况

103、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

104、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

105、4)ps上初＞p2级＞p1级＞ps下初，井口回压预控值修正如下

106、①若ps上≥ps下+δp附，则ps＝ps下+δp附；

107、②若ps上＜ps下+δp附，则ps＝0.5*(ps上+ps下)；

108、5)ps上初＞p2级≥ps下初＞p1级，井口回压预控值修正如下：

109、①钻进工况

110、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

111、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

112、②循环工况

113、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

114、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

115、③控压起下钻工况

116、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

117、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

118、④需要控压的其它工况

119、a、若ps下+δp附≤ps上，则ps＝ps下+δp附；

120、b、若ps下+δp附＞ps上，则ps＝0.5*(ps上+ps下)；

121、6)ps上初＞ps下初＞p2级＞p1级，井口回压预控值修正如下：

122、ps＝p2级。

123、所述p漏＝p孔＞pb0情况时，井口回压预控值修正如下：

124、a、p2级＝p1级

125、1)p2级≥ps上初＝ps下初，井口回压预控值修正如下

126、①若p2级≥ps下+δp附，则ps＝ps下+δp附；

127、②若p2级＜ps下+δp附，则ps＝0.5*(p2级+ps下)；

128、2)ps上初＝ps下初＞p2级，井口回压预控值修正如下

129、①钻进工况

130、a、若ps下初＞pkdb，ps下初≤pklb，则ps＝pkdb；

131、b、若ps下初＞pkdb，ps下初＞pklb，则ps＝pkdb；

132、②循环工况

133、若ps下初＞pklb，则ps＝pklb；

134、③控压起下钻工况

135、若ps下初＞pktb，则ps＝pktb；

136、④需要控压的其它工况

137、若ps下初＞pkom，则ps＝pkom；

138、b、p2级＞p1级

139、1)p2级＞p1级≥ps上初＝ps下初，井口回压预控值修正如下：

140、①若p1级≥ps下+δp附，则ps＝ps下+δp附；

141、②若p1级＜ps下+δp附，则ps＝0.5*(p1级+ps下)；

142、2)p2级≥ps上初＝ps下初＞p1级，井口回压预控值修正如下：

143、①钻进工况

144、a、若ps下+δp附≤p2级，则ps＝ps下+δp附；

145、b、若ps下+δp附＞p2级，则ps＝0.5*(p2级+ps下)；

146、②循环工况

147、a、若ps下+δp附≤p2级，则ps＝ps下+δp附；

148、b、若ps下+δp附＞p2级，则ps＝0.5*(p2级+ps下)；

149、③控压起下钻工况

150、a、若ps下+δp附≤p2级，则ps＝ps下+δp附；

151、b、若ps下+δp附＞p2级，则ps＝0.5*(p2级+ps下)；

152、④需要控压的其它工况

153、a、若ps下+δp附≤p2级，则ps＝ps下+δp附；

154、b、若ps下+δp附＞p2级，则ps＝0.5*(p2级+ps下)；

155、3)ps上初＝ps下初＞p2级＞p1级，井口回压预控值修正如下：

156、ps＝p2级。

157、所述p漏＞pb0＞p孔情况时，井口回压预控值修正如下：

158、a、p2级＝p1级

159、1)p2级≥ps上初＞ps下初，井口回压预控值修正如下：

160、①若ps上≥ps下初+δp附，则ps＝ps下初+δp附；

161、②若ps上＜ps下初+δp附，则ps＝0.5*(ps上+ps下初)；

162、2)ps上初＞p2级≥ps下初，井口回压预控值修正如下：

163、①若ps上≥ps下初+δp附，则ps＝ps下初+δp附；

164、②若ps上＜ps下初+δp附，则ps＝0.5*(ps上+ps下初)；

165、b、p2级＞p1级

166、1)p2级＞p1级≥ps上初＞ps下初，井口回压预控值修正如下：

167、①若ps上≥ps下初+δp附，则ps＝ps下初+δp附；

168、②若ps上＜ps下初+δp附，则ps＝0.5*(ps上+ps下初)；

169、2)p2级≥ps上初＞p1级≥ps下初，井口回压预控值修正如下：

170、①若ps上≥ps下初+δp附，则ps＝ps下初+δp附；

171、②若ps上＜ps下初+δp附，则ps＝0.5*(ps上+ps下初)；

172、3)ps上初＞p2级＞p1级＞ps下初，井口回压预控值修正如下：

173、①若ps上≥ps下初+δp附，则ps＝ps下初+δp附；

174、②若ps上＜ps下初+δp附，则ps＝0.5*(ps上+ps下初)。

175、所述pb0＞p漏＞p孔情况时，井口回压预控值修正如下：

176、②若p1级≥ps下初+δp附，则ps＝ps下初+δp附；

177、②若p1级＜ps下初+δp附，则ps＝0.5*p1级。

178、所述pb0＞p漏＝p孔情况时，井口回压预控值修正如下：

179、②若p1级≥ps下初+δp附，则ps＝ps下初+δp附；

180、②若p1级＜ps下初+δp附，则ps＝0.5*p1级。

181、所述p漏＝pb0＞p孔情况时，井口回压预控值修正如下：

182、②若p1级≥ps下初+δp附，则ps＝ps下初+δp附；

183、②若p1级＜ps下初+δp附，则ps＝0.5*p1级。

184、所述p漏＝p孔＝pb0情况时，井口回压预控值修正如下：

185、a、p2级＝p1级，井口回压预控值修正如下：

186、①若p2级≥δp附，则ps＝δp附；

187、②若p2级＜δp附，则ps＝0.5*p2级；

188、b、p2级＞p1级，井口回压预控值修正如下：

189、①若p1级≥δp附，则ps＝δp附；

190、②若p1级＜δp附，则ps＝0.5*p1级。

191、所述p漏＞p孔＝pb0情况时，井口回压预控值修正如下：

192、a、p2级＝p1级，井口回压预控值修正如下：

193、①若ps上≥δp附，则ps＝δp附；

194、②若ps上＜δp附，则ps＝0.5*ps上；

195、2)ps上初＞p2级≥ps下初，井口回压预控值修正如下：

196、①若ps上≥δp附，则ps＝δp附；

197、②若ps上＜δp附，则ps＝0.5*ps上；

198、b、p2级＞p1级

199、1)p2级＞p1级≥ps上初＞ps下初，井口回压预控值修正如下：

200、①若ps上≥δp附，则ps＝δp附；

201、②若ps上＜δp附，则ps＝0.5*ps上；

202、2)p2级≥ps上初＞p1级≥ps下初，井口回压预控值修正如下：

203、①若ps上≥δp附，则ps＝ps下+δp附；

204、②若ps上＜δp附，则ps＝0.5*ps上；

205、3)ps上初＞p2级＞p1级＞ps下初，井口回压预控值修正如下：

206、①若ps上≥δp附，则ps＝δp附；

207、②若ps上＜δp附，则ps＝0.5*ps上。

208、所述其它修正情况时，井口回压预控值修正如下：

209、若计算结果ps＜0，则ps＝0。

210、上述式中，pb0：不控压条件下的井底压力，pb0＝p静+p环耗，为自动输入，mpa；

211、p静：环空静液柱压力，为自动输入，mpa；

212、p环耗：环空循环压耗，为自动输入，mpa；

213、p漏：地层漏失压力，为自动输入，mpa；

214、p孔：地层孔隙压力，为自动输入，mpa；

215、ps上初：预测控压上限，优化前控压窗口上限，为自动输入，mpa；

216、ps下初：预测控压下限，优化前控压窗口上限，为自动输入，mpa；

217、ps上：推荐控压上限，优化后控压窗口上限，为自动输入，mpa；

218、ps下：推荐控压下限，优化后控压窗口上限，为自动输入，mpa；

219、ps：推荐控压值，为自动输入，mpa；

220、pdp：控压设备许可工作压力，为手动输入，mpa；

221、δp附：控压安全附加值，为手动输入，mpa；

222、p1级：推荐控压值1级报警，为手动输入，pkdb(钻进)、pklb(循环)、pktb(起下钻)、pk0b(其它)，mpa；

223、p2级：推荐控压值2级报警，为手动输入；pkdb(钻进)、pklb(循环)、pktb(起下钻)、pkom(其它)，mpa；

224、ρ漏：地层漏失压力梯度，根据井眼轨迹输入对应垂深处的漏失压力梯度，g/cm3；

225、ρ孔：地层孔隙压力梯度，根据井眼轨迹输入对应垂深处的孔隙压力梯度。

226、本发明还提供了一种计算机设备，包括处理器、输入设备、输出设备和存储器，处理器、输入设备、输出设备和存储器相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，执行如本发明所述的方法。

227、本发明还提供了一种计算机可读存储介质，其中，上述计算机可读存储介质存储有用于电子数据交换的计算机程序，其中，上述计算机程序使得计算机执行如本发明中所描述的部分或全部步骤。

228、与现有技术相比，本发明所带来的有益的技术效果表现在：

229、1、本发明基于大钩重量、大钩速度、钻压、转速、排量、立管压力等实时工程录井参数，建立了自动判断控压钻井作业工况的分支判断树，利用阈值法进行多参数阈值判断，实现起钻、下钻、滑动钻进、旋转钻井等22种钻井工况(含控压条件下11种)的实时识别，为后续井筒压力实时监测和计算提供基础参数。

230、2、本发明提出一种基于在线监测的井筒压力控制智能决策方法是基于致密油注水开发区的控压钻井技术研发，在理论方法方面建立了控压钻井正常钻进、油侵、水侵三种情况下的井筒流动理论模型，解决了控压钻井作业参数的实时模拟和数据实时反馈，基于方法的数据模拟精度≥95％。

231、3、本发明，实现了井口回压预控值的实时计算。基于环空多相流物理模型，考虑井筒压力一地层溢流耦合的影响，建立了控压钻井期间环空瞬态多相流动计算模型、自动判断优选钻井液流变模式，利用环空多相流模型计算井筒压力，并实时预测地层流体侵入后井筒压力的变化，反推出平衡地层压力所需要的井口回压预控值。

232、4、本发明，实现井口回压预控值修正和井筒压力控制智能决策。针对溢漏同存等窄安全密度窗口地层，基于环空多相流物理模型，考虑井筒压力-地层溢流耦合的影响，建立了控压钻井期间环空瞬态多相流动计算模型。结合地层安全作业压力窗口、控压限制条件(标准、设备能力等)等多个因素，建立了8类井况、29种控压限定关系下的控压范围、井口回压预控值和决策建议推荐优选机制，实现保持安全钻井的控压窗口范围、井口回压预控值的计算、确定和推荐，克服了现有控压系统在无井底压力测量数据条件下的实时分析计算能力有限，缺乏作业井况智能分析、参数自动优化等技术，不能根据工程参数、井况的变化进行准确的控压作业实时分析计算与控制决策推荐的技术缺陷。

233、5、本发明，可以更全面、准确地指导致密油注水开发区控压钻井现场操作，配合控压钻井工艺技术要求，实现实时、快速、安全的井口压力控制作业。本发明可在控压钻井过程中，准确识别水侵或油侵，自动判断工况，确定环空流型，预测地层流体侵入后井筒压力的变化，确定合理的井口压力和控压施工建议，实现全自动控压作业。进一步解决致密油注水开发区储层出水控制难、窄安全密度窗口地层钻进复杂等技术难题，为致密油注水开发区以及其他需要应用控压钻井解决的复杂区域、风险井的安全钻井提供技术支持。

234、6、本发明，基于在线监测的井筒压力控制智能决策方法，能够开发1套基于在线监测的井筒压力控制智能决策计算机程序，应用在自动控压钻井设备上，提高井口套压准确性和及时控制水平，可对现有设备以及未来设备进行升级改造，助力自动控压及精细控压技术规模化应用。

235、7、本发明，通过配备基于在线监测的井筒压力控制智能决策方法，可提高作业效率进而降低井控风险。

236、8、本发明，基于在线监测的井筒压力控制智能决策方法中的控压钻井复杂判定可以为复杂与风险快速判定和报警提供依据，实现事故预防和快速处置。