发明内容
本发明所要解决的技术问题是提供一种控压钻井井身结构设计方法,针对目前控压钻井井身结构设计需依赖人工经验确定裸眼井段最优的控压参数及井筒ECD剖面而导致的工作量大、设计缺乏科学依据、结果可靠性难以保障的问题,通过引入人工蜂群智能算法,高效处理分析不同控压参数条件下的井筒ECD剖面样本,准确优选裸眼井段最佳控压参数及井筒ECD剖面,为控压钻井井身结构设计提供科学依据。
为了解决上述问题,本发明公开了一种基于人工蜂群算法的控压钻井井身结构设计方法,该方法包括以下步骤:
定义井筒ECD剖面为井筒ECD沿井深的分布曲线。
(一).地质必封点确定:依据地质资料,确定必封层位;
(二).安全压力窗口确定:根据钻前地层孔隙压力、地层破裂压力以及地层坍塌压力的预测结果,依据区域井身结构设计系数,确定安全压力窗口,设定井深H处的安全压力窗口的下限为ρ 下(H),上限为ρ 上(H);
(三).控压参数可调范围确定:
控压参数对井筒ECD剖面有直接影响,根据不同控压钻井方式的控压特点及实际工况确定一个或多个控压参数K i (i=1,2,……),以及各控压参数的可调范围:
K L≤K i ≤K H
其中,K i 为第i个控压参数,K L为控压参数的下限,K H为控压参数的上限;
不同的控压钻井方式的控压参数也不一样;另一方面,不同的控压设备的控压参数可调范围也不一样。比如:井底恒压控压钻井方式,控压参数主要有钻井液密度、井口回压。CML双梯度控压钻井方式,控压参数主要有泥浆泵安装位置、泥浆液面深度、排量。
(四).不同控压条件下的井筒ECD计算:
在步骤三中的各控压参数的可调范围内,利用井筒ECD计算模型,可计算不同控压条件下的井筒ECD剖面,其中所述的ECD计算模型指的是ECD的算法,也可以称为井筒压力(场)计算、钻井水力学计算、循环当量密度计算、当量循环密度计算,为成熟的现有技术。
控压参数组合方案X j 条件下的井筒ECD为井深H的函数:
ECD (H) Xj =f (K i , H) (1)
其中,ECD (H) Xj 为控压参数组合方案X j 条件下的井筒ECD;K i 为控压参数组合方案Xj条件下各控压参数。控压参数组合方案Xj为1个或多个控压参数分别取一个值所组成的组合方案。
(五).基于人工蜂群算法的裸眼井段最优控压参数及套管下深确定,其中所述的最优控压参数指裸眼井段最大下深对应的控压参数。
在步骤二所确定的安全压力窗口[ρ 下(H), ρ 上(H)]内,在一定的控压参数组合方案Xj条件下,井筒ECD剖面的最大安全井深H1, Xj 可表示为:
H1,Xj=f -1(ECD (H) Xj ) (2)
其中,ECD (H) Xj 为控压参数组合方案Xj条件下的井筒ECD,它是井深H的函数;X j 为第j个控压参数的组合方案,j=1,2,……。
对于某一控压参数组合方案Xj,裸眼井段内井深H i 处的井筒ECD剖面满足:ρ 下(H i )≤ ECD (H i ) Xj ≤ ρ 上(H i )。
其中ECD (H i ) Xj 为某一控压参数组合方案Xj,裸眼井段内井深H i 处的井筒ECD剖面。
则,裸眼井段的套管最大下深H1为:
H1=max (H1,Xj)=max ( f -1(ECD (H) Xj ))=max (f -1(f (K i , H))= H1,Xj0 (3)
其中,Xj0为套管最大下深对应的控压参数组合方案。
则,该裸眼井段的最大套管下深为H1所对应的控压参数组合方案Xj0。
以上为裸眼井段套管下深及最优控压参数的基本计算方法,但由于控压参数的种类往往多于2个,且各控压参数的组合方案较多(精度要求越高,各控压参数的计算样本越多,组合方案越多),不同参数组合方案下的井筒ECD剖面差异也较大,极大降低了计算效率。为此,引入人工蜂群算法,对套管最大下深H1及其对应的最优控压参数组合方案进行高效寻优。
人工蜂群算法基本步骤描述如下:
S1、 初始化种群:确定蜂群总数
SN、引领蜂与跟随蜂数量(一般同时设为
SN/2)、最大迭代次数
MCN与控制参数
Limit、确定目标函数H=
f -1(
f (
K i ,
H))、各参数搜索范围{
K L≤
K i ≤
K H;
ρ 下(
H i ) ≤ ECD(
H i )
Xj ≤
ρ 上(
H i ),
H i ∈[
H 0 ,H end ]}、并在搜索范围内随机生成初始解
,所述的初始解即人工蜂群算法中的侦查蜂,解指控压参数组合方案。引领蜂和跟随蜂指前文各控压钻井参数组合方案;最大迭代次数和控制参数需要根据实际计算需求进行设定,一般最大迭代次数设为10000,控制参数设定为20。H
0为安全钻井密度窗口最小深度,H
end为安全钻井密度窗口最大深度,其中H
0与H
end可通过安全压力窗口[
ρ 下(
H),
ρ 上(
H)]得出,为本领域的公知常识。
S2、 计算并评估每个初始解的适应度,本发明因为是最大值优化问题,所以适应度函数fitness设为目标函数H=f -1(f (K i , H))本身,将目标函数的解即适应度值的最大输出结果作为本算法的优化目标。
S3、 引领蜂利用公式(4)对邻域进行搜索,得到新的解,按照贪婪选择原则,如果新的解的适应度大于原先解的适应度,则用新的解更新原解,否则,保持原解不变;随着优化步骤的迭代进行,初始解会随着适应度最优搜索解的出现被替换,而迭代的反复进行,重新搜索的解若适应度高于原有解,则将新解对原有解进行替换。当所有引领蜂完成邻域搜索后,根据式(5)计算概率Pj。
式中,K
i(Xj)指第j个蜂群个体第i维的参数,即第j个控压参数组合方案X
j的第i个控压参数的值;同理,K
i(Xj´)指第j´个蜂群个体第i维的参数,即第j´个控压参数组合方案X
j´的第i个控压参数的值,
为在X
j附近产生的一个新的解,
为控压参数组合方案
中第i个控压参数的值,具体的,通过公式(4)求取新的解的第i个控压参数的值,为了控制收敛速度,对每一个蜂群个体每一次迭代仅随机选取一个控压参数进行搜索,更换X
j中第i个控压参数的值为公式(4)求取的值
后得到新的解
,若得到的新的解的拟合度高于原解,则将新解对原有解进行替换;
,
且
,
,
D为问题的维数,即指控压参数的种类;
为[-1,1]之间的随机数;
为第
j个解的适应度值;P
j为跟随蜂选择第j个引领蜂的概率。
S4、跟随蜂利用算得的概率Pj并基于轮盘赌原则选择引领蜂,当跟随蜂完成引领蜂的选择后,利用式(4)对邻域搜索,同样按照贪婪选择原则选择适应度高的解。
S5、判断是否有要放弃的解,如果某引领蜂在Limit次迭代后,结果都没有改变,即最优的解都没有发生改变,则由引领蜂变为侦察蜂,并按照公式(6)随机产生一个新的解取代原解。
式中,
,
与
分别表示第i个控压参数的组合方案的上下顶界,因为每个解中包含一个或多个控压参数,K
imax-K
imin为第i个控压参数的上下顶界求差值后得到的控压参数变化幅度,即通过公式6中将每个控压参数进行分别运算后得到新的解的每个控压参数K
i(Xj)后组成新的解。
S6、 记录到目前为止的蜂群的最优解,即适应度最高的解。
S7、 如果当前迭代次数大于最大次数MCN,则迭代结束,算法结束;否则转向步骤S2,Cycle=Cycle+1,即迭代次数+1。本步骤指整个蜂群的整体循环迭代一次;步骤S5指蜂群中的某一引领蜂个体循环迭代的次数。
发明内容中提供的效果仅仅是实施例的效果,而不是发明所有的全部效果,上述技术方案中的一个技术方案具有如下优点或有益效果:
针对目前控压钻井井身结构设计需依赖人工经验确定裸眼井段最优的控压参数及井筒ECD剖面而导致的工作量大、设计缺乏科学依据、结果可靠性难以保障的问题,通过引入人工蜂群智能算法,高效处理分析不同控压参数条件下的井筒ECD剖面样本,准确优选裸眼井段最佳控压参数及井筒ECD剖面,不但减轻了设计人员的计算工作量,而且提高了设计过程的科学性和设计结果准确性。
具体实施方式
A井地层孔隙压力和地层破裂压力信息如图2所示,目的井深6600m。该地区地质必封点主要为200m处松散砾石层和1700m处易漏失地层,井身结构设计系数如下:抽吸压力系数为0.036 g/cm3,激动压力系数为0.036 g/cm3,破裂压力安全系数为0.024 g/cm3,井涌允量为0.06 g/cm3,压差卡钻系数为22 MPa。该井上部常规地层(1700m以上)采用常规钻井方式钻进,下部复杂压力系统地层采用井底恒压控压钻井方式钻进,控压参数主要为钻井液密度和井口回压,其中,A井井口控压设备的最大控压值为5 MPa。
地质必封点确定。依据A井资料,确定必封层位为200m和1700m。
安全压力窗口确定。根据已有的地层孔隙压力和地层破裂压力剖面,依据该地区的井身结构设计系数,可以确定A井的安全压力窗口,如图2所示。
控压参数可调范围确定。
由于采用井底恒压控压钻井方式,控压参数为钻井液密度和井口回压。
钻井液密度的可调范围需根据安全压力窗口分开次进行确定,以裸眼井段安全压力窗口下限的最小值为钻井液密度的可调范围下限,裸眼井段安全压力窗口上限的最小值为钻井液密度的可调范围上限。
根据该井的资料,1700m以上为常规钻井方式,1700m(即第一层技术套管)以下开始采用控压钻井方式。以1700m以下裸眼井段为例,该裸眼井段安全压力窗口下限的最小值为1.18 g/cm3,安全压力窗口上限的最小值为1.47 g/cm3,则该裸眼井段钻井液密度的可调范围为:1.18~1.47 g/cm3。
井口回压的可调范围与控压装备的性能有关,则A井的井口回压可调范围为0~5MPa。
不同控压条件下的井筒ECD计算。
在步骤3中设计的各控压参数的可调范围内,利用井筒ECD计算模型,可计算不同控压条件下的井筒ECD剖面。
控压参数组合方案Xj条件下的井筒ECD为井深H的函数:
ECD (H) Xj =f (ρ d, P 井口, H)
其中,ECD (H) Xj 为控压参数组合方案Xj条件下的井筒ECD;ρ d为静态钻井液密度,1.18~1.47 g/cm3;P 井口为井口回压,0~5 MPa。
基于人工蜂群算法的裸眼井段最优控压参数及套管下深确定。
在步骤2所确定的安全压力窗口[ρ 下(H), ρ 上(H)]内(图2所示,安全压力窗口下限为ρ 下(H),安全压力窗口上限为ρ 上(H)),利用人工蜂群算法对裸眼井段的套管最大下深H1及其对应的最优控压参数组合方案进行高效寻优。
算法步骤如下:
步骤(1)初始化蜂群:蜂群总数SN为50;其中引领蜂与跟随蜂数量均为25;各蜂群个体维数为2,即指控压参数的种类为2,即钻井液密度ρ d和井口回压P 井口;最大迭代次数MCN为1000;控制参数Limit为20;当前迭代次数Cycle为1;目标函数为H=f -1(f (ρ d, P 井口, H)),各参数搜索范围为{ρ d(1.18~1.47 g/cm3)、P 井口(0~5 MPa)、ρ 下(H i ) ≤ ECD(H i ) Xj ≤ ρ 上(H i ),H i ∈[1700,H Xj ]}。
步骤(2):计算并评估每个初始解的适应度,适应度函数fitness设为目标函数H=f -1(f (ρ d, P 井口, H))。
步骤(3):引领蜂利用公式(4)对各维参数邻域进行搜索,得到新的解,按照贪婪选择原则,对原解进行更新。当所有引领蜂完成邻域搜索后,根据式(5)计算概率Pi。
步骤(4):跟随蜂利用算得的概率
并基于轮盘赌原则选择引领蜂,当跟随蜂完成引领蜂的选择后,利用式(4)对邻域搜索,同样按照贪婪选择原则选择适应度高的解。
步骤(5)判断是否有要放弃的解:如果某引领蜂在20次迭代后,最优的解都没有改变,则由引领蜂变为侦察蜂,并按照公式(6)随机产生一个新的解取代原解。
步骤(6):记录到目前为止的蜂群的最优解。
步骤(7):如果当前迭代次数大于最大次数1000,则迭代结束,算法结束;否则转向步骤(2),Cycle=Cycle+1。
以1700m以下裸眼井段为例,通过人工蜂群算法对控压参数进行组合寻优(如图3所示)发现该开次最大下深稳定在3685m附近,该方案下的最优控压参数为:ρ d=1.44 g/cm3,P 井口=1.3 MPa,该控压参数条件下的裸眼井段ECD如图2所示,其对应的套管最大下深为3685m,即第二层技术套管的设计下深。
同样地,可以分别得到第三层技术套管的设计下深为5260m,对应的控压参数为:ρ d=1.51 g/cm3,P 井口=5 MPa;生产套管下深为6600m,对应的控压参数为:ρ d=1.45 g/cm3,P 井口=2 MPa。
井身结构设计结果及各开次控压参数如表1所示。
表1 A井井身结构设计结果及各开次控压参数