1.一种基于井径仪旋转轨迹与测量数据的地应力预测方法,其特征在于,所述方法包括如下步骤:步骤S1:通过井径仪测量数据,分析井径仪在上提测量过程中的旋转轨迹,找寻目的层段中井眼形状为椭圆形的井段,并记录下椭圆形井眼的长轴测量值 与短轴测量值 ;步骤S2:将目的层段密度测井与声波测井数据代入波动方程,求取目的层段的动态杨氏模量 与动态泊松比 ;再将动态杨氏模量 与动态泊松比 代入目的层段的动静态转换经验公式中,求取目的层段的静态杨氏模量 与静态泊松比 ;步骤S3:根据目的层段的钻井施工数据,获取作用于井壁上的泥浆压力 P m ;通过对目的层段密度测井数据进行积分,求取上覆岩层压力 v ;步骤S4:获取目的层段压裂施工数据,根据压裂施工数据中的停泵压力 P t ,计算获取最小水平主应力 h ;步骤S5:将步骤S1-S4中的参数代入平面应变状态下井壁应力-应变公式中,通过公式计算求得两向应力差 和最大水平主应力 H ,实现目的层段的地应力预测;其中,在步骤S2中,求取的目的层段的动态杨氏模量 E 动 与动态泊松比 分别为:上式中, 为纵波波速; 为横波波速; 为动态杨氏模量; 为动态泊松比; 为岩石密度;在步骤S2中,所述目的层段的动静态转换经验公式的获取方式为:通过测井数据获取目的层段的动态杨氏模量 与动态泊松比 ,通过室内岩心实验获取目的层段的静态杨氏模量 与静态泊松比 ,然后通过数学拟合获取目的层段的动静态转换经验公式。
2.根据权利要求1所述的方法,其特征在于,在步骤S2中,将目的层段密度测井与声波测井数据代入波动方程,所述波动方程为:上式中, 为纵波波速; 为横波波速; 为动态杨氏模量; 为动态泊松比; 为岩石密度。
3.根据权利要求1所述的方法,其特征在于,在步骤S4中,所述最小水平主应力 h 的计算公式为:上式中, h 为最小水平主应力; P t 为停泵压力; 为岩石密度; 为重力加速度; 为井筒中静液柱高度。
4.根据权利要求1所述的方法,其特征在于,在步骤S5中,所述井壁应力-应变公式为:上式中, 为长轴测量值; 为短轴测量值; r 为初始孔眼半径; 为静态杨氏模量; 为静态泊松比; 为两向应力差; v 为上覆岩层压力; P m 为泥浆压力; h 为最小水平主应力。
5.根据权利要求1所述的方法,其特征在于,在步骤S5中,将所述井壁应力-应变公式化简,通过下述公式计算得到两向应力差 :上式中, 为长轴测量值; 为短轴测量值; r 为初始孔眼半径; 为静态杨氏模量; 为静态泊松比; 为两向应力差; v 为上覆岩层压力; P m 为泥浆压力; h 为最小水平主应力。
6.根据权利要求5所述的方法,其特征在于,在步骤S5中,所述最大水平主应力 的计算公式如下:上式中, 为最大水平主应力; h 为最小水平主应力; 为两向应力差。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被计算机执行时,使所述计算机执行如权利要求1-6中任一项所述的方法。
8.一种电子设备,所述电子设备包括处理器、存储器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行程序时,实现如权利要求1-6任一项所述的方法。