1.一种复杂地层钻井井壁失稳风险定量评价方法，其特征在于，包括：步骤S10、根据测井数据或室内试验数据，确定地质力学参数和岩石力学参数；

步骤S20、选择多种常见概率分布函数对地质力学参数和岩石力学参数进行分布函数拟合，获得地质力学参数、岩石力学参数的分布拟合结果；

步骤S30、根据K‑S检验对地质力学参数、岩石力学参数的分布拟合结果进行计算，获得K‑S检验值并排序；

步骤S40、根据Kendall相关分析，对地质力学参数和岩石力学参数进行相关性分析获得相关系数矩阵；

步骤S50、基于坍塌压力和破裂压力模型，建立考虑参数相关性影响的井壁失稳风险可靠度评价模型；

步骤S51、根据坍塌压力模型，建立井壁坍塌失稳功能函数，进一步，建立井壁坍塌失稳的可靠度和失效概率模型；

目前，井壁稳定分析中，坍塌压力计算模型可表示：式中，pc为坍塌压力；pp为孔隙压力；σH为最大水平主应力；σh为最小水平主应力；C0为内聚力；为内摩擦角；A为中间系数；α为Biot系数，0.85；η为应力非线性修正系数，取值为η＝0.95；

对于井壁坍塌失稳，根据可靠度理论，将影响井壁坍塌失稳的因素分为两类综合量，即载荷随机变量Qc和抗力随机变量Rc，载荷随机变量Qc表示坍塌压力，受到水平最大主应力、水平最小主应力、孔隙压力、内聚力、内摩擦角的影响；抗力随机变量Rc表示井筒压力pi；

因此，井壁坍塌失稳功能函数Zc可表示为：假定坍塌压力Qc和井筒压力Rc相互独立，即fQc(Qc)和fRc(Rc)为两个独立的随机变量分布函数，由于描述功能函数的基本变量为随机变量，可以认为是连续函数，因此，功能函数Zc为概率分布函数，井壁稳定概率Pfc与井壁坍塌失稳概率Prc和为1，则井壁坍塌失稳的可靠度和失效概率分别为：

式中：Pfc为井壁稳定概率；Prc为井壁坍塌失稳概率；Zc为功能函数；fRc为抗力随机变量分布函数；fQc为载荷随机变量分布函数；Qc为载荷随机变量；

步骤S52、根据破裂压力模型，建立井壁破裂失稳功能函数，进一步，建立井壁破裂失稳的可靠度和失效概率模型；

目前，井壁稳定分析中，破裂压力计算模型可表示：pf＝3σh‑σH‑αpp+St式中：pf为破裂压力，MPa；pp为孔隙压力；σH为最大水平主应力；σh为最小水平主应力；St为抗拉强度；α为Biot系数，0.85；

步骤S60、根据步骤S30获得的K‑S检验值、步骤S40获得的相关系数矩阵，进行蒙特卡洛模拟，并带入井壁失稳风险可靠度评价模型中计算井壁失稳风险。

2.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法，其特征在于，所述地质力学参数包括最大水平主应力、最小水平主应力、孔隙压力，所述岩石力学参数包括内聚力、内摩擦角、抗张强度。

3.根据权利要求1所述的一种复杂地层钻井井壁失稳风险定量评价方法，其特征在于，所述步骤S60的具体过程为：

步骤S61、对地质力学参数、岩石力学参数的分布拟合结果进行排序，将K‑S检验值最小的作为最佳分布函数，对地质力学参数和岩石力学参数的输入分布进行设定；

步骤S62、对地质力学参数和岩石力学参数的分布进行相关系数设定，如果出现一致性错误时，须对无效矩阵进行调整，使调整后的有效矩阵最接近无效矩阵；

步骤S63、采样蒙特卡洛方法生成10000组符合要求的随机数，将其代入井壁失稳模型计算井壁失稳风险。

4.根据权利要求3所述的一种复杂地层钻井井壁失稳风险定量评价方法，其特征在于，所述步骤S62中的调整过程如下：①找到最小特征值E0；②将‑E0与单位矩阵I相乘，再加到相关系数矩阵C上，使得特征值进行调整，使最小特征值等于C’＝C－E0I。