1.一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,包括:步骤S10、根据测井数据或室内试验数据,确定地质力学参数和岩石力学参数;
步骤S20、选择多种常见概率分布函数对地质力学参数和岩石力学参数进行分布函数拟合,获得地质力学参数、岩石力学参数的分布拟合结果;
步骤S30、根据K‑S检验对地质力学参数、岩石力学参数的分布拟合结果进行计算,获得K‑S检验值并排序;
步骤S40、根据Kendall相关分析,对地质力学参数和岩石力学参数进行相关性分析获得相关系数矩阵;
步骤S50、基于坍塌压力和破裂压力模型,建立考虑参数相关性影响的井壁失稳风险可靠度评价模型;
步骤S60、根据步骤S30获得的K‑S检验值、步骤S40获得的相关系数矩阵,进行蒙特卡洛模拟,并带入井壁失稳风险可靠度评价模型中计算井壁失稳风险。
2.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述地质力学参数包括最大水平主应力、最小水平主应力、孔隙压力,所述岩石力学参数包括内聚力、内摩擦角、抗张强度。
3.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述室内实验数据的室内试验样本不低于50个。
4.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述步骤S30中K‑S检验值的计算公式为:Dn=max|Sn(x)‑F(x)|式中:Dn为K‑S检验值;Sn(x)为Ai(xi)的累积分布函数;Fn(x)为理论上的累积分布函数;
kn为小于或等于x的xi的个数;n为事件数;f(y)为理论上的概率分布函数。
5.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述步骤S40中的计算公式为:
其中:
式中:τ(x,y)为变量x与y的Kendall相关系数;xj为变量x的第j个观测值;yj为变量y的第j个观测值;ti为变量x的第i组结点数量;μi为变量y的第i组结点数量;Ω为函数自变量。
6.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述步骤S50中的井壁失稳风险可靠度评价模型如下:式中:Pfc为井壁稳定概率;Prc为井壁失稳概率;Zc为功能函数;fRc为抗力随机变量分布函数;fQc为载荷随机变量分布函数;Qc为载荷随机变量。
7.根据权利要求6所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述步骤S60的具体过程为:
步骤S61、对地质力学参数、岩石力学参数的分布拟合结果进行排序,将K‑S检验值最小的作为最佳分布函数,对地质力学参数和岩石力学参数的输入分布进行设定;
步骤S62、对地质力学参数和岩石力学参数的分布进行相关系数设定,如果出现一致性错误时,须对无效矩阵进行调整,使调整后的有效矩阵最接近无效矩阵;
步骤S63、采样蒙特卡洛方法生成10000组符合要求的随机数,将其代入井壁失稳模型计算井壁失稳风险。
8.根据权利要求7所述的一种复杂地层钻井井壁失稳风险定量评价方法,其特征在于,所述步骤S62中的调整过程如下:①找到最小特征值E0;②将‑E0与单位矩阵(I)相乘,再加到相关系数矩阵(C)上,使得特征值进行调整,使最小特征值等于C’=C-E0I。